Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] HydroBase/trunk/ (Rev. 62)
by rhaas＠tapir.caltech.edu 15 May '13

15 May '13

User: rhaas Date: 2013/05/15 03:46 PM Modified: /trunk/ param.ccl Log: extent ID keywords by "read from file" if given, we do nothing and assume the user uses IOUtil's filereader capability to read in ID from disk. File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +6 -0 =================================================================== --- trunk/param.ccl 2013-05-10 15:42:45 UTC (rev 61) +++ trunk/param.ccl 2013-05-15 20:46:57 UTC (rev 62) @@ -34,6 +34,7 @@ { "none" :: "inactive" "one" :: "initially set to one" + "read from file" :: "do nothing, filereader will provide data" } "none" KEYWORD Y_e_evolution_method "Evolution method for Y_e" @@ -57,18 +58,21 @@ { "none" :: "inactive" "zero" :: "initially set to zero" + "read from file" :: "do nothing, filereader will provide data" } "none" KEYWORD initial_Avec "Initial value for Avec" { "none" :: "inactive" "zero" :: "initially set to zero" + "read from file" :: "do nothing, filereader will provide data" } "none" KEYWORD initial_Aphi "Initial value for Aphi" { "none" :: "inactive" "zero" :: "initially set to zero" + "read from file" :: "do nothing, filereader will provide data" } "none" KEYWORD Bvec_evolution_method "Evolution method for Bvec" @@ -87,12 +91,14 @@ { "none" :: "inactive" "zero" :: "initially set to zero" + "read from file" :: "do nothing, filereader will provide data" } "none" KEYWORD initial_entropy "Initial value for entropy" { "none" :: "inactive" "zero" :: "initially set to zero" + "read from file" :: "do nothing, filereader will provide data" } "none"

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 1011)
by knarf＠cct.lsu.edu 14 May '13

14 May '13

User: knarf Date: 2013/05/14 11:04 PM Modified: /about/members/ members.txt Log: one more member File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +1 -0 =================================================================== --- about/members/members.txt 2013-05-14 15:23:34 UTC (rev 1010) +++ about/members/members.txt 2013-05-15 04:04:02 UTC (rev 1011) @@ -75,6 +75,7 @@ Jennifer Seiler https://www.camk.edu.pl/en/ Nicolaus Copernicus Astronomical Center (NCAC) + Antonios Manousakis Bhupendra Prakash Mishra http://www.ncsu.edu/ North Carolina State University

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 1010)
by knarf＠cct.lsu.edu 14 May '13

14 May '13

User: knarf Date: 2013/05/14 10:23 AM Added: /release-info/ spine_Cactus-gcc4.7__1_16.log, spine_Cactus-gcc4.7__2_8.log, spine_Cactus-intel13.163__1_16.log, spine_Cactus-intel13.163__2_8.log Modified: /release-info/ spine-gcc4.7__1_16.log, spine-gcc4.7__2_8.log, spine-intel13.163__1_16.log, spine-intel13.163__2_8.log Log: results on spine, including the Cactus release thornlist File Changes: Directory: /release-info/ ========================= File [modified]: spine-gcc4.7__1_16.log Delta lines: <Unavailable> File [modified]: spine-gcc4.7__2_8.log Delta lines: <Unavailable> File [modified]: spine-intel13.163__1_16.log Delta lines: <Unavailable> File [modified]: spine-intel13.163__2_8.log Delta lines: <Unavailable> File [added]: spine_Cactus-gcc4.7__1_16.log Delta lines: <Unavailable> File [added]: spine_Cactus-gcc4.7__2_8.log Delta lines: <Unavailable> File [added]: spine_Cactus-intel13.163__1_16.log Delta lines: <Unavailable> File [added]: spine_Cactus-intel13.163__2_8.log Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro_InitData/ (Rev. 213)
by schnetter＠cct.lsu.edu 14 May '13

14 May '13

