I'm getting ready to start working on a C/C++ project that will be building and solving large tri-diagonal, block tri-diagonal and triangular matrices. I know there are a lot of libraries available ...

We introduce a new sparse sliced inverse regression estimator called Cholesky matrix penalization, and its adaptive version, for achieving sparsity when estimating the dimensions of a central subspace ...