diff options
| author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2001-09-15 18:30:34 +0000 |
|---|---|---|
| committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2001-09-15 18:30:34 +0000 |
| commit | 7994e293c56120881769ab7ba9b3aa836523cbc8 (patch) | |
| tree | f28885d2c4cb90d2ef0c39e25f8e54106616a5be /src/patch/test_patch_system.cc | |
| parent | 86a9bf47894d2134e54807447dbf9d139f9bc4e2 (diff) | |
oops, needed to change the C++ code to actually *use* the new Maple output!
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@347 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/patch/test_patch_system.cc')
| -rw-r--r-- | src/patch/test_patch_system.cc | 41 |
1 files changed, 26 insertions, 15 deletions
diff --git a/src/patch/test_patch_system.cc b/src/patch/test_patch_system.cc index b8e985a..8c27dbc 100644 --- a/src/patch/test_patch_system.cc +++ b/src/patch/test_patch_system.cc @@ -645,30 +645,41 @@ if (which_derivs & which_deriv_fn) if (which_derivs & which_deriv_rho) then sum += deriv_weight_rho - * - cos(rho)*exp(sin(rho))*sinh(-1.0+pow(cos(sigma),2.0)) - / cosh(-1.0+pow(cos(sigma),2.0)); + * 0.138E1*cos(0.138E1*rho)*exp(sin(0.138E1*rho)) + *tanh(0.17+0.83*pow(sin(sigma),2.0)); if (which_derivs & which_deriv_sigma) then sum += deriv_weight_sigma - * 2.0*exp(sin(rho))*sin(sigma)*cos(sigma) - / pow(cosh(-1.0+pow(cos(sigma),2.0)),2.0); + * 0.166E1*exp(sin(0.138E1*rho)) + *(1.0-pow(tanh(0.17+0.83*pow(sin(sigma),2.0)),2.0)) + *sin(sigma)*cos(sigma); if (which_derivs & which_deriv_rho_rho) then sum += deriv_weight_rho_rho - * exp(sin(rho))*sinh(-1.0+pow(cos(sigma),2.0)) - * (sin(rho)-pow(cos(rho),2.0)) - / cosh(-1.0+pow(cos(sigma),2.0)); + * ( + -0.19044E1*sin(0.138E1*rho)*exp(sin(0.138E1*rho)) + *tanh(0.17+0.83*pow(sin(sigma),2.0)) + +0.19044E1*pow(cos(0.138E1*rho),2.0)*exp(sin(0.138E1*rho)) + *tanh(0.17+0.83*pow(sin(sigma),2.0)) + ); if (which_derivs & which_deriv_rho_sigma) then sum += deriv_weight_rho_sigma - * 2.0*cos(rho)*exp(sin(rho))*sin(sigma)*cos(sigma) - / pow(cosh(-1.0+pow(cos(sigma),2.0)),2.0); + * 0.22908E1*cos(0.138E1*rho)*exp(sin(0.138E1*rho)) + *(1.0-pow(tanh(0.17+0.83*pow(sin(sigma),2.0)),2.0)) + *sin(sigma)*cos(sigma); if (which_derivs & which_deriv_sigma_sigma) then sum += deriv_weight_sigma_sigma - * -2.0*exp(sin(rho)) - * ( - 4.0*sinh(-1.0+pow(cos(sigma),2.0))*pow(cos(sigma),2.0) - + 4.0*sinh(-1.0+pow(cos(sigma),2.0))*pow(cos(sigma),4.0) - - 2.0*pow(cos(sigma),2.0)*cosh(-1.0+pow(cos(sigma),2.0)) - + cosh(-1.0+pow(cos(sigma),2.0)) ) - / pow(cosh(-1.0+pow(cos(sigma),2.0)),3.0); + * ( + -0.55112E1*exp(sin(0.138E1*rho)) + *tanh(0.17+0.83*pow(sin(sigma),2.0)) + *(1.0-pow(tanh(0.17+0.83*pow(sin(sigma),2.0)),2.0)) + *pow(sin(sigma),2.0)*pow(cos(sigma),2.0) + +0.166E1*exp(sin(0.138E1*rho)) + *(1.0-pow(tanh(0.17+0.83*pow(sin(sigma),2.0)),2.0)) + *pow(cos(sigma),2.0) + -0.166E1*exp(sin(0.138E1*rho)) + *(1.0-pow(tanh(0.17+0.83*pow(sin(sigma),2.0)),2.0)) + *pow(sin(sigma),2.0) + ); return sum; } |