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

Timestamp:

11/16/13 13:36:32 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: done with SSA 3d

File:

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

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

@@ -1446 +1446 @@
 
 }/*}}}*/
+/*FUNCTION Penta::GetVerticesCoordinatesBase(IssmDouble** pxyz_list){{{*/
+void Penta::GetVerticesCoordinatesBase(IssmDouble** pxyz_list){
+
+        IssmDouble* xyz_list = xNew<IssmDouble>(NUMVERTICES2D*3);
+        ::GetVerticesCoordinates(xyz_list,this->vertices,NUMVERTICES2D);
+
+        /*Assign output pointer*/
+        *pxyz_list = xyz_list;
+
+}/*}}}*/
 /*FUNCTION Penta::GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype) {{{*/
 void Penta::GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype){
@@ -2432 +2442 @@
 }
 /*}}}*/
+/*FUNCTION Penta::JacobianDeterminantBase{{{*/
+void Penta::JacobianDeterminantBase(IssmDouble* pJdet,IssmDouble* xyz_list_base,Gauss* gauss){
+
+        _assert_(gauss->Enum()==GaussPentaEnum);
+        this->GetTriaJacobianDeterminant(pJdet,xyz_list_base,(GaussPenta*)gauss);
+
+}
+/*}}}*/
 /*FUNCTION Penta::NoIceInElement {{{*/
 bool   Penta::NoIceInElement(){
@@ -2508 +2526 @@
 Gauss* Penta::NewGauss(int order){
         return new GaussPenta(order,order);
+}
+/*}}}*/
+/*FUNCTION Penta::NewGaussBase(int order){{{*/
+Gauss* Penta::NewGaussBase(int order){
+        return new GaussPenta(0,1,2,order);
 }
 /*}}}*/