[cig-commits] r14077 - cs/cigma/trunk/src

[luis at geodynamics.org]

Author: luis
Date: 2009-02-18 08:14:25 -0800 (Wed, 18 Feb 2009)
New Revision: 14077

Modified:
   cs/cigma/trunk/src/Cell.cpp
Log:
Fixed jacobian. Z(i) should be zero for 2d cells.

Modified: cs/cigma/trunk/src/Cell.cpp
===================================================================
--- cs/cigma/trunk/src/Cell.cpp	2009-02-18 16:14:24 UTC (rev 14076)
+++ cs/cigma/trunk/src/Cell.cpp	2009-02-18 16:14:25 UTC (rev 14077)
@@ -234,22 +234,22 @@
 
                 jac[0][0] += X(i) * s[0];
                 jac[0][1] += Y(i) * s[0];
-                jac[0][2] += Z(i) * s[0];
+                //jac[0][2] += Z(i) * s[0];
 
                 jac[1][0] += X(i) * s[1];
                 jac[1][1] += Y(i) * s[1];
-                jac[1][2] += Z(i) * s[1];
+                //jac[1][2] += Z(i) * s[1];
             }
             {
                 double a[3], b[3], c[3];
 
                 a[0]= X(1) - X(0);
                 a[1]= Y(1) - Y(0);
-                a[2]= Z(1) - Z(0);
+                a[2]= 0.0; //Z(1) - Z(0);
 
                 b[0]= X(2) - X(0);
                 b[1]= Y(2) - Y(0);
-                b[2]= Z(2) - Z(0);
+                b[2]= 0.0; //Z(2) - Z(0);
 
                 prodve(a, b, c);