Hi,
Congrats on the new release.
I found an error in PDL::LinearAlgebra.
pdl> use PDL::LinearAlgebra
pdl> $a=(1+sequence(2,2))*i()
pdl> p $a
[
[ i 2i]
[3i 4i]
]
pdl> p $a->mdet
0 # instead of 2
I believe the problem lies in the code that distinguishes which lapack
LU factorization routine to use. There is a complex routine cgetrf,
called when the object is a PDL::Complex, and a real routine geterf,
called when the object is a PDL. However, native complex types are of
the PDL type. I made a pull request. Nevertheless, I believe that this
kind of errors might appear elsewhere.
Regards,
Luis
--
o
W. Luis Mochán, | tel:(52)(777)329-1734 /<(*)
Instituto de Ciencias Físicas, UNAM | fax:(52)(777)317-5388 `>/ /\
Av. Universidad s/n CP 62210 | (*)/\/ \
Cuernavaca, Morelos, México | moc...@fis.unam.mx /\_/\__/
GPG: 791EB9EB, C949 3F81 6D9B 1191 9A16 C2DF 5F0A C52B 791E B9EB
_______________________________________________
pdl-general mailing list
pdl-general@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pdl-general
