Re: sqrt function, NaN and complex number

[Dmitri A. Sergatskov]

I can reproduce that with Fedora Core 3. I wonder if this is problem with
glibc (glibc-2.3.3-74).

Dmitri.

Per Persson wrote:
> I can't trigger that bug on OS X 10.3.6 with octave 2.1.64, see below. 
> Problem with Mandrake?
> /Per
>
>  > M
> M =
>
>   1 + 1i  1 + 1i
>   1 + 1i  1 + 1i
>
>  > sqrt(M/0)
> warning: division by zero
> ans =
>
>           NaN +         NaNi          NaN +         NaNi
>           NaN +         NaNi          NaN +         NaNi
>
> On Dec 10, 2004, at 19:00, Samir Sharshar wrote:
>
> > Hi,
> >
> > I'm currently use an octave version 2.1.57 on a Mac OSX 10.3.6 and a 
> > 2.1.64 compiled from sources on a mandrake 10.1.
> >
> > I've a couple of functions which return a matrix M of complex numbers 
> > with NaN in it. Display the content of M give NaN + (or -) NaNi in 
> > case of NaN. I need to calculate the square root of M. On 2.1.57, sqrt 
> > function works perfectly and returns NaN + (or -) NaNi in case of NaN. 
> > But on 2.1.64 octave hangs without error or returning to prompt.
> >
> > Any idea ?
> >
> > Samir.
> >
> >
> >
> > -------------------------------------------------------------
> > Octave is freely available under the terms of the GNU GPL.
> >
> > Octave's home on the web:  http://www.octave.org
> > How to fund new projects:  http://www.octave.org/funding.html
> > Subscription information:  http://www.octave.org/archive.html
> > -------------------------------------------------------------
>
>
>
> -------------------------------------------------------------
> Octave is freely available under the terms of the GNU GPL.
>
> Octave's home on the web:  http://www.octave.org
> How to fund new projects:  http://www.octave.org/funding.html
> Subscription information:  http://www.octave.org/archive.html
> -------------------------------------------------------------



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------