Re: [Numpy-discussion] rcond in polyfit

2006-10-15 Thread Charles R Harris
On 10/14/06, A. M. Archibald <[EMAIL PROTECTED]> wrote: On 14/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote: I don't get the impression that the warnings module is much tested; Ihad similar headaches.Turns out to be a rather simple bug (feature?) in warnings, where it shortcircuits the filters

Re: [Numpy-discussion] rcond in polyfit

2006-10-14 Thread Charles R Harris
On 10/14/06, Tim Hochberg <[EMAIL PROTECTED]> wrote: Charles R Harris wrote:>>> On 10/14/06, *A. M. Archibald* <[EMAIL PROTECTED]> [EMAIL PROTECTED] >> wrote:[SNIP]>>> Hmmm, I wonder if we have a dictionary of precisions indexed by dtype> somewhere?Here's some code I stole from somewhere for comput

Re: [Numpy-discussion] rcond in polyfit

2006-10-14 Thread Tim Hochberg
Charles R Harris wrote: > > > On 10/14/06, *A. M. Archibald* <[EMAIL PROTECTED] > > wrote: [SNIP] > > > Hmmm, I wonder if we have a dictionary of precisions indexed by dtype > somewhere? Here's some code I stole from somewhere for computing EPS. It would easy enough to

Re: [Numpy-discussion] rcond in polyfit

2006-10-14 Thread Charles R Harris
On 10/14/06, A. M. Archibald <[EMAIL PROTECTED]> wrote: On 14/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote:>>> On 10/13/06, A. M. Archibald <[EMAIL PROTECTED] > wrote:> > On 13/10/06, Tim Hochberg <[EMAIL PROTECTED]> wrote:> > > Charles R Harris wrote:> Numerical Recipes (http://www.nrbook.co

Re: [Numpy-discussion] rcond in polyfit

2006-10-14 Thread A. M. Archibald
On 14/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > On 10/13/06, A. M. Archibald <[EMAIL PROTECTED]> wrote: > > On 13/10/06, Tim Hochberg <[EMAIL PROTECTED]> wrote: > > > Charles R Harris wrote: > > > > > > On the other hand if error handling is set to 'raise', then a > > > FloatingPoin

Re: [Numpy-discussion] rcond in polyfit

2006-10-14 Thread Charles R Harris
On 10/13/06, A. M. Archibald <[EMAIL PROTECTED]> wrote: On 13/10/06, Tim Hochberg <[EMAIL PROTECTED]> wrote:> Charles R Harris wrote: > On the other hand if error handling is set to 'raise', then a> FloatingPointError is issued. Is a FloatingPointWarning in order to> mirror the FloatingPointError?

Re: [Numpy-discussion] rcond in polyfit

2006-10-13 Thread A. M. Archibald
On 13/10/06, Tim Hochberg <[EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > That sounds good, but how to do it? Should I raise an exception? > Use the warnings framework: > > >>> import warnings > >>> warnings.warn("condition number is BAD") > __main__:1: UserWarning: conditio

Re: [Numpy-discussion] rcond in polyfit

2006-10-13 Thread Tim Hochberg
Charles R Harris wrote: > > > On 10/13/06, *A. M. Archibald* <[EMAIL PROTECTED] > > wrote: > > On 12/10/06, Charles R Harris <[EMAIL PROTECTED] > > wrote: > > Hi all, > > > > I note that polyfit looks like it should work for s

Re: [Numpy-discussion] rcond in polyfit

2006-10-13 Thread A. M. Archibald
On 13/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > > On 10/13/06, A. M. Archibald <[EMAIL PROTECTED]> wrote: > > On 12/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > Hi all, > > > > > > I note that polyfit looks like it should work for single and double, > real > > > and compl

Re: [Numpy-discussion] rcond in polyfit

2006-10-13 Thread Charles R Harris
On 10/13/06, A. M. Archibald <[EMAIL PROTECTED]> wrote: On 12/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote:> Hi all,>> I note that polyfit looks like it should work for single and double, real > and complex, polynomials. On the otherhand, the default rcond needs to> depend on the underlying pr

Re: [Numpy-discussion] rcond in polyfit

2006-10-13 Thread A. M. Archibald
On 12/10/06, Charles R Harris <[EMAIL PROTECTED]> wrote: > Hi all, > > I note that polyfit looks like it should work for single and double, real > and complex, polynomials. On the otherhand, the default rcond needs to > depend on the underlying precision. On the other, other hand, all the svd > com