this was just a quick-n-dirty draft; cg/horizon.c is much better now

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

1 files changed, 0 insertions, 283 deletions

diff --git a/src/gr/H.c b/src/gr/H.c
deleted file mode 100644
index 2fe4da5..0000000
--- a/src/gr/H.c
+++ /dev/null
@@ -1,283 +0,0 @@
-      t1 = 1/r;
-      t3 = Diff(h,rho);
-      t5 = Diff(h,sigma);
-      t7 = zz*t1-X_ud[0][2]*t3-X_ud[1][2]*t5;
-      t8 = g_uu[0][2]*t7;
-      t12 = xx*t1-X_ud[0][0]*t3-X_ud[1][0]*t5;
-      t13 = g_uu[0][1]*t12;
-      t15 = r*r;
-      t17 = 1/t15/r;
-      t22 = Diff(h,rho,rho);
-      t25 = Diff(h,rho,sigma);
-      t30 = Diff(h,sigma,sigma);
-      t32 = xx*yy*t17-X_udd[0][0][1]*t3-X_udd[1][0][1]*t5-X_ud[0][0]*X_ud[0][1]
-*t22-X_ud[1][0]*X_ud[0][1]*t25-X_ud[0][0]*X_ud[1][1]*t25-X_ud[1][0]*X_ud[1][1]*
-t30;
-      t36 = g_uu[1][2]*t7;
-      t40 = yy*t1-X_ud[0][1]*t3-X_ud[1][1]*t5;
-      t41 = g_uu[1][1]*t40;
-      t42 = xx*xx;
-      t43 = zz*zz;
-      t48 = X_ud[0][1]*X_ud[0][1];
-      t53 = X_ud[1][1]*X_ud[1][1];
-      t55 = (t42+t43)*t17-X_udd[0][1][1]*t3-X_udd[1][1][1]*t5-t48*t22-2.0*X_ud
-[1][1]*X_ud[0][1]*t25-t53*t30;
-      t59 = t12*t12;
-      t60 = g_uu[0][1]*t59;
-      t64 = t59*t12;
-      t66 = Diff(g_uu[0][0],yy);
-      t70 = Diff(g_uu[0][0],zz);
-      t74 = Diff(g_uu[0][0],xx);
-      t77 = g_uu[2][2]*g_uu[2][2];
-      t78 = t7*t7;
-      t80 = yy*yy;
-      t85 = X_ud[0][2]*X_ud[0][2];
-      t90 = X_ud[1][2]*X_ud[1][2];
-      t92 = (t42+t80)*t17-X_udd[0][2][2]*t3-X_udd[1][2][2]*t5-t85*t22-2.0*X_ud
-[1][2]*X_ud[0][2]*t25-t90*t30;
-      t94 = g_uu[0][0]*g_uu[0][0];
-      t100 = X_ud[0][0]*X_ud[0][0];
-      t105 = X_ud[1][0]*X_ud[1][0];
-      t107 = (t80+t43)*t17-X_udd[0][0][0]*t3-X_udd[1][0][0]*t5-t100*t22-2.0*
-X_ud[1][0]*X_ud[0][0]*t25-t105*t30;
-      t109 = t78*t7;
-      t111 = Diff(g_uu[2][2],zz);
-      t114 = Diff(g_uu[0][1],yy);
-      t115 = t114*t40;
-      t118 = g_uu[0][2]*t59;
-      t131 = xx*zz*t17-X_udd[0][0][2]*t3-X_udd[1][0][2]*t5-X_ud[0][0]*X_ud[0]
-[2]*t22-X_ud[1][0]*X_ud[0][2]*t25-X_ud[0][0]*X_ud[1][2]*t25-X_ud[1][0]*X_ud[1]
-[2]*t30;
-      t147 = yy*zz*t17-X_udd[0][1][2]*t3-X_udd[1][1][2]*t5-X_ud[0][1]*X_ud[0]
-[2]*t22-X_ud[1][1]*X_ud[0][2]*t25-X_ud[0][1]*X_ud[1][2]*t25-X_ud[1][1]*X_ud[1]
-[2]*t30;
-      t151 = -2.0*t8*t13*t32-2.0*t36*t41*t55-2.0*t60*g_uu[0][0]*t32-g_uu[0][1]*
-t64*t66/2.0-g_uu[0][2]*t64*t70/2.0-g_uu[0][0]*t64*t74/2.0-t77*t78*t92-t94*t59*
-t107-g_uu[2][2]*t109*t111/2.0-t36*t115*t12-2.0*t118*g_uu[0][0]*t131-2.0*t118*
-g_uu[0][1]*t147;
-      t152 = t40*t40;
-      t153 = g_uu[1][1]*t152;
-      t154 = g_uu[0][1]*t32;
-      t157 = g_uu[1][2]*t152;
-      t164 = g_uu[1][2]*t78;
-      t168 = g_uu[2][2]*t78;
-      t169 = g_uu[0][2]*t131;
-      t172 = g_uu[1][2]*t147;
-      t175 = t152*t40;
-      t177 = Diff(g_uu[1][1],yy);
-      t180 = g_uu[0][1]*g_uu[0][1];
-      t184 = Diff(g_uu[1][1],zz);
-      t187 = t13*t55;
-      t190 = g_uu[2][2]*t7;
-      t191 = g_uu[0][0]*t12;
-      t192 = t191*t131;
-      t195 = g_uu[1][2]*t40;
-      t196 = g_uu[0][2]*t12;
-      t197 = t196*t92;
-      t200 = -2.0*t153*t154-2.0*t157*g_uu[0][1]*t131-2.0*t157*g_uu[1][1]*t147
--2.0*t164*g_uu[0][2]*t32-2.0*t168*t169-2.0*t168*t172-g_uu[1][1]*t175*t177/2.0-
-t180*t59*t55-g_uu[1][2]*t175*t184/2.0-2.0*t36*t187-2.0*t190*t192-2.0*t195*t197;
-      t202 = t196*t147;
-      t205 = t13*t147;
-      t208 = g_uu[0][2]*g_uu[0][2];
-      t217 = g_uu[1][2]*g_uu[1][2];
-      t222 = g_uu[0][1]*t40;
-      t223 = t74*t59;
-      t226 = t66*t59;
-      t229 = t70*t59;
-      t238 = g_uu[0][0]*t59;
-      t239 = Diff(g_uu[0][1],xx);
-      t240 = t239*t40;
-      t242 = -2.0*t36*t202-2.0*t190*t205-2.0*t208*t7*t12*t131-2.0*t180*t40*t12*
-t32-2.0*t217*t7*t40*t147-t222*t223/2.0-t41*t226/2.0-t195*t229/2.0-t8*t223/2.0-
-t36*t226/2.0-t190*t229/2.0-t238*t240;
-      t244 = Diff(g_uu[0][1],zz);
-      t245 = t244*t40;
-      t247 = g_uu[0][1]*t152;
-      t250 = t114*t12;
-      t254 = Diff(g_uu[0][2],xx);
-      t255 = t254*t7;
-      t257 = Diff(g_uu[0][2],yy);
-      t258 = t257*t7;
-      t260 = Diff(g_uu[0][2],zz);
-      t261 = t260*t7;
-      t263 = g_uu[0][2]*t78;
-      t268 = t260*t12;
-      t274 = -t60*t115-t118*t245-t247*t239*t12-t153*t250-t157*t244*t12-t238*
-t255-t60*t258-t118*t261-t263*t254*t12-t164*t257*t12-t168*t268-2.0*t190*t197-
-t208*t78*t107;
-      t280 = Diff(g_uu[2][2],xx);
-      t284 = Diff(g_uu[1][1],xx);
-      t291 = g_uu[1][1]*g_uu[1][1];
-      t295 = Diff(g_uu[2][2],yy);
-      t300 = Diff(g_uu[1][2],xx);
-      t301 = t300*t7;
-      t304 = Diff(g_uu[1][2],zz);
-      t305 = t304*t7;
-      t308 = Diff(g_uu[1][2],yy);
-      t309 = t308*t7;
-      t312 = t177*t152;
-      t315 = -t217*t78*t55-g_uu[0][2]*t109*t280/2.0-g_uu[0][1]*t175*t284/2.0-
-t208*t59*t92-t180*t152*t107-t291*t152*t55-g_uu[1][2]*t109*t295/2.0-t217*t152*
-t92-t191*t301*t40-t196*t305*t40-t13*t309*t40-t13*t312/2.0;
-      t316 = t280*t78;
-      t319 = t295*t78;
-      t322 = t111*t78;
-      t331 = t184*t152;
-      t334 = t284*t152;
-      t344 = -t191*t316/2.0-t13*t319/2.0-t196*t322/2.0-t222*t316/2.0-t41*t319/
-2.0-t195*t322/2.0-t196*t331/2.0-t8*t334/2.0-t36*t312/2.0-t190*t331/2.0-t247*
-t301-t153*t309-t157*t305;
-      t350 = t304*t40;
-      t374 = -t263*t300*t40-t164*t308*t40-t168*t350-t191*t334/2.0-t190*t245*t12
--t195*t261*t12-t8*t240*t12-2.0*t36*t222*t32-2.0*t190*t222*t131-t222*t255*t12
--2.0*t8*t41*t32-2.0*t190*t195*t92;
-      t380 = t191*t107;
-      t383 = t191*t32;
-      t407 = -2.0*t8*t195*t131-2.0*t41*t187-2.0*t222*t380-2.0*t41*t383-2.0*t8*
-t380-2.0*t222*t196*t131-2.0*t195*t192-2.0*t195*t205-2.0*t190*t41*t147-2.0*t8*
-t222*t107-t41*t258*t12-2.0*t36*t383-2.0*t41*t202;
-      t417 = t238+2.0*t222*t12+2.0*t8*t12+t153+2.0*t36*t40+t168;
-      t418 = sqrt(t417);
-      t428 = t74*t12+t250+t268+t240+t177*t40+t350+t255+t309+t111*t7+g_uu[0][0]*
-t107+2.0*t154+2.0*t169;
-      t450 = g_uu[1][1]*t55+2.0*t172+g_uu[2][2]*t92+g_uu[0][0]*
-partial_d_ln_sqrt_g[0]*t12+g_uu[0][1]*partial_d_ln_sqrt_g[1]*t12+g_uu[0][2]*
-partial_d_ln_sqrt_g[2]*t12+g_uu[0][1]*partial_d_ln_sqrt_g[0]*t40+g_uu[1][1]*
-partial_d_ln_sqrt_g[1]*t40+g_uu[1][2]*partial_d_ln_sqrt_g[2]*t40+g_uu[0][2]*
-partial_d_ln_sqrt_g[0]*t7+g_uu[1][2]*partial_d_ln_sqrt_g[1]*t7+g_uu[2][2]*
-partial_d_ln_sqrt_g[2]*t7;
-      t469 = (t151+t200+t242+t274+t315+t344+t374+t407)/t418/t417+(t428+t450)/
-t418+(K_uu[0][0]*t59+2.0*K_uu[0][1]*t40*t12+2.0*K_uu[0][2]*t7*t12+K_uu[1][1]*
-t152+2.0*K_uu[1][2]*t7*t40+K_uu[2][2]*t78)/t417-K;
-      t44 = Diff(h,sigma);
-      t62 = 1/r;
-      t45 = Diff(h,rho);
-      t56 = X_ud[1][0];
-      t57 = X_ud[0][0];
-      t16 = xx*t62-t57*t45-t56*t44;
-      t53 = X_ud[0][2];
-      t52 = X_ud[1][2];
-      t14 = zz*t62-t53*t45-t52*t44;
-      t111 = t16*t14;
-      t39 = Diff(h,rho,sigma);
-      t109 = -2.0*t39;
-      t46 = g_uu[2][2];
-      t47 = g_uu[0][1];
-      t48 = g_uu[0][0];
-      t49 = g_uu[0][2];
-      t50 = g_uu[1][2];
-      t8 = t16*t16;
-      t51 = g_uu[1][1];
-      t11 = t14*t14;
-      t54 = X_ud[1][1];
-      t55 = X_ud[0][1];
-      t15 = yy*t62-t55*t45-t54*t44;
-      t12 = t15*t15;
-      t76 = 2.0*t15;
-      t77 = 2.0*t111;
-      t1 = t48*t8+t46*t11+t49*t77+t51*t12+(t47*t16+t50*t14)*t76;
-      t108 = 1/pow(t1,1.0/2.0);
-      t107 = -2.0*t47;
-      t35 = Diff(t48,zz);
-      t106 = -t35/2.0;
-      t105 = -t50/2.0;
-      t22 = Diff(t46,xx);
-      t104 = -t22/2.0;
-      t37 = Diff(t48,xx);
-      t103 = -t37/2.0;
-      t27 = Diff(t51,yy);
-      t102 = -t27/2.0;
-      t20 = Diff(t46,zz);
-      t101 = -t20/2.0;
-      t100 = -Diff(t51,zz)/2.0;
-      t99 = -t47/2.0;
-      t98 = -Diff(t51,xx)/2.0;
-      t97 = -t51/2.0;
-      t38 = Diff(h,sigma,sigma);
-      t40 = Diff(h,rho,rho);
-      t88 = t40*t53;
-      t59 = yy*yy;
-      t60 = xx*xx;
-      t61 = t62/(r*r);
-      t5 = (t60+t59)*t61-X_udd[0][2][2]*t45-X_udd[1][2][2]*t44-t52*t52*t38+(t52
-*t109-t88)*t53;
-      t96 = t46*t5;
-      t90 = xx*t61;
-      t83 = t40*t55;
-      t84 = t38*t54;
-      t3 = -yy*t90-X_udd[0][0][1]*t45-X_udd[1][0][1]*t44-t57*t83-t56*t84+(-t56*
-t55-t57*t54)*t39;
-      t95 = t47*t3;
-      t58 = zz*zz;
-      t89 = t38*t56;
-      t6 = (t59+t58)*t61-X_udd[0][0][0]*t45-X_udd[1][0][0]*t44-t57*t57*t40+(t57
-*t109-t89)*t56;
-      t94 = t48*t6;
-      t2 = -zz*t90-X_udd[0][0][2]*t45-X_udd[1][0][2]*t44-t57*t88-t52*t89+(-t56*
-t53-t57*t52)*t39;
-      t93 = t49*t2;
-      t4 = -yy*zz*t61-X_udd[0][1][2]*t45-X_udd[1][1][2]*t44-t53*t83-t52*t84+(-
-t54*t53-t55*t52)*t39;
-      t92 = t50*t4;
-      t7 = (t60+t58)*t61-X_udd[0][1][1]*t45-X_udd[1][1][1]*t44-t55*t55*t40+(t55
-*t109-t84)*t54;
-      t91 = t51*t7;
-      t24 = Diff(t50,yy);
-      t31 = Diff(t49,xx);
-      t87 = -t31-t24;
-      t29 = Diff(t49,zz);
-      t33 = Diff(t47,yy);
-      t86 = -t33-t29;
-      t23 = Diff(t50,zz);
-      t34 = Diff(t47,xx);
-      t85 = -t34-t23;
-      t82 = -2.0*t95;
-      t81 = -2.0*t93;
-      t80 = -2.0*t92;
-      t79 = -2.0*t3*t49;
-      t78 = -2.0*t48*t50;
-      t75 = 2.0*t47;
-      t74 = 2.0*t49;
-      t73 = -t95-t96;
-      t72 = -t92-t94;
-      t71 = -t91-t93;
-      t70 = t82+t81-t94;
-      t43 = t49*t49;
-      t42 = t50*t50;
-      t41 = t47*t47;
-      t36 = Diff(t48,yy);
-      t32 = Diff(t47,zz);
-      t30 = Diff(t49,yy);
-      t25 = Diff(t50,xx);
-      t21 = Diff(t46,yy);
-      t19 = partial_d_ln_sqrt_g[0];
-      t18 = partial_d_ln_sqrt_g[1];
-      t17 = partial_d_ln_sqrt_g[2];
-      MapleGenVar1 = (t96+t91+2.0*t92+(t20+t49*t19+t50*t18+t46*t17-t87)*t14+(
-t37+t47*t18+t49*t17+t48*t19-t86)*t16+(t50*t17+t27+t47*t19+t51*t18-t85)*t15-t70)
-*t108;
-      MapleGenVar3 = -K;
-      MapleGenVar9 = (t3*t78+(-t50*t7-t46*t4)*t75+2.0*(-t46*t48-t43)*t2+(t72+
-t73)*t74)*t111+((t51*t79+(-t46*t2-t49*t6)*t75+2.0*(-t42-t46*t51)*t4+2.0*(t71+
-t73)*t50)*t14+(t2*t78+(-t51*t4-t50*t5)*t74+2.0*(-t51*t48-t41)*t3+(t71+t72)*t75+
-(-t51*t30-t46*t32-t48*t25+t86*t50+t85*t49+t87*t47)*t14)*t16)*t15;
-      MapleGenVar10 = MapleGenVar9+((t48*t98-t50*t32-t51*t33+t49*t100+(-t34+
-t102)*t47)*t16+(t50*t100+t27*t97+t47*t98)*t15-t42*t5+t50*t2*t107-t41*t6+(t80+
-t82-t91)*t51+(-t47*t25-t51*t24+t49*t98+t46*t100+(-t23+t102)*t50)*t14)*t12;
-      MapleGenVar8 = MapleGenVar10+((t21*t105+t49*t104+t46*t101)*t14+(t48*t104-
-t50*t30-t46*t29+t21*t99+(t101-t31)*t49)*t16+t50*t79-t42*t7-t43*t6+(t81+t80-t96)
-*t46+(t21*t97-t46*t23+t22*t99-t49*t25+(t101-t24)*t50)*t15)*t11+((t48*t103+t36*
-t99+t49*t106)*t16-t43*t5+t49*t4*t107-t41*t7+t70*t48+(t46*t106+t36*t105-t47*t30-
-t48*t31+(t103-t29)*t49)*t14+(-t48*t34+t35*t105-t49*t32+t36*t97+(t103-t33)*t47)*
-t15)*t8;
-      MapleGenVar9 = t108;
-      MapleGenVar7 = MapleGenVar8*MapleGenVar9;
-      MapleGenVar8 = K_uu[0][0]*t8+K_uu[2][2]*t11;
-      MapleGenVar6 = MapleGenVar7+MapleGenVar8;
-      MapleGenVar5 = MapleGenVar6+K_uu[0][2]*t77+K_uu[1][1]*t12+(K_uu[1][2]*t14
-+K_uu[0][1]*t16)*t76;
-      MapleGenVar6 = 1/t1;
-      MapleGenVar4 = MapleGenVar5*MapleGenVar6;
-      MapleGenVar2 = MapleGenVar3+MapleGenVar4;
-      t9 = MapleGenVar1+MapleGenVar2;