[Octave-bug-tracker] [bug #35779] Incorrect results of logm

[Marco Caliari]

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

Hi all.

For the case [0,-i;i,0] the fix is easy: just replace

if (real_eig)
  s = real (s);
endif

with

if (real_eig && isreal(A))
  s = real (s);
endif

in logm.m.
For the case [0,-1;1,0] the problem is in zrsf2csf.f (and in crsf2csf.f):
around line 48 there is a off-by-one error. You should replace 

c(j+1) = 1

with

if (j+1 < n) then
  c(j+1) = 1
end if


Now I'm trying to build from hg sources in order to make an hg patch.

Marco

    _______________________________________________________

Reply to this item at:

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

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