Changeset 23674 for issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

Timestamp:

01/31/19 14:21:20 (6 years ago)

Author:

seroussi

Message:

NEW: extended distance computation to Penta elements

File:

• issm/trunk-jpl/src/c/classes/Elements/Penta.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

@@ -647 +647 @@
 
 }/*}}}*/
+void       Penta::CreateDistanceInputFromSegmentlist(IssmDouble* distances,int distanceenum){/*{{{*/
+
+        /*Get current field and vertex coordinates*/
+        IssmDouble ls[NUMVERTICES],distance;
+        GetInputListOnVertices(&ls[0],distanceenum);
+
+        /*Get distance from list of segments and reset ls*/
+        for(int j=0;j<NUMVERTICES;j++){
+                distance=distances[this->vertices[j]->Lid()];
+                if(xIsNan<IssmDouble>(distance)) _error_("NaN found in vector");
+                if(xIsInf<IssmDouble>(distance)) _error_("Inf found in vector");
+
+                /*FIXME: do we really need this?*/
+                if(distanceenum==MaskIceLevelsetEnum) if(distance>10000) distance=10000;
+                if(ls[j]>0){
+                        ls[j] = distance;
+                }
+                else{
+                        ls[j] = - distance;
+                }
+        }
+
+        /*Update Levelset*/
+        this->inputs->AddInput(new PentaInput(distanceenum,&ls[0],P1Enum));
+}
+/*}}}*/
 void       Penta::ElementResponse(IssmDouble* presponse,int response_enum){/*{{{*/