[SciPy-user] Scipy's orth is more reliable

[Nils Wagner]

Arnd Baecker schrieb:
> 
> just out of curiosity - what to you mean by "more reliable"?
> And to quote you
>   "A small example would be appreciated."
> (sorry, I could not resist ;-)
> 
> Arnd
> 
This is somehow lengthy to explain, but I will try.
I have computed a spectral projector for a defective matrix.A
where the Jordan curve \Gamma is a circle for simplicity. 
The circle encloses the double eigenvalue, which is defective.
Then I have orthonormalized the projector.
The associated invariant subspace displays the geometric multiplicity
of the eigenvalue.
If I apply scipy's orth functions, I am able to show the defective
nature of that eigenvalue.
However with MLab's orth function, the result is polluted by rounding
errors or whatever.


Nils

> On Fri, 19 Sep 2003, Nils Wagner wrote:
> 
> > Hi all,
> >
> > It seems to me, that scipy's orth module is more reliable than
> > the one, which comes with MLab2.py. .
> >
> > http://pdilib.sourceforge.net/MLab2.py
> >
> > Nils
> >
> > _______________________________________________
> > SciPy-user mailing list
> > SciPy-user at scipy.net
> > http://www.scipy.net/mailman/listinfo/scipy-user
> >
> 
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-user