带状对角矩阵的压缩及乘法运算算法实现

算法名称：带状对角矩阵乘法

 

带状对角矩阵定义：

       带状对角系统紧靠在对角线左边（下边）有m₁≥0个非零元素，紧靠其右边（上边）有m₂≥0个非零元素。当然，这仅在m₁和m₂<<N时才有意义。这种情况下，用LU分解求解线性系统可以比通用N×N的情况完成的更快，占用空间更少。

       用a_ij精确定义带状对角矩阵，表示为

a_ij = 0 当 j > i + m₂ 或 i > j + m₁

带状对角矩阵用一种称为压缩形式的方式进行存储及操作，把矩阵顺时针旋转45°，是非零元素排成一个细长的矩阵，它有m₁+1+m₂列和N行。例如，

其中N=7，m1=2，m2=1，存储在压缩形式为7×4矩阵

其中x表示压缩格式中浪费的空间，它们在运算中不被引用，因而其值可随意。

注：原矩阵对角线元素现在第m₁列中，对角线左下的元素在其左侧，对角线右上的元素在其右侧。

 

运行示例：

The origin left matrix (before compression):

3.0 1.0 0.0 0.0 0.0 0.0 0.0

4.0 1.0 5.0 0.0 0.0 0.0 0.0

9.0 2.0 6.0 5.0 0.0 0.0 0.0

0.0 3.0 5.0 8.0 9.0 0.0 0.0

0.0 0.0 7.0 9.0 3.0 2.0 0.0

0.0 0.0 0.0 3.0 8.0 4.0 6.0

0.0 0.0 0.0 0.0 2.0 4.0 4.0

----------------------------------

The new left matrix (after compression):

0.0 0.0 3.0 1.0

0.0 4.0 1.0 5.0

9.0 2.0 6.0 5.0

3.0 5.0 8.0 9.0

7.0 9.0 3.0 2.0

3.0 8.0 4.0 6.0

2.0 4.0 4.0 0.0

----------------------------------

The right matrix:

1.0

2.0

3.0

4.0

5.0

6.0

7.0

----------------------------------

Solution matix:

5.0

21.0

51.0

98.0

84.0

118.0

62.0

 示例程序：

package  com.nc4nr.chapter02.banmul;

public   class  BanMul  ... {

    double<