Should numpy.sqrt(-1) return 1j rather than nan?

[Travis Oliphant]

Tim Hochberg wrote:

>With python 2.5 out now, perhaps it's time to come up with a with 
>statement context manager. Something like:
>
>    from __future__ import with_statement
>    import numpy
>
>    class errstate(object):
>        def __init__(self, **kwargs):
>            self.kwargs = kwargs
>        def __enter__(self):
>            self.oldstate = numpy.seterr(**self.kwargs)
>        def __exit__(self, *exc_info):
>            numpy.seterr(**self.oldstate)
>           
>    a = numpy.arange(10)
>    a/a # ignores divide by zero
>    with errstate(divide='raise'):
>        a/a # raise exception on divide by zer
>    # Would ignore divide by zero again if we got here.
>
>-tim
>
>  
>

This looks great.  I think most people aren't aware of the with 
statement and what it can do (I'm only aware because of your posts, for 
example). 

So, what needs to be added to your example in order to just add it to 
numpy?

-Travis


-------------------------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642