/// <summary>
/// Reduces matrix to row-echelon (REF/Gauss) or reduced row-echelon (RREF/Gauss-Jordan) form and solves for augmented columns.
/// Returns the matrix solution and outputs the full matrix (reduced matrix along with the solution.
/// Accepts the number of augmeted columns. If the number specified is null, the default number specified in the matrix is used.
/// </summary>
/// <returns>
/// Returns the matrix solution and outputs the full matrix (reduced matrix along with the solution.
/// </returns>
/// <remarks>
/// The default value for <seealso cref="Elsheimy.Components.Linears.Matrix.AugmentedColumnCount"/> is used.
/// </remarks>
public virtual Matrix Solve(int?augmentedColCount, out Matrix fullMatrix)
{
int augmentedCols = augmentedColCount ?? this.AugmentedColumnCount;
if (augmentedCols <= 0)
{
throw new InvalidOperationException(Properties.Resources.Exception_NoAugmentedColumns);
}
var result = MatrixFunctions.Eliminate(this.InnerMatrix, MatrixReductionForm.ReducedRowEchelonForm, augmentedCols);
var state = result.SolutionState;
fullMatrix = new Matrix(result.FullMatrix, augmentedCols);
if (result.Solution == null)
{
return(null);
}
return(new Matrix(result.Solution));
}
/// <summary>
/// Reduces matrix to row-echelon (REF/Gauss) or reduced row-echelon (RREF/Gauss-Jordan) form.
/// Accepts the number of augmeted columns. If the number specified is null, the default number specified in the matrix is used.
/// </summary>
/// <remarks>
/// If <param name="augmentedColCount">augmentedColCount</param> is null, <seealso cref="Elsheimy.Components.Linears.Matrix.AugmentedColumnCount"/> is used.
/// </remarks>
public virtual Matrix Reduce(MatrixReductionForm form, int?augmentedColCount = null)
{
int augmentedCols = augmentedColCount ?? this.AugmentedColumnCount;
return(new Matrix(MatrixFunctions.Eliminate(this.InnerMatrix, form, augmentedCols).FullMatrix));
}
