Hello,

I'm working on a research project involving numerical simulations and I'm interested in using Eigen due to the excellent code readability it provides. Since my problem involves finite differences for PDEs I will naturally be dealing with banded matrices. I know that Eigen has a banded matrix class but, if I understood correctly, operations and solvers can currently only be used via BLAS/LAPACK. Is that correct? If so, are there plans to implement these directly in Eigen in the near future? If not, what is the procedure to use BLAS/LAPACK routines?

Thanks in advance,
(and keep up the good work!)
Paulo Alcino

Hi Paulo,

yes there is plan to provide all kind of products and solvers on BandMatrix, though I cannot say when they will be available. In the meantime, you can call BLAS/LAPACK routines using m.coeffs().data(); to get a pointer to the first elements, m.coeffs().outerStride(); for the leading dimension, m.supers()/m.subs() for the number of bands, etc.