[cig-commits] [commit] master: add another Assert (441e405)

[cig_noreply at geodynamics.org]

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

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

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

commit 441e405c808fdff14ccde87c6129fb1839624311
Author: Timo Heister <timo.heister at gmail.com>
Date:   Mon May 19 16:39:26 2014 -0400

    add another Assert


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

441e405c808fdff14ccde87c6129fb1839624311
 source/simulator/nullspace.cc | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/source/simulator/nullspace.cc b/source/simulator/nullspace.cc
index 122a145..2c7ee73 100644
--- a/source/simulator/nullspace.cc
+++ b/source/simulator/nullspace.cc
@@ -123,6 +123,9 @@ namespace aspect
 
     if (funcs.size()>0)
       {
+        Assert(introspection.block_indices.velocities != introspection.block_indices.pressure,
+            ExcNotImplemented());
+
         net_rotations_translations.resize(funcs.size());
         for (unsigned int i=0; i<funcs.size(); ++i)
           {
@@ -149,7 +152,10 @@ namespace aspect
     if (parameters.nullspace_removal & NullspaceRemoval::net_rotation ||
         parameters.nullspace_removal & NullspaceRemoval::net_translation)
       {
-        for (unsigned int i=0; i<net_rotations_translations.size(); ++i)
+        Assert(introspection.block_indices.velocities != introspection.block_indices.pressure,
+            ExcNotImplemented());
+
+       for (unsigned int i=0; i<net_rotations_translations.size(); ++i)
           {
             // compute the magnitude of the solution vector in direction
             // of this null space vector and subtract the corresponding multiple