Re: [Pdl-general] [Pdl-devel] inconsistent determinant

2023-02-15 Thread Ed .
Hi Luis, Any unexpected behaviour is worth asking about, so thank you (as so often)! In particular, anything that breaks an expectation of being able to chain operations in an OO style is bad. I’ve just pushed what I believe is a fix on GitHub, so please give that a go if you like. I’m holding

Re: [Pdl-general] [Pdl-devel] inconsistent determinant

2023-02-14 Thread Luis Mochan
Hi Ed, Thanks for looking at this. It was not too important, but it was surprising, as I was chaining pdl methods (...->det->approx(0) to detetct singular matrices) and I was getting somewhat strange error messages as det was sometimes an ndarray and other times a scalar. Regards, Luis On Sat, Feb

Re: [Pdl-general] [Pdl-devel] inconsistent determinant

2023-02-10 Thread Ed .
Hi Luis, The key (if you’ll forgive the pun) is the behaviour of the “det” function, which is (at this writing) defined in Basic/MatrixOps/matrixops.pd as: sub det { my($x) = shift; my($opt) = shift; $opt = {} unless defined($opt); my($lu,$perm,$par); if(exists ($opt->{lu}) and (ref $o