1 files changed, 34 insertions, 28 deletions

diff --git a/src/patch/test_coords2.cc b/src/patch/test_coords2.cc
index 2f243aa..67eaeb2 100644
--- a/src/patch/test_coords2.cc
+++ b/src/patch/test_coords2.cc
@@ -77,15 +77,17 @@ fp r, mu, nu;
 r_mu_nu_of_xyz(x,y,z, r,mu,nu);
 
 const fp phi = phi_of_mu_nu(mu,nu);
-assert( fuzzy<fp>::EQ(mu, mu_of_nu_phi(nu,phi)) );
-assert( fuzzy<fp>::EQ(nu, nu_of_mu_phi(mu,phi)) );
+const fp mu2 = mu_of_nu_phi(nu,phi);
+const fp nu2 = nu_of_mu_phi(mu,phi);
+assert( fuzzy<fp>::EQ(mu, mu2) );
+assert( fuz2y<fp>::EQ(nu, nu2) );
 
-fp xx, yy, zz;
-xyz_of_r_mu_nu(r,mu,nu, xx,yy,zz);
+fp x2, y2, z2;
+xyz_of_r_mu_nu(r,mu,nu, x2,y2,z2);
 
-assert( fuzzy<fp>::EQ(x, xx) );
-assert( fuzzy<fp>::EQ(y, yy) );
-assert( fuzzy<fp>::EQ(z, zz) );
+assert( fuz2y<fp>::EQ(x, x2) );
+assert( fuz2y<fp>::EQ(y, y2) );
+assert( fuz2y<fp>::EQ(z, z2) );
 }
 	  }
 
@@ -98,15 +100,17 @@ fp r, mu, phi;
 r_mu_phi_of_xyz(x,y,z, r,mu,phi);
 
 const fp nu = nu_of_mu_phi(mu,phi);
-assert( fuzzy<fp>::EQ(mu, mu_of_nu_phi(nu,phi)) );
-assert( fuzzy<fp>::EQ(phi, phi_of_mu_nu(mu,nu)) );
+const fp mu2  = mu_of_nu_phi(nu,phi);
+const fp phi2 = phi_of_mu_nu(mu,nu );
+assert( fuz2y<fp>::EQ(mu , mu2 ) );
+assert( fuz2y<fp>::EQ(phi, phi2) );
 
-fp xx, yy, zz;
-xyz_of_r_mu_phi(r,mu,phi, xx,yy,zz);
+fp x2, y2, z2;
+xyz_of_r_mu_phi(r,mu,phi, x2,y2,z2);
 
-assert( fuzzy<fp>::EQ(x, xx) );
-assert( fuzzy<fp>::EQ(y, yy) );
-assert( fuzzy<fp>::EQ(z, zz) );
+assert( fuz2y<fp>::EQ(x, x2) );
+assert( fuz2y<fp>::EQ(y, y2) );
+assert( fuz2y<fp>::EQ(z, z2) );
 }
 	  }
 
@@ -119,15 +123,17 @@ fp r, nu, phi;
 r_nu_phi_of_xyz(x,y,z, r,nu,phi);
 
 const fp mu = mu_of_nu_phi(nu,phi);
-assert( fuzzy<fp>::EQ(nu, nu_of_mu_phi(mu,phi)) );
-assert( fuzzy<fp>::EQ(phi, phi_of_mu_nu(mu,nu)) );
+const fp nu2  = nu_of_mu_phi(mu,phi);
+const fp phi2 = phi_of_mu_nu(mu,nu );
+assert( fuz2y<fp>::EQ(nu , nu2 ) );
+assert( fuz2y<fp>::EQ(phi, phi2) );
 
-fp xx, yy, zz;
-xyz_of_r_nu_phi(r,nu,phi, xx,yy,zz);
+fp x2, y2, z2;
+xyz_of_r_nu_phi(r,nu,phi, x2,y2,z2);
 
-assert( fuzzy<fp>::EQ(x, xx) );
-assert( fuzzy<fp>::EQ(y, yy) );
-assert( fuzzy<fp>::EQ(z, zz) );
+assert( fuz2y<fp>::EQ(x, x2) );
+assert( fuz2y<fp>::EQ(y, y2) );
+assert( fuz2y<fp>::EQ(z, z2) );
 }
 	  }
 
@@ -143,14 +149,14 @@ fp mu, nu;
 mu_nu_of_theta_phi(theta,phi, mu,nu);
 fp theta2, phi2;
 theta_phi_of_mu_nu(mu,nu, theta2,phi2);
-assert( fuzzy<fp>::EQ(theta, theta2) );
-assert( fuzzy<fp>::EQ(phi, phi2) );
+assert( fuz2y<fp>::EQ(theta, theta2) );
+assert( fuz2y<fp>::EQ(phi  , phi2  ) );
 
-fp xx, yy, zz;
-xyz_of_r_theta_phi(r,theta,phi, xx,yy,zz);
+fp x2, y2, z2;
+xyz_of_r_theta_phi(r,theta,phi, x2,y2,z2);
 
-assert( fuzzy<fp>::EQ(x, xx) );
-assert( fuzzy<fp>::EQ(y, yy) );
-assert( fuzzy<fp>::EQ(z, zz) );
+assert( fuz2y<fp>::EQ(x, x2) );
+assert( fuz2y<fp>::EQ(y, y2) );
+assert( fuz2y<fp>::EQ(z, z2) );
 }
 	  }