User: eschnett Date: 2013/05/14 10:04 AM Modified: /branches/ET_2013_05/ configuration.ccl Log: Add missing dependency on thorn Fortran Directory Changes: Directory: /branches/ET_2013_05/ ================================ Property changes on: branches/ET_2013_05 ___________________________________________________________________ Directory: /svn:mergeinfo/ ========================== + /trunk:212 File Changes: Directory: /branches/ET_2013_05/ ================================ File [modified]: configuration.ccl Delta lines: +1 -2 =================================================================== --- branches/ET_2013_05/configuration.ccl 2013-05-14 04:26:17 UTC (rev 212) +++ branches/ET_2013_05/configuration.ccl 2013-05-14 15:04:50 UTC (rev 213) @@ -1,2 +1 @@ -REQUIRES EOS_Omni - +REQUIRES EOS_Omni Fortran

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 212)
by schnetter＠cct.lsu.edu 13 May '13

13 May '13

User: eschnett Date: 2013/05/13 11:26 PM Modified: /trunk/ configuration.ccl Log: Add missing dependency on thorn Fortran File Changes: Directory: /trunk/ ================== File [modified]: configuration.ccl Delta lines: +1 -2 =================================================================== --- trunk/configuration.ccl 2013-05-10 15:43:12 UTC (rev 211) +++ trunk/configuration.ccl 2013-05-14 04:26:17 UTC (rev 212) @@ -1,2 +1 @@ -REQUIRES EOS_Omni - +REQUIRES EOS_Omni Fortran

1 0

[svn:einsteintoolkit] GRHydro/ (Rev. 519)
by rhaas＠tapir.caltech.edu 13 May '13

13 May '13

User: rhaas Date: 2013/05/13 01:58 PM Modified: /branches/ET_2013_05/ /branches/ET_2013_05/src/ GRHydro_SourceAM.F90 Log: GRHydro: do not declare grid scalar as OpenMP PRIVATE instead leave it SHARED since it is a global read-only variable. Directory Changes: Directory: /branches/ET_2013_05/ ================================ Property changes on: branches/ET_2013_05 ___________________________________________________________________ Directory: /svn:mergeinfo/ ========================== - /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510,512,515 + /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510,512,515-516 File Changes: Directory: /branches/ET_2013_05/src/ ==================================== File [modified]: GRHydro_SourceAM.F90 Delta lines: +1 -1 =================================================================== --- branches/ET_2013_05/src/GRHydro_SourceAM.F90 2013-05-13 18:54:53 UTC (rev 518) +++ branches/ET_2013_05/src/GRHydro_SourceAM.F90 2013-05-13 18:58:17 UTC (rev 519) @@ -260,7 +260,7 @@ !$OMP tau_source,sx_source,sy_source,sz_source,& !$OMP dx_det_bydet,dy_det_bydet,dz_det_bydet,& !$OMP gdg_x,gdg_y,gdg_z,& - !$OMP Avcx_source,Avcy_source,Avcz_source,evolve_Lorenz_gge,& + !$OMP Avcx_source,Avcy_source,Avcz_source,& !$OMP uxx, uxy, uxz, uyy, uyz, uzz,& !$OMP dx_gxx, dx_gxy, dx_gxz, dx_gyy, dx_gyz, dx_gzz,& !$OMP dy_gxx, dy_gxy, dy_gxz, dy_gyy, dy_gyz, dy_gzz,&

1 0

[svn:einsteintoolkit] GRHydro/ (Rev. 518)
by rhaas＠tapir.caltech.edu 13 May '13

13 May '13

User: rhaas Date: 2013/05/13 01:54 PM Modified: /branches/ET_2013_05/ /branches/ET_2013_05/src/ GRHydro_Eigenproblem.F90 Log: GRHydro: fix use of uninitialized variable in eigenproblem_leftright this is only used when doing flux split FD which is not the default. Directory Changes: Directory: /branches/ET_2013_05/ ================================ Property changes on: branches/ET_2013_05 ___________________________________________________________________ Directory: /svn:mergeinfo/ ========================== - /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510,512 + /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510,512,515 File Changes: Directory: /branches/ET_2013_05/src/ ==================================== File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +3 -0 =================================================================== --- branches/ET_2013_05/src/GRHydro_Eigenproblem.F90 2013-05-13 18:11:19 UTC (rev 517) +++ branches/ET_2013_05/src/GRHydro_Eigenproblem.F90 2013-05-13 18:54:53 UTC (rev 518) @@ -618,6 +618,9 @@ call EOS_Omni_cs2(handle,keytemp,GRHydro_eos_rf_prec,n,& rho,eps,xtemp,xye,cs2,keyerr,anyerr) + call EOS_Omni_DPressByDEps(handle,keytemp,GRHydro_eos_rf_prec,n,& + rho,eps,xtemp,xye,dpdeps,keyerr,anyerr) + ! if (cs2<0) cs2=0 ! this does not modify the roe crashing problem with shocktube enthalpy = one + eps + press / rho

