Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] TwoPunctures/branches/momentum_constraint/ (Rev. 113)
by roland.haas＠physics.gatech.edu 19 Jul '10

19 Jul '10

User: rhaas Date: 2010/07/19 08:46 AM Added: /branches/momentum_constraint/test/ bhns_momentum_fd.par /branches/momentum_constraint/test/bhns_momentum_fd/ alp.d.asc, alp.maximum.asc, alp.minimum.asc, alp.norm1.asc, alp.norm2.asc, alp.x.asc, alp.y.asc, alp.z.asc, betax.d.asc, betax.maximum.asc, betax.minimum.asc, betax.norm1.asc, betax.norm2.asc, betax.x.asc, betax.y.asc, betax.z.asc, betay.d.asc, betay.maximum.asc, betay.minimum.asc, betay.norm1.asc, betay.norm2.asc, betay.x.asc, betay.y.asc, betay.z.asc, betaz.d.asc, betaz.maximum.asc, betaz.minimum.asc, betaz.norm1.asc, betaz.norm2.asc, betaz.x.asc, betaz.y.asc, betaz.z.asc, dens.d.asc, dens.maximum.asc, dens.minimum.asc, dens.norm1.asc, dens.norm2.asc, dens.x.asc, dens.y.asc, dens.z.asc, eps.d.asc, eps.maximum.asc, eps.minimum.asc, eps.norm1.asc, eps.norm2.asc, eps.x.asc, eps.y.asc, eps.z.asc, gxx.d.asc, gxx.maximum.asc, gxx.minimum.asc, gxx.norm1.asc, gxx.norm2.asc, gxx.x.asc, gxx.y.asc, gxx.z.asc, gxy.d.asc, gxy.maximum.asc, gxy.minimum.asc, gxy.norm1.asc, gxy.norm2.asc, gxy.x.asc, gxy.y.asc, gxy.z.asc, gxz.d.asc, gxz.maximum.asc, gxz.minimum.asc, gxz.norm1.asc, gxz.norm2.asc, gxz.x.asc, gxz.y.asc, gxz.z.asc, gyy.d.asc, gyy.maximum.asc, gyy.minimum.asc, gyy.norm1.asc, gyy.norm2.asc, gyy.x.asc, gyy.y.asc, gyy.z.asc, gyz.d.asc, gyz.maximum.asc, gyz.minimum.asc, gyz.norm1.asc, gyz.norm2.asc, gyz.x.asc, gyz.y.asc, gyz.z.asc, gzz.d.asc, gzz.maximum.asc, gzz.minimum.asc, gzz.norm1.asc, gzz.norm2.asc, gzz.x.asc, gzz.y.asc, gzz.z.asc, kxx.d.asc, kxx.maximum.asc, kxx.minimum.asc, kxx.norm1.asc, kxx.norm2.asc, kxx.x.asc, kxx.y.asc, kxx.z.asc, kxy.d.asc, kxy.maximum.asc, kxy.minimum.asc, kxy.norm1.asc, kxy.norm2.asc, kxy.x.asc, kxy.y.asc, kxy.z.asc, kxz.d.asc, kxz.maximum.asc, kxz.minimum.asc, kxz.norm1.asc, kxz.norm2.asc, kxz.x.asc, kxz.y.asc, kxz.z.asc, kyy.d.asc, kyy.maximum.asc, kyy.minimum.asc, kyy.norm1.asc, kyy.norm2.asc, kyy.x.asc, kyy.y.asc, kyy.z.asc, kyz.d.asc, kyz.maximum.asc, kyz.minimum.asc, kyz.norm1.asc, kyz.norm2.asc, kyz.x.asc, kyz.y.asc, kyz.z.asc, kzz.d.asc, kzz.maximum.asc, kzz.minimum.asc, kzz.norm1.asc, kzz.norm2.asc, kzz.x.asc, kzz.y.asc, kzz.z.asc, press.d.asc, press.maximum.asc, press.minimum.asc, press.norm1.asc, press.norm2.asc, press.x.asc, press.y.asc, press.z.asc, puncture_u.d.asc, puncture_u.maximum.asc, puncture_u.minimum.asc, puncture_u.norm1.asc, puncture_u.norm2.asc, puncture_u.x.asc, puncture_u.y.asc, puncture_u.z.asc, rho.d.asc, rho.maximum.asc, rho.minimum.asc, rho.norm1.asc, rho.norm2.asc, rho.x.asc, rho.y.asc, rho.z.asc, scon[0].d.asc, scon[0].maximum.asc, scon[0].minimum.asc, scon[0].norm1.asc, scon[0].norm2.asc, scon[0].x.asc, scon[0].y.asc, scon[0].z.asc, scon[1].d.asc, scon[1].maximum.asc, scon[1].minimum.asc, scon[1].norm1.asc, scon[1].norm2.asc, scon[1].x.asc, scon[1].y.asc, scon[1].z.asc, scon[2].d.asc, scon[2].maximum.asc, scon[2].minimum.asc, scon[2].norm1.asc, scon[2].norm2.asc, scon[2].x.asc, scon[2].y.asc, scon[2].z.asc, tau.d.asc, tau.maximum.asc, tau.minimum.asc, tau.norm1.asc, tau.norm2.asc, tau.x.asc, tau.y.asc, tau.z.asc, vel[0].d.asc, vel[0].maximum.asc, vel[0].minimum.asc, vel[0].norm1.asc, vel[0].norm2.asc, vel[0].x.asc, vel[0].y.asc, vel[0].z.asc, vel[1].d.asc, vel[1].maximum.asc, vel[1].minimum.asc, vel[1].norm1.asc, vel[1].norm2.asc, vel[1].x.asc, vel[1].y.asc, vel[1].z.asc, vel[2].d.asc, vel[2].maximum.asc, vel[2].minimum.asc, vel[2].norm1.asc, vel[2].norm2.asc, vel[2].x.asc, vel[2].y.asc, vel[2].z.asc, w_lorentz.d.asc, w_lorentz.maximum.asc, w_lorentz.minimum.asc, w_lorentz.norm1.asc, w_lorentz.norm2.asc, w_lorentz.x.asc, w_lorentz.y.asc, w_lorentz.z.asc Modified: /branches/momentum_constraint/ param.ccl /branches/momentum_constraint/src/ Source.c, TwoPunctures.c, TwoPunctures.h Log: TwoPunctures: add option to compute Aij using finite differencing Adds an option use_finite_differencing_for_Aij (defaults to no) which computes the extrinsic curvature via finite differencing when putting it onto the Cactus grid. Gives smoother results than using the spectral values directly. Test suite (and example parameter file) are included in this commit. Currently terribly inefficient since it evaluates the spectral representation nine times for each Cactus grid point. File Changes: Directory: /branches/momentum_constraint/test/bhns_momentum_fd/ =============================================================== File [added]: alp.d.asc Delta lines: <Unavailable> File [added]: alp.maximum.asc Delta lines: <Unavailable> File [added]: alp.minimum.asc Delta lines: <Unavailable> File [added]: alp.norm1.asc Delta lines: <Unavailable> File [added]: alp.norm2.asc Delta lines: <Unavailable> File [added]: alp.x.asc Delta lines: <Unavailable> File [added]: alp.y.asc Delta lines: <Unavailable> File [added]: alp.z.asc Delta lines: <Unavailable> File [added]: betax.d.asc Delta lines: <Unavailable> File [added]: betax.maximum.asc Delta lines: <Unavailable> File [added]: betax.minimum.asc Delta lines: <Unavailable> File [added]: betax.norm1.asc Delta lines: <Unavailable> File [added]: betax.norm2.asc Delta lines: <Unavailable> File [added]: betax.x.asc Delta lines: <Unavailable> File [added]: betax.y.asc Delta lines: <Unavailable> File [added]: betax.z.asc Delta lines: <Unavailable> File [added]: betay.d.asc Delta lines: <Unavailable> File [added]: betay.maximum.asc Delta lines: <Unavailable> File [added]: betay.minimum.asc Delta lines: <Unavailable> File [added]: betay.norm1.asc Delta lines: <Unavailable> File [added]: betay.norm2.asc Delta lines: <Unavailable> File [added]: betay.x.asc Delta lines: <Unavailable> File [added]: betay.y.asc Delta lines: <Unavailable> File [added]: betay.z.asc Delta lines: <Unavailable> File [added]: betaz.d.asc Delta lines: <Unavailable> File [added]: betaz.maximum.asc Delta lines: <Unavailable> File [added]: betaz.minimum.asc Delta lines: <Unavailable> File [added]: betaz.norm1.asc Delta lines: <Unavailable> File [added]: betaz.norm2.asc Delta lines: <Unavailable> File [added]: betaz.x.asc Delta lines: <Unavailable> File [added]: betaz.y.asc Delta lines: <Unavailable> File [added]: betaz.z.asc Delta lines: <Unavailable> File [added]: dens.d.asc Delta lines: <Unavailable> File [added]: dens.maximum.asc Delta lines: <Unavailable> File [added]: dens.minimum.asc Delta lines: <Unavailable> File [added]: dens.norm1.asc Delta lines: <Unavailable> File [added]: dens.norm2.asc Delta lines: <Unavailable> File [added]: dens.x.asc Delta lines: <Unavailable> File [added]: dens.y.asc Delta lines: <Unavailable> File [added]: dens.z.asc Delta lines: <Unavailable> File [added]: eps.d.asc Delta lines: <Unavailable> File [added]: eps.maximum.asc Delta lines: <Unavailable> File [added]: eps.minimum.asc Delta lines: <Unavailable> File [added]: eps.norm1.asc Delta lines: <Unavailable> File [added]: eps.norm2.asc Delta lines: <Unavailable> File [added]: eps.x.asc Delta lines: <Unavailable> File [added]: eps.y.asc Delta lines: <Unavailable> File [added]: eps.z.asc Delta lines: <Unavailable> File [added]: gxx.d.asc Delta lines: <Unavailable> File [added]: gxx.maximum.asc Delta lines: <Unavailable> File [added]: gxx.minimum.asc Delta lines: <Unavailable> File [added]: gxx.norm1.asc Delta lines: <Unavailable> File [added]: gxx.norm2.asc Delta lines: <Unavailable> File [added]: gxx.x.asc Delta lines: <Unavailable> File [added]: gxx.y.asc Delta lines: <Unavailable> File [added]: gxx.z.asc Delta lines: <Unavailable> File [added]: gxy.d.asc Delta lines: <Unavailable> File [added]: gxy.maximum.asc Delta lines: <Unavailable> File [added]: gxy.minimum.asc Delta lines: <Unavailable> File [added]: gxy.norm1.asc Delta lines: <Unavailable> File [added]: gxy.norm2.asc Delta lines: <Unavailable> File [added]: gxy.x.asc Delta lines: <Unavailable> File [added]: gxy.y.asc Delta lines: <Unavailable> File [added]: gxy.z.asc Delta lines: <Unavailable> File [added]: gxz.d.asc Delta lines: <Unavailable> File [added]: gxz.maximum.asc Delta lines: <Unavailable> File [added]: gxz.minimum.asc Delta lines: <Unavailable> File [added]: gxz.norm1.asc Delta lines: <Unavailable> File [added]: gxz.norm2.asc Delta lines: <Unavailable> File [added]: gxz.x.asc Delta lines: <Unavailable> File [added]: gxz.y.asc Delta lines: <Unavailable> File [added]: gxz.z.asc Delta lines: <Unavailable> File [added]: gyy.d.asc Delta lines: <Unavailable> File [added]: gyy.maximum.asc Delta lines: <Unavailable> File [added]: gyy.minimum.asc Delta lines: <Unavailable> File [added]: gyy.norm1.asc Delta lines: <Unavailable> File [added]: gyy.norm2.asc Delta lines: <Unavailable> File [added]: gyy.x.asc Delta lines: <Unavailable> File [added]: gyy.y.asc Delta lines: <Unavailable> File [added]: gyy.z.asc Delta lines: <Unavailable> File [added]: gyz.d.asc Delta lines: <Unavailable> File [added]: gyz.maximum.asc Delta lines: <Unavailable> File [added]: gyz.minimum.asc Delta lines: <Unavailable> File [added]: gyz.norm1.asc Delta lines: <Unavailable> File [added]: gyz.norm2.asc Delta lines: <Unavailable> File [added]: gyz.x.asc Delta lines: <Unavailable> File [added]: gyz.y.asc Delta lines: <Unavailable> File [added]: gyz.z.asc Delta lines: <Unavailable> File [added]: gzz.d.asc Delta lines: <Unavailable> File [added]: gzz.maximum.asc Delta lines: <Unavailable> File [added]: gzz.minimum.asc Delta lines: <Unavailable> File [added]: gzz.norm1.asc Delta lines: <Unavailable> File [added]: gzz.norm2.asc Delta lines: <Unavailable> File [added]: gzz.x.asc Delta lines: <Unavailable> File [added]: gzz.y.asc Delta lines: <Unavailable> File [added]: gzz.z.asc Delta lines: <Unavailable> File [added]: kxx.d.asc Delta lines: <Unavailable> File [added]: kxx.maximum.asc Delta lines: <Unavailable> File [added]: kxx.minimum.asc Delta lines: <Unavailable> File [added]: kxx.norm1.asc Delta lines: <Unavailable> File [added]: kxx.norm2.asc Delta lines: <Unavailable> File [added]: kxx.x.asc Delta lines: <Unavailable> File [added]: kxx.y.asc Delta lines: <Unavailable> File [added]: kxx.z.asc Delta lines: <Unavailable> File [added]: kxy.d.asc Delta lines: <Unavailable> File [added]: kxy.maximum.asc Delta lines: <Unavailable> File [added]: kxy.minimum.asc Delta lines: <Unavailable> File [added]: kxy.norm1.asc Delta lines: <Unavailable> File [added]: kxy.norm2.asc Delta lines: <Unavailable> File [added]: kxy.x.asc Delta lines: <Unavailable> File [added]: kxy.y.asc Delta lines: <Unavailable> File [added]: kxy.z.asc Delta lines: <Unavailable> File [added]: kxz.d.asc Delta lines: <Unavailable> File [added]: kxz.maximum.asc Delta lines: <Unavailable> File [added]: kxz.minimum.asc Delta lines: <Unavailable> File [added]: kxz.norm1.asc Delta lines: <Unavailable> File [added]: kxz.norm2.asc Delta lines: <Unavailable> File [added]: kxz.x.asc Delta lines: <Unavailable> File [added]: kxz.y.asc Delta lines: <Unavailable> File [added]: kxz.z.asc Delta lines: <Unavailable> File [added]: kyy.d.asc Delta lines: <Unavailable> File [added]: kyy.maximum.asc Delta lines: <Unavailable> File [added]: kyy.minimum.asc Delta lines: <Unavailable> File [added]: kyy.norm1.asc Delta lines: <Unavailable> File [added]: kyy.norm2.asc Delta lines: <Unavailable> File [added]: kyy.x.asc Delta lines: <Unavailable> File [added]: kyy.y.asc Delta lines: <Unavailable> File [added]: kyy.z.asc Delta lines: <Unavailable> File [added]: kyz.d.asc Delta lines: <Unavailable> File [added]: kyz.maximum.asc Delta lines: <Unavailable> File [added]: kyz.minimum.asc Delta lines: <Unavailable> File [added]: kyz.norm1.asc Delta lines: <Unavailable> File [added]: kyz.norm2.asc Delta lines: <Unavailable> File [added]: kyz.x.asc Delta lines: <Unavailable> File [added]: kyz.y.asc Delta lines: <Unavailable> File [added]: kyz.z.asc Delta lines: <Unavailable> File [added]: kzz.d.asc Delta lines: <Unavailable> File [added]: kzz.maximum.asc Delta lines: <Unavailable> File [added]: kzz.minimum.asc Delta lines: <Unavailable> File [added]: kzz.norm1.asc Delta lines: <Unavailable> File [added]: kzz.norm2.asc Delta lines: <Unavailable> File [added]: kzz.x.asc Delta lines: <Unavailable> File [added]: kzz.y.asc Delta lines: <Unavailable> File [added]: kzz.z.asc Delta lines: <Unavailable> File [added]: press.d.asc Delta lines: <Unavailable> File [added]: press.maximum.asc Delta lines: <Unavailable> File [added]: press.minimum.asc Delta lines: <Unavailable> File [added]: press.norm1.asc Delta lines: <Unavailable> File [added]: press.norm2.asc Delta lines: <Unavailable> File [added]: press.x.asc Delta lines: <Unavailable> File [added]: press.y.asc Delta lines: <Unavailable> File [added]: press.z.asc Delta lines: <Unavailable> File [added]: puncture_u.d.asc Delta lines: <Unavailable> File [added]: puncture_u.maximum.asc Delta lines: <Unavailable> File [added]: puncture_u.minimum.asc Delta lines: <Unavailable> File [added]: puncture_u.norm1.asc Delta lines: <Unavailable> File [added]: puncture_u.norm2.asc Delta lines: <Unavailable> File [added]: puncture_u.x.asc Delta lines: <Unavailable> File [added]: puncture_u.y.asc Delta lines: <Unavailable> File [added]: puncture_u.z.asc Delta lines: <Unavailable> File [added]: rho.d.asc Delta lines: <Unavailable> File [added]: rho.maximum.asc Delta lines: <Unavailable> File [added]: rho.minimum.asc Delta lines: <Unavailable> File [added]: rho.norm1.asc Delta lines: <Unavailable> File [added]: rho.norm2.asc Delta lines: <Unavailable> File [added]: rho.x.asc Delta lines: <Unavailable> File [added]: rho.y.asc Delta lines: <Unavailable> File [added]: rho.z.asc Delta lines: <Unavailable> File [added]: scon[0].d.asc Delta lines: <Unavailable> File [added]: scon[0].maximum.asc Delta lines: <Unavailable> File [added]: scon[0].minimum.asc Delta lines: <Unavailable> File [added]: scon[0].norm1.asc Delta lines: <Unavailable> File [added]: scon[0].norm2.asc Delta lines: <Unavailable> File [added]: scon[0].x.asc Delta lines: <Unavailable> File [added]: scon[0].y.asc Delta lines: <Unavailable> File [added]: scon[0].z.asc Delta lines: <Unavailable> File [added]: scon[1].d.asc Delta lines: <Unavailable> File [added]: scon[1].maximum.asc Delta lines: <Unavailable> File [added]: scon[1].minimum.asc Delta lines: <Unavailable> File [added]: scon[1].norm1.asc Delta lines: <Unavailable> File [added]: scon[1].norm2.asc Delta lines: <Unavailable> File [added]: scon[1].x.asc Delta lines: <Unavailable> File [added]: scon[1].y.asc Delta lines: <Unavailable> File [added]: scon[1].z.asc Delta lines: <Unavailable> File [added]: scon[2].d.asc Delta lines: <Unavailable> File [added]: scon[2].maximum.asc Delta lines: <Unavailable> File [added]: scon[2].minimum.asc Delta lines: <Unavailable> File [added]: scon[2].norm1.asc Delta lines: <Unavailable> File [added]: scon[2].norm2.asc Delta lines: <Unavailable> File [added]: scon[2].x.asc Delta lines: <Unavailable> File [added]: scon[2].y.asc Delta lines: <Unavailable> File [added]: scon[2].z.asc Delta lines: <Unavailable> File [added]: tau.d.asc Delta lines: <Unavailable> File [added]: tau.maximum.asc Delta lines: <Unavailable> File [added]: tau.minimum.asc Delta lines: <Unavailable> File [added]: tau.norm1.asc Delta lines: <Unavailable> File [added]: tau.norm2.asc Delta lines: <Unavailable> File [added]: tau.x.asc Delta lines: <Unavailable> File [added]: tau.y.asc Delta lines: <Unavailable> File [added]: tau.z.asc Delta lines: <Unavailable> File [added]: vel[0].d.asc Delta lines: <Unavailable> File [added]: vel[0].maximum.asc Delta lines: <Unavailable> File [added]: vel[0].minimum.asc Delta lines: <Unavailable> File [added]: vel[0].norm1.asc Delta lines: <Unavailable> File [added]: vel[0].norm2.asc Delta lines: <Unavailable> File [added]: vel[0].x.asc Delta lines: <Unavailable> File [added]: vel[0].y.asc Delta lines: <Unavailable> File [added]: vel[0].z.asc Delta lines: <Unavailable> File [added]: vel[1].d.asc Delta lines: <Unavailable> File [added]: vel[1].maximum.asc Delta lines: <Unavailable> File [added]: vel[1].minimum.asc Delta lines: <Unavailable> File [added]: vel[1].norm1.asc Delta lines: <Unavailable> File [added]: vel[1].norm2.asc Delta lines: <Unavailable> File [added]: vel[1].x.asc Delta lines: <Unavailable> File [added]: vel[1].y.asc Delta lines: <Unavailable> File [added]: vel[1].z.asc Delta lines: <Unavailable> File [added]: vel[2].d.asc Delta lines: <Unavailable> File [added]: vel[2].maximum.asc Delta lines: <Unavailable> File [added]: vel[2].minimum.asc Delta lines: <Unavailable> File [added]: vel[2].norm1.asc Delta lines: <Unavailable> File [added]: vel[2].norm2.asc Delta lines: <Unavailable> File [added]: vel[2].x.asc Delta lines: <Unavailable> File [added]: vel[2].y.asc Delta lines: <Unavailable> File [added]: vel[2].z.asc Delta lines: <Unavailable> File [added]: w_lorentz.d.asc Delta lines: <Unavailable> File [added]: w_lorentz.maximum.asc Delta lines: <Unavailable> File [added]: w_lorentz.minimum.asc Delta lines: <Unavailable> File [added]: w_lorentz.norm1.asc Delta lines: <Unavailable> File [added]: w_lorentz.norm2.asc Delta lines: <Unavailable> File [added]: w_lorentz.x.asc Delta lines: <Unavailable> File [added]: w_lorentz.y.asc Delta lines: <Unavailable> File [added]: w_lorentz.z.asc Delta lines: <Unavailable> Directory: /branches/momentum_constraint/ ========================================= File [modified]: param.ccl Delta lines: <Unavailable> Directory: /branches/momentum_constraint/src/ ============================================= File [modified]: Source.c Delta lines: <Unavailable> File [modified]: TwoPunctures.c Delta lines: <Unavailable> File [modified]: TwoPunctures.h Delta lines: <Unavailable> Directory: /branches/momentum_constraint/test/ ============================================== File [added]: bhns_momentum_fd.par Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] www/ (Rev. 316)
by gallen＠cct.lsu.edu 15 Jul '10

