Changeset 23962 for issm/trunk-jpl/src/c/classes/Elements/TriaRef.cpp

Timestamp:

05/30/19 13:56:00 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: working on channels

File:

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

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

@@ -328 +328 @@
         /*Clean up*/
         xDelete<IssmDouble>(triabasis);
+}
+/*}}}*/
+void TriaRef::GetSegmentNodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list_tria,Gauss* gauss,int index1,int index2,int finiteelement){/*{{{*/
+        /*This routine returns the values of the nodal functions  at the gaussian point.*/
+
+        _assert_(index1>=0 && index1<3);
+        _assert_(index2>=0 && index2<3);
+
+        /*Fetch number of nodes for this finite element*/
+        int numnodes = this->NumberofNodes(finiteelement);
+
+        /*Get nodal functions*/
+        IssmDouble* dtriabasis=xNew<IssmDouble>(2*numnodes);
+        GetNodalFunctionsDerivatives(dtriabasis,xyz_list_tria,gauss,finiteelement);
+
+        switch(finiteelement){
+                case P1Enum: case P1DGEnum:
+                        dbasis[2*0+0] = dtriabasis[numnodes*0+index1];
+                        dbasis[2*0+1] = dtriabasis[numnodes*1+index1];
+                        dbasis[2*1+0] = dtriabasis[numnodes*0+index2];
+                        dbasis[2*1+1] = dtriabasis[numnodes*1+index2];
+                        break;
+                default:
+                        _error_("Element type "<<EnumToStringx(finiteelement)<<" not supported yet");
+        }
+
+        /*Clean up*/
+        xDelete<IssmDouble>(dtriabasis);
 }
 /*}}}*/