public void ResetCalculationState() { var criterium = new ResidualStopCriterium(1e-3, 10); Assert.IsNotNull(criterium, "There should be a criterium"); var solution = new DenseVector(new[] {0.001, 0.001, 0.002}); var source = new DenseVector(new[] {0.001, 0.001, 0.002}); var residual = new DenseVector(new[] {1.000, 1.000, 2.001}); var status = criterium.DetermineStatus(5, solution, source, residual); Assert.AreEqual(IterationStatus.Continue, status, "Should be running"); criterium.Reset(); Assert.AreEqual(IterationStatus.Continue, criterium.Status, "Should not have started"); }
public void DetermineStatusWithSourceNaN() { var criterium = new ResidualStopCriterium(1e-3, 10); Assert.IsNotNull(criterium, "There should be a criterium"); var solution = new DenseVector(new[] {1.0, 1.0, 2.0}); var source = new DenseVector(new[] {1.0, 1.0, double.NaN}); var residual = new DenseVector(new[] {1000.0, 1000.0, 2001.0}); var status = criterium.DetermineStatus(5, solution, source, residual); Assert.AreEqual(IterationStatus.Diverged, status, "Should be diverged"); }
public void DetermineStatusWithNonMatchingSourceVectorThrowsArgumentException() { var criterium = new ResidualStopCriterium(1e-8, 50); Assert.IsNotNull(criterium, "There should be a criterium"); Assert.Throws<ArgumentException>(() => criterium.DetermineStatus( 1, DenseVector.Create(3, i => 4), DenseVector.Create(4, i => 4), DenseVector.Create(3, i => 4))); }
public void DetermineStatusWithIllegalIterationNumberThrowsArgumentOutOfRangeException() { var criterium = new ResidualStopCriterium(1e-8, 50); Assert.IsNotNull(criterium, "There should be a criterium"); Assert.Throws<ArgumentOutOfRangeException>(() => criterium.DetermineStatus( -1, DenseVector.Create(3, i => 4), DenseVector.Create(3, i => 5), DenseVector.Create(3, i => 6))); }
public void DetermineStatusWithConvergenceAtFirstIteration() { var criterium = new ResidualStopCriterium(); Assert.IsNotNull(criterium, "There should be a criterium"); var solution = new DenseVector(new[] {1.0, 1.0, 1.0}); var source = new DenseVector(new[] {1.0, 1.0, 1.0}); var residual = new DenseVector(new[] {0.0, 0.0, 0.0}); var status = criterium.DetermineStatus(0, solution, source, residual); Assert.AreEqual(IterationStatus.Converged, status, "Should be done"); }
public void DetermineStatus() { var criterium = new ResidualStopCriterium(1e-3, 10); Assert.IsNotNull(criterium, "There should be a criterium"); // the solution vector isn't actually being used so ... var solution = new DenseVector(new[] {double.NaN, double.NaN, double.NaN}); // Set the source values var source = new DenseVector(new[] {1.000, 1.000, 2.001}); // Set the residual values var residual = new DenseVector(new[] {0.001, 0.001, 0.002}); var status = criterium.DetermineStatus(5, solution, source, residual); Assert.AreEqual(IterationStatus.Continue, status, "Should still be running"); var status2 = criterium.DetermineStatus(16, solution, source, residual); Assert.AreEqual(IterationStatus.Converged, status2, "Should be done"); }