source: issm/oecreview/Archive/16554-17801/ISSM-17140-17141.diff@ 17802

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

@@ -5 +5 @@
 #include "../modules/modules.h"
 #include "../solutionsequences/solutionsequences.h"
 
-//#define FSANALYTICAL 1
+//#define FSANALYTICAL 2
 
 /*Model processing*/
 int  StressbalanceAnalysis::DofsPerNode(int** pdoftype,int meshtype,int approximation){/*{{{*/
@@ -2761 +2761 @@
                 y_coord=element->GetYcoord(gauss);
                 z_coord=element->GetZcoord(gauss);
 
-                forcex=fx1(x_coord,y_coord,z_coord);
-                forcey=fy1(x_coord,y_coord,z_coord);
-                forcez=fz1(x_coord,y_coord,z_coord);
+                #if FSANALYTICAL == 1
+                        forcex=fx1(x_coord,y_coord,z_coord);
+                        forcey=fy1(x_coord,y_coord,z_coord);
+                        forcez=fz1(x_coord,y_coord,z_coord);
+                #elseif FSANALYTICAL == 2
+                        forcex=fx2(x_coord,y_coord,z_coord);
+                        forcey=fy2(x_coord,y_coord,z_coord);
+                        forcez=fz2(x_coord,y_coord,z_coord);
+                #else 
+                        _error_("FS analytical not implemented yet");
+                #endif
 
                 for(i=0;i<vnumnodes;i++){
                         pe->values[i*dim+0] += forcex *Jdet*gauss->weight*vbasis[i];
@@ -3154 +3162 @@
          *                                      [ phi ]
          *
          *      In 3d
-         *         Bvi=[ dh/dx          0             0      ]
-         *                                      [   0           dh/dy           0      ]
-         *                                      [   0             0           dh/dz    ]
-         *                                      [ 1/2*dh/dy    1/2*dh/dx        0      ]
-         *                                      [ 1/2*dh/dz       0         1/2*dh/dx  ]
-         *                                      [   0          1/2*dh/dz    1/2*dh/dy  ]
-         *                                      [   0             0             0      ]
-         *                                      [ dh/dx         dh/dy         dh/dz    ]
+         *         Bvi=[ dh/dx     0        0    ]
+         *                                      [   0      dh/dy      0    ]
+         *                                      [   0        0      dh/dz  ]
+         *                                      [ dh/dy    dh/dx      0    ]
+         *                                      [ dh/dz      0      dh/dx  ]
+         *                                      [   0      dh/dz    dh/dy  ]
+         *                                      [ dh/dx    dh/dy    dh/dz  ]
+         *                                      [   0        0        0    ]
          *
          *         Bpi=[ 0 ]
          *                                      [ 0 ]
@@ -3169 +3177 @@
          *                                      [ 0 ]
          *                                      [ 0 ]
          *                                      [ 0 ]
+         *                                      [ 0 ]
          *                                      [ h ]
-         *                                      [ 0 ]
          *      where phi is the finiteelement function for node i.
          *      In 3d:
          */