Changeset 3854

Timestamp:

05/19/10 16:09:28 (15 years ago)

Author:

Mathieu Morlighem

Message:

new switch

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (2 diffs)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -1089 +1089 @@
 
                                 /*Get strain rate from velocity: */
-                                inputs->GetStrainRate3d(&epsilon[0],&xyz_list[0][0],gauss_coord,VxEnum,VyEnum,VzEnum,PattynFormulationEnum);
-                                inputs->GetStrainRate3d(&oldepsilon[0],&xyz_list[0][0],gauss_coord,VxOldEnum,VyOldEnum,VzObsEnum,PattynFormulationEnum);
+                                inputs->GetStrainRate3dPattyn(&epsilon[0],&xyz_list[0][0],gauss_coord,VxEnum,VyEnum);
+                                inputs->GetStrainRate3dPattyn(&oldepsilon[0],&xyz_list[0][0],gauss_coord,VxOldEnum,VyOldEnum);
 
                                 /*Get viscosity: */
@@ -4370 +4370 @@
 /*FUNCTION Penta::UpdateInputsFromVector(double* vector, int name, int type);{{{1*/
 void  Penta::UpdateInputsFromVector(double* vector, int name, int type){
-        ISSMERROR(" not supported yet!");
+
+        switch(type){
+
+                case VertexEnum:
+
+                        /*New PentaVertexInpu*/
+                        double values[6];
+
+                        /*Get values on the 6 vertices*/
+                        for (int i=0;i<6;i++){
+                                values[i]=vector[this->nodes[i]->GetVertexDof()];
+                        }
+
+                        /*update input*/
+                        this->inputs->AddInput(new PentaVertexInput(name,values));
+                        return;
+
+                default:
+
+                        ISSMERROR("type %i (%s) not implemented yet",type,EnumAsString(type));
+        }
 }
 /*}}}*/