Commits September 2010

commits@lists.einsteintoolkit.org

7 participants
24 discussions

[svn:einsteintoolkit] www/ (Rev. 323)
by cjordan＠cct.lsu.edu 16 Sep '10

16 Sep '10

User: cjordan Date: 2010/09/16 08:24 PM Modified: /info/documentation/ send-account.php /toolkit/community/ index.php Log: corrected a slight error in the mail list links that caused the mailman interface to not show list listings File Changes: Directory: /toolkit/community/ ============================== File [modified]: index.php Delta lines: +3 -3 =================================================================== --- toolkit/community/index.php 2010-09-02 18:09:01 UTC (rev 322) +++ toolkit/community/index.php 2010-09-17 01:24:07 UTC (rev 323) @@ -7,17 +7,17 @@ <li> <a href="mailto:users@einsteintoolkit.org">users(a)einsteintoolkit.org</a>: Primary list for the Einstein Toolkit including user questions and developer discussions. - [<a href="http://mail.einsteintoolkit.org/mailman/listinfo/users">Subscribe</a>][<a href="http://lists.einsteintoolkit.org/pipermail/users/">Archive</a>] + [<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users">Subscribe</a>][<a href="http://lists.einsteintoolkit.org/pipermail/users/">Archive</a>] <li> <a href="mailto:commits@einsteintoolkit.org">commits(a)einsteintoolkit.org</a>: All commit messages to the Einstein Toolkit SVN are posted to this list. - [<a href="http://mail.einsteintoolkit.org/mailman/listinfo/commits">Subscribe</a>][<a href="http://lists.einsteintoolkit.org/pipermail/commits/">Archive</a>] + [<a href="http://lists.einsteintoolkit.org/mailman/listinfo/commits">Subscribe</a>][<a href="http://lists.einsteintoolkit.org/pipermail/commits/">Archive</a>] <li> <a href="mailto:users@einsteintoolkit.org">announce(a)einsteintoolkit.org</a>: For broadcasting general news and information about the Einstein Toolkit for those not wanting to see the traffic on the users mail list. - [<a href="http://mail.einsteintoolkit.org/mailman/listinfo/announce">Subscribe</a>][<a href="http://lists.einsteintoolkit.org/pipermail/announce/">Archive</a>] + [<a href="http://lists.einsteintoolkit.org/mailman/listinfo/announce">Subscribe</a>][<a href="http://lists.einsteintoolkit.org/pipermail/announce/">Archive</a>] </ul> Directory: /info/documentation/ =============================== File [modified]: send-account.php Delta lines: +2 -2 =================================================================== --- info/documentation/send-account.php 2010-09-02 18:09:01 UTC (rev 322) +++ info/documentation/send-account.php 2010-09-17 01:24:07 UTC (rev 323) @@ -32,7 +32,7 @@ $institution = 'at '.$institution.', '; } -$message = "Einstein Toolkit maintainers: \n\n".$name.", ".$position.$department.$institution."has submitted a request for a test account for the Einstein Toolkit. His or her email address is ".$email.". ".$mentionadvisor." \n\n Thanks,\n Einstein Toolkit Request Bot\n"; +$message = "Einstein Toolkit maintainers: \n\n".$name.", ".$position.$department.$institution."has submitted a request for a test account for the Einstein Toolkit. Her or his email address is ".$email.". ".$mentionadvisor." \n\n Thanks,\n Einstein Toolkit Request Bot\n"; /* PHP form validation: the script checks that the Email field contains a valid email address and the Subject field isn't empty. preg_match performs a regular expression match. It's a very powerful PHP function to validate form fields and other strings - see PHP manual for details. */ if (empty($name)) { @@ -56,4 +56,4 @@ -<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?> \ No newline at end of file +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?>

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 156)
by knarf＠cct.lsu.edu 15 Sep '10

15 Sep '10

