/// <summary>
/// 解方程，求出本次迭代的外方位元素
/// </summary>
private double[,] SolveEqu(Matrix[] aMat, Matrix[] lMat)
{
Matrix mergeA = MergeMatrix(aMat), mergeL = MergeMatrix(lMat);
var AT = MatrixOperation.MatrixTrans(mergeA);
var final = (1 / (AT * mergeA)) * AT * mergeL;
// 结果一定是6*1的矩阵
return(final.Data);
}
public static Matrix operator +(Matrix a, Matrix b)
{
return(MatrixOperation.Add(a, b));
}
public static Matrix operator *(double a, Matrix b)
{
return(MatrixOperation.SimpleMultiply(a, b));
}
public static Matrix operator *(Matrix a, Matrix b)
{
return(MatrixOperation.Multiply(a, b));
}
public static Matrix operator -(Matrix a, Matrix b)
{
return(MatrixOperation.Sub(a, b));
}