aboutsummaryrefslogtreecommitdiff
path: root/src/gr.cg/horizon_function.c
diff options
context:
space:
mode:
authorjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2002-11-21 20:39:51 +0000
committerjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2002-11-21 20:39:51 +0000
commit03bd2968bc2cafe03a4d76c3609526bd23896a0f (patch)
treee53cca7d80221f9157510fdacf982cd4fb504273 /src/gr.cg/horizon_function.c
parent47f5db5fc54dfed7c03e6d5264a25ee02e382347 (diff)
change in terminology/notation:
LHS of apparent horizon equation was formerly called $H$, changed to $\Theta$ because it is in fact precisely the expansion of the surface r = h(angle). This implies renaming a lot of variables & functions & a few parameters (which aren't specified in most par files outside my own for testing this thorn) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@901 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/gr.cg/horizon_function.c')
-rw-r--r--src/gr.cg/horizon_function.c182
1 files changed, 0 insertions, 182 deletions
diff --git a/src/gr.cg/horizon_function.c b/src/gr.cg/horizon_function.c
deleted file mode 100644
index 7e6e717..0000000
--- a/src/gr.cg/horizon_function.c
+++ /dev/null
@@ -1,182 +0,0 @@
-/*
- * inputs = {r, partial_d_ln_sqrt_g, partial_d_g_uu, X_ud, X_udd, g_uu, K_uu, h}
- * outputs = {HA, HB, HC, HD}
- * cost = 134*assignments+3*divisions+5*functions+401*multiplications+173*additions
- */
-fp t1, t2, t4, t5, t7, t8, t10, t11, t12, t14;
-fp t16, t18, t19, t20, t23, t24, t25, t26, t29, t30;
-fp t31, t33, t35, t37, t38, t40, t41, t42, t43, t45;
-fp t46, t47, t48, t50, t51, t52, t53, t55, t58, t59;
-fp t61, t64, t67, t68, t70, t72, t73, t75, t77, t83;
-fp t84, t87, t90, t91, t93, t95, t112, t130, t134, t137;
-fp t139, t143, t147, t150, t152, t159, t164, t166, t168, t176;
-fp t181, t183, t199, t203, t206, t209, t213, t216, t217, t219;
-fp t228, t229, t232, t235, t236, t237, t240, t243, t247, t251;
-fp t254, t259, t260, t262, t263, t269, t272, t273, t276, t279;
-fp t282, t283, t290, t293, t296, t299, t302, t305, t312, t315;
-fp t318, t324, t325, t329, t330, t338, t358, t381, t392, t395;
-fp t396, t397, t408, t411, t431, t440, t444, t447, t450, t465;
- t1 = g_uu_33;
- t2 = 1/r;
- t4 = X_ud_13;
- t5 = PARTIAL_RHO(h);
- t7 = X_ud_23;
- t8 = PARTIAL_SIGMA(h);
- t10 = zz*t2-t4*t5-t7*t8;
- t11 = t1*t10;
- t12 = partial_d_g_uu_311;
- t14 = X_ud_11;
- t16 = X_ud_21;
- t18 = xx*t2-t14*t5-t16*t8;
- t19 = t18*t18;
- t20 = t12*t19;
- t23 = g_uu_23;
- t24 = t23*t10;
- t25 = partial_d_g_uu_211;
- t26 = t25*t19;
- t29 = g_uu_13;
- t30 = t29*t19;
- t31 = partial_d_g_uu_312;
- t33 = X_ud_12;
- t35 = X_ud_22;
- t37 = yy*t2-t33*t5-t35*t8;
- t38 = t31*t37;
- t40 = g_uu_12;
- t41 = t40*t19;
- t42 = partial_d_g_uu_212;
- t43 = t42*t37;
- t45 = g_uu_11;
- t46 = t45*t19;
- t47 = partial_d_g_uu_112;
- t48 = t47*t37;
- t50 = g_uu_22;
- t51 = t37*t37;
- t52 = t50*t51;
- t53 = t42*t18;
- t55 = t40*t51;
- t58 = partial_d_g_uu_113;
- t59 = t58*t10;
- t61 = t23*t51;
- t64 = t29*t10;
- t67 = partial_d_g_uu_213;
- t68 = t67*t10;
- t70 = t45*t45;
- t72 = yy*yy;
- t73 = zz*zz;
- t75 = r*r;
- t77 = 1/t75/r;
- t83 = t14*t14;
- t84 = PARTIAL_RHO_RHO(h);
- t87 = PARTIAL_RHO_SIGMA(h);
- t90 = t16*t16;
- t91 = PARTIAL_SIGMA_SIGMA(h);
- t93 = (t72+t73)*t77-X_udd_111*t5-X_udd_211*t8-t83*t84-2.0*t16*t14*t87-t90
-*t91;
- t95 = RATIONAL(-1.0,2.0)*t11*t20+RATIONAL(-1.0,2.0)*t24*t26-t30*t38-t41*
-t43-t46*t48-t52*t53-t55*t47*t18-t46*t59-t61*t31*t18-t64*t48*t18-t41*t68-t70*t19
-*t93;
- t112 = -xx*yy*t77-X_udd_112*t5-X_udd_212*t8-t14*t33*t84-t16*t33*t87-t14*
-t35*t87-t16*t35*t91;
- t130 = -xx*zz*t77-X_udd_113*t5-X_udd_213*t8-t14*t4*t84-t16*t4*t87-t14*t7*
-t87-t16*t7*t91;
- t134 = t40*t37;
- t137 = t51*t37;
- t139 = partial_d_g_uu_322;
- t143 = partial_d_g_uu_122;
- t147 = partial_d_g_uu_222;
- t150 = t40*t40;
- t152 = xx*xx;
- t159 = t33*t33;
- t164 = t35*t35;
- t166 = (t152+t73)*t77-X_udd_122*t5-X_udd_222*t8-t159*t84-2.0*t35*t33*t87-
-t164*t91;
- t168 = t29*t29;
- t176 = t4*t4;
- t181 = t7*t7;
- t183 = (t152+t72)*t77-X_udd_133*t5-X_udd_233*t8-t176*t84-2.0*t7*t4*t87-
-t181*t91;
- t199 = -yy*zz*t77-X_udd_123*t5-X_udd_223*t8-t33*t4*t84-t35*t4*t87-t33*t7*
-t87-t35*t7*t91;
- t203 = t40*t112;
- t206 = t50*t37;
- t209 = -t24*t43*t18-2.0*t41*t45*t112-2.0*t30*t45*t130-t134*t59*t18+
-RATIONAL(-1.0,2.0)*t23*t137*t139+RATIONAL(-1.0,2.0)*t40*t137*t143+RATIONAL(-1.0
-,2.0)*t50*t137*t147-t150*t19*t166-t168*t19*t183-2.0*t30*t40*t199-2.0*t52*t203-
-t206*t68*t18;
- t213 = t50*t50;
- t216 = t10*t10;
- t217 = t216*t10;
- t219 = partial_d_g_uu_333;
- t228 = partial_d_g_uu_123;
- t229 = t228*t10;
- t232 = partial_d_g_uu_233;
- t235 = t23*t37;
- t236 = partial_d_g_uu_313;
- t237 = t236*t10;
- t240 = t23*t23;
- t243 = t29*t18;
- t247 = t243*t183;
- t251 = partial_d_g_uu_133;
- t254 = -t150*t51*t93-t213*t51*t166+RATIONAL(-1.0,2.0)*t1*t217*t219-2.0*
-t61*t40*t130-2.0*t61*t50*t199-t55*t229+RATIONAL(-1.0,2.0)*t23*t217*t232-t235*
-t237*t18-t240*t51*t183-2.0*t134*t243*t130-2.0*t235*t247+RATIONAL(-1.0,2.0)*t29*
-t217*t251;
- t259 = partial_d_g_uu_323;
- t260 = t259*t10;
- t262 = partial_d_g_uu_223;
- t263 = t262*t10;
- t269 = t243*t199;
- t272 = t40*t18;
- t273 = t272*t166;
- t276 = t272*t199;
- t279 = t219*t216;
- t282 = t45*t18;
- t283 = t251*t216;
- t290 = t232*t216;
- t293 = t23*t216;
- t296 = -2.0*t150*t37*t18*t112-t61*t260-t52*t263-2.0*t168*t10*t18*t130-2.0
-*t206*t269-2.0*t206*t273-2.0*t235*t276+RATIONAL(-1.0,2.0)*t243*t279+RATIONAL(
--1.0,2.0)*t282*t283+RATIONAL(-1.0,2.0)*t235*t279+RATIONAL(-1.0,2.0)*t134*t283+
-RATIONAL(-1.0,2.0)*t206*t290-t293*t262*t37;
- t299 = t29*t216;
- t302 = t147*t51;
- t305 = t139*t51;
- t312 = t282*t93;
- t315 = t282*t130;
- t318 = t282*t112;
- t324 = t1*t216;
- t325 = t236*t18;
- t329 = -t299*t228*t37+RATIONAL(-1.0,2.0)*t272*t302+RATIONAL(-1.0,2.0)*t11
-*t305+RATIONAL(-1.0,2.0)*t24*t302+RATIONAL(-1.0,2.0)*t272*t290-2.0*t134*t312
--2.0*t235*t315-2.0*t206*t318-t30*t237-2.0*t24*t269-t324*t325-2.0*t11*t247;
- t330 = t259*t37;
- t338 = t143*t51;
- t358 = -t324*t330-2.0*t11*t315-t272*t263*t37-2.0*t11*t276+RATIONAL(-1.0,
-2.0)*t282*t338-t243*t260*t37-2.0*t64*t312-t282*t229*t37-t299*t58*t18-2.0*t24*
-t273-2.0*t24*t318+RATIONAL(-1.0,2.0)*t64*t338-2.0*t64*t235*t130;
- t381 = t29*t130;
- t392 = t23*t199;
- t395 = -2.0*t240*t10*t37*t199-2.0*t11*t235*t183-2.0*t11*t134*t130-2.0*t64
-*t134*t93-2.0*t24*t206*t166-t11*t38*t18-2.0*t293*t29*t112-2.0*t324*t381-2.0*t11
-*t206*t199-2.0*t64*t206*t112-t168*t216*t93-2.0*t324*t392;
- t396 = partial_d_g_uu_111;
- t397 = t396*t19;
- t408 = t1*t1;
- t411 = t19*t18;
- t431 = RATIONAL(-1.0,2.0)*t134*t397+RATIONAL(-1.0,2.0)*t206*t26+RATIONAL(
--1.0,2.0)*t235*t20+RATIONAL(-1.0,2.0)*t64*t397-t240*t216*t166-t408*t216*t183+
-RATIONAL(-1.0,2.0)*t45*t411*t396+RATIONAL(-1.0,2.0)*t40*t411*t25+RATIONAL(-1.0,
-2.0)*t29*t411*t12-2.0*t24*t134*t112+RATIONAL(-1.0,2.0)*t243*t305-t293*t67*t18
--2.0*t64*t272*t112;
- HA = t95+t209+t254+t296+t329+t358+t395+t431;
- t440 = t396*t18+t53+t325+t48+t147*t37+t330+t59+t263+t219*t10+t45*t93+2.0*
-t203+2.0*t381;
- t444 = partial_d_ln_sqrt_g_1;
- t447 = partial_d_ln_sqrt_g_2;
- t450 = partial_d_ln_sqrt_g_3;
- t465 = t50*t166+2.0*t392+t1*t183+t45*t444*t18+t40*t447*t18+t29*t450*t18+
-t40*t444*t37+t50*t447*t37+t23*t450*t37+t29*t444*t10+t23*t447*t10+t1*t450*t10;
- HB = t440+t465;
- HC = K_uu_11*t19+2.0*K_uu_12*t37*t18+2.0*K_uu_13*t10*t18+K_uu_22*t51+2.0*
-K_uu_23*t10*t37+K_uu_33*t216;
- HD = t46+2.0*t134*t18+2.0*t64*t18+t52+2.0*t24*t37+t324;