User: knarf Date: 2010/09/15 01:50 PM Modified: /trunk/ param.ccl, schedule.ccl /trunk/src/ GRHydro_UpdateMask.F90 /trunk/test/ GRHydro_test_tov_ppm.par, GRHydro_test_tov_ppm_ML.par, GRHydro_test_tov_ppm_ML_disable_internal_excision.par, GRHydro_test_tov_ppm_ML_eos_omni.par, GRHydro_test_tov_ppm_disable_internal_excision.par, GRHydro_test_tov_ppm_no_trp.par /trunk/test/GRHydro_test_tov_ppm_eosg/ alp_maximum.xg, alp_norm1.xg, alp_norm2.xg, dens_3D_diagonal.xg, dens_minimum.xg, dens_norm1.xg, dens_norm2.xg, dens_x_[16][16].xg, dens_y_[16][16].xg, dens_z_[16][16].xg, eps_3D_diagonal.xg, eps_minimum.xg, eps_norm1.xg, eps_norm2.xg, eps_x_[16][16].xg, eps_y_[16][16].xg, eps_z_[16][16].xg, gxx_minimum.xg, gxx_norm1.xg, gxx_norm2.xg, gxy_maximum.xg, gxy_minimum.xg, gxy_norm1.xg, gxy_norm2.xg, gxz_maximum.xg, gxz_minimum.xg, gxz_norm1.xg, gxz_norm2.xg, gyy_minimum.xg, gyy_norm1.xg, gyy_norm2.xg, gyz_maximum.xg, gyz_minimum.xg, gyz_norm1.xg, gyz_norm2.xg, gzz_minimum.xg, gzz_norm1.xg, gzz_norm2.xg, ham_3D_diagonal.xg, ham_maximum.xg, ham_norm1.xg, ham_norm2.xg, ham_x_[16][16].xg, ham_y_[16][16].xg, ham_z_[16][16].xg, kxx_minimum.xg, kxx_norm1.xg, kxx_norm2.xg, kxy_maximum.xg, kxy_minimum.xg, kxy_norm1.xg, kxy_norm2.xg, kxz_maximum.xg, kxz_minimum.xg, kxz_norm1.xg, kxz_norm2.xg, kyy_minimum.xg, kyy_norm1.xg, kyy_norm2.xg, kyz_maximum.xg, kyz_minimum.xg, kyz_norm1.xg, kyz_norm2.xg, kzz_minimum.xg, kzz_norm1.xg, kzz_norm2.xg, momx_3D_diagonal.xg, momx_maximum.xg, momx_minimum.xg, momx_norm1.xg, momx_norm2.xg, momx_x_[16][16].xg, momx_y_[16][16].xg, momx_z_[16][16].xg, momy_3D_diagonal.xg, momy_maximum.xg, momy_minimum.xg, momy_norm1.xg, momy_norm2.xg, momy_x_[16][16].xg, momy_y_[16][16].xg, momy_z_[16][16].xg, momz_3D_diagonal.xg, momz_maximum.xg, momz_minimum.xg, momz_norm1.xg, momz_norm2.xg, momz_x_[16][16].xg, momz_y_[16][16].xg, momz_z_[16][16].xg, press_3D_diagonal.xg, press_minimum.xg, press_norm1.xg, press_x_[16][16].xg, press_y_[16][16].xg, press_z_[16][16].xg, rho_3D_diagonal.xg, rho_norm1.xg, rho_norm2.xg, rho_x_[16][16].xg, rho_y_[16][16].xg, rho_z_[16][16].xg, scon[0]_3D_diagonal.xg, scon[0]_norm1.xg, scon[0]_norm2.xg, scon[0]_x_[16][16].xg, scon[0]_y_[16][16].xg, scon[0]_z_[16][16].xg, scon[1]_3D_diagonal.xg, scon[1]_norm1.xg, scon[1]_norm2.xg, scon[1]_x_[16][16].xg, scon[1]_y_[16][16].xg, scon[1]_z_[16][16].xg, scon[2]_3D_diagonal.xg, scon[2]_norm1.xg, scon[2]_norm2.xg, scon[2]_x_[16][16].xg, scon[2]_y_[16][16].xg, scon[2]_z_[16][16].xg, tau_3D_diagonal.xg, tau_minimum.xg, tau_norm1.xg, tau_x_[16][16].xg, tau_y_[16][16].xg, tau_z_[16][16].xg, vel[0]_3D_diagonal.xg, vel[0]_maximum.xg, vel[0]_minimum.xg, vel[0]_norm1.xg, vel[0]_norm2.xg, vel[0]_x_[16][16].xg, vel[0]_y_[16][16].xg, vel[0]_z_[16][16].xg, vel[1]_3D_diagonal.xg, vel[1]_maximum.xg, vel[1]_minimum.xg, vel[1]_norm1.xg, vel[1]_norm2.xg, vel[1]_x_[16][16].xg, vel[1]_y_[16][16].xg, vel[1]_z_[16][16].xg, vel[2]_3D_diagonal.xg, vel[2]_maximum.xg, vel[2]_minimum.xg, vel[2]_norm1.xg, vel[2]_norm2.xg, vel[2]_x_[16][16].xg, vel[2]_y_[16][16].xg, vel[2]_z_[16][16].xg, w_lorentz_3D_diagonal.xg, w_lorentz_maximum.xg, w_lorentz_norm1.xg, w_lorentz_norm2.xg, w_lorentz_x_[16][16].xg, w_lorentz_y_[16][16].xg, w_lorentz_z_[16][16].xg /trunk/test/GRHydro_test_tov_ppm_eosg_ML/ alp_maximum.xg, alp_norm1.xg, alp_norm2.xg, dens_3D_diagonal.xg, dens_minimum.xg, dens_norm1.xg, dens_norm2.xg, dens_x_[16][16].xg, dens_y_[16][16].xg, dens_z_[16][16].xg, eps_3D_diagonal.xg, eps_minimum.xg, eps_norm1.xg, eps_norm2.xg, eps_x_[16][16].xg, eps_y_[16][16].xg, eps_z_[16][16].xg, gxx_minimum.xg, gxx_norm1.xg, gxx_norm2.xg, gxy_maximum.xg, gxy_minimum.xg, gxy_norm1.xg, gxy_norm2.xg, gxz_maximum.xg, gxz_minimum.xg, gxz_norm1.xg, gxz_norm2.xg, gyy_minimum.xg, gyy_norm1.xg, gyy_norm2.xg, gyz_maximum.xg, gyz_minimum.xg, gyz_norm1.xg, gyz_norm2.xg, gzz_minimum.xg, gzz_norm1.xg, gzz_norm2.xg, ham_3D_diagonal.xg, ham_maximum.xg, ham_norm1.xg, ham_norm2.xg, ham_x_[16][16].xg, ham_y_[16][16].xg, ham_z_[16][16].xg, kxx_minimum.xg, kxx_norm1.xg, kxx_norm2.xg, kxy_maximum.xg, kxy_minimum.xg, kxy_norm1.xg, kxy_norm2.xg, kxz_maximum.xg, kxz_minimum.xg, kxz_norm1.xg, kxz_norm2.xg, kyy_minimum.xg, kyy_norm1.xg, kyy_norm2.xg, kyz_maximum.xg, kyz_minimum.xg, kyz_norm1.xg, kyz_norm2.xg, kzz_minimum.xg, kzz_norm1.xg, kzz_norm2.xg, momx_3D_diagonal.xg, momx_maximum.xg, momx_minimum.xg, momx_norm1.xg, momx_norm2.xg, momx_x_[16][16].xg, momx_y_[16][16].xg, momx_z_[16][16].xg, momy_3D_diagonal.xg, momy_maximum.xg, momy_minimum.xg, momy_norm1.xg, momy_norm2.xg, momy_x_[16][16].xg, momy_y_[16][16].xg, momy_z_[16][16].xg, momz_3D_diagonal.xg, momz_maximum.xg, momz_minimum.xg, momz_norm1.xg, momz_norm2.xg, momz_x_[16][16].xg, momz_y_[16][16].xg, momz_z_[16][16].xg, press_3D_diagonal.xg, press_minimum.xg, press_norm1.xg, press_x_[16][16].xg, press_y_[16][16].xg, press_z_[16][16].xg, rho_3D_diagonal.xg, rho_norm1.xg, rho_norm2.xg, rho_x_[16][16].xg, rho_y_[16][16].xg, rho_z_[16][16].xg, scon[0]_3D_diagonal.xg, scon[0]_norm1.xg, scon[0]_norm2.xg, scon[0]_x_[16][16].xg, scon[0]_y_[16][16].xg, scon[0]_z_[16][16].xg, scon[1]_3D_diagonal.xg, scon[1]_norm1.xg, scon[1]_norm2.xg, scon[1]_x_[16][16].xg, scon[1]_y_[16][16].xg, scon[1]_z_[16][16].xg, scon[2]_3D_diagonal.xg, scon[2]_norm1.xg, scon[2]_norm2.xg, scon[2]_x_[16][16].xg, scon[2]_y_[16][16].xg, scon[2]_z_[16][16].xg, tau_3D_diagonal.xg, tau_minimum.xg, tau_norm1.xg, tau_x_[16][16].xg, tau_y_[16][16].xg, tau_z_[16][16].xg, vel[0]_3D_diagonal.xg, vel[0]_maximum.xg, vel[0]_minimum.xg, vel[0]_norm1.xg, vel[0]_norm2.xg, vel[0]_x_[16][16].xg, vel[0]_y_[16][16].xg, vel[0]_z_[16][16].xg, vel[1]_3D_diagonal.xg, vel[1]_maximum.xg, vel[1]_minimum.xg, vel[1]_norm1.xg, vel[1]_norm2.xg, vel[1]_x_[16][16].xg, vel[1]_y_[16][16].xg, vel[1]_z_[16][16].xg, vel[2]_3D_diagonal.xg, vel[2]_maximum.xg, vel[2]_minimum.xg, vel[2]_norm1.xg, vel[2]_norm2.xg, vel[2]_x_[16][16].xg, vel[2]_y_[16][16].xg, vel[2]_z_[16][16].xg, w_lorentz_3D_diagonal.xg, w_lorentz_maximum.xg, w_lorentz_norm1.xg, w_lorentz_norm2.xg, w_lorentz_x_[16][16].xg, w_lorentz_y_[16][16].xg, w_lorentz_z_[16][16].xg Log: let GRHydro set the atmosphere of initial data, so that ID thorns don't have to worry about it. This changes the data slitely for EOSG testsuites because this now uses EOSG calls instead of the otherwise fixed calculations of the TOVSolver, thus the change of two testsuites File Changes: Directory: /trunk/test/GRHydro_test_tov_ppm_eosg/ ================================================= File [modified]: alp_maximum.xg Delta lines: <Unavailable> File [modified]: alp_norm1.xg Delta lines: <Unavailable> File [modified]: alp_norm2.xg Delta lines: <Unavailable> File [modified]: dens_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: dens_minimum.xg Delta lines: <Unavailable> File [modified]: dens_norm1.xg Delta lines: <Unavailable> File [modified]: dens_norm2.xg Delta lines: <Unavailable> File [modified]: dens_x_[16][16].xg Delta lines: <Unavailable> File [modified]: dens_y_[16][16].xg Delta lines: <Unavailable> File [modified]: dens_z_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: eps_minimum.xg Delta lines: <Unavailable> File [modified]: eps_norm1.xg Delta lines: <Unavailable> File [modified]: eps_norm2.xg Delta lines: <Unavailable> File [modified]: eps_x_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_y_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_z_[16][16].xg Delta lines: <Unavailable> File [modified]: gxx_minimum.xg Delta lines: <Unavailable> File [modified]: gxx_norm1.xg Delta lines: <Unavailable> File [modified]: gxx_norm2.xg Delta lines: <Unavailable> File [modified]: gxy_maximum.xg Delta lines: <Unavailable> File [modified]: gxy_minimum.xg Delta lines: <Unavailable> File [modified]: gxy_norm1.xg Delta lines: <Unavailable> File [modified]: gxy_norm2.xg Delta lines: <Unavailable> File [modified]: gxz_maximum.xg Delta lines: <Unavailable> File [modified]: gxz_minimum.xg Delta lines: <Unavailable> File [modified]: gxz_norm1.xg Delta lines: <Unavailable> File [modified]: gxz_norm2.xg Delta lines: <Unavailable> File [modified]: gyy_minimum.xg Delta lines: <Unavailable> File [modified]: gyy_norm1.xg Delta lines: <Unavailable> File [modified]: gyy_norm2.xg Delta lines: <Unavailable> File [modified]: gyz_maximum.xg Delta lines: <Unavailable> File [modified]: gyz_minimum.xg Delta lines: <Unavailable> File [modified]: gyz_norm1.xg Delta lines: <Unavailable> File [modified]: gyz_norm2.xg Delta lines: <Unavailable> File [modified]: gzz_minimum.xg Delta lines: <Unavailable> File [modified]: gzz_norm1.xg Delta lines: <Unavailable> File [modified]: gzz_norm2.xg Delta lines: <Unavailable> File [modified]: ham_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: ham_maximum.xg Delta lines: <Unavailable> File [modified]: ham_norm1.xg Delta lines: <Unavailable> File [modified]: ham_norm2.xg Delta lines: <Unavailable> File [modified]: ham_x_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_y_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_z_[16][16].xg Delta lines: <Unavailable> File [modified]: kxx_minimum.xg Delta lines: <Unavailable> File [modified]: kxx_norm1.xg Delta lines: <Unavailable> File [modified]: kxx_norm2.xg Delta lines: <Unavailable> File [modified]: kxy_maximum.xg Delta lines: <Unavailable> File [modified]: kxy_minimum.xg Delta lines: <Unavailable> File [modified]: kxy_norm1.xg Delta lines: <Unavailable> File [modified]: kxy_norm2.xg Delta lines: <Unavailable> File [modified]: kxz_maximum.xg Delta lines: <Unavailable> File [modified]: kxz_minimum.xg Delta lines: <Unavailable> File [modified]: kxz_norm1.xg Delta lines: <Unavailable> File [modified]: kxz_norm2.xg Delta lines: <Unavailable> File [modified]: kyy_minimum.xg Delta lines: <Unavailable> File [modified]: kyy_norm1.xg Delta lines: <Unavailable> File [modified]: kyy_norm2.xg Delta lines: <Unavailable> File [modified]: kyz_maximum.xg Delta lines: <Unavailable> File [modified]: kyz_minimum.xg Delta lines: <Unavailable> File [modified]: kyz_norm1.xg Delta lines: <Unavailable> File [modified]: kyz_norm2.xg Delta lines: <Unavailable> File [modified]: kzz_minimum.xg Delta lines: <Unavailable> File [modified]: kzz_norm1.xg Delta lines: <Unavailable> File [modified]: kzz_norm2.xg Delta lines: <Unavailable> File [modified]: momx_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: momx_maximum.xg Delta lines: <Unavailable> File [modified]: momx_minimum.xg Delta lines: <Unavailable> File [modified]: momx_norm1.xg Delta lines: <Unavailable> File [modified]: momx_norm2.xg Delta lines: <Unavailable> File [modified]: momx_x_[16][16].xg Delta lines: <Unavailable> File [modified]: momx_y_[16][16].xg Delta lines: <Unavailable> File [modified]: momx_z_[16][16].xg Delta lines: <Unavailable> File [modified]: momy_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: momy_maximum.xg Delta lines: <Unavailable> File [modified]: momy_minimum.xg Delta lines: <Unavailable> File [modified]: momy_norm1.xg Delta lines: <Unavailable> File [modified]: momy_norm2.xg Delta lines: <Unavailable> File [modified]: momy_x_[16][16].xg Delta lines: <Unavailable> File [modified]: momy_y_[16][16].xg Delta lines: <Unavailable> File [modified]: momy_z_[16][16].xg Delta lines: <Unavailable> File [modified]: momz_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: momz_maximum.xg Delta lines: <Unavailable> File [modified]: momz_minimum.xg Delta lines: <Unavailable> File [modified]: momz_norm1.xg Delta lines: <Unavailable> File [modified]: momz_norm2.xg Delta lines: <Unavailable> File [modified]: momz_x_[16][16].xg Delta lines: <Unavailable> File [modified]: momz_y_[16][16].xg Delta lines: <Unavailable> File [modified]: momz_z_[16][16].xg Delta lines: <Unavailable> File [modified]: press_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: press_minimum.xg Delta lines: <Unavailable> File [modified]: press_norm1.xg Delta lines: <Unavailable> File [modified]: press_x_[16][16].xg Delta lines: <Unavailable> File [modified]: press_y_[16][16].xg Delta lines: <Unavailable> File [modified]: press_z_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: rho_norm1.xg Delta lines: <Unavailable> File [modified]: rho_norm2.xg Delta lines: <Unavailable> File [modified]: rho_x_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_y_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_z_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[0]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[0]_norm1.xg Delta lines: <Unavailable> File [modified]: scon[0]_norm2.xg Delta lines: <Unavailable> File [modified]: scon[0]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[0]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[0]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[1]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[1]_norm1.xg Delta lines: <Unavailable> File [modified]: scon[1]_norm2.xg Delta lines: <Unavailable> File [modified]: scon[1]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[1]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[1]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[2]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[2]_norm1.xg Delta lines: <Unavailable> File [modified]: scon[2]_norm2.xg Delta lines: <Unavailable> File [modified]: scon[2]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[2]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[2]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: tau_minimum.xg Delta lines: <Unavailable> File [modified]: tau_norm1.xg Delta lines: <Unavailable> File [modified]: tau_x_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_y_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_z_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[0]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[0]_maximum.xg Delta lines: <Unavailable> File [modified]: vel[0]_minimum.xg Delta lines: <Unavailable> File [modified]: vel[0]_norm1.xg Delta lines: <Unavailable> File [modified]: vel[0]_norm2.xg Delta lines: <Unavailable> File [modified]: vel[0]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[0]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[0]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[1]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[1]_maximum.xg Delta lines: <Unavailable> File [modified]: vel[1]_minimum.xg Delta lines: <Unavailable> File [modified]: vel[1]_norm1.xg Delta lines: <Unavailable> File [modified]: vel[1]_norm2.xg Delta lines: <Unavailable> File [modified]: vel[1]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[1]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[1]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[2]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[2]_maximum.xg Delta lines: <Unavailable> File [modified]: vel[2]_minimum.xg Delta lines: <Unavailable> File [modified]: vel[2]_norm1.xg Delta lines: <Unavailable> File [modified]: vel[2]_norm2.xg Delta lines: <Unavailable> File [modified]: vel[2]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[2]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[2]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: w_lorentz_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: w_lorentz_maximum.xg Delta lines: <Unavailable> File [modified]: w_lorentz_norm1.xg Delta lines: <Unavailable> File [modified]: w_lorentz_norm2.xg Delta lines: <Unavailable> File [modified]: w_lorentz_x_[16][16].xg Delta lines: <Unavailable> File [modified]: w_lorentz_y_[16][16].xg Delta lines: <Unavailable> File [modified]: w_lorentz_z_[16][16].xg Delta lines: <Unavailable> Directory: /trunk/test/GRHydro_test_tov_ppm_eosg_ML/ ==================================================== File [modified]: alp_maximum.xg Delta lines: <Unavailable> File [modified]: alp_norm1.xg Delta lines: <Unavailable> File [modified]: alp_norm2.xg Delta lines: <Unavailable> File [modified]: dens_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: dens_minimum.xg Delta lines: <Unavailable> File [modified]: dens_norm1.xg Delta lines: <Unavailable> File [modified]: dens_norm2.xg Delta lines: <Unavailable> File [modified]: dens_x_[16][16].xg Delta lines: <Unavailable> File [modified]: dens_y_[16][16].xg Delta lines: <Unavailable> File [modified]: dens_z_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: eps_minimum.xg Delta lines: <Unavailable> File [modified]: eps_norm1.xg Delta lines: <Unavailable> File [modified]: eps_norm2.xg Delta lines: <Unavailable> File [modified]: eps_x_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_y_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_z_[16][16].xg Delta lines: <Unavailable> File [modified]: gxx_minimum.xg Delta lines: <Unavailable> File [modified]: gxx_norm1.xg Delta lines: <Unavailable> File [modified]: gxx_norm2.xg Delta lines: <Unavailable> File [modified]: gxy_maximum.xg Delta lines: <Unavailable> File [modified]: gxy_minimum.xg Delta lines: <Unavailable> File [modified]: gxy_norm1.xg Delta lines: <Unavailable> File [modified]: gxy_norm2.xg Delta lines: <Unavailable> File [modified]: gxz_maximum.xg Delta lines: <Unavailable> File [modified]: gxz_minimum.xg Delta lines: <Unavailable> File [modified]: gxz_norm1.xg Delta lines: <Unavailable> File [modified]: gxz_norm2.xg Delta lines: <Unavailable> File [modified]: gyy_minimum.xg Delta lines: <Unavailable> File [modified]: gyy_norm1.xg Delta lines: <Unavailable> File [modified]: gyy_norm2.xg Delta lines: <Unavailable> File [modified]: gyz_maximum.xg Delta lines: <Unavailable> File [modified]: gyz_minimum.xg Delta lines: <Unavailable> File [modified]: gyz_norm1.xg Delta lines: <Unavailable> File [modified]: gyz_norm2.xg Delta lines: <Unavailable> File [modified]: gzz_minimum.xg Delta lines: <Unavailable> File [modified]: gzz_norm1.xg Delta lines: <Unavailable> File [modified]: gzz_norm2.xg Delta lines: <Unavailable> File [modified]: ham_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: ham_maximum.xg Delta lines: <Unavailable> File [modified]: ham_norm1.xg Delta lines: <Unavailable> File [modified]: ham_norm2.xg Delta lines: <Unavailable> File [modified]: ham_x_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_y_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_z_[16][16].xg Delta lines: <Unavailable> File [modified]: kxx_minimum.xg Delta lines: <Unavailable> File [modified]: kxx_norm1.xg Delta lines: <Unavailable> File [modified]: kxx_norm2.xg Delta lines: <Unavailable> File [modified]: kxy_maximum.xg Delta lines: <Unavailable> File [modified]: kxy_minimum.xg Delta lines: <Unavailable> File [modified]: kxy_norm1.xg Delta lines: <Unavailable> File [modified]: kxy_norm2.xg Delta lines: <Unavailable> File [modified]: kxz_maximum.xg Delta lines: <Unavailable> File [modified]: kxz_minimum.xg Delta lines: <Unavailable> File [modified]: kxz_norm1.xg Delta lines: <Unavailable> File [modified]: kxz_norm2.xg Delta lines: <Unavailable> File [modified]: kyy_minimum.xg Delta lines: <Unavailable> File [modified]: kyy_norm1.xg Delta lines: <Unavailable> File [modified]: kyy_norm2.xg Delta lines: <Unavailable> File [modified]: kyz_maximum.xg Delta lines: <Unavailable> File [modified]: kyz_minimum.xg Delta lines: <Unavailable> File [modified]: kyz_norm1.xg Delta lines: <Unavailable> File [modified]: kyz_norm2.xg Delta lines: <Unavailable> File [modified]: kzz_minimum.xg Delta lines: <Unavailable> File [modified]: kzz_norm1.xg Delta lines: <Unavailable> File [modified]: kzz_norm2.xg Delta lines: <Unavailable> File [modified]: momx_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: momx_maximum.xg Delta lines: <Unavailable> File [modified]: momx_minimum.xg Delta lines: <Unavailable> File [modified]: momx_norm1.xg Delta lines: <Unavailable> File [modified]: momx_norm2.xg Delta lines: <Unavailable> File [modified]: momx_x_[16][16].xg Delta lines: <Unavailable> File [modified]: momx_y_[16][16].xg Delta lines: <Unavailable> File [modified]: momx_z_[16][16].xg Delta lines: <Unavailable> File [modified]: momy_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: momy_maximum.xg Delta lines: <Unavailable> File [modified]: momy_minimum.xg Delta lines: <Unavailable> File [modified]: momy_norm1.xg Delta lines: <Unavailable> File [modified]: momy_norm2.xg Delta lines: <Unavailable> File [modified]: momy_x_[16][16].xg Delta lines: <Unavailable> File [modified]: momy_y_[16][16].xg Delta lines: <Unavailable> File [modified]: momy_z_[16][16].xg Delta lines: <Unavailable> File [modified]: momz_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: momz_maximum.xg Delta lines: <Unavailable> File [modified]: momz_minimum.xg Delta lines: <Unavailable> File [modified]: momz_norm1.xg Delta lines: <Unavailable> File [modified]: momz_norm2.xg Delta lines: <Unavailable> File [modified]: momz_x_[16][16].xg Delta lines: <Unavailable> File [modified]: momz_y_[16][16].xg Delta lines: <Unavailable> File [modified]: momz_z_[16][16].xg Delta lines: <Unavailable> File [modified]: press_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: press_minimum.xg Delta lines: <Unavailable> File [modified]: press_norm1.xg Delta lines: <Unavailable> File [modified]: press_x_[16][16].xg Delta lines: <Unavailable> File [modified]: press_y_[16][16].xg Delta lines: <Unavailable> File [modified]: press_z_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: rho_norm1.xg Delta lines: <Unavailable> File [modified]: rho_norm2.xg Delta lines: <Unavailable> File [modified]: rho_x_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_y_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_z_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[0]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[0]_norm1.xg Delta lines: <Unavailable> File [modified]: scon[0]_norm2.xg Delta lines: <Unavailable> File [modified]: scon[0]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[0]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[0]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[1]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[1]_norm1.xg Delta lines: <Unavailable> File [modified]: scon[1]_norm2.xg Delta lines: <Unavailable> File [modified]: scon[1]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[1]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[1]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[2]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[2]_norm1.xg Delta lines: <Unavailable> File [modified]: scon[2]_norm2.xg Delta lines: <Unavailable> File [modified]: scon[2]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[2]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: scon[2]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: tau_minimum.xg Delta lines: <Unavailable> File [modified]: tau_norm1.xg Delta lines: <Unavailable> File [modified]: tau_x_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_y_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_z_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[0]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[0]_maximum.xg Delta lines: <Unavailable> File [modified]: vel[0]_minimum.xg Delta lines: <Unavailable> File [modified]: vel[0]_norm1.xg Delta lines: <Unavailable> File [modified]: vel[0]_norm2.xg Delta lines: <Unavailable> File [modified]: vel[0]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[0]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[0]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[1]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[1]_maximum.xg Delta lines: <Unavailable> File [modified]: vel[1]_minimum.xg Delta lines: <Unavailable> File [modified]: vel[1]_norm1.xg Delta lines: <Unavailable> File [modified]: vel[1]_norm2.xg Delta lines: <Unavailable> File [modified]: vel[1]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[1]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[1]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[2]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[2]_maximum.xg Delta lines: <Unavailable> File [modified]: vel[2]_minimum.xg Delta lines: <Unavailable> File [modified]: vel[2]_norm1.xg Delta lines: <Unavailable> File [modified]: vel[2]_norm2.xg Delta lines: <Unavailable> File [modified]: vel[2]_x_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[2]_y_[16][16].xg Delta lines: <Unavailable> File [modified]: vel[2]_z_[16][16].xg Delta lines: <Unavailable> File [modified]: w_lorentz_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: w_lorentz_maximum.xg Delta lines: <Unavailable> File [modified]: w_lorentz_norm1.xg Delta lines: <Unavailable> File [modified]: w_lorentz_norm2.xg Delta lines: <Unavailable> File [modified]: w_lorentz_x_[16][16].xg Delta lines: <Unavailable> File [modified]: w_lorentz_y_[16][16].xg Delta lines: <Unavailable> File [modified]: w_lorentz_z_[16][16].xg Delta lines: <Unavailable> Directory: /trunk/src/ ====================== File [modified]: GRHydro_UpdateMask.F90 Delta lines: <Unavailable> Directory: /trunk/test/ ======================= File [modified]: GRHydro_test_tov_ppm.par Delta lines: <Unavailable> File [modified]: GRHydro_test_tov_ppm_ML.par Delta lines: <Unavailable> File [modified]: GRHydro_test_tov_ppm_ML_disable_internal_excision.par Delta lines: <Unavailable> File [modified]: GRHydro_test_tov_ppm_ML_eos_omni.par Delta lines: <Unavailable> File [modified]: GRHydro_test_tov_ppm_disable_internal_excision.par Delta lines: <Unavailable> File [modified]: GRHydro_test_tov_ppm_no_trp.par Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: <Unavailable> File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 155)
by knarf＠cct.lsu.edu 15 Sep '10

