[Octave-bug-tracker] [bug #36099] diag(1,N,M) produces NaNs

[Twan van Laarhoven]

Follow-up Comment #1, bug #36099 (project octave):

Probable cause seems to be:
+verbose+
return Array<ST> (dim_vector (1, 1), scalar).diag (k);
-verbose-
The diag() function of Array produces a dense matrix instead of DiagMatrix
object. I think an MArray<ST> should be used instead.

That would also fix the issue that a dense matrix is produced, while it
shouldn't be:
+verbose+
octave:1> length(diag([1,2],1e5,1e5))
ans =  100000
octave:2> length(diag([1],1e5,1e5))
error: memory exhausted or requested size too large for range of Octave's
index type -- trying to return to prompt
-verbose-

This doesn't explain why Array::diag or Array::resize messes things up,
though.

    _______________________________________________________

Reply to this item at:

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

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