1 0

[svn:einsteintoolkit] www/about/gallery/ (Rev. 1009)
by knarf＠cct.lsu.edu 13 May '13

13 May '13

User: knarf Date: 2013/05/13 01:12 PM Modified: /about/gallery/ index.php Log: don't show non-ready things File Changes: Directory: /about/gallery/ ========================== File [modified]: index.php Delta lines: +11 -2 =================================================================== --- about/gallery/index.php 2013-05-13 18:06:52 UTC (rev 1008) +++ about/gallery/index.php 2013-05-13 18:12:20 UTC (rev 1009) @@ -1,8 +1,6 @@ <?php $title='Einstein Toolkit Gallery'; include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> -<p><b>This page is currently under construction (12/2012)</b></p> - <p>This page contains images and movies of example simulations using the Einstein Toolkit, either exclusively or combined with external codes. Some (hopefully most) also include the corresponding parameter files @@ -14,11 +12,16 @@ <h3>Using exclusively the Einstein Toolkit</h3> <ul> + + <li> Single, stable neutron star <p>This simulation shows how to evolve a stable, single neutron star. It uses a TOV solution @@ -38,16 +41,22 @@ </table> <a href="tov/rho_max.pdf"><img src="tov/rho_max.png" alt="Maximum density over time"></a> </li> + + </ul> + <?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 517)
by rhaas＠tapir.caltech.edu 13 May '13

13 May '13