15 Sep '10

User: knarf Date: 2010/09/15 01:48 PM Modified: /trunk/src/ make.configuration.deps Log: add lstdc__ to Flags for utils File Changes: Directory: /trunk/src/ ====================== File [modified]: make.configuration.deps Delta lines: +1 -1 =================================================================== --- trunk/src/make.configuration.deps 2010-09-15 13:32:42 UTC (rev 154) +++ trunk/src/make.configuration.deps 2010-09-15 18:48:46 UTC (rev 155) @@ -8,7 +8,7 @@ $(UTIL_DIR)$(DIRSEP)%: $(BUILD_DIR)$(DIRSEP)GRHydro$(DIRSEP)%.o @echo Creating $* in $(UTIL_DIR) from $< if [ ! -d $(UTIL_DIR) ]; then $(MKDIR) $(MKDIRFLAGS) $(UTIL_DIR) ; fi - $(F90) $(CREATEEXE)$(OPTIONSEP)$@ $(DEBUG_LD) $(LDFLAGS) $(EXTRAFLAGS) $(GENERAL_LIBRARIES) $< + $(F90) $(CREATEEXE)$(OPTIONSEP)$@ $(DEBUG_LD) $(LDFLAGS) $(EXTRAFLAGS) $(GENERAL_LIBRARIES) -lstdc++ $< $(BUILD_DIR)$(DIRSEP)GRHydro$(DIRSEP)%.o: $(PACKAGE_DIR)$(DIRSEP)EinsteinEvolve$(DIRSEP)GRHydro$(DIRSEP)src$(DIRSEP)util$(DIRSEP)%.f90

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 154)
by knarf＠cct.lsu.edu 15 Sep '10

