[SciPy-user] Inconsistent standard deviation and variance implementation in scipy vs. scipy.stats

Robert Kern robert.kern@gmail....
Mon Oct 6 21:40:38 CDT 2008

On Mon, Oct 6, 2008 at 21:30, Anne Archibald <peridot.faceted@gmail.com> wrote:
> The case I was concerned about was
> import scipy as sp
> x = sp.cos(2*sp.arccos(y))
> If this is changed to
> import numpy as np
> x = np.cos(2*np.arccos(y))
> it suddenly stops working for values y>1. To keep the same behaviour
> it needs to be
> import scipy as sp
> import numpy as np
> x = np.cos(2*sp.arccos(y))
> This is perhaps all right, but it does mean that users need to pay attention.

Well, if we remove some names from scipy/__init__.py, we should remove
them all. The actual definitions of those extended-domain functions
are actually in numpy.lib.scimath, not scipy. We can make a convenient
module inside numpy that basically does this:

  from numpy import *
  from numpy.lib.scimath import *

Then the transition for people using "import scipy" becomes quite easy.

Robert Kern

"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
  -- Umberto Eco

More information about the SciPy-user mailing list