【分享】解稀疏对称方程组的ICCG法fortran源程序

zhujd2003

求解稀疏对称方程组的ICCG法源程序  
PROGRAM ICCGSOLVER----主程序  
SILLT-----不完全三角分解子程序  
ICCG------ICCG迭代子程序  
SINA------被子程序ICCG调用的子程序,可作三角矩阵(LDLt)-1与相应列矩阵的乘法  
SMI-------被子程序ICCG调用的子程序,可作原系数矩阵与相应列矩阵的乘法  
说明:本程序采用文件INPUT.DAT构成数据模块  
    ND(NND)---整型数组,NND为未知数总数,用来存放紧缩形式的系数矩阵中对角线元素  的序号.  
    NC(NNC)---整型数组,NNC为系数矩阵中非零元素的总数,用来存放紧缩形式的系数矩阵中各元素的列号.  
    AA(NNC)---双精度型数组,用来存放L和D的乘积;  
    SS(NNC)---双精度型数组,用来存放紧缩形式的系数矩阵中的非零元素  
给出一个6*6阶矩阵:  
    原系数矩阵A=[ 1.1   0    0    0    0    0 ]  
                          0   2.5  2.2  7.6   0    0   
                          0   2.2  4.0   0   9.2  5.3  
                          0   7.6   0   1.3   0    0  
                          0    0   9.2   0    3.4  0            
                        [ 0    0   5.3   0     0  1.6 ]
  矩阵A的紧缩存储形式 SS=[ 1.1 2.5 2.2 4.0 7.6 1.3 9.2 3.4 5.3 1.6]  
     ND=[1 2 4 6 8 10]  
     NC=[1 2 2 3 2 4 3 5 3 6]  
     L=[1   0    0    0    0   0]  
          0   1    0    0    0   0  
          0  L32  1    0    0   0  
          0  L42  0    1    0   0  
          0   0   L53  0    1   0  
         [0   0   L63  0    0  1]  
     
      D= [D1   0   0   0   0   0]  
             0  D2   0   0   0   0  
             0   0  D3   0   0   0  
             0   0   0   D4  0   0  
             0   0   0   0   D5  0  
            [0   0   0   0   0  D6]   
  
     LD= [ D1    0       0     0   0   0  
               0   D2       0     0   0   0  
               0  L32D2  D3    0   0   0  
               0  L42D2   0    D4  0   0                          
               0    0   L53D3   0   D5 0         
               0    0   L63D3   0   0  D6]            
  矩阵LD的紧缩存储形式为  
       SS=[D1  D2  L32D2  D3  L42D2  D4  L53D3  D5  L63D3   D6 ]      
  输入数据中的其它符号----双精度型数组,用来存放方程组的右端项
程序见附件:  
       

zhujd2003

这是一本书里的一个程序,可我看不太懂.
想请这里的高手们帮分析!
不胜感激.

zhgao520

这程序也只能解个数少的方程组,我拿它解了一个42个方程组,结果就不对阿,有高手知道要改那儿吗,请指点

homage79

徐士良的《常用算法程序集》里面有相关的文件，可以参考

zhgao520

我这儿没那本书,网上的下下来那么多的程序,没一个说明,也不知是哪个?不知楼上的能不能给传一下,先谢谢了