[Octave-bug-tracker] [bug #35787] Diagonal matrices don't broadcast

[Lachlan Andrew]

Follow-up Comment #4, bug #35787 (project octave):

This is still an issue in version 4.0.0.

It is easy to work around by placing the diagonal matrix in [...]:

> [eye(2)] + [1 0]
ans =

   2   0
   1   1

However, that is not done in some library routines (I was bitten by mvnrnd.m)
and so Matlab-compliant user code that does not rely on automatic broadcasting
can be affected.  The bug has been open for three years; is it time to rewrite
all library routines to work around it?


Another surprising twist is that a slice of a diagonal matrix is still
considered "diagonal".
> a = eye(2);
> b = a(1,:)
b =

Diagonal Matrix

   1   0
> b + [1; 0]
error: operator +: nonconformant arguments (op1 is 1x2, op2 is 2x1)




    _______________________________________________________

Reply to this item at:

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

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