您所在位置：主题：[讨论]求助：CVF下imsl计算行列式
faye100
[专家分：0] 发布于 2010-08-11 15:15:00
请教：
我使用imsl库的det计算行列式，发现矩阵值很大时，算的有问题诶，如下：
[u]real :: a(2,2),value
a=2e10
value=det(a)
print *,'该方阵的行列式值为',value[/u]
该矩阵内各个位置值相同，理论上其行列式应该是0
但，
若取a=2e10，则输出的值为1.813..e13
若取a=2e0，则输出的值为-1.656..e-8
。。。。这让我很疑惑啊，请教各位大牛，这是什么原因？？
PS：我是在做有限元刚度矩阵计算的，理论上都该为0，但。。。甚是疑惑
我使用imsl库的det计算行列式，发现矩阵值很大时，算的有问题诶，如下：
[u]real :: a(2,2),value
a=2e10
value=det(a)
print *,'该方阵的行列式值为',value[/u]
该矩阵内各个位置值相同，理论上其行列式应该是0
但，
若取a=2e10，则输出的值为1.813..e13
若取a=2e0，则输出的值为-1.656..e-8
。。。。这让我很疑惑啊，请教各位大牛，这是什么原因？？
PS：我是在做有限元刚度矩阵计算的，理论上都该为0，但。。。甚是疑惑
