/// <summary>
/// Determina a solução de um sistema linear de equações.
/// </summary>
/// <remarks>
/// A matriz e o vector independente são alterados durante o processo.
/// </remarks>
/// <param name="coefficientsMatrix">A matriz dos coeficientes.</param>
/// <param name="independentVector">O vector independente.</param>
/// <returns>A solução do sistema.</returns>
public LinearSystemSolution <ElementType> Run(
IMathMatrix <ElementType> coefficientsMatrix,
IMathMatrix <ElementType> independentVector)
{
this.condensationAlgorithm.Run(coefficientsMatrix, independentVector);
var result = new LinearSystemSolution <ElementType>();
var matrixLines = coefficientsMatrix.GetLength(0);
var matrixColumns = coefficientsMatrix.GetLength(1);
var currentPivotLine = 0;
var currentPivotColumn = 0;
var lastNonZeroColumn = -1;
var independentSolutionVector = new ArrayMathVector <ElementType>(matrixColumns, this.field.AdditiveUnity);
while (currentPivotLine < matrixLines && currentPivotColumn < matrixColumns)
{
var pivotValue = coefficientsMatrix[currentPivotLine, currentPivotColumn];
if (this.field.IsAdditiveUnity(pivotValue))
{
if (this.field.IsAdditiveUnity(independentSolutionVector[currentPivotLine]))
{
var basisVector = new ArrayMathVector <ElementType>(matrixColumns, this.field.AdditiveUnity);
basisVector[currentPivotColumn] = this.field.AdditiveInverse(
this.field.MultiplicativeUnity);
var i = currentPivotLine - 1;
var j = lastNonZeroColumn;
while (i >= 0 && j >= 0)
{
var linePivotValue = coefficientsMatrix[i, j];
if (this.field.IsMultiplicativeUnity(linePivotValue))
{
basisVector[j] = coefficientsMatrix[i, currentPivotColumn];
--i;
}
--j;
}
result.VectorSpaceBasis.Add(basisVector);
}
else
{
result.Vector = null;
result.VectorSpaceBasis.Clear();
return(result);
}
}
else
{
lastNonZeroColumn = currentPivotColumn;
independentSolutionVector[currentPivotColumn] = independentVector[currentPivotLine, 0];
++currentPivotLine;
}
++currentPivotColumn;
}
result.Vector = independentSolutionVector;
return(result);
}
/// <summary>
/// Constrói uma cópida do vector actual.
/// </summary>
/// <returns>A cópia do vector actual.</returns>
public override IVector <CoeffType> Clone()
{
var result = new ArrayMathVector <CoeffType>();
result.vectorEntries = new CoeffType[this.vectorEntries.Length];
Array.Copy(this.vectorEntries, result.vectorEntries, this.vectorEntries.Length);
return(result);
}
