请教：FORTRAN的imsl库 计算行列式

faye100

请教：
我使用imsl库的det计算行列式，发现矩阵值很大时，算的有问题诶，如下：

real :: a(2,2),value
  a=2e10
  value=det(a)
print *,'该方阵的行列式值为',value

该矩阵内各个位置值相同，理论上其行列式应该是0
但，
若取a=2e10，则输出的值为1.813..e13
若取a=2e0，则输出的值为-1.656..e-8

。。。。这让我很疑惑啊，请教各位大牛，这是什么原因？？

PS：我是在做有限元刚度矩阵计算的，理论上都该为0，但。。。甚是疑惑

pasuka

请教：
我使用imsl库的det计算行列式，发现矩阵值很大时，算的有问题诶，如下：

real :: a(2,2),value
  a=2e10
  value=det(a)
print *,'该方阵的行列式值为',value

该矩阵内各个位置值相同，理论上其行列 ...
faye100 发表于 2010-8-12 09:19

用双精度也是同样的结果？！

faye100

2# pasuka 恩，是的
program  d
use imsl
implicit none

integer :: n,i,j
real(8):: a(2,2)
real(8) det_a

a=1.0d10

det_a=det(a)
print *,'该方阵的行列式值为',det_a
end

计算的值为26973.983057096722，晕