This forum has been archived. All content is frozen. Please use KDE Discuss instead.

Symmetric Band Matrix + Lapack

Tags: None
(comma "," separated)
davidhigh
Registered Member
Posts
18
Karma
0

Symmetric Band Matrix + Lapack

Fri Nov 08, 2013 8:40 am
Hi,

I've seen there is an old Band Matrix implementation, and also read some threads here. It seems that only the storage scheme is supported by now, and I read that it is possible to apply Lapack functionality to the stored matrix, at least for the non-symmetric case.
Some questions:

(i) is this Band Matrix storage scheme still safely usable? Are there extensions planned?
(ii) Is there also a symmetric Band Matrix compatible with Lapack's symmetric storage scheme? Or can it be realized by SymmetricView -- or similar things -- in a way that it is also possible to use the symmetric solvers of Lapack?

Thanks in advance,
David
User avatar
ggael
Moderator
Posts
3447
Karma
19
OS

Re: Symmetric Band Matrix + Lapack

Sun Nov 10, 2013 8:51 am
(I) yes, and yes.
(II) BandMatrix is naturally compatible with symmetric band storage, you simply has to set the number of super diagonals to 0 or get the address of the (0,0) element as a starting point
davidhigh
Registered Member
Posts
18
Karma
0

Re: Symmetric Band Matrix + Lapack

Sun Nov 10, 2013 7:33 pm
Thank you,
David
davidhigh
Registered Member
Posts
18
Karma
0

Re: Symmetric Band Matrix + Lapack

Sun Nov 10, 2013 8:34 pm
could you provide a short example of how to initialize the band matrix? If I include the Skyline module from the "unsupported" folder, I get a bunch full of errors ...
Thanks in advance,
David
davidhigh
Registered Member
Posts
18
Karma
0

Re: Symmetric Band Matrix + Lapack

Sun Nov 10, 2013 8:39 pm
ah, ok,

int N=1000;
Eigen::internal::BandMatrix(double,N,2,2) mat;

seems to work. I forgot the "::internal". But how do I fill it?

mat(0,0)=1.0;

doesn't work...

Best,
David
User avatar
ggael
Moderator
Posts
3447
Karma
19
OS

Re: Symmetric Band Matrix + Lapack

Sun Nov 10, 2013 9:19 pm
you have to fill it through its diagonals.


Bookmarks



Who is online

Registered users: Bing [Bot], Google [Bot], kesang, Sogou [Bot], Yahoo [Bot]