source: issm/oecreview/Archive/16554-17801/ISSM-16774-16775.diff@ 17802

../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

@@ -1171 +1171 @@
         for(i=0;i<numdofFSp;i++) FSvalues[numdofFSv+i]=solution[doflistFSp[i]];
 
         /*Transform solution in Cartesian Space*/
-        element->TransformSolutionCoord(FSvalues,cs_list);
-        element->TransformSolutionCoord(HOvalues,XYEnum);
+        element->TransformSolutionCoord(FSvalues,2*numnodes,cs_list);
+        element->TransformSolutionCoord(HOvalues,numnodes,XYEnum);
 
         /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
         for(i=0;i<numnodes;i++){

../trunk-jpl/src/c/classes/Elements/Element.h

@@ -52 +52 @@
                 virtual IssmDouble GetMaterialParameter(int enum_in)=0;
                 virtual void   TransformSolutionCoord(IssmDouble* values,int transformenum)=0;
                 virtual void   TransformSolutionCoord(IssmDouble* values,int* transformenum_list)=0;
+                virtual void   TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum)=0;
+                virtual void   TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list)=0;
                 virtual Element* GetBasalElement(void)=0;
                 virtual void    GetDofList(int** pdoflist,int approximation_enum,int setenum)=0;
                 virtual void    GetDofListVelocity(int** pdoflist,int setenum)=0;

../trunk-jpl/src/c/classes/Elements/Tria.h

@@ -272 +272 @@
                 void             SurfaceNormal(IssmDouble* surface_normal, IssmDouble xyz_list[3][3]);
                 void           TransformSolutionCoord(IssmDouble* values,int transformenum);
                 void           TransformSolutionCoord(IssmDouble* values,int* transformenum_list);
+                void           TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){_error_("not implemented yet");};
+                void           TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){_error_("not implemented yet");};
                 #ifdef _HAVE_STRESSBALANCE_
                 ElementMatrix* CreateKMatrixStressbalanceSSA(void);
                 ElementMatrix* CreateKMatrixStressbalanceSSAViscous(void);

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

@@ -2015 +2015 @@
                         else if (original_input->ObjectEnum()==DoubleInputEnum)
                          total_integrated_input=new DoubleInput(average_enum_type,0.0);
                         else{
-                                printf("id %i \n",this->id);
-                                printf("enum %i \n",average_enum_type);
                          _error_("object " << EnumToStringx(original_input->ObjectEnum()) << " not supported yet");
                         }
                 }
@@ -3328 +3326 @@
 
         return dt;
 }/*}}}*/
-/*FUNCTION Penta::TransformSolutionCoord{{{*/
+/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int transformenum){{{*/
 void Penta::TransformSolutionCoord(IssmDouble* values,int transformenum){
 
         ::TransformSolutionCoord(values,this->nodes,this->NumberofNodes(),transformenum);
 
 }
 /*}}}*/
-/*FUNCTION Penta::TransformSolutionCoord{{{*/
+/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int* transformenum_list){{{*/
 void Penta::TransformSolutionCoord(IssmDouble* values,int* transformenum_list){
 
         ::TransformSolutionCoord(values,this->nodes,this->NumberofNodes(),transformenum_list);
 
 }
 /*}}}*/
+/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){{{*/
+void Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){
+
+        ::TransformSolutionCoord(values,this->nodes,numnodes,transformenum);
+
+}
+/*}}}*/
+/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){{{*/
+void Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){
+
+        ::TransformSolutionCoord(values,this->nodes,numnodes,transformenum_list);
+
+}
+/*}}}*/
 /*FUNCTION Penta::Update {{{*/
 void Penta::Update(int index,IoModel* iomodel,int analysis_counter,int analysis_type,int finiteelement_type){ 
 

../trunk-jpl/src/c/classes/Elements/Penta.h

@@ -248 +248 @@
                 void             SurfaceNormal(IssmDouble* surface_normal, IssmDouble xyz_list[3][3]);
                 void           TransformSolutionCoord(IssmDouble* values,int transformenum);
                 void           TransformSolutionCoord(IssmDouble* values,int* transformenum_list);
+                void           TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum);      /*Tiling only*/
+                void           TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list);/*Tiling only*/
 
                 #ifdef _HAVE_STRESSBALANCE_
                 ElementMatrix* CreateKMatrixCouplingSSAHO(void);

../trunk-jpl/src/c/classes/Elements/Seg.h

@@ -153 +153 @@
                 IssmDouble  TimeAdapt(){_error_("not implemented yet");};
                 void        TransformSolutionCoord(IssmDouble* values,int transformenum){_error_("not implemented yet");};
                 void        TransformSolutionCoord(IssmDouble* values,int* transformenum_list){_error_("not implemented yet");};
+                void        TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){_error_("not implemented yet");};
+                void        TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){_error_("not implemented yet");};
                 void UpdateConstraintsExtrudeFromBase(){_error_("not implemented");};
                 void UpdateConstraintsExtrudeFromTop(){_error_("not implemented");};
                 void UpdateConstraintsL2ProjectionEPL(){_error_("not implemented");};