included tolerance

cardmech/src/electrophysiology/solvers/LinearImplicitSolver.cpp

@@ -21,7 +21,7 @@ void LinearImplicitSolver::Initialize(Vector &V, double dt) {
 void LinearImplicitSolver::Solve(Vector &V, TimeSeries &timeSeries) {
   Initialize(V, timeSeries.StepSize());
   InitialPlotting(V, V);
-
+  const double tTol =1e-10;
   int step = 1;
   for (double t = timeSeries.NextTStepJ();; t = timeSeries.NextTStepJ(), ++step) {
 
@@ -41,7 +41,7 @@ void LinearImplicitSolver::Solve(Vector &V, TimeSeries &timeSeries) {
 
     vout(10).EndBlock();
 
-    if (t >= timeSeries.LastTStep()) break;
+    if (t + tTol > timeSeries.LastTStep()) break;
 
 
     A.updateActivationTime((*gating)[0], t);
@@ -250,6 +250,7 @@ void ImplicitSolver::Solve(Vector &V, TimeSeries &timeSeries) {
   Initialize(V, timeSeries.StepSize());
   InitialPlotting(V, V);
 
+  const double tTol = 1e-10;
   int step = 1;
   for (double t = timeSeries.NextTStepJ();; t = timeSeries.NextTStepJ(), ++step) {
 
@@ -264,7 +265,7 @@ void ImplicitSolver::Solve(Vector &V, TimeSeries &timeSeries) {
     StepPlotting(step, (*gating)[0]);
     vout(10).EndBlock();
 
-    if (t >= timeSeries.LastTStep()) break;
+    if (t + tTol > timeSeries.LastTStep()) break;
 
 
     A.updateActivationTime((*gating)[0], t);
@@ -312,6 +313,7 @@ void SemiImplicitSolver::Solve(Vector &V, TimeSeries &timeSeries) {
   Initialize(V, timeSeries.StepSize());
 
   InitialPlotting(V, V);
+  const double tTol=1e-10;
   int step = 1;
   for (double t = timeSeries.NextTStepJ();; t = timeSeries.NextTStepJ(), ++step) {
 
@@ -324,7 +326,7 @@ void SemiImplicitSolver::Solve(Vector &V, TimeSeries &timeSeries) {
     StaggeredStep(timeSeries.PreviousTStep(), timeSeries.StepSize());
     StepPlotting(step, (*gating)[0]);
     vout(10).EndBlock();
-    if (t >= timeSeries.LastTStep()) break;
+    if (t + tTol > timeSeries.LastTStep()) break;
     A.updateActivationTime((*gating)[0], t);
   }
   V = (*gating)[0];