Changeset 10393

Timestamp:

11/01/11 07:27:08 (13 years ago)

Author:

Mathieu Morlighem

Message:

Forgot to commit 4d transformation

File:

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

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

@@ -15 +15 @@
         /*Some checks*/
         _assert_(numnodes && nodes);
-        if(dimension!=2 && dimension!=3) _error_("only 2d and 3d coordinate systems are supported");
+        if(dimension!=2 && dimension!=3 && dimension!=4) _error_("only 2d and 3d coordinate systems are supported");
 
         /*Allocate and initialize transform matrix*/
@@ -52 +52 @@
                         transform[(dimension*numnodes)*(i*dimension+2) + i*dimension+2] = coord_system[2][2];
                 }
+                else if(dimension==4){
+                        /*Only the first 3 coordinates are changed (x,y,z), leave the others (P) unchanged*/
+                        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+0) + i*dimension+2] = coord_system[0][2];
+                        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];
+                        transform[(dimension*numnodes)*(i*dimension+1) + i*dimension+2] = coord_system[1][2];
+                        transform[(dimension*numnodes)*(i*dimension+2) + i*dimension+0] = coord_system[2][0];
+                        transform[(dimension*numnodes)*(i*dimension+2) + i*dimension+1] = coord_system[2][1];
+                        transform[(dimension*numnodes)*(i*dimension+2) + i*dimension+2] = coord_system[2][2];
+                        transform[(dimension*numnodes)*(i*dimension+3) + i*dimension+3] = 1.0;
+                }
         }