15 Sep '10

User: knarf Date: 2010/09/15 08:32 AM Modified: /trunk/src/ GRHydro_Minima.cc Log: tabs to spaces File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Minima.cc Delta lines: +36 -36 =================================================================== --- trunk/src/GRHydro_Minima.cc 2010-09-08 19:04:46 UTC (rev 153) +++ trunk/src/GRHydro_Minima.cc 2010-09-15 13:32:42 UTC (rev 154) @@ -38,7 +38,7 @@ @desc After initial data, set up the scalar GRHydro_rho_min used for checking the atmosphere. This is computed taking into account the actual maximum rest-mass density on the grid. - Version for Carpet. + Version for Carpet. @enddesc @calls @calledby @@ -62,41 +62,41 @@ if (flag) // Actually run this reduction only the first time the routine is called. { if (rho_abs_min > 0.0) - { - *GRHydro_rho_min = rho_abs_min; - } + { + *GRHydro_rho_min = rho_abs_min; + } else - { - // Go to global mode + { + // Go to global mode - BEGIN_GLOBAL_MODE (cctkGH) { + BEGIN_GLOBAL_MODE (cctkGH) { // Find the global maximum of rho const int Reduction_Handle = CCTK_ReductionHandle ("maximum"); - + const CCTK_INT input_array_variable_indices= {CCTK_VarIndex("HydroBase::rho")}; - + const int ierr = CCTK_Reduce(cctkGH, - -1, // target processors; -1 -> all - Reduction_Handle, - 1, // number of output variables - CCTK_VARIABLE_REAL, - &max_rho, - 1, // number of variables to be reduced - input_array_variable_indices); - + -1, // target processors; -1 -> all + Reduction_Handle, + 1, // number of output variables + CCTK_VARIABLE_REAL, + &max_rho, + 1, // number of variables to be reduced + input_array_variable_indices); + if (ierr != 0) - { - CCTK_WARN(0, "Failed to compute the global maximum of rho"); - } - + { + CCTK_WARN(0, "Failed to compute the global maximum of rho"); + } + *GRHydro_rho_min = max_rho * rho_rel_min; - + // Go back to local mode - } END_GLOBAL_MODE; - - } + } END_GLOBAL_MODE; + + } // After this has run once, set the flag so that this does not run again flag = false; } // end if (flag) @@ -152,22 +152,22 @@ // Find the global maximum of rho const int Reduction_Handle = CCTK_ReductionHandle("maximum"); - + const CCTK_INT input_array_variable_indices={CCTK_VarIndex("HydroBase::rho")}; const int ierr = CCTK_Reduce(cctkGH, - -1, // target processors; -1 -> all - Reduction_Handle, - 1, // number of output variables - CCTK_VARIABLE_REAL, - &max_rho, - 1, // number of variables to be reduced - input_array_variable_indices); + -1, // target processors; -1 -> all + Reduction_Handle, + 1, // number of output variables + CCTK_VARIABLE_REAL, + &max_rho, + 1, // number of variables to be reduced + input_array_variable_indices); if (ierr != 0) - { - CCTK_WARN(0, "Failed to compute the global maximum of rho"); - } + { + CCTK_WARN(0, "Failed to compute the global maximum of rho"); + } *GRHydro_rho_min = max_rho * rho_rel_min; }

1 0

[svn:einsteintoolkit] TOVSolver/trunk/ (Rev. 119)
by knarf＠cct.lsu.edu 14 Sep '10

14 Sep '10

