论坛 › 编程工具交流 › 使用lapack的dsptrf输出上下三角矩阵计算结果不同问题

nong 发表于 2025-2-17 11:26:46

使用lapack的dsptrf输出上下三角矩阵计算结果不同问题

同一个对称矩阵用lapack的dsptrf计算厚上三角存储与下三角存储不同。   7   4000.00000      0.00000         0.00000       0.00000       1500.00000       0.00000    2000.00000
      0.00000   6685.75977      1735.68018   768.00000      -5000.00000       0.00000       0.00000
      0.00000   1735.68018      3448.23999   924.00000          0.00000    -750.00000   500.00000
      0.00000      768.00000       924.00000    7200.00000          0.00000   -1500.00000    2000.00000
矩阵A
1500.00000    -5000.00000         0.00000       0.00000       5750.00000       0.00000    1500.00000
      0.00000      0.00000      -750.00000   -1500.00000          0.00000    5750.00000   -1500.00000
   2000.00000      0.00000       500.00000    2000.00000       1500.00000   -1500.00000    8000.00000

   矩阵B   0.00000000       0.00000000      -37.0000000      -4.25000000       0.00000000      -12.0000000       8.00000000
求A*X=B    call dspsv(uplo,n,nrhs,ap,ipiv,b,ldb,info)
L
         0         1         2         3         4         5         6         7
   4000.0000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000       0.37500000000000000      0.0000000000000000       0.50000000000000000      6685.7597699999997       0.25960851716333805       0.11487101338072755      -0.74785816003077832      0.0000000000000000      0.0000000000000000      2997.6426322004045       0.24173016857805102       0.43302112529133219      -0.25019660180422049       0.50039320360844097      6936.6163877122535      3.7565723581641050E-002 -0.19010743853479539       0.23605237130216328      876.34049945563072       0.42712148195006694      4.4455343516748964E-002   5151.7842019021118      -0.16112128318997046      5727.4243615039459   

X=
-9.2481364876117224E-003   2.1208594097849018E-002-2.3487786720192479E-002-1.8023544606899295E-003   1.9928814217908649E-002-4.6947614371289238E-003   3.5496623117919584E-003

U
         0         1         2         3         4         5         6         7
   2212.4072924458269       0.85251684486045420      1368.7491062482197       -8.6353629983775451E-0020.47571300813309153      3084.9973674698904       -5.5590296671079252E-002   7.3994136421592915E-002   6.7846637375845634E-002   6450.1309586170764       0.20272393923520168      -0.91671031953902560       -4.7145102147721323E-002-5.8145625982189628E-002   5454.2857142857147      6.8571428571428575E-002   0.0000000000000000       -8.5714285714285715E-002 -0.20571428571428571      5.1428571428571428E-002   5468.7500000000000       0.25000000000000000      0.0000000000000000       0.18750000000000000       0.25000000000000000       0.18750000000000000      -0.18750000000000000      8000.0000000000000   
X=
-9.2481364876117276E-003   2.1208594097849032E-002-2.3487786720192489E-002-1.8023544606899312E-003   1.9928814217908663E-002-4.6947614371289264E-003   3.5496623117919580E-003

查看完整版本: 使用lapack的dsptrf输出上下三角矩阵计算结果不同问题