🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: sci.math.num-analysis
Article #288

Re: Eigenvalue,Eigenvector (part II)

#288
From: Ferhat F. Hatay
Date: Tue, 14 Apr 1992 14:42
73 lines
2802 bytes
I am having difficulties to reach Eurepan dudes (dudesses!) by mail, 
somewhat they keep rejecting my mail. So, please forgive my overuse
of network to communicate the following -deep thoughts- to them. 
I don't claim any liability for my suggestions, anyone is welcomed 
to find any holes in my reasoning. feel free to poke my failing eyes. 
 
In article <17229@sophia.inria.fr> falzon@sophia.inria.fr writes:
>
> Does anybody have the general forms of eigenvalues and eigenvectors of an NxN
> pseudo-tridiagonal matrix of the form :
> 
> For example with N=4;
> 
>
>	[(1+1/a) -1/a       0         0 ]
>  A = 	[-1/a    (1+2/a)  -1/a        0 ]
>	[ 0       -1/a    (1+2/a)  -1/a ]
>	[ 0         0      -1/a  (1+1/a)] 
>
> Thanks in advance.
>
>
>      Fred.
>-----------------------------------------------------------------------------------------

Hmm, The problem adds to its lure when we consider it with 
some boundary conditions. It stops being circulant and
turns to a headache. Assuming that this is a solution
to a second order pde like Helmholtz or Poisson, let me 
suggest the following:

For dirichlet boundary conditions we can fold them 
onto (or into?) a system of order of N-2. The 
resulting -nonsingular- system (matrix) will again
have a circulant structure and a different right 
hand side. 

For neumann boundary conditions: I use my eigenvalue solver!

I still think there is an analytical solution for the problem
mentioned here. If we solve the eigenproblem for a matrix A:

    A = Q D Qt where D is diagonal matrix of eigenvalues.

    Q A Qt = D 

what would be the solution for

    B = A + S ? where S is a diagonal matrix. 

let's try:

   Q B Qt = Q A Qt + Q S Qt
    
   Q B Qt =  D + Q S Qt 

then what? eigenvectors of matrix B are Q same as the original 
matrix A and eigenvalues are given by D + S?

Do you see any holes in the reasoning above? It was rather
spontaneous so I should have probably checked myself.
  
cheers, 

Ferhat


-- 
Ferhat F. Hatay                        || Disclamer: |||||||||||||||||
Aerospace Engineering Sciences         || I speak for myself unless || 
University of Colorado, Boulder        || otherwise indicated.      ||
e-mail: hatay@spot.colorado.edu        |||||||||||||||||||||||||||||||
Technical Details
Message-ID: <1992Apr14.144201.2845@ucsu.Colorado.EDU>
Path: arc1-ca.pugleaf.net!sparky!uunet!boulder!ucsu!spot.Colorado.EDU!hatay
References: <17229@sophia.inria.fr>

Show Full Headers
Complete Article Headers
Newsgroups: sci.math.num-analysis
Path: arc1-ca.pugleaf.net!sparky!uunet!boulder!ucsu!spot.Colorado.EDU!hatay
From: hatay@spot.Colorado.EDU (Ferhat F. Hatay)
Subject: Re: Eigenvalue,Eigenvector (part II)
Message-ID: <1992Apr14.144201.2845@ucsu.Colorado.EDU>
Sender: news@ucsu.Colorado.EDU (USENET News System)
Nntp-Posting-Host: spot.colorado.edu
Organization: University of Colorado, Boulder
References: <17229@sophia.inria.fr>
Date: Tue, 14 Apr 1992 14:42:01 GMT
Lines: 73