User: knarf Date: 2010/09/14 02:45 PM Modified: /trunk/src/ external.inc, tov.c /trunk/test/test_two_av/ dens_3D_diagonal.xg, dens_minimum.xg, dens_norm1.xg, dens_norm2.xg, dens_x_[16][16].xg, dens_y_[16][16].xg, dens_z_[16][16].xg, eps_3D_diagonal.xg, eps_minimum.xg, eps_norm1.xg, eps_norm2.xg, eps_x_[16][16].xg, eps_y_[16][16].xg, eps_z_[16][16].xg, ham_3D_diagonal.xg, ham_maximum.xg, ham_norm1.xg, ham_norm2.xg, ham_x_[16][16].xg, ham_y_[16][16].xg, ham_z_[16][16].xg, press_3D_diagonal.xg, press_minimum.xg, press_norm1.xg, press_norm2.xg, press_x_[16][16].xg, press_y_[16][16].xg, press_z_[16][16].xg, rho_3D_diagonal.xg, rho_minimum.xg, rho_norm1.xg, rho_norm2.xg, rho_x_[16][16].xg, rho_y_[16][16].xg, rho_z_[16][16].xg, tau_3D_diagonal.xg, tau_minimum.xg, tau_norm1.xg, tau_norm2.xg, tau_x_[16][16].xg, tau_y_[16][16].xg, tau_z_[16][16].xg Log: Simplify and correct atmosphere handling - This will be taken out of this thorn entirely later, this is only a cleanup to prepare this - This actually also fixes a bug with the atmosphere of multiple stars using the average method, thus the change to the testsuite Remove some unused function arguments of internal function File Changes: Directory: /trunk/src/ ====================== File [modified]: external.inc Delta lines: <Unavailable> File [modified]: tov.c Delta lines: <Unavailable> Directory: /trunk/test/test_two_av/ =================================== File [modified]: dens_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: dens_minimum.xg Delta lines: <Unavailable> File [modified]: dens_norm1.xg Delta lines: <Unavailable> File [modified]: dens_norm2.xg Delta lines: <Unavailable> File [modified]: dens_x_[16][16].xg Delta lines: <Unavailable> File [modified]: dens_y_[16][16].xg Delta lines: <Unavailable> File [modified]: dens_z_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: eps_minimum.xg Delta lines: <Unavailable> File [modified]: eps_norm1.xg Delta lines: <Unavailable> File [modified]: eps_norm2.xg Delta lines: <Unavailable> File [modified]: eps_x_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_y_[16][16].xg Delta lines: <Unavailable> File [modified]: eps_z_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: ham_maximum.xg Delta lines: <Unavailable> File [modified]: ham_norm1.xg Delta lines: <Unavailable> File [modified]: ham_norm2.xg Delta lines: <Unavailable> File [modified]: ham_x_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_y_[16][16].xg Delta lines: <Unavailable> File [modified]: ham_z_[16][16].xg Delta lines: <Unavailable> File [modified]: press_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: press_minimum.xg Delta lines: <Unavailable> File [modified]: press_norm1.xg Delta lines: <Unavailable> File [modified]: press_norm2.xg Delta lines: <Unavailable> File [modified]: press_x_[16][16].xg Delta lines: <Unavailable> File [modified]: press_y_[16][16].xg Delta lines: <Unavailable> File [modified]: press_z_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: rho_minimum.xg Delta lines: <Unavailable> File [modified]: rho_norm1.xg Delta lines: <Unavailable> File [modified]: rho_norm2.xg Delta lines: <Unavailable> File [modified]: rho_x_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_y_[16][16].xg Delta lines: <Unavailable> File [modified]: rho_z_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: tau_minimum.xg Delta lines: <Unavailable> File [modified]: tau_norm1.xg Delta lines: <Unavailable> File [modified]: tau_norm2.xg Delta lines: <Unavailable> File [modified]: tau_x_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_y_[16][16].xg Delta lines: <Unavailable> File [modified]: tau_z_[16][16].xg Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] TOVSolver/trunk/src/ (Rev. 118)
by knarf＠cct.lsu.edu 14 Sep '10

14 Sep '10

