Hi Luis,

I’ll take a look shortly (in the final stages of getting the multi-C to work). 
Can you confirm whether this worked according to your expectation in 2.057?

Best regards,
Ed

From: Luis Mochan<mailto:moc...@icf.unam.mx>
Sent: 25 September 2021 15:49
To: perldl<mailto:pdl-gene...@lists.sourceforge.net>; 
perldl<mailto:pdl-devel@lists.sourceforge.net>
Subject: [Pdl-devel] native complex and PDL::LinearAlgebra


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-devel mailing list
pdl-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pdl-devel

_______________________________________________
pdl-devel mailing list
pdl-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pdl-devel

Reply via email to