#1186 DynMatrixHelper?
Metadata
Property | Value |
---|---|
Reported by | Arne Morten Kvarving (arne.morten.kvarving@sintef.no) |
Reported at | Oct 3, 2012 16:09 |
Type | Discussion |
Version | 2.2 |
Operating System | Unspecified / All |
Last edited by | Markus Blatt (markus@dr-blatt.de) |
Last edited at | Nov 21, 2012 16:20 |
Closed by | Markus Blatt (markus@dr-blatt.de) |
Closed at | Nov 21, 2012 16:20 |
Closed in version | Unknown |
Resolution | Implemented |
Comment | in rev, 7061 |
Description
Hi,
I have had the need for eigenvalues of matrices, and the available utils are a bit lacking.
i have written some code structured on fmatrixev.[hh/cpp] and now i would like to contribute this back to the project.
since i don't know how you want the shed painted, i ask how and if before i waste time painting it up.
sources attached for reference. there is two sets of files; dynmatrixev.hh/cpp does complex eigenvalues for non-symmetric dynamic matrices, while fmatrixev_ext.hh/cc implements the same for fmatrices. of course lapack is used to do the heavy lifting.
please advice (whether it's wanted, and if so how to do it).
if this goes in, i can also contribute (basic but fully working) code to calculate gauss-lobatto-legendre quadratures if this is of interest (this is what i needed it for). currently dune needs alglib to do this, and alglib is marked as deprecated.