diff --git a/src/elasticity/solvers/IterativePressureSolver.cpp b/src/elasticity/solvers/IterativePressureSolver.cpp
index 9fa9b21eb82b72c27668cd16a165ac27435a2d9d..0359df8f9305e0838b25c14e1a84ee419db30b69 100644
--- a/src/elasticity/solvers/IterativePressureSolver.cpp
+++ b/src/elasticity/solvers/IterativePressureSolver.cpp
@@ -33,7 +33,7 @@ bool IterativePressureSolver::Method(IElasticity &assemble, Vector &u) {
   while (!iteration.IsFinished()) {
     int n = iteration.Step();
 
-    if (n > 1) {
+    if (n > 0) {
       uOld = uNew;  // u^{n-1}
     }
 
@@ -55,8 +55,8 @@ bool IterativePressureSolver::Method(IElasticity &assemble, Vector &u) {
         return Method(assemble, u);
       return false;
     }
-    // Extrapolation ab n = 2
-    if (useExtrapolation && n >= 2) {
+    // Extrapolation ab n = 1
+    if (useExtrapolation && n >= 1) {
       deltaU = uNew - uOld;
       uNew += dampingExtrapolation * deltaU;
     }