Fortran Coder

标题: 使用lapack的dsptrf输出上下三角矩阵计算结果不同问题 [打印本页]

作者: nong 时间: 5 天前
标题: 使用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-002  0.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

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)