User: rhaas Date: 2013/05/13 01:11 PM Removed: /trunk/src/ GRHydro_InterfacesAM.h Modified: /trunk/src/ GRHydro_Con2PrimAM.F90 Log: GRHydro: remove duplicate (and incorrect) GRHydro_InterfacesAM.h sync up GRHydro_Con2PrimAM with GRHydro_Con2PrimM File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2PrimAM.F90 Delta lines: +35 -6 =================================================================== --- trunk/src/GRHydro_Con2PrimAM.F90 2013-05-13 17:54:52 UTC (rev 516) +++ trunk/src/GRHydro_Con2PrimAM.F90 2013-05-13 18:11:19 UTC (rev 517) @@ -13,7 +13,6 @@ #include "cctk_Arguments.h" #include "cctk_Functions.h" #include "SpaceMask.h" -#include "GRHydro_InterfacesAM.h" #include "GRHydro_Macros.h" #define ITER_TOL (1.0e-8) @@ -349,6 +348,18 @@ !of the comoving B-field, b^{\mu} b_{\mu}. It is overwritten !in this routine, but we may need to find a better notation !avoid future confusions. + if(GRHydro_eos_handle .eq. 1 .or. GRHydro_eos_handle .eq. 2) then + + call GRHydro_Con2PrimM_ptold(GRHydro_eos_handle, local_gam(1), dens(i,j,k), & + scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3), tau(i,j,k), & + Bconsx_tmp,Bconsy_tmp,Bconsz_tmp,rho_tmp, & + velx_tmp,vely_tmp,velz_tmp,eps_tmp,press_tmp, & + Bvecx_tmp,Bvecy_tmp,Bvecz_tmp,b2, w_lorentz_tmp,& + g11(i,j,k),g12(i,j,k),g13(i,j,k),g22(i,j,k),g23(i,j,k),g33(i,j,k), & + uxx,uxy,uxz,uyy,uyz,uzz,det, & + epsnegative,GRHydro_C2P_failed(i,j,k)) + + else call GRHydro_Con2PrimM_pt2(GRHydro_eos_handle, keytemp, GRHydro_eos_rf_prec, GRHydro_perc_ptol, local_gam(1), dens(i,j,k), & scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3), tau(i,j,k), & Bconsx_tmp, Bconsy_tmp, Bconsz_tmp,xye(1), & @@ -358,9 +369,16 @@ g22(i,j,k),g23(i,j,k),g33(i,j,k), & uxx,uxy,uxz,uyy,uyz,uzz,det, & epsnegative,GRHydro_C2P_failed(i,j,k)) + endif if(evolve_entropy.ne.0) then + if(GRHydro_C2P_failed(i,j,k).ne.0) then + !Use previous time step for rho: entropy(i,j,k) = entropycons(i,j,k)/dens(i,j,k)*rho(i,j,k) + else + !Use the current correct value of rho returned by con2prim: + entropy(i,j,k) = entropycons(i,j,k)/dens(i,j,k)*rho_tmp + endif endif if(GRHydro_C2P_failed(i,j,k).ne.0) then @@ -449,6 +467,19 @@ Bconsy_tmp = sdet*Bvecy_tmp Bconsz_tmp = sdet*Bvecz_tmp + if(evolve_entropy.ne.0) then + call GRHydro_Con2PrimM_ptee(GRHydro_eos_handle, keytemp, & + GRHydro_eos_rf_prec, local_gam(1), dens(i,j,k), & + scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3), tau(i,j,k), & + Bconsx_tmp,Bconsy_tmp,Bconsz_tmp, & + entropycons(i,j,k), xye(1), & + xtemp(1),rho_tmp,velx_tmp,vely_tmp,velz_tmp,& + eps_tmp,press_tmp,Bvecx_tmp,Bvecy_tmp,Bvecz_tmp,b2,& + w_lorentz_tmp,g11(i,j,k),g12(i,j,k),g13(i,j,k),& + g22(i,j,k),g23(i,j,k),g33(i,j,k), & + uxx,uxy,uxz,uyy,uyz,uzz,det, & + epsnegative,GRHydro_C2P_failed(i,j,k)) + else call GRHydro_Con2PrimM_Polytype_pt(GRHydro_eos_handle, local_pgam, & dens(i,j,k),scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3), sc, & Bconsx_tmp, Bconsy_tmp, Bconsz_tmp, rho_tmp,& @@ -457,6 +488,7 @@ g11(i,j,k),g12(i,j,k),g13(i,j,k),g22(i,j,k),g23(i,j,k),g33(i,j,k), & uxx,uxy,uxz,uyy,uyz,uzz,det, & epsnegative,GRHydro_C2P_failed(i,j,k)) + end if rho(i,j,k) = rho_tmp press(i,j,k) = press_tmp @@ -682,12 +714,8 @@ !!$ tau does need to take into account the existing B-field !!$ with w_lorentz=1, we find tau = sqrtdet*(rho (1+eps+b^2/2)) - dens [Press drops out] - tau(i,j,k) = sdet * (rho(i,j,k)*(1.0+eps(i,j,k)+b2/2.0)) - dens(i,j,k) + tau(i,j,k) = sdet * (rho(i,j,k)*eps(i,j,k)+b2/2.0) - if(tau(i,j,k).le.sdet*b2*0.5d0)then - tau(i,j,k) = GRHydro_tau_min + sdet*b2*0.5d0 - endif - end if end do @@ -870,6 +898,7 @@ if(evolve_Y_e.ne.0) then Y_e(i,j,k) = Y_e_con(i,j,k) / dens(i,j,k) endif + Bconsx_tmp = sqrt(avg_detl)*Bvecxminus(i,j,k) Bconsy_tmp = sqrt(avg_detl)*Bvecyminus(i,j,k) Bconsz_tmp = sqrt(avg_detl)*Bveczminus(i,j,k) File [removed]: GRHydro_InterfacesAM.h Delta lines: +0 -95 =================================================================== --- trunk/src/GRHydro_InterfacesAM.h 2013-05-13 17:54:52 UTC (rev 516) +++ trunk/src/GRHydro_InterfacesAM.h 2013-05-13 18:11:19 UTC (rev 517) @@ -1,95 +0,0 @@ -module Con2PrimAM_fortran_interfaces - implicit none - - interface - - subroutine GRHydro_Con2PrimM_pt( handle, keytemp, prec, & - local_gam, dens, & - sx, sy, sz, & - tau, & - Bconsx, Bconsy, Bconsz, & - xye, xtemp, & - rho, & - velx, vely, velz,& - epsilon, pressure,& - Bx, By, Bz, & - bsq,& - w_lorentz, & - gxx, gxy, gxz, & - gyy, gyz, gzz, & - uxx, uxy, uxz,& - uyy, uyz, uzz,& - det,& - epsnegative, & - retval) - - implicit none - CCTK_INT handle - CCTK_INT keytemp - CCTK_REAL prec - CCTK_REAL local_gam - CCTK_REAL dens - CCTK_REAL sx, sy, sz - CCTK_REAL tau - CCTK_REAL Bconsx, Bconsy, Bconsz - CCTK_REAL xye - CCTK_REAL xtemp - CCTK_REAL rho - CCTK_REAL velx, vely, velz - CCTK_REAL epsilon, pressure - CCTK_REAL Bx, By, Bz - CCTK_REAL bsq - CCTK_REAL w_lorentz - CCTK_REAL gxx, gxy, gxz - CCTK_REAL gyy, gyz, gzz - CCTK_REAL uxx, uxy, uxz - CCTK_REAL uyy, uyz, uzz - CCTK_REAL det - CCTK_INT epsnegative - CCTK_REAL retval - end subroutine GRHydro_Con2PrimM_pt - - subroutine GRHydro_Con2PrimM_Polytype_pt( handle, local_gam,& - dens, & - sx, sy, sz, & - sc, & - Bconsx, Bconsy, Bconsz, & - rho, & - velx, vely, velz,& - epsilon, pressure,& - Bx, By, Bz, & - bsq,& - w_lorentz, & - gxx, gxy, gxz, & - gyy, gyz, gzz, & - uxx, uxy, uxz,& - uyy, uyz, uzz,& - det,& - epsnegative, & - retval) - - implicit none - CCTK_INT handle - CCTK_REAL local_gam - CCTK_REAL dens - CCTK_REAL sx, sy, sz - CCTK_REAL sc - CCTK_REAL Bconsx, Bconsy, Bconsz - CCTK_REAL rho - CCTK_REAL velx, vely, velz - CCTK_REAL epsilon, pressure - CCTK_REAL Bx, By, Bz - CCTK_REAL bsq - CCTK_REAL w_lorentz - CCTK_REAL gxx, gxy, gxz - CCTK_REAL gyy, gyz, gzz - CCTK_REAL uxx, uxy, uxz - CCTK_REAL uyy, uyz, uzz - CCTK_REAL det - CCTK_INT epsnegative - CCTK_REAL retval - end subroutine GRHydro_Con2PrimM_Polytype_pt - - end interface - -end module Con2PrimAM_fortran_interfaces

