change test fn for FD tests slightly,

and more importantly, use the *SAME* test fn for both C++ and Maple code!
(before C++ had sin(4.0*rho) while Maple had sin(rho) :( :( :( )


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

1 files changed, 63 insertions, 47 deletions

diff --git a/src/patch/test_patch_system.out b/src/patch/test_patch_system.out
index afb6f86..be3ea52 100644
--- a/src/patch/test_patch_system.out
+++ b/src/patch/test_patch_system.out
@@ -4,69 +4,85 @@
  <____ ____>  Waterloo Maple Inc.
       |       Type ? for help.
 # test_patch_system.maple -- test function for  test_fd_grid
-# $Id: test_patch_system.maple,v 1.2 2001/07/18 16:20:05 jthorn Exp $
+# $Id: test_patch_system.maple,v 1.3 2001/09/09 14:04:04 jthorn Exp $
 > 
 # test function
-> fn := exp(sin(rho)) * tanh(sin(sigma)^2);
-                                                        2
-                     fn := exp(sin(rho)) tanh(sin(sigma) )
+> fn := exp(sin(1.38*rho)) * tanh(0.17+0.83*sin(sigma)^2);
+                                                                2
+              fn := exp(sin(1.38 rho)) tanh(.17 + .83 sin(sigma) )
 
 > 
 # 1st derivatives
 > simplify(diff(fn, rho));		codegen[C](%);
-                                                             2
-                  cos(rho) exp(sin(rho)) sinh(-1 + cos(sigma) )
-                - ---------------------------------------------
-                                                 2
-                             cosh(-1 + cos(sigma) )
-
-      t0 = -cos(rho)*exp(sin(rho))*sinh(-1.0+pow(cos(sigma),2.0))/cosh(-1.0+pow
-(cos(sigma),2.0));
+-1.380000000000000 cos(1.380000000000000 rho) exp(sin(1.380000000000000 rho))
+
+                                           2    /
+    sinh(-1. + .8300000000000000 cos(sigma) )  /
+                                              /
+
+                                           2
+    cosh(-1. + .8300000000000000 cos(sigma) )
+
+      t0 = -0.138E1*cos(0.138E1*rho)*exp(sin(0.138E1*rho))*sinh(-1.0+0.83*pow(
+cos(sigma),2.0))/cosh(-1.0+0.83*pow(cos(sigma),2.0));
 > simplify(diff(fn, sigma));		codegen[C](%);
-                       exp(sin(rho)) sin(sigma) cos(sigma)
-                     2 -----------------------------------
-                                                 2 2
-                             cosh(-1 + cos(sigma) )
+                      exp(sin(1.380000000000000 rho)) sin(sigma) cos(sigma)
+    1.660000000000000 -----------------------------------------------------
+                                                                  2 2
+                           cosh(-1. + .8300000000000000 cos(sigma) )
 
-      t0 = 2.0*exp(sin(rho))*sin(sigma)*cos(sigma)/pow(cosh(-1.0+pow(cos(sigma)
-,2.0)),2.0);
+      t0 = 0.166E1*exp(sin(0.138E1*rho))*sin(sigma)*cos(sigma)/pow(cosh(-1.0+
+0.83*pow(cos(sigma),2.0)),2.0);
 > 
 # 2nd derivatives
 > simplify(diff(fn, rho,rho));		codegen[C](%);
-                                            2                      2
-          exp(sin(rho)) sinh(-1 + cos(sigma) ) (sin(rho) - cos(rho) )
-          -----------------------------------------------------------
-                                                2
-                            cosh(-1 + cos(sigma) )
-
-      t0 = 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));
+bytes used=2000136, alloc=1834672, time=0.15
+1.904400000000000 exp(sin(1.380000000000000 rho))
+
+                                           2
+    sinh(-1. + .8300000000000000 cos(sigma) )
+
+                                                               2    /
+    (sin(1.380000000000000 rho) - 1. cos(1.380000000000000 rho) )  /
+                                                                  /
+
+                                           2
+    cosh(-1. + .8300000000000000 cos(sigma) )
+
+      t0 = 0.19044E1*exp(sin(0.138E1*rho))*sinh(-1.0+0.83*pow(cos(sigma),2.0))*
+(sin(0.138E1*rho)-1.0*pow(cos(0.138E1*rho),2.0))/cosh(-1.0+0.83*pow(cos(sigma),
+2.0));
 > simplify(diff(fn, rho,sigma));		codegen[C](%);
-                   cos(rho) exp(sin(rho)) sin(sigma) cos(sigma)
-                 2 --------------------------------------------
-                                                 2 2
-                             cosh(-1 + cos(sigma) )
+2.290800000000000 cos(1.380000000000000 rho) exp(sin(1.380000000000000 rho))
+
+                                                                 2 2
+    sin(sigma) cos(sigma)/cosh(-1. + .8300000000000000 cos(sigma) )
 
-      t0 = 2.0*cos(rho)*exp(sin(rho))*sin(sigma)*cos(sigma)/pow(cosh(-1.0+pow(
-cos(sigma),2.0)),2.0);
+      t0 = 0.22908E1*cos(0.138E1*rho)*exp(sin(0.138E1*rho))*sin(sigma)*cos(
+sigma)/pow(cosh(-1.0+0.83*pow(cos(sigma),2.0)),2.0);
 > simplify(diff(fn, sigma,sigma));	codegen[C](%);
-bytes used=2000632, alloc=1834672, time=0.14
-                                         2            2
--2 exp(sin(rho)) (-4 sinh(-1 + cos(sigma) ) cos(sigma)
+-.06640000000000000 exp(sin(1.380000000000000 rho)) (
+
+                                                2            2
+    -83. sinh(-1. + .8300000000000000 cos(sigma) ) cos(sigma)
+
+                                                  2            4
+     + 83. sinh(-1. + .8300000000000000 cos(sigma) ) cos(sigma)
 
-                             2            4
-     + 4 sinh(-1 + cos(sigma) ) cos(sigma)
+                     2                                        2
+     - 50. cos(sigma)  cosh(-1. + .8300000000000000 cos(sigma) )
 
-                   2                     2                        2     /
-     - 2 cos(sigma)  cosh(-1 + cos(sigma) ) + cosh(-1 + cos(sigma) ))  /
-                                                                      /
+                                                  2     /
+     + 25. cosh(-1. + .8300000000000000 cos(sigma) ))  /
+                                                      /
 
-                        2 3
-    cosh(-1 + cos(sigma) )
+                                           2 3
+    cosh(-1. + .8300000000000000 cos(sigma) )
 
-      t0 = -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);
+      t0 = -0.664E-1*exp(sin(0.138E1*rho))*(-83.0*sinh(-1.0+0.83*pow(cos(sigma)
+,2.0))*pow(cos(sigma),2.0)+83.0*sinh(-1.0+0.83*pow(cos(sigma),2.0))*pow(cos(
+sigma),4.0)-50.0*pow(cos(sigma),2.0)*cosh(-1.0+0.83*pow(cos(sigma),2.0))+25.0*
+cosh(-1.0+0.83*pow(cos(sigma),2.0)))/pow(cosh(-1.0+0.83*pow(cos(sigma),2.0)),
+3.0);
 > quit
-bytes used=2313224, alloc=1834672, time=0.17
+bytes used=3055016, alloc=1965720, time=0.24