public void InvertMat() { CSqMatrix mat = new CSqMatrix(new Complex[, ] { { 11, 12, 4 }, { 0, 9, 30 }, { 40, 5, 3 } }); double n1 = (mat.Invert() * mat - SqMatrix.E(3)).CubeNorm; n1.Show(); double n2 = (mat.InvertSum * mat - SqMatrix.E(3)).CubeNorm; n2.Show(); Assert.IsTrue(n1 <= 1e-6 && n2 <= 1e-6); }
/// <summary> /// Вывести показатели качества обращения матриц /// </summary> private void NevaskShow() { $"Невязка D1: {(D1 * D1m - SqMatrix.E(D1.RowCount)).OctNorn}".Show(); $"Невязка D2: {(D2 * D2m - SqMatrix.E(D1.RowCount)).OctNorn}".Show(); }