private static void TestFactorization() { //TestSettings.RunMultiproviderTest(providers, delegate () { // var skyline = SkylineMatrix.CreateFromArrays(SparsePosDef10by10.Order, SparsePosDef10by10.SkylineValues, SparsePosDef10by10.SkylineDiagOffsets, true, true); var expectedU = Matrix.CreateFromArray(SparsePosDef10by10.CholeskyU); CholeskySkyline factorization = skyline.FactorCholesky(false); TriangularUpper computedU = factorization.GetFactorU(); comparer.AssertEqual(expectedU, computedU); //}); }
private static void TestSystemSolution(/*LinearAlgebraProviderChoice providers*/) { //TestSettings.RunMultiproviderTest(providers, delegate () { // var skyline = SkylineMatrix.CreateFromArrays(SparsePosDef10by10.Order, SparsePosDef10by10.SkylineValues, SparsePosDef10by10.SkylineDiagOffsets, true, true); var b = Vector.CreateFromArray(SparsePosDef10by10.Rhs); var xExpected = Vector.CreateFromArray(SparsePosDef10by10.Lhs); CholeskySkyline factor = skyline.FactorCholesky(false); Vector xComputed = factor.SolveLinearSystem(b); comparer.AssertEqual(xExpected, xComputed); //}); }