User: knarf Date: 2010/09/14 01:58 PM Modified: /trunk/src/ external.inc, utils.inc Log: parallelize some parts with openmp File Changes: Directory: /trunk/src/ ====================== File [modified]: external.inc Delta lines: +21 -28 =================================================================== --- trunk/src/external.inc 2010-07-13 15:47:07 UTC (rev 117) +++ trunk/src/external.inc 2010-09-14 18:58:03 UTC (rev 118) @@ -50,10 +50,6 @@ const CCTK_REAL *z ) { DECLARE_CCTK_PARAMETERS - CCTK_REAL r_to_star, press, phi, r; - CCTK_REAL rho, eps, gxx, w_lorentz_2; - CCTK_INT star_i=0, star=0; - CCTK_INT i; assert(TOV_Surface!=0); assert(TOV_R_Surface!=0); @@ -69,8 +65,12 @@ /* TOV_debug_input_points(size,x,y,z); */ /* loop over all points we got */ - for (i=0; i<size; i++) + #pragma omp parallel for + for (CCTK_INT i=0; i<size; i++) { + CCTK_REAL r_to_star, press, phi, r; + CCTK_REAL rho, eps, gxx, w_lorentz_2; + CCTK_INT star_i=0, star=0; /* calculate the distance to the star */ star=get_nearest_star(x,y,z,i); star_i = star * TOV_Num_Radial; @@ -165,8 +165,7 @@ w_lorentz_2 = 1./(1. - gxx*TOV_Velocity_x[star]*TOV_Velocity_x[star] - gxx*TOV_Velocity_y[star]*TOV_Velocity_y[star] - gxx*TOV_Velocity_z[star]*TOV_Velocity_z[star]); - rho_eps = w_lorentz_2 * rho + - w_lorentz_2*TOV_Gamma[star]/(TOV_Gamma[star]-1.) * press; + rho_eps = w_lorentz_2 * (rho + TOV_Gamma[star]/(TOV_Gamma[star]-1.) * press); switch(dir) { case 0: source[i]=psip*rho_eps*gxx*TOV_Velocity_x[star]; break; @@ -365,7 +364,7 @@ *vel1_p, *vel2_p, *scon0_p, *scon1_p, *scon2_p; CCTK_REAL *rho_p_p, *press_p_p, *eps_p_p, *dens_p_p, *tau_p_p, *w_lorentz_p_p, *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 *source, *mom_source[3]; CCTK_REAL psip, my_psi4, sqrt_det, vx, vy, vz, v_2, w_lorentz_2, D_h_w; CCTK_REAL vlowx, vlowy, vlowz; @@ -417,10 +416,6 @@ scon1_p_p =(CCTK_REAL*)CCTK_VarDataPtr(cctkGH, 2, "GRHydro::scon[1]"); scon2_p_p =(CCTK_REAL*)CCTK_VarDataPtr(cctkGH, 2, "GRHydro::scon[2]"); - /* get u */ - u=calloc(size, sizeof(CCTK_REAL)); - TOV_Set_Initial_Guess_for_u(cctkGH, size, u, x, y, z); - /* get the old source */ source=calloc(size, sizeof(CCTK_REAL)); TOV_Set_Rho_ADM(cctkGH, size, source, x, y, z); @@ -434,30 +429,29 @@ if (debug) { debugfile=fopen("tovdebug.dat", "w"); + fprintf(debugfile, "# 1:x 2:y 3:psi4 4:source 5-7:mom_source_?, 8,9: ham-source, 10,11: momx-source\n"); } /* loop over all points we got */ for (i3D=0; i3D<size; i3D++) { if (source[i3D]>0.0) { - if (debug && (fabs(z[i3D])<1.e-15)) - { - fprintf(debugfile, - "%.15g %.15g %.15g %.15g %.15g %.15g ", - x[i3D], y[i3D], source[i3D], - mom_source[0][i3D], mom_source[1][i3D], mom_source[1][i3D]); - } if (psi != NULL) my_psi4 = pow(psi[i3D], 4.0)*gxx[i3D]; else my_psi4 = gxx[i3D]; - /* we do not need u, we need gxx here, which will be stored in u! */ - u[i3D] = pow(u[i3D]+1.0, 4.0); - /* We also store the original J^i in mom_source[i] */ - psip = pow(u[i3D], TOV_Momentum_Psi_Power/4.); + /* We store the original J^i in mom_source[i] */ + psip = pow(my_psi4, TOV_Momentum_Psi_Power/4.); mom_source[0][i3D] /= psip; mom_source[1][i3D] /= psip; mom_source[2][i3D] /= psip; + if (debug && (fabs(z[i3D])<1.e-15)) + { + fprintf(debugfile, + "%.15g %.15g %.15g %.15g %.15g %.15g %.15g ", + x[i3D], y[i3D], my_psi4, source[i3D], + mom_source[0][i3D], mom_source[1][i3D], mom_source[1][i3D]); + } if (rho[i3D] < 1.e-10) rho[i3D] = 1.e-10; @@ -523,12 +517,12 @@ (fabs(rhonew - rhoold) > 1.e-12) ); if (count==1) - printf("Fast rescale! %g %g %g %g\n", u[i3D], my_psi4, - mom_source[0][i3D]/u[i3D]/ - (source[i3D]/u[i3D]/u[i3D]+TOV_K[0]* + printf("Fast rescale! %g %g %g %g\n", my_psi4, my_psi4, + mom_source[0][i3D]/my_psi4/ + (source[i3D]/my_psi4/my_psi4+TOV_K[0]* pow(rhonew, TOV_Gamma[0])), mom_source[0][i3D]/my_psi4/ - (source[i3D]/u[i3D]/my_psi4+TOV_K[0]* + (source[i3D]/my_psi4/my_psi4+TOV_K[0]* pow(rhonew, TOV_Gamma[0]))); #else newton_raphson(&vx,&vy,&vz,&rhoold,&rhonew,&w_lorentz_2,&v_2, @@ -614,7 +608,6 @@ TOV_Copy(size, w_lorentz_p, w_lorentz); } free(source); - free(u); return 1; } File [modified]: utils.inc Delta lines: +3 -5 =================================================================== --- trunk/src/utils.inc 2010-07-13 15:47:07 UTC (rev 117) +++ trunk/src/utils.inc 2010-09-14 18:58:03 UTC (rev 118) @@ -84,9 +84,7 @@ void TOV_Copy(CCTK_INT size, CCTK_REAL *var_p, CCTK_REAL *var) { - for(; size; ) - { - --size; - var_p[size] = var[size]; - } +#pragma omp parallel for + for(int i=0; i<size; i++) + var_p[i] = var[i]; }

1 0

[svn:einsteintoolkit] WeylScal4/trunk/test/ (Rev. 62)
by ian.hinder＠aei.mpg.de 13 Sep '10

13 Sep '10

User: hinder Date: 2010/09/13 06:19 AM Modified: /trunk/test/ teukolsky.par, teukolskyID.par Log: Update test suite to remove obsolete parameter File Changes: Directory: /trunk/test/ ======================= File [modified]: teukolsky.par Delta lines: +0 -3 =================================================================== --- trunk/test/teukolsky.par 2010-09-13 11:18:30 UTC (rev 61) +++ trunk/test/teukolsky.par 2010-09-13 11:19:14 UTC (rev 62) @@ -88,9 +88,6 @@ ML_BSSN_O2::ML_shift_bound = "flat" ML_BSSN_O2::ML_dtshift_bound = "flat" -WeylScal4::Psi4r_group_bound = "flat" -WeylScal4::Psi4i_group_bound = "flat" - # Output IO::out_dir = $parfile IO::out_fileinfo = "none" File [modified]: teukolskyID.par Delta lines: +0 -3 =================================================================== --- trunk/test/teukolskyID.par 2010-09-13 11:18:30 UTC (rev 61) +++ trunk/test/teukolskyID.par 2010-09-13 11:19:14 UTC (rev 62) @@ -73,9 +73,6 @@ # Evolution ADMBase::evolution_method = "static" -WeylScal4::Psi4r_group_bound = "flat" -WeylScal4::Psi4i_group_bound = "flat" - # Output IO::out_dir = $parfile IO::out_fileinfo = "none"

1 0

[svn:einsteintoolkit] WeylScal4/trunk/ (Rev. 61)
by ian.hinder＠aei.mpg.de 13 Sep '10

13 Sep '10

User: hinder Date: 2010/09/13 06:18 AM Modified: /trunk/ interface.ccl, param.ccl, schedule.ccl /trunk/src/ Boundaries.c, psis_calc_2nd.c, psis_calc_4th.c Log: Update automatically generated code File Changes: Directory: /trunk/src/ ====================== File [modified]: Boundaries.c Delta lines: +0 -208 =================================================================== --- trunk/src/Boundaries.c 2010-09-13 11:16:27 UTC (rev 60) +++ trunk/src/Boundaries.c 2010-09-13 11:18:30 UTC (rev 61) @@ -31,219 +31,11 @@ DECLARE_CCTK_PARAMETERS; CCTK_INT ierr = 0; - - if (CCTK_EQUALS(Psi4i_group_bound, "none" ) || - CCTK_EQUALS(Psi4i_group_bound, "static") || - CCTK_EQUALS(Psi4i_group_bound, "flat" ) || - CCTK_EQUALS(Psi4i_group_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4i_group", Psi4i_group_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4i_group_bound BC for WeylScal4::Psi4i_group!"); - } - - if (CCTK_EQUALS(Psi4r_group_bound, "none" ) || - CCTK_EQUALS(Psi4r_group_bound, "static") || - CCTK_EQUALS(Psi4r_group_bound, "flat" ) || - CCTK_EQUALS(Psi4r_group_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4r_group", Psi4r_group_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4r_group_bound BC for WeylScal4::Psi4r_group!"); - } - - if (CCTK_EQUALS(Psi4i_bound, "none" ) || - CCTK_EQUALS(Psi4i_bound, "static") || - CCTK_EQUALS(Psi4i_bound, "flat" ) || - CCTK_EQUALS(Psi4i_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4i", Psi4i_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4i_bound BC for WeylScal4::Psi4i!"); - } - - if (CCTK_EQUALS(Psi4r_bound, "none" ) || - CCTK_EQUALS(Psi4r_bound, "static") || - CCTK_EQUALS(Psi4r_bound, "flat" ) || - CCTK_EQUALS(Psi4r_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4r", Psi4r_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4r_bound BC for WeylScal4::Psi4r!"); - } - - if (CCTK_EQUALS(Psi4i_group_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4i_group_bound = -1; - if (handle_Psi4i_group_bound < 0) handle_Psi4i_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_group_bound , Psi4i_group_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4i_group_bound ,Psi4i_group_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_group_bound, - "WeylScal4::Psi4i_group", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4i_group!"); - - } - - if (CCTK_EQUALS(Psi4r_group_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4r_group_bound = -1; - if (handle_Psi4r_group_bound < 0) handle_Psi4r_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_group_bound , Psi4r_group_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4r_group_bound ,Psi4r_group_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_group_bound, - "WeylScal4::Psi4r_group", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4r_group!"); - - } - - if (CCTK_EQUALS(Psi4i_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4i_bound = -1; - if (handle_Psi4i_bound < 0) handle_Psi4i_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_bound , Psi4i_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4i_bound ,Psi4i_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_bound, - "WeylScal4::Psi4i", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4i!"); - - } - - if (CCTK_EQUALS(Psi4r_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4r_bound = -1; - if (handle_Psi4r_bound < 0) handle_Psi4r_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_bound , Psi4r_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4r_bound ,Psi4r_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_bound, - "WeylScal4::Psi4r", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4r!"); - - } - - if (CCTK_EQUALS(Psi4i_group_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4i_group_bound = -1; - if (handle_Psi4i_group_bound < 0) handle_Psi4i_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_group_bound ,Psi4i_group_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_group_bound, - "WeylScal4::Psi4i_group", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for WeylScal4::Psi4i_group!"); - - } - - if (CCTK_EQUALS(Psi4r_group_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4r_group_bound = -1; - if (handle_Psi4r_group_bound < 0) handle_Psi4r_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_group_bound ,Psi4r_group_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_group_bound, - "WeylScal4::Psi4r_group", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for WeylScal4::Psi4r_group!"); - - } - - if (CCTK_EQUALS(Psi4i_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4i_bound = -1; - if (handle_Psi4i_bound < 0) handle_Psi4i_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_bound ,Psi4i_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_bound, - "WeylScal4::Psi4i", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for WeylScal4::Psi4i!"); - - } - - if (CCTK_EQUALS(Psi4r_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4r_bound = -1; - if (handle_Psi4r_bound < 0) handle_Psi4r_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_bound ,Psi4r_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_bound, - "WeylScal4::Psi4r", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for WeylScal4::Psi4r!"); - - } return; } /* template for entries in parameter file: -#$bound$#WeylScal4::Psi4i_group_bound = "skip" -#$bound$#WeylScal4::Psi4i_group_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4i_group_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4i_group_bound_scalar = 0.0 - -#$bound$#WeylScal4::Psi4r_group_bound = "skip" -#$bound$#WeylScal4::Psi4r_group_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4r_group_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4r_group_bound_scalar = 0.0 - -#$bound$#WeylScal4::Psi4i_bound = "skip" -#$bound$#WeylScal4::Psi4i_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4i_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4i_bound_scalar = 0.0 - -#$bound$#WeylScal4::Psi4r_bound = "skip" -#$bound$#WeylScal4::Psi4r_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4r_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4r_bound_scalar = 0.0 - */ File [modified]: psis_calc_2nd.c Delta lines: +15 -0 =================================================================== --- trunk/src/psis_calc_2nd.c 2010-09-13 11:16:27 UTC (rev 60) +++ trunk/src/psis_calc_2nd.c 2010-09-13 11:18:30 UTC (rev 61) @@ -19,6 +19,21 @@ #define CUB(x) ((x) * (x) * (x)) #define QAD(x) ((x) * (x) * (x) * (x)) +void psis_calc_2nd_SelectBCs(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + CCTK_INT ierr = 0; + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4i_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4i_group."); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4r_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4r_group."); + return; +} + void psis_calc_2nd_Body(cGH const * restrict const cctkGH, int const dir, int const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], int const min[3], int const max[3], int const n_subblock_gfs, CCTK_REAL * restrict const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS; File [modified]: psis_calc_4th.c Delta lines: +15 -0 =================================================================== --- trunk/src/psis_calc_4th.c 2010-09-13 11:16:27 UTC (rev 60) +++ trunk/src/psis_calc_4th.c 2010-09-13 11:18:30 UTC (rev 61) @@ -19,6 +19,21 @@ #define CUB(x) ((x) * (x) * (x)) #define QAD(x) ((x) * (x) * (x) * (x)) +void psis_calc_4th_SelectBCs(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + CCTK_INT ierr = 0; + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4i_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4i_group."); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4r_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4r_group."); + return; +} + void psis_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, int const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], int const min[3], int const max[3], int const n_subblock_gfs, CCTK_REAL * restrict const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS; Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +1 -1 =================================================================== --- trunk/interface.ccl 2010-09-13 11:16:27 UTC (rev 60) +++ trunk/interface.ccl 2010-09-13 11:18:30 UTC (rev 61) @@ -2,7 +2,7 @@ implements: WeylScal4 -inherits: admbase Grid GenericFD Boundary +inherits: admbase methodoflines Grid GenericFD Boundary File [modified]: param.ccl Delta lines: +10 -130 =================================================================== --- trunk/param.ccl 2010-09-13 11:16:27 UTC (rev 60) +++ trunk/param.ccl 2010-09-13 11:18:30 UTC (rev 61) @@ -52,25 +52,31 @@ } "2nd" restricted: +CCTK_INT WeylScal4_MaxNumEvolvedVars "Number of evolved variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars +{ + 0:0 :: "Number of evolved variables used by this thorn" +} 0 + +restricted: CCTK_INT timelevels "Number of active timelevels" { 0:3 :: "" } 3 restricted: -CCTK_INT psis_calc_2nd_calc_every "psis_calc_2nd_calc_every" +CCTK_INT rhs_timelevels "Number of active RHS timelevels" { - *:* :: "" + 0:3 :: "" } 1 restricted: -CCTK_INT psis_calc_4th_calc_every "psis_calc_4th_calc_every" +CCTK_INT psis_calc_2nd_calc_every "psis_calc_2nd_calc_every" { *:* :: "" } 1 restricted: -CCTK_INT dummycalc_calc_every "dummycalc_calc_every" +CCTK_INT psis_calc_4th_calc_every "psis_calc_4th_calc_every" { *:* :: "" } 1 @@ -87,129 +93,3 @@ *:* :: "" } 0 -restricted: -CCTK_INT dummycalc_calc_offset "dummycalc_calc_offset" -{ - *:* :: "" -} 0 - -private: -KEYWORD Psi4i_bound "Boundary condition to implement" -{ - "flat" :: "Flat boundary condition" - "none" :: "No boundary condition" - "static" :: "Boundaries held fixed" - "radiative" :: "Radiation boundary condition" - "scalar" :: "Dirichlet boundary condition" - "newrad" :: "Improved radiative boundary condition" - "skip" :: "skip boundary condition code" -} "skip" - -private: -KEYWORD Psi4r_bound "Boundary condition to implement" -{ - "flat" :: "Flat boundary condition" - "none" :: "No boundary condition" - "static" :: "Boundaries held fixed" - "radiative" :: "Radiation boundary condition" - "scalar" :: "Dirichlet boundary condition" - "newrad" :: "Improved radiative boundary condition" - "skip" :: "skip boundary condition code" -} "skip" - -private: -KEYWORD Psi4i_group_bound "Boundary condition to implement" -{ - "flat" :: "Flat boundary condition" - "none" :: "No boundary condition" - "static" :: "Boundaries held fixed" - "radiative" :: "Radiation boundary condition" - "scalar" :: "Dirichlet boundary condition" - "newrad" :: "Improved radiative boundary condition" - "skip" :: "skip boundary condition code" -} "skip" - -private: -KEYWORD Psi4r_group_bound "Boundary condition to implement" -{ - "flat" :: "Flat boundary condition" - "none" :: "No boundary condition" - "static" :: "Boundaries held fixed" - "radiative" :: "Radiation boundary condition" - "scalar" :: "Dirichlet boundary condition" - "newrad" :: "Improved radiative boundary condition" - "skip" :: "skip boundary condition code" -} "skip" - -private: -CCTK_REAL Psi4i_bound_speed "characteristic speed at boundary" -{ - "0:*" :: "outgoing characteristic speed > 0" -} 1. - -private: -CCTK_REAL Psi4r_bound_speed "characteristic speed at boundary" -{ - "0:*" :: "outgoing characteristic speed > 0" -} 1. - -private: -CCTK_REAL Psi4i_group_bound_speed "characteristic speed at boundary" -{ - "0:*" :: "outgoing characteristic speed > 0" -} 1. - -private: -CCTK_REAL Psi4r_group_bound_speed "characteristic speed at boundary" -{ - "0:*" :: "outgoing characteristic speed > 0" -} 1. - -private: -CCTK_REAL Psi4i_bound_limit "limit value for r -> infinity" -{ - "*:*" :: "value of limit value is unrestricted" -} 0. - -private: -CCTK_REAL Psi4r_bound_limit "limit value for r -> infinity" -{ - "*:*" :: "value of limit value is unrestricted" -} 0. - -private: -CCTK_REAL Psi4i_group_bound_limit "limit value for r -> infinity" -{ - "*:*" :: "value of limit value is unrestricted" -} 0. - -private: -CCTK_REAL Psi4r_group_bound_limit "limit value for r -> infinity" -{ - "*:*" :: "value of limit value is unrestricted" -} 0. - -private: -CCTK_REAL Psi4i_bound_scalar "Dirichlet boundary value" -{ - "*:*" :: "unrestricted" -} 0. - -private: -CCTK_REAL Psi4r_bound_scalar "Dirichlet boundary value" -{ - "*:*" :: "unrestricted" -} 0. - -private: -CCTK_REAL Psi4i_group_bound_scalar "Dirichlet boundary value" -{ - "*:*" :: "unrestricted" -} 0. - -private: -CCTK_REAL Psi4r_group_bound_scalar "Dirichlet boundary value" -{ - "*:*" :: "unrestricted" -} 0. - File [modified]: schedule.ccl Delta lines: +93 -50 =================================================================== --- trunk/schedule.ccl 2010-09-13 11:16:27 UTC (rev 60) +++ trunk/schedule.ccl 2010-09-13 11:18:30 UTC (rev 61) @@ -1,31 +1,9 @@ # File produced by Kranc -if (timelevels == 1) -{ - STORAGE: Psi4i_group[1] -} -if (timelevels == 2) -{ - STORAGE: Psi4i_group[2] -} -if (timelevels == 3) -{ - STORAGE: Psi4i_group[3] -} +STORAGE: Psi4i_group[3] -if (timelevels == 1) -{ - STORAGE: Psi4r_group[1] -} -if (timelevels == 2) -{ - STORAGE: Psi4r_group[2] -} -if (timelevels == 3) -{ - STORAGE: Psi4r_group[3] -} +STORAGE: Psi4r_group[3] schedule WeylScal4_Startup at STARTUP { @@ -46,61 +24,126 @@ } "register symmetries" -schedule group WeylScal4_Calculate in MoL_PseudoEvolution +if (CCTK_EQUALS(fd_order, "2nd")) { -} "Calculate the Newman-Penrose quantities" + schedule group psis_calc_2nd_group in MoL_PseudoEvolution + { + # no language specified + } "psis_calc_2nd" +} +schedule psis_calc_2nd in psis_calc_2nd_group +{ + LANG: C +} "psis_calc_2nd" + +schedule psis_calc_2nd_SelectBCs in psis_calc_2nd_bc_group +{ + LANG: C + OPTIONS: level + SYNC: Psi4i_group + SYNC: Psi4r_group +} "psis_calc_2nd_SelectBCs" + +schedule group ApplyBCs as psis_calc_2nd_ApplyBCs in psis_calc_2nd_bc_group after psis_calc_2nd_SelectBCs +{ + # no language specified +} "Apply BCs for groups set in psis_calc_2nd" + + if (CCTK_EQUALS(fd_order, "2nd")) { - schedule psis_calc_2nd in WeylScal4_Calculate as calc_np + schedule group psis_calc_2nd_bc_group in psis_calc_2nd_group after psis_calc_2nd { - LANG: C - SYNC: Psi4i_group - SYNC: Psi4r_group + # no language specified } "psis_calc_2nd" } +if (CCTK_EQUALS(fd_order, "2nd")) +{ + schedule group psis_calc_2nd_bc_group at CCTK_POSTRESTRICT + { + # no language specified + } "psis_calc_2nd" +} + + +if (CCTK_EQUALS(fd_order, "2nd")) +{ + schedule group psis_calc_2nd_bc_group at CCTK_POSTRESTRICTINITIAL + { + # no language specified + } "psis_calc_2nd" +} + + if (CCTK_EQUALS(fd_order, "4th")) { - schedule psis_calc_4th in WeylScal4_Calculate as calc_np + schedule group psis_calc_4th_group in MoL_PseudoEvolution { - LANG: C - SYNC: Psi4i_group - SYNC: Psi4r_group + # no language specified } "psis_calc_4th" } -schedule WeylScal4_SelectBoundConds in WeylScal4_Boundaries +schedule psis_calc_4th in psis_calc_4th_group { LANG: C +} "psis_calc_4th" + +schedule psis_calc_4th_SelectBCs in psis_calc_4th_bc_group +{ + LANG: C OPTIONS: level SYNC: Psi4i_group SYNC: Psi4r_group -} "select boundary conditions" +} "psis_calc_4th_SelectBCs" -schedule WeylScal4_CheckBoundaries at BASEGRID +schedule group ApplyBCs as psis_calc_4th_ApplyBCs in psis_calc_4th_bc_group after psis_calc_4th_SelectBCs { - LANG: C - OPTIONS: meta -} "check boundaries treatment" + # no language specified +} "Apply BCs for groups set in psis_calc_4th" -schedule group ApplyBCs as WeylScal4_ApplyBCs in WeylScal4_Boundaries after WeylScal4_SelectBoundConds + +if (CCTK_EQUALS(fd_order, "4th")) { - # no language specified -} "Apply boundary conditions controlled by thorn Boundary" + schedule group psis_calc_4th_bc_group in psis_calc_4th_group after psis_calc_4th + { + # no language specified + } "psis_calc_4th" +} -schedule group WeylScal4_Boundaries in WeylScal4_Calculate after calc_np + +if (CCTK_EQUALS(fd_order, "4th")) { - # no language specified -} "Apply boundary conditions controlled by thorn Boundary" + schedule group psis_calc_4th_bc_group at CCTK_POSTRESTRICT + { + # no language specified + } "psis_calc_4th" +} -schedule group WeylScal4_Boundaries in CCTK_POSTRESTRICTINITIAL + +if (CCTK_EQUALS(fd_order, "4th")) { - # no language specified -} "Apply boundary conditions controlled by thorn Boundary" + schedule group psis_calc_4th_bc_group at CCTK_POSTRESTRICTINITIAL + { + # no language specified + } "psis_calc_4th" +} -schedule group WeylScal4_Boundaries in CCTK_POSTRESTRICT +schedule WeylScal4_SelectBoundConds in MoL_PostStep { + LANG: C + OPTIONS: level +} "select boundary conditions" + +schedule WeylScal4_CheckBoundaries at BASEGRID +{ + LANG: C + OPTIONS: meta +} "check boundaries treatment" + +schedule group ApplyBCs as WeylScal4_ApplyBCs in MoL_PostStep after WeylScal4_SelectBoundConds +{ # no language specified } "Apply boundary conditions controlled by thorn Boundary"

