public LU(double[][] A)
{
this.m_LU = LU.new_copy(A);
this.m_pivot = new int[A.Length];
LU.factor(this.m_LU, this.m_pivot);
}
public virtual double[] solve(double[] b)
{
double[] b1 = LU.new_copy(b);
LU.solve(this.m_LU, this.m_pivot, b1);
return(b1);
}