[eigen] banded matrices in Eigen

[Deaglan O Halligan <dhalliga@xxxxxxxxxx>]

Hi,

I'm a grad student working on Trilinos, a collection of algorithms for large-scale scientific computation from Sandia Labs, and a related fluid dynamics application code called Tramonto. We're implementing algorithms that use David Bailey's qdlib data types, and this work led us to look at Eigen for doing a few necessary dense computations. I'm currently working on adding some Eigen code to Trilinos as a third-party-library (TPL).

So far I've added Eigen's householderQR and PartialPivLU to Trilinos. I'm also interested in adding the symmetric eigensolver and a banded LU factorization. I tried using Eigen's unsupported Skyline module but ran into problems, and I'm not sure what the status of that module is. I then realized that there is a BandMatrix class that uses the banded storage format, and there's been some mention on forum threads about more functionality being added for the BandMatrix class. It seems like the BandMatrix class is where future development will go on and I should ignore the Skyline module? I'm very interested in knowing what the status of BandMatrix is and how I can help with further development (in particular implementing LU and products) if that's needed.

Thanks,
Deaglan Halligan