Changeset 10388

Timestamp:

10/31/11 16:17:54 (13 years ago)

Author:

Mathieu Morlighem

Message:

Forgot to divide by norm

File:

• issm/trunk/src/c/shared/Elements/CoordinateSystemTransform.cpp (modified) (2 diffs)

issm/trunk/src/c/shared/Elements/CoordinateSystemTransform.cpp

@@ -3 +3 @@
  */ 
 #include "./elements.h"
+#include <math.h>
 
 void CoordinateSystemTransform(double** ptransform,Node** nodes,int numnodes,int dimension){
 
         int     i;
+        double  x_norm,y_norm;
         double *transform         = NULL;
         double *values            = NULL;
@@ -31 +33 @@
                 nodes[i]->GetCoordinateSystem(&coord_system[0][0]);
                 if(dimension==2){
-                        transform[(dimension*numnodes)*(i*dimension+0) + i*dimension+0] = coord_system[0][0];
-                        transform[(dimension*numnodes)*(i*dimension+0) + i*dimension+1] = coord_system[0][1];
-                        transform[(dimension*numnodes)*(i*dimension+1) + i*dimension+0] = coord_system[1][0];
-                        transform[(dimension*numnodes)*(i*dimension+1) + i*dimension+1] = coord_system[1][1];
+                        /*We remove the z component, we need to renormalize x and y*/
+                        x_norm = sqrt( coord_system[0][0]*coord_system[0][0] + coord_system[1][0]*coord_system[1][0]);
+                        y_norm = sqrt( coord_system[0][1]*coord_system[0][1] + coord_system[1][1]*coord_system[1][1]);
+                        transform[(dimension*numnodes)*(i*dimension+0) + i*dimension+0] = coord_system[0][0]/x_norm;
+                        transform[(dimension*numnodes)*(i*dimension+0) + i*dimension+1] = coord_system[0][1]/y_norm;
+                        transform[(dimension*numnodes)*(i*dimension+1) + i*dimension+0] = coord_system[1][0]/x_norm;
+                        transform[(dimension*numnodes)*(i*dimension+1) + i*dimension+1] = coord_system[1][1]/y_norm;
                 }
                 else if(dimension==3){