Changeset 15974

Timestamp:

08/27/13 13:46:28 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: xyz_list_tria was not well created, the z coordinate was never updated

File:

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

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

@@ -4664 +4664 @@
         /*Retrieve all inputs and parameters*/
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);
         rho_ice=matpar->GetRhoIce();
@@ -4916 +4916 @@
         /*Retrieve all inputs and parameters*/
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);
         rho_ice=matpar->GetRhoIce();
@@ -6211 +6211 @@
         GradientIndexing(&vertexpidlist[0],control_index);
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         Input* adjointx_input=inputs->GetInput(AdjointxEnum);               _assert_(adjointx_input);
         Input* adjointy_input=inputs->GetInput(AdjointyEnum);               _assert_(adjointy_input);
@@ -6283 +6283 @@
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         GradientIndexing(&vertexpidlist[0],control_index);
         Input* drag_input    =inputs->GetInput(FrictionCoefficientEnum); _assert_(drag_input);
@@ -7193 +7193 @@
         /*retrieve inputs :*/
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);
         Input* vx_input=inputs->GetInput(VxEnum);           _assert_(vx_input);
@@ -8048 +8048 @@
         /*Retrieve all inputs and parameters*/
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<2;j++) xyz_list_tria[i][j]=xyz_list[i][j];
+        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
         parameters->FindParam(&migration_style,GroundinglineMigrationEnum);
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);