本文整理汇总了C#中DotNetMatrix.GeneralMatrix.LUD方法的典型用法代码示例。如果您正苦于以下问题:C# GeneralMatrix.LUD方法的具体用法?C# GeneralMatrix.LUD怎么用?C# GeneralMatrix.LUD使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类DotNetMatrix.GeneralMatrix
的用法示例。
在下文中一共展示了GeneralMatrix.LUD方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C#代码示例。
示例1: LUDecomposition
public void LUDecomposition()
{
GeneralMatrix A = new GeneralMatrix(columnwise, 4);
int n = A.ColumnDimension;
A = A.GetMatrix(0, n - 1, 0, n - 1);
A.SetElement(0, 0, 0.0);
LUDecomposition LU = A.LUD();
Assert.IsTrue(GeneralTests.Check(A.GetMatrix(LU.Pivot, 0, n - 1), LU.L.Multiply(LU.U)));
}
示例2: Main
//.........这里部分代码省略.........
}
catch (System.SystemException e)
{
errorCount = try_failure(errorCount, "SingularValueDecomposition...", "incorrect singular value decomposition calculation");
System.Console.Out.WriteLine(e.Message);
}
DEF = new GeneralMatrix(rankdef);
try
{
check(DEF.Rank(), System.Math.Min(DEF.RowDimension, DEF.ColumnDimension) - 1);
try_success("Rank()...", "");
}
catch (System.SystemException e)
{
errorCount = try_failure(errorCount, "Rank()...", "incorrect Rank calculation");
System.Console.Out.WriteLine(e.Message);
}
B = new GeneralMatrix(condmat);
SVD = B.SVD();
double[] singularvalues = SVD.SingularValues;
try
{
check(B.Condition(), singularvalues[0] / singularvalues[System.Math.Min(B.RowDimension, B.ColumnDimension) - 1]);
try_success("Condition()...", "");
}
catch (System.SystemException e)
{
errorCount = try_failure(errorCount, "Condition()...", "incorrect condition number calculation");
System.Console.Out.WriteLine(e.Message);
}
int n = A.ColumnDimension;
A = A.GetMatrix(0, n - 1, 0, n - 1);
A.SetElement(0, 0, 0.0);
LUDecomposition LU = A.LUD();
try
{
check(A.GetMatrix(LU.Pivot, 0, n - 1), LU.L.Multiply(LU.U));
try_success("LUDecomposition...", "");
}
catch (System.SystemException e)
{
errorCount = try_failure(errorCount, "LUDecomposition...", "incorrect LU decomposition calculation");
System.Console.Out.WriteLine(e.Message);
}
X = A.Inverse();
try
{
check(A.Multiply(X), GeneralMatrix.Identity(3, 3));
try_success("Inverse()...", "");
}
catch (System.SystemException e)
{
errorCount = try_failure(errorCount, "Inverse()...", "incorrect Inverse calculation");
System.Console.Out.WriteLine(e.Message);
}
O = new GeneralMatrix(SUB.RowDimension, 1, 1.0);
SOL = new GeneralMatrix(sqSolution);
SQ = SUB.GetMatrix(0, SUB.RowDimension - 1, 0, SUB.RowDimension - 1);
try
{
check(SQ.Solve(SOL), O);
try_success("Solve()...", "");
}
catch (System.ArgumentException e1)
{
errorCount = try_failure(errorCount, "Solve()...", e1.Message);