MathGroup Archive 1999

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Eigensystems


Tony Harker <ahh at baker.phys.ucl.ac.uk> wrote

>  Can anybody recommend an efficient method for solving the modified
>eigensystem
>       A x = lambda B x
>numerically in Mathematica?
>

If Binv,the inverse of B exists, left-multiply through resulting in
(Binv.A).x = lambda x. This is then soluble by eigensystems. If B is not
invertable, then do the analogous thing with Ainv. In this case the
eigenvalue is 1/lambdq. If neither inverse exists, you have a bigger
problem.

    ==

Hugh Walker
Gnarly Oaks




  • Prev by Date: Re: single precision complex bug in g77 - was Testing g77 with LAPACK 3.0
  • Next by Date: Hypergeometric2F1[a,b,c,1] in v4
  • Previous by thread: Eigensystems
  • Next by thread: Re: Eigensystems