revisions to test program -- not finished yet

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

1 files changed, 50 insertions, 47 deletions

diff --git a/src/patch/test_patch_system.out b/src/patch/test_patch_system.out
index 197877b..5dda17b 100644
--- a/src/patch/test_patch_system.out
+++ b/src/patch/test_patch_system.out
@@ -4,65 +4,68 @@
  <____ ____>  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/06/20 13:12:23 jthorn Exp $
+# $Id: test_patch_system.maple,v 1.1 2001/07/18 16:16:53 jthorn Exp $
 > 
 # test function
-> fn1(x,y) := exp(sin(x)) * tanh(sin(y)^2);
-                                                         2
-                     fn1(x, y) := exp(sin(x)) tanh(sin(y) )
+> fn := proc(x,y) exp(sin(x)) * tanh(sin(y)^2) end;
+              fn := proc(x, y) exp(sin(x))*tanh(sin(y)^2) end proc
 
 > 
 # 1st derivatives
-> simplify(diff(fn1(x,y), x));		codegen[C](%);
-                                                         2
-                      cos(x) exp(sin(x)) sinh(-1 + cos(y) )
-                    - -------------------------------------
-                                               2
-                               cosh(-1 + cos(y) )
+> simplify(diff(fn(rho,sigma), rho));		codegen[C](%);
+                                                             2
+                  cos(rho) exp(sin(rho)) sinh(-1 + cos(sigma) )
+                - ---------------------------------------------
+                                                 2
+                             cosh(-1 + cos(sigma) )
 
-      t0 = -cos(x)*exp(sin(x))*sinh(-1.0+pow(cos(y),2.0))/cosh(-1.0+pow(cos(y),
-2.0));
-> simplify(diff(fn1(x,y), y));		codegen[C](%);
-                            exp(sin(x)) sin(y) cos(y)
-                          2 -------------------------
-                                               2 2
-                               cosh(-1 + cos(y) )
+      t0 = -cos(rho)*exp(sin(rho))*sinh(-1.0+pow(cos(sigma),2.0))/cosh(-1.0+pow
+(cos(sigma),2.0));
+> simplify(diff(fn(rho,sigma), sigma));		codegen[C](%);
+                       exp(sin(rho)) sin(sigma) cos(sigma)
+                     2 -----------------------------------
+                                                 2 2
+                             cosh(-1 + cos(sigma) )
 
-      t0 = 2.0*exp(sin(x))*sin(y)*cos(y)/pow(cosh(-1.0+pow(cos(y),2.0)),2.0);
+      t0 = 2.0*exp(sin(rho))*sin(sigma)*cos(sigma)/pow(cosh(-1.0+pow(cos(sigma)
+,2.0)),2.0);
 > 
 # 2nd derivatives
-> simplify(diff(fn1(x,y), x,x));		codegen[C](%);
-                                           2                  2
-               exp(sin(x)) sinh(-1 + cos(y) ) (sin(x) - cos(x) )
-               -------------------------------------------------
-                                              2
-                              cosh(-1 + cos(y) )
+> simplify(diff(fn(rho,sigma), 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(x))*sinh(-1.0+pow(cos(y),2.0))*(sin(x)-pow(cos(x),2.0))/cosh
-(-1.0+pow(cos(y),2.0));
-> simplify(diff(fn1(x,y), x,y));		codegen[C](%);
-bytes used=2000760, alloc=1703624, time=0.28
-                         cos(x) exp(sin(x)) sin(y) cos(y)
-                       2 --------------------------------
-                                               2 2
-                               cosh(-1 + cos(y) )
+      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));
+> simplify(diff(fn(rho,sigma), rho,sigma));	codegen[C](%);
+bytes used=2000576, alloc=1703624, time=0.28
+                   cos(rho) exp(sin(rho)) sin(sigma) cos(sigma)
+                 2 --------------------------------------------
+                                                 2 2
+                             cosh(-1 + cos(sigma) )
 
-      t0 = 2.0*cos(x)*exp(sin(x))*sin(y)*cos(y)/pow(cosh(-1.0+pow(cos(y),2.0)),
-2.0);
-> simplify(diff(fn1(x,y), y,y));		codegen[C](%);
-                                   2        2                     2        4
--2 exp(sin(x)) (-4 sinh(-1 + cos(y) ) cos(y)  + 4 sinh(-1 + cos(y) ) cos(y)
+      t0 = 2.0*cos(rho)*exp(sin(rho))*sin(sigma)*cos(sigma)/pow(cosh(-1.0+pow(
+cos(sigma),2.0)),2.0);
+> simplify(diff(fn(rho,sigma), sigma,sigma));	codegen[C](%);
+                                         2            2
+-2 exp(sin(rho)) (-4 sinh(-1 + cos(sigma) ) cos(sigma)
 
-               2                 2                    2     /
-     - 2 cos(y)  cosh(-1 + cos(y) ) + cosh(-1 + cos(y) ))  /
-                                                          /
+                             2            4
+     + 4 sinh(-1 + cos(sigma) ) cos(sigma)
 
-                    2 3
-    cosh(-1 + cos(y) )
+                   2                     2                        2     /
+     - 2 cos(sigma)  cosh(-1 + cos(sigma) ) + cosh(-1 + cos(sigma) ))  /
+                                                                      /
 
-      t0 = -2.0*exp(sin(x))*(-4.0*sinh(-1.0+pow(cos(y),2.0))*pow(cos(y),2.0)+
-4.0*sinh(-1.0+pow(cos(y),2.0))*pow(cos(y),4.0)-2.0*pow(cos(y),2.0)*cosh(-1.0+
-pow(cos(y),2.0))+cosh(-1.0+pow(cos(y),2.0)))/pow(cosh(-1.0+pow(cos(y),2.0)),3.0
-);
+                        2 3
+    cosh(-1 + 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);
 > quit
-bytes used=2657856, alloc=1703624, time=0.41
+bytes used=2655448, alloc=1703624, time=0.39