Mon May 13, 2013 1:41 pm
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 ?
Tue May 14, 2013 3:19 pm
We just add a support for that in the devel branch
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