[Numpy-discussion] Problem in LinearAlgebra? A solution?

Konrad Hinsen hinsen at cnrs-orleans.fr
Tue Nov 4 02:40:06 EST 2003


On Tuesday 04 November 2003 12:19, Nadav Horesh wrote:
> The condition of the matrix is about 1.0E7 and its dimensions are
> 10000x36: This is not a stable linear system, at least not for a simple

You can cut down the number of data points to much less (I use 400) and the 
problem persists.

>  But the solution to the polynomial fit turns to be much simpler:
> In the "fitPolynomial" function the 5th and the 4th lines before the end
> are commented. These lines uses the "generalized_inverse" procedure to
> solve the set of equations. just uncomment these lines and comment the

That is exactly what I recommended Rob to do as well.

Konrad.
-- 
-------------------------------------------------------------------------------
Konrad Hinsen                            | E-Mail: hinsen at cnrs-orleans.fr
Centre de Biophysique Moleculaire (CNRS) | Tel.: +33-2.38.25.56.24
Rue Charles Sadron                       | Fax:  +33-2.38.63.15.17
45071 Orleans Cedex 2                    | Deutsch/Esperanto/English/
France                                   | Nederlands/Francais
-------------------------------------------------------------------------------





More information about the NumPy-Discussion mailing list