diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-11-21 20:39:51 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-11-21 20:39:51 +0000 |
commit | 03bd2968bc2cafe03a4d76c3609526bd23896a0f (patch) | |
tree | e53cca7d80221f9157510fdacf982cd4fb504273 /src/gr.cg/horizon_function.c | |
parent | 47f5db5fc54dfed7c03e6d5264a25ee02e382347 (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.c | 182 |
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; |