解稀疏对称方程组的iccg法fortran源程序.txt

解稀疏对称方程组的ICCG法fortran源程序 有算例

解稀疏对称方程组的ICCG法fortran源程序


求解稀疏对称方程组的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 ]      
  输入数据中的其它符号----双精度型数组,用来存放方程组的右端项 
程序见附件