[Octave-bug-tracker] [bug #36437] chol2inv can't handle sparse matrix

[Nir Krakauer]

URL:
  <http://savannah.gnu.org/bugs/?36437>

                 Summary: chol2inv can't handle sparse matrix
                 Project: GNU Octave
            Submitted by: nir_krakauer
            Submitted on: Wed 09 May 2012 01:02:15 AM GMT
                Category: None
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: None
                  Status: None
             Assigned to: None
         Originator Name: Nir Krakauer
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 3.6.1
        Operating System: GNU/Linux

    _______________________________________________________

Details:


chol2inv seems not to work when the input is a sparse upper triangular matrix.
For example:

n = 10;
X = spdiags(ones(n, 1) * [1 2 1], [-1 0 1], n, n);
Y = chol(X);
Z = chol2inv(Y)
error: chol2inv: spchol2inv requires triangular matrix

On the other hand, 
Z = chol2inv(full(Y))
returns the expected result




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?36437>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/