Changeset 18240

Timestamp:

07/14/14 11:12:10 (11 years ago)

Author:

seroussi

Message:

NEW: analytical solution for Stokes test

File:

• issm/trunk-jpl/src/c/shared/FSanalyticals/fsanalyticals.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/shared/FSanalyticals/fsanalyticals.cpp

@@ -45 +45 @@
                 case 10:
                         return 4*pow(PI, 2)*pow(p, 2)*(cos(PI*p*x) - 1)*sin(PI*p*y)*sin(PI*p*z) + 2*pow(PI, 2)*pow(p, 2)*sin(PI*p*y)*sin(PI*p*z)*cos(PI*p*x) + PI*q*pow(z, a)*sin(PI*q*y)*cos(PI*q*x);
+                case 12:
+                        z=y;
+                        return 32.0/5.0*pow(PI, 2)*sin(8*PI*z);
                 case 21:
                         z=y;
@@ -87 +90 @@
                 case 10:
                         return -2*pow(PI, 2)*pow(p, 2)*(cos(PI*p*y) - 1)*sin(PI*p*x)*sin(PI*p*z) - pow(PI, 2)*pow(p, 2)*sin(PI*p*x)*sin(PI*p*z)*cos(PI*p*y) + PI*q*pow(z, a)*sin(PI*q*x)*cos(PI*q*y);
+                case 12:
+                        z=y;
+                        return 32.0/5.0*pow(PI, 2)*sin(8*PI*x) - 100;
                 case 21:
                         z=y;
@@ -106 +112 @@
 
         switch(testid){
-                case 1: case 2: case 7: case 8: case 11: case 18: case 21:
+                case 1: case 2: case 7: case 8: case 11:  case 12: case 18: case 21:
                         return 0.;
                 case 3: 
@@ -132 +138 @@
 
         switch(testid){
-                case 1: case 2: case 3: case 4: case 5: case 6: case 7: case 8: case 9: case 10: case 21:
+                case 1: case 2: case 3: case 4: case 5: case 6: case 7: case 8: case 9: case 10: case 12: case 21:
                         return 0.;
                 case 11: