diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-04-21 12:26:58 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-04-21 12:26:58 +0000 |
commit | 13c293dd852f79fc22610838a95671ae4109e4c7 (patch) | |
tree | 4703d53ee66070b9f052a23ea1ce25e7636c7c43 /src/patch | |
parent | 54b313aa1e685327a397e32792a2570402ea004a (diff) |
remove unused fn args, eg phi = phi(x,y), so there's no need for a z
argument, and no need to compute d/dz or d^2/dwdz for any w
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@537 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/patch')
-rw-r--r-- | src/patch/coords.cc | 46 | ||||
-rw-r--r-- | src/patch/coords.hh | 42 |
2 files changed, 32 insertions, 56 deletions
diff --git a/src/patch/coords.cc b/src/patch/coords.cc index 0def972..3f48507 100644 --- a/src/patch/coords.cc +++ b/src/patch/coords.cc @@ -260,17 +260,14 @@ fp phi_of_xy(fp x, fp y) { return arctan_xy(x,y); } // partial {mu,nu,phi} / partial {x,y,z} // as computed by the maple file "coord_derivs.{maple,out}" in this directory // -fp partial_mu_wrt_x(fp x, fp y, fp z) { return 0.0; } -fp partial_mu_wrt_y(fp x, fp y, fp z) { return z / (y*y + z*z); } -fp partial_mu_wrt_z(fp x, fp y, fp z) { return -y / (y*y + z*z); } +fp partial_mu_wrt_y(fp y, fp z) { return z / (y*y + z*z); } +fp partial_mu_wrt_z(fp y, fp z) { return -y / (y*y + z*z); } -fp partial_nu_wrt_x(fp x, fp y, fp z) { return z / (x*x + z*z); } -fp partial_nu_wrt_y(fp x, fp y, fp z) { return 0.0; } -fp partial_nu_wrt_z(fp x, fp y, fp z) { return -x / (x*x + z*z); } +fp partial_nu_wrt_x(fp x, fp z) { return z / (x*x + z*z); } +fp partial_nu_wrt_z(fp x, fp z) { return -x / (x*x + z*z); } -fp partial_phi_wrt_x(fp x, fp y, fp z) { return -y / (x*x + y*y); } -fp partial_phi_wrt_y(fp x, fp y, fp z) { return x / (x*x + y*y); } -fp partial_phi_wrt_z(fp x, fp y, fp z) { return 0.0; } +fp partial_phi_wrt_x(fp x, fp y) { return -y / (x*x + y*y); } +fp partial_phi_wrt_y(fp x, fp y) { return x / (x*x + y*y); } //****************************************************************************** @@ -279,26 +276,17 @@ fp partial_phi_wrt_z(fp x, fp y, fp z) { return 0.0; } // partial {mu,nu,phi} / partial {xx,xy,xz,yy,yz,zz} // as computed by the maple file "coord_derivs.{maple,out}" in this directory // -fp partial2_mu_wrt_xx(fp x, fp y, fp z) {return 0.0;} -fp partial2_mu_wrt_xy(fp x, fp y, fp z) {return 0.0;} -fp partial2_mu_wrt_xz(fp x, fp y, fp z) {return 0.0;} -fp partial2_mu_wrt_yy(fp x, fp y, fp z) {return -2.0*y*z / pow2(y*y + z*z);} -fp partial2_mu_wrt_yz(fp x, fp y, fp z) {return (y*y - z*z) / pow2(y*y + z*z);} -fp partial2_mu_wrt_zz(fp x, fp y, fp z) {return 2.0*y*z / pow2(y*y + z*z);} - -fp partial2_nu_wrt_xx(fp x, fp y, fp z) {return -2.0*x*z / pow2(x*x + z*z);} -fp partial2_nu_wrt_xy(fp x, fp y, fp z) {return 0.0;} -fp partial2_nu_wrt_xz(fp x, fp y, fp z) {return (x*x - z*z) / pow2(x*x + z*z);} -fp partial2_nu_wrt_yy(fp x, fp y, fp z) {return 0.0;} -fp partial2_nu_wrt_yz(fp x, fp y, fp z) {return 0.0;} -fp partial2_nu_wrt_zz(fp x, fp y, fp z) {return 2.0*x*z / pow2(x*x + z*z);} - -fp partial2_phi_wrt_xx(fp x, fp y, fp z) {return 2.0*x*y / pow2(x*x + y*y);} -fp partial2_phi_wrt_xy(fp x, fp y, fp z) {return (y*y - x*x) / pow2(x*x + y*y);} -fp partial2_phi_wrt_xz(fp x, fp y, fp z) {return 0.0;} -fp partial2_phi_wrt_yy(fp x, fp y, fp z) {return -2.0*x*y / pow2(x*x + y*y);} -fp partial2_phi_wrt_yz(fp x, fp y, fp z) {return 0.0;} -fp partial2_phi_wrt_zz(fp x, fp y, fp z) {return 0.0;} +fp partial2_mu_wrt_yy(fp y, fp z) {return -2.0*y*z / pow2(y*y + z*z);} +fp partial2_mu_wrt_yz(fp y, fp z) {return (y*y - z*z) / pow2(y*y + z*z);} +fp partial2_mu_wrt_zz(fp y, fp z) {return 2.0*y*z / pow2(y*y + z*z);} + +fp partial2_nu_wrt_xx(fp x, fp z) {return -2.0*x*z / pow2(x*x + z*z);} +fp partial2_nu_wrt_xz(fp x, fp z) {return (x*x - z*z) / pow2(x*x + z*z);} +fp partial2_nu_wrt_zz(fp x, fp z) {return 2.0*x*z / pow2(x*x + z*z);} + +fp partial2_phi_wrt_xx(fp x, fp y) {return 2.0*x*y / pow2(x*x + y*y);} +fp partial2_phi_wrt_xy(fp x, fp y) {return (y*y - x*x) / pow2(x*x + y*y);} +fp partial2_phi_wrt_yy(fp x, fp y) {return -2.0*x*y / pow2(x*x + y*y);} //****************************************************************************** //****************************************************************************** diff --git a/src/patch/coords.hh b/src/patch/coords.hh index 12a6baa..b2c460a 100644 --- a/src/patch/coords.hh +++ b/src/patch/coords.hh @@ -138,35 +138,23 @@ fp nu_of_mu_phi(fp mu, fp phi); fp mu_of_nu_phi(fp nu, fp phi); // partial {mu,nu,phi} / partial {x,y,z} -fp partial_mu_wrt_x(fp x, fp y, fp z); -fp partial_mu_wrt_y(fp x, fp y, fp z); -fp partial_mu_wrt_z(fp x, fp y, fp z); -fp partial_nu_wrt_x(fp x, fp y, fp z); -fp partial_nu_wrt_y(fp x, fp y, fp z); -fp partial_nu_wrt_z(fp x, fp y, fp z); -fp partial_phi_wrt_x(fp x, fp y, fp z); -fp partial_phi_wrt_y(fp x, fp y, fp z); -fp partial_phi_wrt_z(fp x, fp y, fp z); +fp partial_mu_wrt_y(fp y, fp z); +fp partial_mu_wrt_z(fp y, fp z); +fp partial_nu_wrt_x(fp x, fp z); +fp partial_nu_wrt_z(fp x, fp z); +fp partial_phi_wrt_x(fp x, fp y); +fp partial_phi_wrt_y(fp x, fp y); // partial^2 {mu,nu,phi} / partial {x,y,z}{x,y,z} -fp partial2_mu_wrt_xx(fp x, fp y, fp z); -fp partial2_mu_wrt_xy(fp x, fp y, fp z); -fp partial2_mu_wrt_xz(fp x, fp y, fp z); -fp partial2_mu_wrt_yy(fp x, fp y, fp z); -fp partial2_mu_wrt_yz(fp x, fp y, fp z); -fp partial2_mu_wrt_zz(fp x, fp y, fp z); -fp partial2_nu_wrt_xx(fp x, fp y, fp z); -fp partial2_nu_wrt_xy(fp x, fp y, fp z); -fp partial2_nu_wrt_xz(fp x, fp y, fp z); -fp partial2_nu_wrt_yy(fp x, fp y, fp z); -fp partial2_nu_wrt_yz(fp x, fp y, fp z); -fp partial2_nu_wrt_zz(fp x, fp y, fp z); -fp partial2_phi_wrt_xx(fp x, fp y, fp z); -fp partial2_phi_wrt_xy(fp x, fp y, fp z); -fp partial2_phi_wrt_xz(fp x, fp y, fp z); -fp partial2_phi_wrt_yy(fp x, fp y, fp z); -fp partial2_phi_wrt_yz(fp x, fp y, fp z); -fp partial2_phi_wrt_zz(fp x, fp y, fp z); +fp partial2_mu_wrt_yy(fp y, fp z); +fp partial2_mu_wrt_yz(fp y, fp z); +fp partial2_mu_wrt_zz(fp y, fp z); +fp partial2_nu_wrt_xx(fp x, fp z); +fp partial2_nu_wrt_yy(fp x, fp z); +fp partial2_nu_wrt_zz(fp x, fp z); +fp partial2_phi_wrt_xx(fp x, fp y); +fp partial2_phi_wrt_xy(fp x, fp y); +fp partial2_phi_wrt_yy(fp x, fp y); // usual polar spherical (r,theta,phi) <--> (x,y,z) void xyz_of_r_theta_phi(fp r, fp theta, fp phi, fp& x, fp& y, fp& z); |