主题:Fortran求复数矩阵行列式的det函数
通过IMSl库里的det函数求复数矩阵的行列式总是得到真实值(已经通过matlab验证好几个复数矩阵的行列式)的共轭。以下面这个简单的例子为例,det函数求的值为2i,但是真实值是-2i(手算都可以了)。哪位高手比较了解库函数det,帮我解决一下这个问题:
program main
use imsl
implicit none
integer n
complex :: u(2,2),udet
u(1,1)=(1.,0.)
u(2,1)=(1.,2.)
u(1,2)=(2.,0.)
u(2,2)=(2.,2.)
udet=det(u)
write(*,*)udet
end
program main
use imsl
implicit none
integer n
complex :: u(2,2),udet
u(1,1)=(1.,0.)
u(2,1)=(1.,2.)
u(1,2)=(2.,0.)
u(2,2)=(2.,2.)
udet=det(u)
write(*,*)udet
end