[cig-commits] commit 1994 by buerg to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1994

Boundary conditions should be satisfied in every iteration step.

U   trunk/aspire/source/simulator/solver.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1994&peg=1994

Diff:
Modified: trunk/aspire/source/simulator/solver.cc
===================================================================
--- trunk/aspire/source/simulator/solver.cc	2013-10-25 16:30:17 UTC (rev 1993)
+++ trunk/aspire/source/simulator/solver.cc	2013-10-25 17:34:45 UTC (rev 1994)
@@ -288,15 +288,15 @@
       
         solver.solve (system_matrix_velocity, distributed_stokes_solution, system_rhs_velocity, preconditioner);
       }    
-    // distribute hanging node and
-    // other constraints
-    current_constraints_velocity.distribute (distributed_stokes_solution);
-
     // then copy back the solution from the temporary (non-ghosted) vector
     // into the ghosted one with all solution components
     solution_velocity.block (0).equ (1.0 - parameters.damping, old_solution_velocity.block (0), parameters.damping, distributed_stokes_solution.block (0));
     solution_velocity.block (1).equ (1.0 - parameters.damping, old_solution_velocity.block (1), parameters.damping, distributed_stokes_solution.block (1));
 
+    // distribute hanging node and
+    // other constraints
+    current_constraints_velocity.distribute (solution_velocity);
+
     // now rescale the pressure back to real physical units
     normalize_pressure (solution_velocity);