Re: [Numpy-discussion] rcond in polyfit

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

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

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

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

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

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

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

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

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

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

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

[Numpy-discussion] rcond in polyfit

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 computations are done with dgelsd or zgelsd, i.e., double precision. Eve