MathGroup Archive 1998

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

Search the Archive

Re: Numerical Determinants

  • To: mathgroup at smc.vnet.net
  • Subject: [mg13001] Re: Numerical Determinants
  • From: hello at there.com (Rod Pinna)
  • Date: Tue, 30 Jun 1998 00:26:01 -0400
  • Organization: UWA
  • References: <6n4o6m$nin@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Thanks for the reply. I eventually realised that I should be  using the
eigenvalue function, however, it turns out that for the  matricies I'm
examining, B is non-invertible. To get around this,  I've been using
the following:

            1
--------------------------
Eigenvalues[Inverse[A]].B]

This seems to be giving me the expected answers for A.x=p B.x, However,
I'm not convinced by the mathematics of it...

Thanks again though.

In article <6n4o6m$nin at smc.vnet.net>, Hugh Walker  <hwalker at gvtc.com>
wrote:

>Yes, I have a suggestion.  Consider the eigen-problem -a.x = p b.x.
>Non-trivial eigenvectors exit provided Det[a + b p] = 0, which is your
>situation.  Assuming Inverse[b] exists, the eigen values of m.x = p x
>where m = -Inverse[b].a
>
>Eigenvalues[-Inverse[b].a]
>
>will coincide with your p's.
>
>I tested cases in which a and b are each random 10x10 matrices. The
>results were the same, and timings for the two approaches were
>
>Det = 9.0 sec
>eigen = 0.0 sec
>
>I hope this helps.
>
>
>Hugh Walker
>Professsor Emeritus
>University of Houston
>
>Gnarly Oaks
>
>

Rod Pinna
(rpinnaX at XcivilX.uwa.edu.au  Remove the X for email)



  • Prev by Date: Re: (a-b)/(c-d)=!=(b-a)/(d-c) ???
  • Next by Date: Re: Converting from StandarForm to Text
  • Previous by thread: Re: Numerical Determinants
  • Next by thread: Find Max of "Concave" List