本文整理汇总了C#中IMatrix.Submatrix方法的典型用法代码示例。如果您正苦于以下问题：C# IMatrix.Submatrix方法的具体用法？C# IMatrix.Submatrix怎么用？C# IMatrix.Submatrix使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类IMatrix的用法示例。
在下文中一共展示了IMatrix.Submatrix方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C#代码示例。
示例1: Solve
public IMatrix Solve(IMatrix b)
{
if (b.Rows != _lu.Rows) throw new ArgumentException("Invalid matrix dimensions.");
if (!IsNonSingular) throw new InvalidOperationException("Matrix is singular");
// Copy right hand side with pivoting
int count = b.Columns;
IMatrix x = b.Submatrix(_pivotVector, 0, count - 1);
//int rows = LU.Rows;
int columns = _lu.Columns;
double[][] lu = _lu.Array;
// Solve L*Y = B(piv,:)
for (int k = 0; k < columns; k++)
{
for (int i = k + 1; i < columns; i++)
{
for (int j = 0; j < count; j++)
{
x[i, j] -= x[k, j]*lu[i][k];
}
}
}
// Solve U*X = Y;
for (int k = columns - 1; k >= 0; k--)
{
for (int j = 0; j < count; j++)
{
x[k, j] /= lu[k][k];
}
for (int i = 0; i < k; i++)
{
for (int j = 0; j < count; j++)
{
x[i, j] -= x[k, j]*lu[i][k];
}
}
}
return x;
}