- add partial_[xyz]_wrt_r to partial_xyz_wrt_r_rho_sigma() & friends

- add   dh/drho and dh/dsigma terms to Jacobian computation for
  surface integrals ==> now works to very high accuracy even
  when Schw BH is offset from origin


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@731 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 15 insertions, 12 deletions

diff --git a/src/patch/coords.hh b/src/patch/coords.hh
index 787c203..50bd989 100644
--- a/src/patch/coords.hh
+++ b/src/patch/coords.hh
@@ -147,18 +147,21 @@ fp nu_of_mu_phi(fp mu, fp phi);
 fp mu_of_nu_phi(fp nu, fp phi);
 
 // partial {x,y,z} / partial {mu,nu,phi}
-void partial_xyz_wrt_mu_nu(fp r, fp mu, fp nu,
-			   fp& partial_x_wrt_mu, fp& partial_x_wrt_nu,
-			   fp& partial_y_wrt_mu, fp& partial_y_wrt_nu,
-			   fp& partial_z_wrt_mu, fp& partial_z_wrt_nu);
-void partial_xyz_wrt_mu_phi(fp r, fp mu, fp phi,
-			    fp& partial_x_wrt_mu, fp& partial_x_wrt_phi,
-			    fp& partial_y_wrt_mu, fp& partial_y_wrt_phi,
-			    fp& partial_z_wrt_mu, fp& partial_z_wrt_phi);
-void partial_xyz_wrt_nu_phi(fp r, fp nu, fp phi,
-			    fp& partial_x_wrt_nu, fp& partial_x_wrt_phi,
-			    fp& partial_y_wrt_nu, fp& partial_y_wrt_phi,
-			    fp& partial_z_wrt_nu, fp& partial_z_wrt_phi);
+void partial_xyz_wrt_r_mu_nu
+	(fp r, fp mu, fp nu,
+	 fp& partial_x_wrt_r, fp& partial_x_wrt_mu, fp& partial_x_wrt_nu,
+	 fp& partial_y_wrt_r, fp& partial_y_wrt_mu, fp& partial_y_wrt_nu,
+	 fp& partial_y_wrt_r, fp& partial_z_wrt_mu, fp& partial_z_wrt_nu);
+void partial_xyz_wrt_r_mu_phi
+	(fp r, fp mu, fp phi,
+	 fp& partial_x_wrt_r, fp& partial_x_wrt_mu, fp& partial_x_wrt_phi,
+	 fp& partial_y_wrt_r, fp& partial_y_wrt_mu, fp& partial_y_wrt_phi,
+	 fp& partial_z_wrt_r, fp& partial_z_wrt_mu, fp& partial_z_wrt_phi);
+void partial_xyz_wrt_r_nu_phi
+	(fp r, fp nu, fp phi,
+	 fp& partial_x_wrt_r, fp& partial_x_wrt_nu, fp& partial_x_wrt_phi,
+	 fp& partial_y_wrt_r, fp& partial_y_wrt_nu, fp& partial_y_wrt_phi,
+	 fp& partial_z_wrt_r, fp& partial_z_wrt_nu, fp& partial_z_wrt_phi);
 
 // partial {mu,nu,phi} / partial {x,y,z}
 fp partial_mu_wrt_y(fp y, fp z);