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

[cig_noreply at geodynamics.org]

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

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

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

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

    add an assert


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

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

diff --git a/source/simulator/assembly.cc b/source/simulator/assembly.cc
index 3d4442d..fe674a2 100644
--- a/source/simulator/assembly.cc
+++ b/source/simulator/assembly.cc
@@ -1568,6 +1568,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) /