cholinv/spcholinv return incorrect results

[Thomas Weber]

On Fri, Dec 12, 2008 at 12:08:25PM -0600, Gabriele Pannocchia wrote:
> Hi,
> 
> the following code shows for me a bug in cholinv when the argument is a
> sparse matrix.
> 
> A = [2,0.2;0.2,1];
> issymmetric(A)
> min(eig(A))
> Ainv = inv(A);
> Ainv1 = cholinv(A);
> Ainv2 = inv(sparse(A));
> Ainv3 = cholinv(sparse(A));
> Ainv4 = spcholinv(sparse(A));
> err1 = norm(Ainv-Ainv1)
> err2 = norm(Ainv-Ainv2)
> err3 = norm(Ainv-Ainv3)
> err4 = norm(Ainv-Ainv4)
> 
> I am using:
> 
> ----------------------------------------------------------------------
> GNU Octave Version 3.0.1
> GNU Octave License: GNU General Public License
> Operating System: Linux 2.6.24-1-686 #1 SMP Thu May 8 02:16:39 UTC 2008
> i686
> ----------------------------------------------------------------------

Oh, Debian Lenny :)

> Can anyone reproduce this?
Yes, happens with 3.0.3 as well. No idea what's going wrong there. To
the list: this version of Octave use suitesparse 3.1.0.

	Thomas