Re: Inverse Matrix Function appears a bit wonky

[Geraint Paul Bevan]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Dmitri A. Sergatskov wrote:

|
| How come I do not get a warning?
| (Octave 2.1.60 / ATLAS 2.7.8)
|  octave:35> m = [1 2 3 4; 2 2 3 3; 1 2 1 2; 3 2 2 1];
| octave:36> det (m)
| ans =  -1.3323e-15

Perhaps because the answer is greater than your machine's precision?
FWIW, eps is 2.22044604925031E-16 on mine.

I am running 2.1.60, Debian package octave2.1 2.1.60-1 which is linked
against:

address@hidden:~$ ldd /usr/bin/octave
~                liboctinterp.so =>
/usr/lib/octave-2.1.60/liboctinterp.so (0x40018000)
~        liboctave.so => /usr/lib/octave-2.1.60/liboctave.so (0x407dd000)
~        libcruft.so => /usr/lib/octave-2.1.60/libcruft.so (0x40e0e000)
~        liblapack.so.3 => /usr/lib/atlas/liblapack.so.3 (0x40e7c000)
~        libblas.so.3 => /usr/lib/atlas/libblas.so.3 (0x414c8000)
~        libfftw3.so.3 => /usr/lib/libfftw3.so.3 (0x41809000)
~        libreadline.so.4 => /lib/libreadline.so.4 (0x418c1000)
~        libncurses.so.5 => /lib/libncurses.so.5 (0x418ed000)
~        libdl.so.2 => /lib/tls/libdl.so.2 (0x4192d000)
~        libhdf5-1.6.2.so.0 => /usr/lib/libhdf5-1.6.2.so.0 (0x41930000)
~        libz.so.1 => /usr/lib/libz.so.1 (0x41a3a000)
~        libg2c.so.0 => /usr/lib/libg2c.so.0 (0x41a4b000)
~        libm.so.6 => /lib/tls/libm.so.6 (0x41a69000)
~        libgcc_s.so.1 => /lib/libgcc_s.so.1 (0x41a8c000)
~        libstdc++.so.5 => /usr/lib/libstdc++.so.5 (0x41a96000)
~        libc.so.6 => /lib/tls/libc.so.6 (0x41b4f000)
~        /lib/ld-linux.so.2 => /lib/ld-linux.so.2 (0x40000000)
~        libpthread.so.0 => /lib/tls/libpthread.so.0 (0x41c8a000)

- --
Geraint Bevan
http://homepage.ntlworld.com/geraint.bevan

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.4 (GNU/Linux)

iEYEARECAAYFAkGJT3QACgkQcXV3N50QmNPsXgCaAj/YPPRW48RI2JW9u5iGDofV
XNYAoIFbBAJZ6c6ZK+SbsBLkgNn87Rw3
=fGV5
-----END PGP SIGNATURE-----



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------