2008/5/12 Pierre GM <[EMAIL PROTECTED]>: > I fixed the power function in numpy.ma following Anne's suggestion: compute > first, mask the problems afterwards. It's a quick and dirty fix that crashes > if the user has set its error system to raise an exception on invalid > (np.seterr(invalid='raise')), but it works otherwise and keeps subclasses > (such as TimeSeries). > I will have to modify the .__pow__ method so that ma.power is called: right > now, a**b calls ndarray(a).__pow__(b), which may yield NaNs and Infs. > What should I do with oldnumeric.ma.power ? Try to fix it the same way, or > leave the bug ? I'm not that enthusiastic to have to debug the old package, > but if it's part of the job...
We should leave the oldnumeric.ma package alone. Even if its `power` is broken, some packages may depend on it. We'll provide it in 1.2 for backward compatibility, and get rid of it after 1.3. Regards Stéfan _______________________________________________ Numpy-discussion mailing list [email protected] http://projects.scipy.org/mailman/listinfo/numpy-discussion
