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

[dealii.demon at gmail.com]

Revision 2031

Make initial solution solve the Stokes problem.

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


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

Diff:
Modified: trunk/aspire/source/simulator/solver.cc
===================================================================
--- trunk/aspire/source/simulator/solver.cc	2013-11-21 18:05:40 UTC (rev 2030)
+++ trunk/aspire/source/simulator/solver.cc	2013-11-21 18:19:58 UTC (rev 2031)
@@ -290,8 +290,14 @@
       }    
     // 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));
+    if (time_step == 0)
+      solution_velocity = distributed_stokes_solution;
+    
+    else
+    {
+      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