Changeset 24270

Timestamp:

10/24/19 10:09:13 (5 years ago)

Author:

tsantos

Message:

BUG: do not compute complete adjoint if there is no ice on the element

File:

• issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp

@@ -56 +56 @@
 ElementMatrix* AdjointHorizAnalysis::CreateKMatrixFS(Element* element){/*{{{*/
 
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
+
         /*Intermediaries */
         bool        incomplete_adjoint;
@@ -144 +147 @@
 }/*}}}*/
 ElementMatrix* AdjointHorizAnalysis::CreateKMatrixHO(Element* element){/*{{{*/
+        
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
 
         /*Intermediaries */