1 0

[svn:einsteintoolkit] WeylScal4/trunk/m/ (Rev. 60)
by ian.hinder＠aei.mpg.de 13 Sep '10

13 Sep '10

User: hinder Date: 2010/09/13 06:16 AM Modified: /trunk/m/ Makefile, WeylScal4.m, WeylScal4.patch Log: Schedule calculation automatically including boundary conditions Newer versions of Kranc default to scheduling calculations in MoL_PseudoEvolution and applying the boundary conditions if no Schedule string is supplied. Here WeylScal4 is updated to use this feature, which removes the need for manual patching. File Changes: Directory: /trunk/m/ ==================== File [modified]: Makefile Delta lines: +0 -1 =================================================================== --- trunk/m/Makefile 2010-09-13 11:13:03 UTC (rev 59) +++ trunk/m/Makefile 2010-09-13 11:16:27 UTC (rev 60) @@ -8,7 +8,6 @@ WeylScal4.out: WeylScal4.m rm -rf WeylScal4 ./runmath.sh $^ - patch -p1 -d WeylScal4 <WeylScal4.patch for thorn in WeylScal4; do \ ./copy-if-changed.sh $$thorn/configuration.ccl ../../$$thorn/configuration.ccl && \ ./copy-if-changed.sh $$thorn/interface.ccl ../../$$thorn/interface.ccl && \ File [modified]: WeylScal4.m Delta lines: +2 -26 =================================================================== --- trunk/m/WeylScal4.m 2010-09-13 11:13:03 UTC (rev 59) +++ trunk/m/WeylScal4.m 2010-09-13 11:16:27 UTC (rev 60) @@ -124,7 +124,6 @@ PsisCalc[fdOrder_, PD_] := { Name -> "psis_calc_" <> fdOrder, - Schedule -> {"in WeylScal4_Calculate as calc_np"}, Where -> Interior, ConditionalOnKeyword -> {"fd_order", fdOrder}, Shorthands -> shorthands, @@ -270,28 +269,6 @@ } }; -(* Dummy calculation routine to fool Kranc into thinking that the Weyl scalars - are evolved quantities. This way we get code to register boundary conditions. - Unfortunately it also registers with MoL so we have to remove that code by hand. -*) -DummyCalc := -{ - Name -> "dummycalc", - Schedule -> {"in Nowhere_Nevertoberun"}, - Equations -> - { - dot[Psi4r] -> 0, - dot[Psi4i] -> 0 (*, - dot[Psi3r] -> 0, - dot[Psi3i] -> 0, - dot[Psi2r] -> 0, - dot[Psi2i] -> 0, - dot[Psi1r] -> 0, - dot[Psi1i] -> 0, - dot[Psi0r] -> 0, - dot[Psi0i] -> 0 *) - } -} (**************************************************************************** Construct the thorn @@ -312,8 +289,7 @@ calculations = { PsisCalc["2nd", PDstandard2nd], - PsisCalc["4th", PDstandard4th], - DummyCalc + PsisCalc["4th", PDstandard4th] }; CreateKrancThornTT[groups, ".", "WeylScal4", @@ -322,4 +298,4 @@ PartialDerivatives -> derivatives, KeywordParameters -> keywordParameters, RealParameters -> realParameters, - InheritedImplementations -> {"admbase"}]; + InheritedImplementations -> {"admbase", "methodoflines"}]; File [modified]: WeylScal4.patch Delta lines: +0 -181 =================================================================== --- trunk/m/WeylScal4.patch 2010-09-13 11:13:03 UTC (rev 59) +++ trunk/m/WeylScal4.patch 2010-09-13 11:16:27 UTC (rev 60) @@ -1,181 +0,0 @@ -diff -r -c a/interface.ccl b/interface.ccl -*** a/interface.ccl 2010-05-03 23:48:16.810592882 -0400 ---- b/interface.ccl 2010-05-03 23:51:37.206592882 -0400 -*************** -*** 34,48 **** - { - Psi4r - } "Psi4r_group" -- -- public: -- CCTK_REAL Psi4i_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1 tensorparity=-1' -- { -- Psi4irhs -- } "Psi4i_grouprhs" -- -- public: -- CCTK_REAL Psi4r_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1 tensorparity=1' -- { -- Psi4rrhs -- } "Psi4r_grouprhs" ---- 34,36 ---- -diff -r -c a/param.ccl b/param.ccl -*** a/param.ccl 2010-05-03 23:48:16.811592882 -0400 ---- b/param.ccl 2010-05-03 23:51:37.207592882 -0400 -*************** -*** 52,75 **** - } "2nd" - - restricted: -- CCTK_INT WeylScal4_MaxNumEvolvedVars "Number of evolved variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars -- { -- 2:2 :: "Number of evolved variables used by this thorn" -- } 2 -- -- restricted: - CCTK_INT timelevels "Number of active timelevels" - { - 0:3 :: "" - } 3 - - restricted: -- CCTK_INT rhs_timelevels "Number of active RHS timelevels" -- { -- 0:3 :: "" -- } 1 -- -- restricted: - CCTK_INT psis_calc_2nd_calc_every "psis_calc_2nd_calc_every" - { - *:* :: "" ---- 52,63 ---- -diff -r -c a/schedule.ccl b/schedule.ccl -*** a/schedule.ccl 2010-05-03 23:48:16.811592882 -0400 ---- b/schedule.ccl 2010-05-03 23:51:37.206592882 -0400 -*************** -*** 27,58 **** - STORAGE: Psi4r_group[3] - } - -- if (rhs_timelevels == 1) -- { -- STORAGE: Psi4i_grouprhs[1] -- } -- if (rhs_timelevels == 2) -- { -- STORAGE: Psi4i_grouprhs[2] -- } -- if (rhs_timelevels == 3) -- { -- STORAGE: Psi4i_grouprhs[3] -- } -- -- if (rhs_timelevels == 1) -- { -- STORAGE: Psi4r_grouprhs[1] -- } -- if (rhs_timelevels == 2) -- { -- STORAGE: Psi4r_grouprhs[2] -- } -- if (rhs_timelevels == 3) -- { -- STORAGE: Psi4r_grouprhs[3] -- } -- - schedule WeylScal4_Startup at STARTUP - { - LANG: C ---- 27,32 ---- -*************** -*** 72,77 **** ---- 46,55 ---- - } "register symmetries" - - -+ schedule group WeylScal4_Calculate in MoL_PseudoEvolution -+ { -+ } "Calculate the Newman-Penrose quantities" -+ - if (CCTK_EQUALS(fd_order, "2nd")) - { - schedule psis_calc_2nd in WeylScal4_Calculate as calc_np -*************** -*** 93,104 **** - } "psis_calc_4th" - } - -! schedule dummycalc in Nowhere_Nevertoberun -! { -! LANG: C -! } "dummycalc" -! -! schedule WeylScal4_SelectBoundConds in MoL_PostStep - { - LANG: C - OPTIONS: level ---- 70,76 ---- - } "psis_calc_4th" - } - -! schedule WeylScal4_SelectBoundConds in WeylScal4_Boundaries - { - LANG: C - OPTIONS: level -*************** -*** 112,118 **** - OPTIONS: meta - } "check boundaries treatment" - -! schedule group ApplyBCs as WeylScal4_ApplyBCs in MoL_PostStep after WeylScal4_SelectBoundConds - { - # no language specified - } "Apply boundary conditions controlled by thorn Boundary" ---- 84,105 ---- - OPTIONS: meta - } "check boundaries treatment" - -! schedule group ApplyBCs as WeylScal4_ApplyBCs in WeylScal4_Boundaries after WeylScal4_SelectBoundConds -! { -! # no language specified -! } "Apply boundary conditions controlled by thorn Boundary" -! -! schedule group WeylScal4_Boundaries in WeylScal4_Calculate after calc_np -! { -! # no language specified -! } "Apply boundary conditions controlled by thorn Boundary" -! -! schedule group WeylScal4_Boundaries in CCTK_POSTRESTRICTINITIAL -! { -! # no language specified -! } "Apply boundary conditions controlled by thorn Boundary" -! -! schedule group WeylScal4_Boundaries in CCTK_POSTRESTRICT - { - # no language specified - } "Apply boundary conditions controlled by thorn Boundary" -diff -r -c a/src/RegisterMoL.c b/src/RegisterMoL.c -*** a/src/RegisterMoL.c 2010-05-03 23:48:16.811592882 -0400 ---- b/src/RegisterMoL.c 2010-05-03 23:51:37.206592882 -0400 -*************** -*** 12,18 **** - CCTK_INT ierr = 0; - - /* Register all the evolved grid functions with MoL */ -- ierr += MoLRegisterEvolved(CCTK_VarIndex("WeylScal4::Psi4i"), CCTK_VarIndex("WeylScal4::Psi4irhs")); -- ierr += MoLRegisterEvolved(CCTK_VarIndex("WeylScal4::Psi4r"), CCTK_VarIndex("WeylScal4::Psi4rrhs")); - return; - } ---- 12,16 ---- -diff -r -c a/src/make.code.defn b/src/make.code.defn -*** a/src/make.code.defn 2010-05-03 23:48:16.811592882 -0400 ---- b/src/make.code.defn 2010-05-03 23:51:37.206592882 -0400 -*************** -*** 1,3 **** - # File produced by Kranc - -! SRCS = Startup.c RegisterMoL.c RegisterSymmetries.c psis_calc_2nd.c psis_calc_4th.c dummycalc.c Boundaries.c ---- 1,3 ---- - # File produced by Kranc - -! SRCS = Startup.c RegisterMoL.c RegisterSymmetries.c psis_calc_2nd.c psis_calc_4th.c Boundaries.c

