解稀疏對稱方程組的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 ]      
  輸入數據中的其它符號----雙精度型數組,用來存放方程組的右端項 
程序見附件