Changeset 17388

Timestamp:

03/05/14 15:19:14 (11 years ago)

Author:

seroussi

Message:

NEW: another analytical solution

File:

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

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

@@ -40 +40 @@
                 case 8: 
                         return 1.0;
+                case 9:
+                        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) + 2*x*z;
 
                 case 101: 
@@ -72 +74 @@
                 case 8: 
                         return 1.0;
+                case 9: 
+                        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) + 2*y*z;
 
                 case 101: 
@@ -95 +99 @@
                 case 6: 
                         return - 2. * PI*PI * p*p * (cos(p*PI*z)-1) * sin(p*PI*x) * sin(p*PI*y)-PI*PI * p*p *  sin(p*PI*x)    * sin(p*PI*y) * cos(p*PI*z)+ q * PI * sin(q*PI*x) * sin(q*PI*y) * cos(q*PI*z);
+                case 9:
+                        return -2*pow(PI, 3)*pow(p, 3)*(cos(PI*p*z)/(PI*p) - 1/(PI*p))*sin(PI*p*x)*sin(PI*p*y) - pow(PI, 2)*pow(p, 2)*sin(PI*p*x)*sin(PI*p*y)*cos(PI*p*z) + pow(x, 2) + pow(y, 2);
                 default:
                         _error_("FS analytical solution"<<testid<<" not implemented yet");