Re: Sparse + complex arithmetic inconsistency

[Rafael Laboissiere]

* David Bateman <address@hidden> [2007-03-30 14:13]:

> Michael Goffioul wrote:
> >
> > I already reported this kind of problem some months ago. This is
> > due to the complex arthmetic implementation and it also appears
> > under MSVC. See
> > http://www.nabble.com/Sparse-matrix-problem-tf2844890.html#a8149128
> >
> > You can even illustrate the problem without the use of sparse
> > matrices, like:
> >
> > [1+i 1+i] ./ [0 i]
> >
> Following the bug report trail from this thread leads to
> 
> [snip]
> 
> No resolution from the GCC people is proposed, and all of the above bugs
> are open. Hopefully, one day the gcc/g++ complex classes will be
> mathematically correct :-(

Does the gcc/g++ bug have something to do with the fact that the two
following lines produce different results in Octave (but shouldn't):

    octave2.9:1> ([1+i 1+i] ./ [0 i]) (1)
    ans = NaN - NaNi
    octave2.9:2> (1+i) / 0
    warning: division by zero
    ans = Inf + Infi

It looks like a problem in Octave.

-- 
Rafael