[cig-commits] [commit] master: Add some commentary. (d97cebf)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/6b7c137fa297caa534f0c498bb594d05ed012a95...dd3794be780d4879f029fcef1efb957c61665109

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

commit d97cebfba4ff75aedf96be0d023ff44ebf741e93
Author: Wolfgang Bangerth <bangerth at math.tamu.edu>
Date:   Sun Jun 22 06:52:56 2014 -0500

    Add some commentary.


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

d97cebfba4ff75aedf96be0d023ff44ebf741e93
 source/simulator/core.cc | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/source/simulator/core.cc b/source/simulator/core.cc
index 1409a48..a7bb1ac 100644
--- a/source/simulator/core.cc
+++ b/source/simulator/core.cc
@@ -1230,6 +1230,17 @@ namespace aspect
       // transfer the data previously stored into the vectors indexed by
       // system_tmp. then ensure that the interpolated solution satisfies
       // hanging node constraints
+      //
+      // note that the 'constraints' variable contains hanging node constraints
+      // and constraints from periodic boundary conditions (as well as from
+      // zero and tangential velocity boundary conditions), but not from
+      // non-homogeneous boundary conditions. the latter are added not in setup_dofs(),
+      // which we call above, but added to 'current_constraints' in start_timestep(),
+      // which we do not want to call here.
+      //
+      // however, what we have should be sufficient: we have everything that
+      // is necessary to make the solution vectors *conforming* on the current
+      // mesh.
       system_trans.interpolate (system_tmp);
 
       constraints.distribute (distributed_system);