[Octave-bug-tracker] [bug #47415] out of memory with factorized sparse matrix using mldivide

[Carlo de Falco]

Follow-up Comment #3, bug #47415 (project octave):

Intrestingly enough, I tried separating the solve line into 
different steps and it works:


>> load sparse_matrix_bug_data.binary.gz 
>> [jacL, jacU, jacP, jacQ, jacR] = lu (jac);
>> delta = - jacQ * (jacU \ (jacL \ (jacP * (jacR \ res))));
error: out of memory or dimension too large for Octave's index type
>> a = jacR \ res;
>> b = jacP \ a;
>> c = jacL \ b;
>> d = jacU \ c;
>> e = jacQ * d;
>> delta = -e;


maybe the error lies in the parser somewhere?



    _______________________________________________________

Reply to this item at:

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

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