Hi all,

I am using Eigen 3.2.0-beta1.
I want to solve
L X = B
U Y = C
for [L,U] = SparseLU(A), independently.

I found SparseLU::matrixL() so I can solve L X = B,
however, I could not find SparseLU::matrixU().
Is there any method to solve U Y = C ?

Thanks.

We just add a support for that in the devel branch
https://bitbucket.org/eigen/eigen/commits/6b984a1fb368/
Changeset: 6b984a1fb368
User: dnuentsa
Date: 2013-05-14 17:15:23
Summary: Add support in SparseLU to solve with L and U factors independently

Note that you will need to apply the rows and columns permutations as well.
You can find the typical steps in the _solve method of SparseLU.h