15 Jul '10

User: gallen Date: 2010/07/15 06:45 PM Added: /about/ funding.php Log: Just prototyping a page for comments Directory Changes: Directory: /svn:executable/ =========================== + * File Changes: Directory: /about/ ================== File [added]: funding.php Delta lines: +23 -0 =================================================================== --- about/funding.php (rev 0) +++ about/funding.php 2010-07-15 23:45:09 UTC (rev 316) @@ -0,0 +1,23 @@ +<?php $title='Funding'; +include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> + +The development and support of the Einstein Toolkit components and packages have been funded from multiple sources, some of which are listed below. Please + write to the <a href="mailto:maintainers@einsteintoolkit.org">maintainers</a> to add an item to this page. + +<ul> + +<li> NSF #0904015: <a href="">Collaborative Research: Community Infrastructure for General Relativity MHD (CIGR)</a> +Core development and support of the community software for general relativistic MHD</li> + +<li> NSF #1005165: <a href="">Interdisciplinary Research Experience in Computational Sciences</a> + Development of GetComponents for code assembly</li> + +<li> NSF #0932251: <a href="http://www.nsf.gov/awardsearch/showAward.do?AwardNumber=0932251">TeraGrid Extension: Bridging to XD</a> +Development of SimFactory into a science gateway for TeraGrid</li> + + + +</ul> + + +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?> Property changes on: about/funding.php ___________________________________________________________________

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 315)
by gallen＠cct.lsu.edu 15 Jul '10

