public static SafeMatrix SolveLinearSystem(SafeMatrix A, SafeMatrix X)
{
double[] a = A.GetElements();
double[] x = X.GetElements();
double[] y = new double[A.ColumnCount * X.ColumnCount];
SolveLinearSystem(a, A.RowCount, A.ColumnCount, x, X.RowCount, X.ColumnCount, y);
SafeMatrix rv = new SafeMatrix(X.RowCount, X.ColumnCount);
rv.SetElements(y);
return rv;
}
/// <summary>
/// Warning: This method is not thread safe!
/// </summary>
public static SafeMatrix SolveLinearSystemFast(SafeMatrix A, SafeMatrix X)
{
if (A.RowCount > 35 * 35)
throw new InvalidOperationException("Not a PSF fitting linear system.");
if (s_NumVariables != A.ColumnCount)
{
LinearSystemFastInitialiseSolution(A.ColumnCount, 35 * 35);
s_NumVariables = A.ColumnCount;
}
double[] a = A.GetElements();
double[] x = X.GetElements();
double[] y = new double[A.ColumnCount * X.ColumnCount];
SolveLinearSystemFast(a, x, A.RowCount, y);
SafeMatrix rv = new SafeMatrix(X.RowCount, X.ColumnCount);
rv.SetElements(y);
return rv;
}
