add patch::partial_xyz_wrt_rho_sigma()

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

1 files changed, 53 insertions, 0 deletions

diff --git a/src/patch/patch.hh b/src/patch/patch.hh
index 950a7f6..025aecd 100644
--- a/src/patch/patch.hh
+++ b/src/patch/patch.hh
@@ -218,6 +218,14 @@ public:
 					 fp& xcos, fp& ycos, fp& zcos)
 		const = 0;
 
+	// partial (x,y,z) / partial (rho,sigma)
+	virtual void partial_xyz_wrt_rho_sigma
+		(fp r, fp rho, fp sigma,
+		 fp& partial_x_wrt_rho, fp& partial_x_wrt_sigma,
+		 fp& partial_y_wrt_rho, fp& partial_y_wrt_sigma,
+		 fp& partial_z_wrt_rho, fp& partial_z_wrt_sigma)
+		const = 0;
+
 	// partial (rho,sigma) / partial (x,y,z)
 	virtual fp partial_rho_wrt_x(fp x, fp y, fp z) const = 0;
 	virtual fp partial_rho_wrt_y(fp x, fp y, fp z) const = 0;
@@ -604,6 +612,21 @@ public:
 		const
 		{ local_coords::xyzcos_of_mu_nu(rho,sigma, xcos,ycos,zcos); }
 
+	// partial (x,y,z) / partial (rho,sigma)
+	void partial_xyz_wrt_rho_sigma
+		(fp r, fp rho, fp sigma,
+		 fp& partial_x_wrt_rho, fp& partial_x_wrt_sigma,
+		 fp& partial_y_wrt_rho, fp& partial_y_wrt_sigma,
+		 fp& partial_z_wrt_rho, fp& partial_z_wrt_sigma)
+		const
+		{
+		local_coords::
+		 partial_xyz_wrt_mu_nu(r, rho, sigma,
+				       partial_x_wrt_rho, partial_x_wrt_sigma,
+				       partial_y_wrt_rho, partial_y_wrt_sigma,
+				       partial_z_wrt_rho, partial_z_wrt_sigma);
+		}
+
 	// partial (rho,sigma) / partial (x,y,z)
 	fp partial_rho_wrt_x(fp x, fp y, fp z) const { return 0.0; }
 	fp partial_rho_wrt_y(fp x, fp y, fp z) const
@@ -723,6 +746,21 @@ public:
 		const
 		{ local_coords::xyzcos_of_nu_phi(rho,sigma, xcos,ycos,zcos); }
 
+	// partial (x,y,z) / partial (rho,sigma)
+	void partial_xyz_wrt_rho_sigma
+		(fp r, fp rho, fp sigma,
+		 fp& partial_x_wrt_rho, fp& partial_x_wrt_sigma,
+		 fp& partial_y_wrt_rho, fp& partial_y_wrt_sigma,
+		 fp& partial_z_wrt_rho, fp& partial_z_wrt_sigma)
+		const
+		{
+		local_coords::
+		 partial_xyz_wrt_nu_phi(r, rho, sigma,
+					partial_x_wrt_rho, partial_x_wrt_sigma,
+					partial_y_wrt_rho, partial_y_wrt_sigma,
+					partial_z_wrt_rho, partial_z_wrt_sigma);
+		}
+
 	// partial (rho,sigma) / partial (x,y,z)
 	fp partial_rho_wrt_x(fp x, fp y, fp z) const
 		{ return local_coords::partial_nu_wrt_x(x,z); }
@@ -844,6 +882,21 @@ public:
 		const
 		{ local_coords::xyzcos_of_mu_phi(rho,sigma, xcos,ycos,zcos); }
 
+	// partial (x,y,z) / partial (rho,sigma)
+	void partial_xyz_wrt_rho_sigma
+		(fp r, fp rho, fp sigma,
+		 fp& partial_x_wrt_rho, fp& partial_x_wrt_sigma,
+		 fp& partial_y_wrt_rho, fp& partial_y_wrt_sigma,
+		 fp& partial_z_wrt_rho, fp& partial_z_wrt_sigma)
+		const
+		{
+		local_coords::
+		 partial_xyz_wrt_mu_phi(r, rho, sigma,
+					partial_x_wrt_rho, partial_x_wrt_sigma,
+					partial_y_wrt_rho, partial_y_wrt_sigma,
+					partial_z_wrt_rho, partial_z_wrt_sigma);
+		}
+
 	// partial (rho,sigma) / partial (x,y,z)
 	fp partial_rho_wrt_x(fp x, fp y, fp z) const { return 0.0; }
 	fp partial_rho_wrt_y(fp x, fp y, fp z) const