,
Ed
From: Luis Mochan<mailto:moc...@icf.unam.mx>
Sent: 14 February 2023 19:05
To: Ed .<mailto:ej...@hotmail.com>
Cc: perldl<mailto:pdl-gene...@lists.sourceforge.net>;
perldl<mailto:pdl-devel@lists.sourceforge.net>
Subject: Re: [Pdl-devel] inconsistent determinant
Hi E
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
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
