[cig-commits] [commit] master: fix the residuum of the iterated IMPES solver (9a9fdd2)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/aed4d3ee69f4a126db2602d3d0bb6fcc4a49633a...4812dd65891608db0aa2a1b96dd6c083d000fd89

>---------------------------------------------------------------

commit 9a9fdd2d3e161d9760a7b738270d9bfda7b02158
Author: Juliane Dannberg <dannberg at gfz-potsdam.de>
Date:   Tue Nov 11 16:29:02 2014 +0100

    fix the residuum of the iterated IMPES solver


>---------------------------------------------------------------

9a9fdd2d3e161d9760a7b738270d9bfda7b02158
 source/simulator/solver.cc | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/source/simulator/solver.cc b/source/simulator/solver.cc
index dac5a13..40576c1 100644
--- a/source/simulator/solver.cc
+++ b/source/simulator/solver.cc
@@ -391,6 +391,9 @@ namespace aspect
       introspection.index_sets.system_partitioning[block_idx],
       mpi_communicator);
 
+    // solve the linear system:
+    current_constraints.set_zero(distributed_solution);
+
     // Compute the residual before we solve and return this at the end.
     // This is used in the nonlinear solver.
     const double initial_residual = system_matrix.block(block_idx,block_idx).residual
@@ -398,8 +401,6 @@ namespace aspect
                                      distributed_solution.block(block_idx),
                                      system_rhs.block(block_idx));
 
-    // solve the linear system:
-    current_constraints.set_zero(distributed_solution);
     try
       {
         solver.solve (system_matrix.block(block_idx,block_idx),