Changeset 9035

Timestamp:

07/18/11 09:02:33 (14 years ago)

Author:

seroussi

Message:

added GetStrainRate for double inputs

Location:

issm/trunk/src/c/objects/Inputs

Files:

• DoubleInput.cpp (modified) (1 diff)
• DoubleInput.h (modified) (1 diff)

issm/trunk/src/c/objects/Inputs/DoubleInput.cpp

@@ -194 +194 @@
 void DoubleInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussPenta* gauss){_error_(" not supported yet!");}
 /*}}}*/
+/*FUNCTION DoubleInput::GetVxStrainRate2d(double* epsilonvx,double* xyz_list, GaussTria* gauss){{{1*/
+void DoubleInput::GetVxStrainRate2d(double* epsilonvx,double* xyz_list, GaussTria* gauss){
+        /*Epsilon is zero as vx is constant over the element*/
+        for(int i=0;i<3;i++) epsilonvx[i]=0;
+}
+/*}}}*/
+/*FUNCTION DoubleInput::GetVyStrainRate2d(double* epsilonvy,double* xyz_list, GaussTria* gauss){{{1*/
+void DoubleInput::GetVyStrainRate2d(double* epsilonvy,double* xyz_list, GaussTria* gauss){
+        /*Epsilon is zero as vy is constant over the element*/
+        for(int i=0;i<3;i++) epsilonvy[i]=0;
+}
+/*}}}*/
+/*FUNCTION DoubleInput::GetVxStrainRate3d(double* epsilonvx,double* xyz_list, GaussPenta* gauss){{{1*/
+void DoubleInput::GetVxStrainRate3d(double* epsilonvx,double* xyz_list, GaussPenta* gauss){
+        /*Epsilon is zero as vx is constant over the element*/
+        for(int i=0;i<6;i++) epsilonvx[i]=0;
+}
+/*}}}*/
+/*FUNCTION DoubleInput::GetVyStrainRate3d(double* epsilonvy,double* xyz_list, GaussPenta* gauss){{{1*/
+void DoubleInput::GetVyStrainRate3d(double* epsilonvy,double* xyz_list, GaussPenta* gauss){
+        /*Epsilon is zero as vy is constant over the element*/
+        for(int i=0;i<6;i++) epsilonvy[i]=0;
+}
+/*}}}*/
+/*FUNCTION DoubleInput::GetVzStrainRate3d(double* epsilonvz,double* xyz_list, GaussPenta* gauss){{{1*/
+void DoubleInput::GetVzStrainRate3d(double* epsilonvz,double* xyz_list, GaussPenta* gauss){
+        /*Epsilon is zero as vz is constant over the element*/
+        for(int i=0;i<6;i++) epsilonvz[i]=0;
+}
+/*}}}*/
+/*FUNCTION DoubleInput::GetVxStrainRate3dPattyn(double* epsilonvx,double* xyz_list, GaussPenta* gauss){{{1*/
+void DoubleInput::GetVxStrainRate3dPattyn(double* epsilonvx,double* xyz_list, GaussPenta* gauss){
+        /*Epsilon is zero as vx is constant over the element*/
+        for(int i=0;i<5;i++) epsilonvx[i]=0;
+}
+/*}}}*/
+/*FUNCTION DoubleInput::GetVyStrainRate3dPattyn(double* epsilonvy,double* xyz_list, GaussPenta* gauss){{{1*/
+void DoubleInput::GetVyStrainRate3dPattyn(double* epsilonvy,double* xyz_list, GaussPenta* gauss){
+        /*Epsilon is zero as vy is constant over the element*/
+        for(int i=0;i<5;i++) epsilonvy[i]=0;
+}
+/*}}}*/
 /*FUNCTION DoubleInput::ChangeEnum{{{1*/
 void DoubleInput::ChangeEnum(int newenumtype){

issm/trunk/src/c/objects/Inputs/DoubleInput.h

@@ -57 +57 @@
                 void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussPenta* gauss);
                 void GetParameterAverage(double* pvalue);
-                void GetVxStrainRate2d(double* epsilonvx,double* xyz_list, GaussTria* gauss){_error_("not implemented yet");};
-                void GetVyStrainRate2d(double* epsilonvy,double* xyz_list, GaussTria* gauss){_error_("not implemented yet");};
-                void GetVxStrainRate3d(double* epsilonvx,double* xyz_list, GaussPenta* gauss){_error_("not implemented yet");};
-                void GetVyStrainRate3d(double* epsilonvy,double* xyz_list, GaussPenta* gauss){_error_("not implemented yet");};
-                void GetVzStrainRate3d(double* epsilonvz,double* xyz_list, GaussPenta* gauss){_error_("not implemented yet");};
-                void GetVxStrainRate3dPattyn(double* epsilonvx,double* xyz_list, GaussPenta* gauss){_error_("not implemented yet");};
-                void GetVyStrainRate3dPattyn(double* epsilonvy,double* xyz_list, GaussPenta* gauss){_error_("not implemented yet");};
+                void GetVxStrainRate2d(double* epsilonvx,double* xyz_list, GaussTria* gauss);
+                void GetVyStrainRate2d(double* epsilonvy,double* xyz_list, GaussTria* gauss);
+                void GetVxStrainRate3d(double* epsilonvx,double* xyz_list, GaussPenta* gauss);
+                void GetVyStrainRate3d(double* epsilonvy,double* xyz_list, GaussPenta* gauss);
+                void GetVzStrainRate3d(double* epsilonvz,double* xyz_list, GaussPenta* gauss);
+                void GetVxStrainRate3dPattyn(double* epsilonvx,double* xyz_list, GaussPenta* gauss);
+                void GetVyStrainRate3dPattyn(double* epsilonvy,double* xyz_list, GaussPenta* gauss);
                 void ChangeEnum(int newenumtype);
                 void SquareMin(double* psquaremin, bool process_units,Parameters* parameters);