[Numpy-discussion] Unexpected warnings in where

[T. C. <al...@an...>]

Hi,

I cannot understand why I do receive these warnings:

Warning: Encountered invalid numeric result(s)  in sqrt
Warning: Encountered invalid numeric result(s)  in divide

after

>>> newtimes =3D where(logical_and(RV<0,discr>0), (-RV-sqrt(discr))/VV, far=
_future)


RV, discr, RR, VV are NxN arrays of reals (Float64), while
far_future=3D1e20 (scalar).

1. sqrt(<0) should be impossible, since it's explicitly ruled out that
discr<0 in the where condition.

2. /0 should also be impossible, since VV=3D=3D0 is True only for the
diagonal...

I've done some testing and where(logical_and(RV<0, discr>0), discr, 666)
doesn't show any negative numbers, neither does
where(logical_and(RV<0,discr>0),VV,666) show any zeroes...


=09
Please.. what am I doing wrong here?


--=20
=C1lvaro Tejero Cantero

http://alqua.org    --	documentos libres
			free documents