public void CanSolveForRandomVectorWhenResultVectorGivenUsingThinQR(int order) { var matrixA = new UserDefinedMatrix(Matrix <Complex> .Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorQR = matrixA.QR(QRMethod.Thin); var vectorb = new UserDefinedVector(Vector <Complex> .Build.Random(order, 1).ToArray()); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorQR.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { AssertHelpers.AlmostEqual(vectorb[i], matrixBReconstruct[i], 10); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorAndSymmetricMatrixWhenResultVectorGiven([Values(1, 2, 5, 10, 50, 100)] int order) { var A = new UserDefinedMatrix(Matrix <Complex32> .Build.RandomPositiveDefinite(order, 1).ToArray()); MatrixHelpers.ForceConjugateSymmetric(A); var ACopy = A.Clone(); var evd = A.Evd(); var b = new UserDefinedVector(Vector <Complex32> .Build.Random(order, 1).ToArray()); var bCopy = b.Clone(); var x = new UserDefinedVector(order); evd.Solve(b, x); var bReconstruct = A * x; // Check the reconstruction. AssertHelpers.AlmostEqual(b, bReconstruct, 2); // Make sure A/B didn't change. AssertHelpers.AlmostEqual(ACopy, A, 14); AssertHelpers.AlmostEqual(bCopy, b, 14); }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = new UserDefinedMatrix(Matrix <Complex> .Build.RandomPositiveDefinite(order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var chol = matrixA.Cholesky(); var b = new UserDefinedVector(Vector <Complex> .Build.Random(order, 1).ToArray()); var matrixBCopy = b.Clone(); var x = new UserDefinedVector(order); chol.Solve(b, x); Assert.AreEqual(b.Count, x.Count); var matrixBReconstruct = matrixA * x; // Check the reconstruction. for (var i = 0; i < order; i++) { AssertHelpers.AlmostEqual(b[i], matrixBReconstruct[i], 10); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < order; i++) { Assert.AreEqual(matrixBCopy[i], b[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = MatrixLoader.GenerateRandomUserDefinedMatrix(order, order); var matrixACopy = matrixA.Clone(); var factorGramSchmidt = matrixA.GramSchmidt(); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(order); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorGramSchmidt.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var bReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreApproximatelyEqual(vectorb[i], bReconstruct[i], 1.0e-11); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = new UserDefinedMatrix(Matrix <double> .Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorLU = matrixA.LU(); var vectorb = new UserDefinedVector(Vector <double> .Build.Random(order, 1).ToArray()); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorLU.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1.0e-11); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = MatrixLoader.GenerateRandomPositiveDefiniteUserDefinedMatrix(order); var matrixACopy = matrixA.Clone(); var chol = matrixA.Cholesky(); var b = MatrixLoader.GenerateRandomUserDefinedVector(order); var matrixBCopy = b.Clone(); var x = new UserDefinedVector(order); chol.Solve(b, x); Assert.AreEqual(b.Count, x.Count); var matrixBReconstruct = matrixA * x; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreEqual(b[i], matrixBReconstruct[i], 1.0e-11); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < order; i++) { Assert.AreEqual(matrixBCopy[i], b[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = MatrixLoader.GenerateRandomUserDefinedMatrix(order, order); var matrixACopy = matrixA.Clone(); var factorGramSchmidt = matrixA.GramSchmidt(); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(order); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorGramSchmidt.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i].Real, matrixBReconstruct[i].Real, 1e-3f); Assert.AreEqual(vectorb[i].Imaginary, matrixBReconstruct[i].Imaginary, 1e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGivenUsingThinQR(int order) { var matrixA = new UserDefinedMatrix(Matrix<double>.Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorQR = matrixA.QR(QRMethod.Thin); var vectorb = new UserDefinedVector(Vector<double>.Build.Random(order, 1).ToArray()); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorQR.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1.0e-11); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven([Values(1, 2, 5, 10, 50, 100)] int order) { var matrixA = MatrixLoader.GenerateRandomUserDefinedMatrix(order, order); var matrixACopy = matrixA.Clone(); var factorLU = matrixA.LU(); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(order); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorLU.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { AssertHelpers.AlmostEqual(vectorb[i], matrixBReconstruct[i], 9); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int row, int column) { var matrixA = new UserDefinedMatrix(Matrix<float>.Build.Random(row, column, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorSvd = matrixA.Svd(); var vectorb = new UserDefinedVector(Vector<float>.Build.Random(row, 1).ToArray()); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(column); factorSvd.Solve(vectorb, resultx); var matrixBReconstruct = matrixA*resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1e-4); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = new UserDefinedMatrix(Matrix<Complex>.Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorGramSchmidt = matrixA.GramSchmidt(); var vectorb = new UserDefinedVector(Vector<Complex>.Build.Random(order, 1).ToArray()); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorGramSchmidt.Solve(vectorb, resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { AssertHelpers.AlmostEqual(vectorb[i], matrixBReconstruct[i], 10); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVector(int order) { var matrixA = new UserDefinedMatrix(Matrix<Complex32>.Build.RandomPositiveDefinite(order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var chol = matrixA.Cholesky(); var b = new UserDefinedVector(Vector<Complex32>.Build.Random(order, 1).ToArray()); var x = chol.Solve(b); Assert.AreEqual(b.Count, x.Count); var matrixBReconstruct = matrixA * x; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreEqual(b[i].Real, matrixBReconstruct[i].Real, 1e-3f); Assert.AreEqual(b[i].Imaginary, matrixBReconstruct[i].Imaginary, 1e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int row, int column) { var matrixA = MatrixLoader.GenerateRandomUserDefinedMatrix(row, column); var matrixACopy = matrixA.Clone(); var factorSvd = matrixA.Svd(true); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(row); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(column); factorSvd.Solve(vectorb,resultx); var bReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreApproximatelyEqual(vectorb[i].Real, bReconstruct[i].Real, 1e-3f); Assert.AreApproximatelyEqual(vectorb[i].Imaginary, bReconstruct[i].Imaginary, 1e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorAndSymmetricMatrixWhenResultVectorGiven(int order) { var matrixA = new UserDefinedMatrix(Matrix<double>.Build.RandomPositiveDefinite(order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorEvd = matrixA.Evd(); var vectorb = new UserDefinedVector(Vector<double>.Build.Random(order, 1).ToArray()); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorEvd.Solve(vectorb, resultx); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1.0e-10); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorAndSymmetricMatrix([Values(1, 2, 5, 10, 50, 100)] int order) { var matrixA = new UserDefinedMatrix(Matrix<Complex32>.Build.RandomPositiveDefinite(order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorEvd = matrixA.Evd(); var vectorb = new UserDefinedVector(Vector<Complex32>.Build.Random(order, 1).ToArray()); var resultx = factorEvd.Solve(vectorb); Assert.AreEqual(matrixA.ColumnCount, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i].Real, matrixBReconstruct[i].Real, 1e-3f); Assert.AreEqual(vectorb[i].Imaginary, matrixBReconstruct[i].Imaginary, 1e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } }
public void CanSolveForRandomVectorUsingThinQR(int order) { var matrixA = new UserDefinedMatrix(Matrix<Complex>.Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorQR = matrixA.QR(QRMethod.Thin); var vectorb = new UserDefinedVector(Vector<Complex>.Build.Random(order, 1).ToArray()); var resultx = factorQR.Solve(vectorb); Assert.AreEqual(matrixA.ColumnCount, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < order; i++) { AssertHelpers.AlmostEqual(vectorb[i], matrixBReconstruct[i], 10); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } }
public void CanSolveForRandomVector(int row, int column) { var matrixA = new UserDefinedMatrix(Matrix<Complex>.Build.Random(row, column, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorSvd = matrixA.Svd(); var vectorb = new UserDefinedVector(Vector<Complex>.Build.Random(row, 1).ToArray()); var resultx = factorSvd.Solve(vectorb); Assert.AreEqual(matrixA.ColumnCount, resultx.Count); var matrixBReconstruct = matrixA*resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { AssertHelpers.AlmostEqual(vectorb[i], matrixBReconstruct[i], 10); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } }
public void CanSolveForRandomVectorWhenResultVectorGiven([Values(1, 2, 5, 9, 50, 90)] int row, [Values(1, 2, 5, 10, 50, 100)] int column) { var matrixA = MatrixLoader.GenerateRandomUserDefinedMatrix(row, column); var matrixACopy = matrixA.Clone(); var factorSvd = matrixA.Svd(true); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(row); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(column); factorSvd.Solve(vectorb, resultx); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1.0e-11); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVector(int order) { var matrixA = new UserDefinedMatrix(Matrix<Complex32>.Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorGramSchmidt = matrixA.GramSchmidt(); var vectorb = new UserDefinedVector(Vector<Complex32>.Build.Random(order, 1).ToArray()); var resultx = factorGramSchmidt.Solve(vectorb); Assert.AreEqual(matrixA.ColumnCount, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreEqual(vectorb[i].Real, matrixBReconstruct[i].Real, 1e-3f); Assert.AreEqual(vectorb[i].Imaginary, matrixBReconstruct[i].Imaginary, 1e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } }
public void CanSolveForRandomVector(int order) { var matrixA = new UserDefinedMatrix(Matrix<float>.Build.Random(order, order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var factorLU = matrixA.LU(); var vectorb = new UserDefinedVector(Vector<float>.Build.Random(order, 1).ToArray()); var resultx = factorLU.Solve(vectorb); Assert.AreEqual(matrixA.ColumnCount, resultx.Count); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1e-4); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } }
public void CanSolveForRandomVectorAndSymmetricMatrixWhenResultVectorGiven([Values(1, 2, 5, 10, 50, 100)] int order) { var A = new UserDefinedMatrix(Matrix<Complex>.Build.RandomPositiveDefinite(order, 1).ToArray()); MatrixHelpers.ForceHermitian(A); var ACopy = A.Clone(); var evd = A.Evd(Symmetricity.Hermitian); var b = new UserDefinedVector(Vector<Complex>.Build.Random(order, 1).ToArray()); var bCopy = b.Clone(); var x = new UserDefinedVector(order); evd.Solve(b, x); var bReconstruct = A * x; // Check the reconstruction. AssertHelpers.AlmostEqual(b, bReconstruct, 9); // Make sure A/B didn't change. AssertHelpers.AlmostEqual(ACopy, A, 14); AssertHelpers.AlmostEqual(bCopy, b, 14); }
public void CanSolveForRandomVectorWhenResultVectorGiven([Values(1, 2, 5, 10, 50, 100)] int order) { var matrixA = MatrixLoader.GenerateRandomPositiveDefiniteHermitianUserDefinedMatrix(order); var matrixACopy = matrixA.Clone(); var chol = matrixA.Cholesky(); var b = MatrixLoader.GenerateRandomUserDefinedVector(order); var matrixBCopy = b.Clone(); var x = new UserDefinedVector(order); chol.Solve(b, x); Assert.AreEqual(b.Count, x.Count); var matrixBReconstruct = matrixA * x; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreEqual(b[i].Real, matrixBReconstruct[i].Real, 1e-3f); Assert.AreEqual(b[i].Imaginary, matrixBReconstruct[i].Imaginary, 1e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < order; i++) { Assert.AreEqual(matrixBCopy[i], b[i]); } }
public void SolveVectorIfVectorsNotComputedThrowsInvalidOperationException() { var matrixA = new UserDefinedMatrix(Matrix<float>.Build.Random(10, 10, 1).ToArray()); var factorSvd = matrixA.Svd(false); var vectorb = new UserDefinedVector(Vector<float>.Build.Random(10, 1).ToArray()); Assert.That(() => factorSvd.Solve(vectorb), Throws.InvalidOperationException); }
public void CanSolveForRandomVectorAndSymmetricMatrix([Values(1, 2, 5, 10, 50, 100)] int order) { var A = new UserDefinedMatrix(Matrix<float>.Build.RandomPositiveDefinite(order, 1).ToArray()); MatrixHelpers.ForceSymmetric(A); var ACopy = A.Clone(); var evd = A.Evd(); var b = new UserDefinedVector(Vector<float>.Build.Random(order, 1).ToArray()); var bCopy = b.Clone(); var x = evd.Solve(b); var bReconstruct = A * x; // Check the reconstruction. AssertHelpers.AlmostEqual(b, bReconstruct, -1); // Make sure A/B didn't change. AssertHelpers.AlmostEqual(ACopy, A, 14); AssertHelpers.AlmostEqual(bCopy, b, 14); }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = MatrixLoader.GenerateRandomUserDefinedMatrix(order, order); var matrixACopy = matrixA.Clone(); var factorQR = matrixA.QR(); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(order); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorQR.Solve(vectorb,resultx); Assert.AreEqual(vectorb.Count, resultx.Count); var bReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreApproximatelyEqual(vectorb[i], bReconstruct[i], 1.0e-11); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorAndSymmetricMatrixWhenResultVectorGiven([Values(1, 2, 5, 10, 50, 100)] int order) { var matrixA = MatrixLoader.GenerateRandomPositiveDefiniteUserDefinedMatrix(order); var matrixACopy = matrixA.Clone(); var factorEvd = matrixA.Evd(); var vectorb = MatrixLoader.GenerateRandomUserDefinedVector(order); var vectorbCopy = vectorb.Clone(); var resultx = new UserDefinedVector(order); factorEvd.Solve(vectorb, resultx); var matrixBReconstruct = matrixA * resultx; // Check the reconstruction. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorb[i], matrixBReconstruct[i], 1e-3); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < vectorb.Count; i++) { Assert.AreEqual(vectorbCopy[i], vectorb[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = new UserDefinedMatrix(Matrix<float>.Build.RandomPositiveDefinite(order, 1).ToArray()); var matrixACopy = matrixA.Clone(); var chol = matrixA.Cholesky(); var b = new UserDefinedVector(Vector<float>.Build.Random(order, 1).ToArray()); var matrixBCopy = b.Clone(); var x = new UserDefinedVector(order); chol.Solve(b, x); Assert.AreEqual(b.Count, x.Count); var matrixBReconstruct = matrixA * x; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreEqual(b[i], matrixBReconstruct[i], 0.5); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < order; i++) { Assert.AreEqual(matrixBCopy[i], b[i]); } }
public void CanSolveForRandomVectorWhenResultVectorGiven(int order) { var matrixA = MatrixLoader.GenerateRandomPositiveDefiniteUserDefinedMatrix(order); var matrixACopy = matrixA.Clone(); var chol = matrixA.Cholesky(); var b = MatrixLoader.GenerateRandomUserDefinedVector(order); var bCopy = b.Clone(); var x = new UserDefinedVector(order); chol.Solve(b, x); Assert.AreEqual(b.Count, x.Count); var bReconstruct = matrixA * x; // Check the reconstruction. for (var i = 0; i < order; i++) { Assert.AreApproximatelyEqual(b[i], bReconstruct[i], 10e-3f); } // Make sure A didn't change. for (var i = 0; i < matrixA.RowCount; i++) { for (var j = 0; j < matrixA.ColumnCount; j++) { Assert.AreEqual(matrixACopy[i, j], matrixA[i, j]); } } // Make sure b didn't change. for (var i = 0; i < order; i++) { Assert.AreEqual(bCopy[i], b[i]); } }