|
From: | Sebastien Loisel |
Subject: | Re: Good eigenvalue/eigenvector test problems for eigs? |
Date: | Fri, 22 Sep 2006 11:38:59 -0400 |
Does any one have a simple to create set of matrices of variable order
for symmetric, non-symmetric and complex problems for me to test my eigs
code against. Basically, it seems that the easy to create sparse
matrices that I can think of are all ill-conditioned and result in
incorrect results in many case. What I currently have is
## Symmetric real PD
A = spdiags([ones(n,1),4*ones(n,1),ones(n,1)],[-2,0,2],n,n);
## Non-symmetric
A = spdiags([ones(n,1),4*ones(n,1),-ones(n,1)],[-2,0,2],n,n);
## Complex
A = spdiags([ones(n,1),4*ones(n,1),1i*ones(n,1)],[-2,0,2],n,n);
which as I say are ill-conditioned. So I'm looking for something as easy
to express as the above.. Apart from this issue, I appear to have a
relative complete eigs function now working (more info later)...
D.
--
David Bateman address@hidden
Motorola Labs - Paris +33 1 69 35 48 04 (Ph)
Parc Les Algorithmes, Commune de St Aubin +33 6 72 01 06 33 (Mob)
91193 Gif-Sur-Yvette FRANCE +33 1 69 35 77 01 (Fax)
The information contained in this communication has been classified as:
[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary
[Prev in Thread] | Current Thread | [Next in Thread] |