public void SolveLU(Matrix B, Matrix X)
{
if (_LU != null)
{
if (B.RowCount == _N & X.RowCount == _N & B.ColumnCount == X.ColumnCount)
{
for (int j = 0; j < B.ColumnCount; j++)
{
Vector _X = new Vector(_N);
SolveLU(B.Column(j), _X);
for (int i = 0; i < _N; i++)
{
X[i, j] = _X[i];
}
}
}
else
{
throw new Exception("wrong sizes");
}
}
}
