[cig-commits] [commit] master: add another Assert (65aeaa2)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/76d4275352ef2cae5de9a073acd1c03a92c2670c...4f3d06fd1f3754419813db37ec9ef7f0f6f3cb15

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

commit 65aeaa2cd340e0573472e0b8862bee35bb046b10
Author: Timo Heister <timo.heister at gmail.com>
Date:   Mon May 19 20:05:59 2014 -0400

    add another Assert


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

65aeaa2cd340e0573472e0b8862bee35bb046b10
 source/simulator/nullspace.cc | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/source/simulator/nullspace.cc b/source/simulator/nullspace.cc
index aafd273..2091c35 100644
--- a/source/simulator/nullspace.cc
+++ b/source/simulator/nullspace.cc
@@ -188,6 +188,9 @@ namespace aspect
   Simulator<2>::remove_net_angular_momentum( LinearAlgebra::BlockVector &relevant_dst,
                                              LinearAlgebra::BlockVector &tmp_distributed_stokes )
   {
+    Assert(introspection.block_indices.velocities != introspection.block_indices.pressure,
+        ExcNotImplemented());
+
     // compute and remove angular momentum from velocity field, by computing
     // int rho V \cdot r_orth = omega  * int rho x^2
     const unsigned int dim=2;