The program at the bottom of this message returns the following runtime warning:
python test.py
test.py:5: RuntimeWarning: invalid value encountered in divide
return np.where(x==0., 1., np.sin(x)/x)
The function works correctly returning
x = np.array([ 0., 1., 2., 3., 4., 5., 6.,
Jaime,
If you are going to work on this, you should also take a look at the recent
thread
http://mail.scipy.org/pipermail/numpy-discussion/2013-February/065649.html,
which is about the weighting function, which is in a confused state in the
current version of polyfit. By the way, Numerical Recipe
it is clear. I
would be happy to contribute both to improving the documentation and
software.
David
If
On Wed, Feb 27, 2013 at 12:47 PM, Pauli Virtanen wrote:
> 27.02.2013 16:40, David Pine kirjoitti:
> [clip]
> > 2. I am sorry but I don't understand your response. The
eduction and Error Analysis for the
Physical Sciences" by Bevington, both standard works.
Dave
On Wed, Feb 27, 2013 at 9:03 AM, Charles R Harris wrote:
>
>
> On Wed, Feb 27, 2013 at 6:46 AM, David Pine wrote:
>
>> As of NumPy v1.7, numpy.polyfit includes an option for
he references given in the poyfit code make no mention of it. I think it
would be much better to return the standard definition of the covariance
matrix (see, for example, the book Numerical Recipes).
David Pine
___
NumPy-Discussion mailing list
NumPy
