[cig-commits] [commit] master: add an assert (5eeb1dd)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/426fbfcad3568313f61490f77127c5e1c97a005d...a871ef437febf4fd67ce6f1ff7349b8a57987c8a

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

commit 5eeb1dd797dccf93be572569fd875b2e87647d3b
Author: Juliane Dannberg <dannberg at gfz-potsdam.de>
Date:   Mon May 19 21:37:20 2014 +0200

    add an assert


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

5eeb1dd797dccf93be572569fd875b2e87647d3b
 source/simulator/assembly.cc | 1 +
 1 file changed, 1 insertion(+)

diff --git a/source/simulator/assembly.cc b/source/simulator/assembly.cc
index bd6fe5b..1242e6a 100644
--- a/source/simulator/assembly.cc
+++ b/source/simulator/assembly.cc
@@ -1576,6 +1576,7 @@ namespace aspect
             *
             (density_c_P + latent_heat_LHS);
 
+        AssertThrow(density_c_P + latent_heat_LHS, ExcMessage("mass matrix must be positive"));
 
         const Tensor<1,dim> current_u = scratch.current_velocity_values[q];
         const double factor = (use_bdf2_scheme)? ((2*time_step + old_time_step) /