1 0

[svn:einsteintoolkit] WeylScal4/trunk/m/ (Rev. 59)
by ian.hinder＠aei.mpg.de 13 Sep '10

13 Sep '10

User: hinder Date: 2010/09/13 06:13 AM Modified: /trunk/m/ runmath.sh Log: runmath.sh: Find Kranc correctly in new Einstein Toolkit checkout File Changes: Directory: /trunk/m/ ==================== File [modified]: runmath.sh Delta lines: +3 -1 =================================================================== --- trunk/m/runmath.sh 2010-07-12 14:50:05 UTC (rev 58) +++ trunk/m/runmath.sh 2010-09-13 11:13:03 UTC (rev 59) @@ -17,7 +17,7 @@ KRANCPATH="${KRANCPATH:+$KRANCPATH}" if [ -z "$KRANCPATH" ] ; then - KRANCPATHS=(. ../../../../arrangements/.. ../../../arrangements/.. $HOME) + KRANCPATHS=(. ../../../../arrangements/../repos ../../../arrangements/../repos $HOME) # look for Kranc in a number of likely locations # choice 2 and 3 are $CCTK_HOME if the thorn resides in arrangements or # git-repos @@ -55,6 +55,8 @@ echo -e "Could not find Kranc. Expect trouble unless you have made sure Mathematica\nwill find Kranc. You might also consider setting \$KRANCPATH." fi +echo "MATHPATH = $MATHPATH" + error=$(basename $script .m).err output=$(basename $script .m).out

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits September 2010