1 0

[svn:einsteintoolkit] www/about/gallery/ (Rev. 1008)
by knarf＠cct.lsu.edu 13 May '13

13 May '13

User: knarf Date: 2013/05/13 01:06 PM Modified: /about/gallery/ tov.par /about/gallery/tov/ rho_max.pdf, rho_max.png Log: update par-file and plot File Changes: Directory: /about/gallery/ ========================== File [modified]: tov.par Delta lines: +7 -7 =================================================================== --- about/gallery/tov.par 2013-05-11 06:54:09 UTC (rev 1007) +++ about/gallery/tov.par 2013-05-13 18:06:52 UTC (rev 1008) @@ -192,15 +192,15 @@ # Parameters for initial star TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::TOV_Gamma[0] = 2.75 -TOVSolver::TOV_K[0] = 30000 +TOVSolver::TOV_Gamma[0] = 2 +TOVSolver::TOV_K[0] = 100 # Set equation of state for evolution -EOS_Omni::poly_gamma = 2.75 -EOS_Omni::poly_gamma_ini = 2.75 -EOS_Omni::poly_k = 30000 -EOS_Omni::gl_gamma = 2.75 -EOS_Omni::gl_k = 30000 +EOS_Omni::poly_gamma = 2 +EOS_Omni::poly_gamma_ini = 2 +EOS_Omni::poly_k = 100 +EOS_Omni::gl_gamma = 2 +EOS_Omni::gl_k = 100 # I/O Directory: /about/gallery/tov/ ============================== File [modified]: rho_max.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) File [modified]: rho_max.png Delta lines: +0 -0 =================================================================== (Binary files differ)

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits