Changeset 17295
- Timestamp:
- 02/18/14 14:50:58 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/shared/FSanalyticals/FSanalytical3.cpp
r17294 r17295 46 46 IssmDouble fx5(IssmDouble x,IssmDouble y,IssmDouble z){ 47 47 IssmDouble functionx; 48 intp=2;49 intq=2;48 IssmDouble p=2; 49 IssmDouble q=2; 50 50 51 51 functionx = 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*sin(PI*q*y)*sin(PI*q*z)*cos(PI*q*x); … … 55 55 IssmDouble fy5(IssmDouble x,IssmDouble y,IssmDouble z){ 56 56 IssmDouble functiony; 57 intp=2;58 intq=2;57 IssmDouble p=2; 58 IssmDouble q=2; 59 59 60 60 functiony = -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*sin(PI*q*x)*sin(PI*q*z)*cos(PI*q*y); … … 64 64 IssmDouble fz5(IssmDouble x,IssmDouble y,IssmDouble z){ 65 65 IssmDouble functionz; 66 intp=2;67 intq=2;66 IssmDouble p=2; 67 IssmDouble q=2; 68 68 69 69 functionz = -2*pow(PI, 2)*pow(p, 2)*(cos(PI*p*z) - 1)*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) + PI*q*sin(PI*q*x)*sin(PI*q*y)*cos(PI*q*z);
Note:
See TracChangeset
for help on using the changeset viewer.