15 Jul '10

User: gallen Date: 2010/07/15 04:58 PM Modified: /about/members/ index.php Log: Billy Vazquez and Mew Bing Wan File Changes: Directory: /about/members/ ========================== File [modified]: index.php Delta lines: +6 -0 =================================================================== --- about/members/index.php 2010-07-10 14:10:23 UTC (rev 314) +++ about/members/index.php 2010-07-15 21:58:54 UTC (rev 315) @@ -34,6 +34,11 @@ <li>Bela Szilagyi</li> </ul> + <li>Chinese Academy of Sciences</li> + <ul> + <li>Mew Bing Wan</li> + </ul> + <li><a href="http://www.cct.lsu.edu/numerical/" class="nostyle">Louisiana State University</a></li> <ul> <li>Ernazar Abdikamalov</li> @@ -86,6 +91,7 @@ <li>Scott Noble</li> <li>Marcelo Ponce</li> <li>David Saroff</li> + <li>Billy Vazquez</li> <li>Yosef Zlochower</li> </ul>

1 0

[svn:einsteintoolkit] TOVSolver/trunk/src/ (Rev. 117)
by knarf＠cct.lsu.edu 13 Jul '10

13 Jul '10

User: knarf Date: 2010/07/13 10:47 AM Modified: /trunk/src/ external.inc, tov.c Log: fix compiler warnings File Changes: Directory: /trunk/src/ ====================== File [modified]: external.inc Delta lines: +6 -7 =================================================================== --- trunk/src/external.inc 2010-07-13 15:41:27 UTC (rev 116) +++ trunk/src/external.inc 2010-07-13 15:47:07 UTC (rev 117) @@ -367,10 +367,10 @@ *vel0_p_p, *vel1_p_p, *vel2_p_p, *scon0_p_p, *scon1_p_p, *scon2_p_p; CCTK_REAL *u, *source, *mom_source[3]; CCTK_REAL psip, my_psi4, sqrt_det, vx, vy, vz, v_2, - w_lorentz_2, d_w_lorentz_2, D_h_w; + w_lorentz_2, D_h_w; CCTK_REAL vlowx, vlowy, vlowz; - CCTK_INT i, type; - CCTK_REAL rhoold, rhonew, f, df; + CCTK_INT type; + CCTK_REAL rhoold, rhonew; int i3D; static int debug = 1; FILE *debugfile = NULL; @@ -440,7 +440,6 @@ { if (source[i3D]>0.0) { - int count=0; if (debug && (fabs(z[i3D])<1.e-15)) { fprintf(debugfile, @@ -464,8 +463,10 @@ rho[i3D] = 1.e-10; rhoold = rhonew = rho[i3D]; #ifdef OLD_NR + int count=0; do { + CCTK_REAL f, df; count++; vx = mom_source[0][i3D]/my_psi4/ (source[i3D]/my_psi4/my_psi4+TOV_K[0] * pow(rhonew, TOV_Gamma[0])); @@ -498,7 +499,7 @@ TOV_K[0] * TOV_Gamma[0] * pow(rhonew, TOV_Gamma[0]-1.); */ /* d_w_lorentz_2/drhonew */ - d_w_lorentz_2 = + CCTK_REAL d_w_lorentz_2 = -2 * w_lorentz_2*w_lorentz_2 * v_2 * TOV_K[0]*TOV_Gamma[0]*pow(rhonew, TOV_Gamma[0]-1) / (source[i3D]/my_psi4/my_psi4 + TOV_K[0]*pow(rhonew, TOV_Gamma[0])); @@ -647,8 +648,6 @@ void TOVSolverC_export_local_variables(CCTK_REAL **exported_TOV_press_1d, CCTK_REAL **exported_TOV_m_1d, CCTK_REAL **exported_TOV_phi_1d, CCTK_REAL **exported_TOV_rbar_1d, CCTK_REAL **exported_TOV_r_1d) { - int i; - *exported_TOV_press_1d = TOV_press_1d; *exported_TOV_m_1d = TOV_m_1d; *exported_TOV_phi_1d = TOV_phi_1d; File [modified]: tov.c Delta lines: +4 -4 =================================================================== --- trunk/src/tov.c 2010-07-13 15:41:27 UTC (rev 116) +++ trunk/src/tov.c 2010-07-13 15:47:07 UTC (rev 117) @@ -149,7 +149,7 @@ r_minus_two_m = r - 2.0 * m; - if ((r==0.0) && (m==0.0)) + if ((r<=0.0) && (m<=0.0)) { source_data[1] = 0.0; source_data[2] = 0.0; @@ -219,7 +219,7 @@ star_i = star * TOV_Num_Radial; /* check for parameters */ - if ((TOV_Rho_Central[star]==0.0) && (GRHydro_rho_central>0.0)) + if ((TOV_Rho_Central[star]<=0.0) && (GRHydro_rho_central>0.0)) rho_central=GRHydro_rho_central; else rho_central=TOV_Rho_Central[star]; @@ -283,7 +283,7 @@ old_data[0] = TOV_press_1d[i]; old_data[1] = TOV_m_1d[i]; old_data[2] = TOV_phi_1d[i]; - if (TOV_rbar_1d[i]==TOV_r_1d[i]) + if (fabs(TOV_rbar_1d[i] - TOV_r_1d[i]) < LOCAL_TINY) old_data[3] = 0.0; else old_data[3] = log(TOV_rbar_1d[i] / TOV_r_1d[i]); @@ -377,7 +377,7 @@ CCTK_VInfo(CCTK_THORNSTRING, "Information about the TOVs used:"); CCTK_VInfo("", "TOV radius mass bary_mass mass(g) cent.rho rho(cgi) K K(cgi) Gamma"); for (i=0; i<TOV_Num_TOVs; i++) - if (TOV_Gamma[i]==2.0) + if (fabs(TOV_Gamma[i] - 2.0) < LOCAL_TINY) CCTK_VInfo(""," %d %8g %8g %8g %8.3g %8g %8.3g %8g %8.3g %8g", i+1, TOV_R_Surface[i], TOV_m_1d[(i+1)*TOV_Num_Radial-1],

1 0

[svn:einsteintoolkit] TOVSolver/trunk/src/ (Rev. 116)
by knarf＠cct.lsu.edu 13 Jul '10

13 Jul '10

User: knarf Date: 2010/07/13 10:41 AM Modified: /trunk/src/ external.inc Log: combine error messages File Changes: Directory: /trunk/src/ ====================== File [modified]: external.inc Delta lines: +3 -3 =================================================================== --- trunk/src/external.inc 2010-07-12 14:50:50 UTC (rev 115) +++ trunk/src/external.inc 2010-07-13 15:41:27 UTC (rev 116) @@ -308,11 +308,11 @@ if (count > 100) { CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, \ - "Failed to converge IVP rescaling\nat (%.3g %.3g %.3g)\n"\ - " orig_rho:%g rho:%g, v^2=%g, rel_diff=%g", \ + "Failed to converge IVP rescaling, trying slow bisection method\n"\ + " at (%.3g %.3g %.3g)\n"\ + " orig_rho:%g rho:%g, v^2=%g, rel_diff=%g", \ x, y, z, rho_orig, *rhonew, *v_2, fabs(*rhonew-*rhoold)/fabs(*rhonew)); - CCTK_WARN(1, "trying slow bisection method"); bisection(vx,vy,vz,rhoold,rhonew,w_lorentz_2,v_2, Gamma, K, source, mom_source_x, mom_source_y, mom_source_z,

1 0

[svn:einsteintoolkit] TGRtensor/tags/ET_2010_06_v0/ (Rev. 54)
by knarf＠cct.lsu.edu 12 Jul '10

12 Jul '10