Commits November 2012

commits@lists.einsteintoolkit.org

9 participants
90 discussions

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 434)
by rhaas＠tapir.caltech.edu 08 Nov '12

08 Nov '12

User: rhaas Date: 2012/11/08 07:53 PM Modified: /trunk/src/ GRHydro_PPMM.F90 Log: GRhydro: add recovstruct_Wv to magnetized PPM routine From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_PPMM.F90 Delta lines: +110 -48 =================================================================== --- trunk/src/GRHydro_PPMM.F90 2012-11-09 01:53:08 UTC (rev 433) +++ trunk/src/GRHydro_PPMM.F90 2012-11-09 01:53:11 UTC (rev 434) @@ -79,7 +79,25 @@ logical :: cond + CCTK_REAL, dimension(nx) :: vx,vy,vz + + + !!$ reconstruct w^i = w_lorentz*v^i to ensure slower than light + !!$ speeds? + if (reconstruct_Wv.ne.0) then + ! all velocity like quantities are now w_lorentz*velocity. We will convert + ! back to ordinary velocity at the end, using w_lorentz = sqrt(1 + g_{ij} + ! w^i w^j). + vx = w_lorentz*velx + vy = w_lorentz*vely + vz = w_lorentz*velz + else + vx = velx + vy = vely + vz = velz + end if + !!$ Initially, all the Riemann problems will be trivial trivial_rp = .true. @@ -89,9 +107,9 @@ do i = 2, nx - 1 drho(i) = 0.5d0 * (rho(i+1) - rho(i-1)) - dvelx(i) = 0.5d0 * (velx(i+1) - velx(i-1)) - dvely(i) = 0.5d0 * (vely(i+1) - vely(i-1)) - dvelz(i) = 0.5d0 * (velz(i+1) - velz(i-1)) + dvelx(i) = 0.5d0 * (vx(i+1) - vx(i-1)) + dvely(i) = 0.5d0 * (vy(i+1) - vy(i-1)) + dvelz(i) = 0.5d0 * (vz(i+1) - vz(i-1)) dBvcx(i) = 0.5d0 * (Bvcx(i+1) - Bvcx(i-1)) dBvcy(i) = 0.5d0 * (Bvcy(i+1) - Bvcy(i-1)) dBvcz(i) = 0.5d0 * (Bvcz(i+1) - Bvcz(i-1)) @@ -143,13 +161,13 @@ rhoplus(i) = 0.5d0 * (rho(i) + rho(i+1)) + & (dmrho(i) - dmrho(i+1)) / 6.d0 rhominus(i+1) = rhoplus(i) - velxplus(i) = 0.5d0 * (velx(i) + velx(i+1)) + & + velxplus(i) = 0.5d0 * (vx(i) + vx(i+1)) + & (dmvelx(i) - dmvelx(i+1)) / 6.d0 velxminus(i+1) = velxplus(i) - velyplus(i) = 0.5d0 * (vely(i) + vely(i+1)) + & + velyplus(i) = 0.5d0 * (vy(i) + vy(i+1)) + & (dmvely(i) - dmvely(i+1)) / 6.d0 velyminus(i+1) = velyplus(i) - velzplus(i) = 0.5d0 * (velz(i) + velz(i+1)) + & + velzplus(i) = 0.5d0 * (vz(i) + vz(i+1)) + & (dmvelz(i) - dmvelz(i+1)) / 6.d0 velzminus(i+1) = velzplus(i) @@ -371,12 +389,12 @@ end if rhoplus(i) = flatten * rhoplus(i) + (1.d0 - flatten) * rho(i) rhominus(i) = flatten * rhominus(i) + (1.d0 - flatten) * rho(i) - velxplus(i) = flatten * velxplus(i) + (1.d0 - flatten) * velx(i) - velxminus(i) = flatten * velxminus(i) + (1.d0 - flatten) * velx(i) - velyplus(i) = flatten * velyplus(i) + (1.d0 - flatten) * vely(i) - velyminus(i) = flatten * velyminus(i) + (1.d0 - flatten) * vely(i) - velzplus(i) = flatten * velzplus(i) + (1.d0 - flatten) * velz(i) - velzminus(i) = flatten * velzminus(i) + (1.d0 - flatten) * velz(i) + velxplus(i) = flatten * velxplus(i) + (1.d0 - flatten) * vx(i) + velxminus(i) = flatten * velxminus(i) + (1.d0 - flatten) * vx(i) + velyplus(i) = flatten * velyplus(i) + (1.d0 - flatten) * vy(i) + velyminus(i) = flatten * velyminus(i) + (1.d0 - flatten) * vy(i) + velzplus(i) = flatten * velzplus(i) + (1.d0 - flatten) * vz(i) + velzminus(i) = flatten * velzminus(i) + (1.d0 - flatten) * vz(i) Bvcxplus(i) = flatten * Bvcxplus(i) + (1.d0 - flatten) * Bvcx(i) Bvcxminus(i) = flatten * Bvcxminus(i) + (1.d0 - flatten) * Bvcx(i) Bvcyplus(i) = flatten * Bvcyplus(i) + (1.d0 - flatten) * Bvcy(i) @@ -402,12 +420,12 @@ end if rhoplus(i) = flatten * rhoplus(i) + (1.d0 - flatten) * rho(i) rhominus(i) = flatten * rhominus(i) + (1.d0 - flatten) * rho(i) - velxplus(i) = flatten * velxplus(i) + (1.d0 - flatten) * velx(i) - velxminus(i) = flatten * velxminus(i) + (1.d0 - flatten) * velx(i) - velyplus(i) = flatten * velyplus(i) + (1.d0 - flatten) * vely(i) - velyminus(i) = flatten * velyminus(i) + (1.d0 - flatten) * vely(i) - velzplus(i) = flatten * velzplus(i) + (1.d0 - flatten) * velz(i) - velzminus(i) = flatten * velzminus(i) + (1.d0 - flatten) * velz(i) + velxplus(i) = flatten * velxplus(i) + (1.d0 - flatten) * vx(i) + velxminus(i) = flatten * velxminus(i) + (1.d0 - flatten) * vx(i) + velyplus(i) = flatten * velyplus(i) + (1.d0 - flatten) * vy(i) + velyminus(i) = flatten * velyminus(i) + (1.d0 - flatten) * vy(i) + velzplus(i) = flatten * velzplus(i) + (1.d0 - flatten) * vz(i) + velzminus(i) = flatten * velzminus(i) + (1.d0 - flatten) * vz(i) Bvcxplus(i) = flatten * Bvcxplus(i) + (1.d0 - flatten) * Bvcx(i) Bvcxminus(i) = flatten * Bvcxminus(i) + (1.d0 - flatten) * Bvcx(i) Bvcyplus(i) = flatten * Bvcyplus(i) + (1.d0 - flatten) * Bvcy(i) @@ -494,12 +512,12 @@ if (cond) then rhominus(i)=rho(i) rhoplus(i)=rho(i) - velxminus(i)=velx(i) - velxplus(i)=velx(i) - velyminus(i)=vely(i) - velyplus(i)=vely(i) - velzminus(i)=velz(i) - velzplus(i)=velz(i) + velxminus(i)=vx(i) + velxplus(i)=vx(i) + velyminus(i)=vy(i) + velyplus(i)=vy(i) + velzminus(i)=vz(i) + velzplus(i)=vz(i) Bvcxminus(i)=Bvcx(i) Bvcxplus(i)=Bvcx(i) Bvcyminus(i)=Bvcy(i) @@ -508,12 +526,12 @@ Bvczplus(i)=Bvcz(i) rhominus(i-1)=rho(i) rhoplus(i-1)=rho(i) - velxminus(i-1)=velx(i) - velxplus(i-1)=velx(i) - velyminus(i-1)=vely(i) - velyplus(i-1)=vely(i) - velzminus(i-1)=velz(i) - velzplus(i-1)=velz(i) + velxminus(i-1)=vx(i) + velxplus(i-1)=vx(i) + velyminus(i-1)=vy(i) + velyplus(i-1)=vy(i) + velzminus(i-1)=vz(i) + velzplus(i-1)=vz(i) Bvcxminus(i-1)=Bvcx(i) Bvcxplus(i-1)=Bvcx(i) Bvcyminus(i-1)=Bvcy(i) @@ -539,9 +557,9 @@ end if if (cond) then call PPM_TVD(rho(i-1), rho(i), rho(i+1), rhominus(i), rhoplus(i)) - call PPM_TVD(velx(i-1), velx(i), velx(i+1), velxminus(i), velxplus(i)) - call PPM_TVD(vely(i-1), vely(i), vely(i+1), velyminus(i), velyplus(i)) - call PPM_TVD(velz(i-1), velz(i), velz(i+1), velzminus(i), velzplus(i)) + call PPM_TVD(vx(i-1), vx(i), vx(i+1), velxminus(i), velxplus(i)) + call PPM_TVD(vy(i-1), vy(i), vy(i+1), velyminus(i), velyplus(i)) + call PPM_TVD(vz(i-1), vz(i), vz(i+1), velzminus(i), velzplus(i)) call PPM_TVD(Bvcx(i-1), Bvcx(i), Bvcx(i+1), Bvcxminus(i), Bvcxplus(i)) call PPM_TVD(Bvcy(i-1), Bvcy(i), Bvcy(i+1), Bvcyminus(i), Bvcyplus(i)) call PPM_TVD(Bvcz(i-1), Bvcz(i), Bvcz(i+1), Bvczminus(i), Bvczplus(i)) @@ -560,12 +578,12 @@ if (cond) then rhominus(i)=rho(i) rhoplus(i)=rho(i) - velxminus(i)=velx(i) - velxplus(i)=velx(i) - velyminus(i)=vely(i) - velyplus(i)=vely(i) - velzminus(i)=velz(i) - velzplus(i)=velz(i) + velxminus(i)=vx(i) + velxplus(i)=vx(i) + velyminus(i)=vy(i) + velyplus(i)=vy(i) + velzminus(i)=vz(i) + velzplus(i)=vz(i) Bvcxminus(i)=Bvcx(i) Bvcxplus(i)=Bvcx(i) Bvcyminus(i)=Bvcy(i) @@ -574,12 +592,12 @@ Bvczplus(i)=Bvcz(i) rhominus(i+1)=rho(i) rhoplus(i+1)=rho(i) - velxminus(i+1)=velx(i) - velxplus(i+1)=velx(i) - velyminus(i+1)=vely(i) - velyplus(i+1)=vely(i) - velzminus(i+1)=velz(i) - velzplus(i+1)=velz(i) + velxminus(i+1)=vx(i) + velxplus(i+1)=vx(i) + velyminus(i+1)=vy(i) + velyplus(i+1)=vy(i) + velzminus(i+1)=vz(i) + velzplus(i+1)=vz(i) Bvcxminus(i+1)=Bvcx(i) Bvcxplus(i+1)=Bvcx(i) Bvcyminus(i+1)=Bvcy(i) @@ -605,9 +623,9 @@ end if if (cond) then call PPM_TVD(rho(i-1), rho(i), rho(i+1), rhominus(i), rhoplus(i)) - call PPM_TVD(velx(i-1), velx(i), velx(i+1), velxminus(i), velxplus(i)) - call PPM_TVD(vely(i-1), vely(i), vely(i+1), velyminus(i), velyplus(i)) - call PPM_TVD(velz(i-1), velz(i), velz(i+1), velzminus(i), velzplus(i)) + call PPM_TVD(vx(i-1), vx(i), vx(i+1), velxminus(i), velxplus(i)) + call PPM_TVD(vy(i-1), vy(i), vy(i+1), velyminus(i), velyplus(i)) + call PPM_TVD(vz(i-1), vz(i), vz(i+1), velzminus(i), velzplus(i)) call PPM_TVD(Bvcx(i-1), Bvcx(i), Bvcx(i+1), Bvcxminus(i), Bvcxplus(i)) call PPM_TVD(Bvcy(i-1), Bvcy(i), Bvcy(i+1), Bvcyminus(i), Bvcyplus(i)) call PPM_TVD(Bvcz(i-1), Bvcz(i), Bvcz(i+1), Bvczminus(i), Bvczplus(i)) @@ -621,6 +639,50 @@ end if end if end do + + !!$ transform back to vel if Wv was used! + if (reconstruct_Wv.ne.0) then + do i = grhydro_stencil, nx - grhydro_stencil + 1 + ! divide out the Loretnz factor obtained from w_lorentz = + ! sqrt(1+g_{ij} w^i w^j) for both the + ! plus and minus quantities this should by construction ensure + ! that any Lorentz factor calculated from them later on is + ! physical (ie. > 1.d0) + agxx = 0.5d0*( gxx(i) + gxx(i-1) ) + agxy = 0.5d0*( gxy(i) + gxy(i-1) ) + agxz = 0.5d0*( gxz(i) + gxz(i-1) ) + agyy = 0.5d0*( gyy(i) + gyy(i-1) ) + agyz = 0.5d0*( gyz(i) + gyz(i-1) ) + agzz = 0.5d0*( gzz(i) + gzz(i-1) ) + w = sqrt( 1.d0 + agxx*velxminus(i)*velxminus(i) + & + agyy*velyminus(i)*velyminus(i) & + + agzz*velzminus(i)*velzminus(i) + & + 2.d0*agxy*velxminus(i)*velyminus(i) & + + 2.d0*agxz*velxminus(i)*velzminus(i) + & + 2.d0*agyz*velyminus(i)*velzminus(i) ) + velxminus(i) = velxminus(i)/w + velyminus(i) = velyminus(i)/w + velzminus(i) = velzminus(i)/w + + agxx = 0.5d0*( gxx(i) + gxx(i+1) ) + agxy = 0.5d0*( gxy(i) + gxy(i+1) ) + agxz = 0.5d0*( gxz(i) + gxz(i+1) ) + agyy = 0.5d0*( gyy(i) + gyy(i+1) ) + agyz = 0.5d0*( gyz(i) + gyz(i+1) ) + agzz = 0.5d0*( gzz(i) + gzz(i+1) ) + w = sqrt( 1.d0 + agxx*velxplus(i)*velxplus(i) + & + agyy*velyplus(i)*velyplus(i) & + + agzz*velzplus(i)*velzplus(i) + & + 2.d0*agxy*velxplus(i)*velyplus(i) & + + 2.d0*agxz*velxplus(i)*velzplus(i) + & + 2.d0*agyz*velyplus(i)*velzplus(i) ) + velxplus(i) = velxplus(i)/w + velyplus(i) = velyplus(i)/w + velzplus(i) = velzplus(i)/w + end do + end if + + return end subroutine SimplePPM_1dM

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 433)
by rhaas＠tapir.caltech.edu 08 Nov '12

08 Nov '12

User: rhaas Date: 2012/11/08 07:53 PM Modified: /trunk/src/ GRHydro_SourceM.F90, GRHydro_TmunuM.F90 Log: GRHydro: make sure to pass distinct dummy arguments to Fortran routines since Fortran does not allow aliased arguments From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_SourceM.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_SourceM.F90 2012-11-09 01:53:05 UTC (rev 432) +++ trunk/src/GRHydro_SourceM.F90 2012-11-09 01:53:08 UTC (rev 433) @@ -97,7 +97,7 @@ CCTK_REAL :: dx_det_bydet, dy_det_bydet, dz_det_bydet CCTK_REAL :: gdg_x, gdg_y, gdg_z !! g^{ik} d_k g_{ij} - CCTK_REAL :: Bvecxlow,Bvecylow,Bveczlow,bdotv,b2,dum,bxlow,bylow,bzlow + CCTK_REAL :: Bvecxlow,Bvecylow,Bveczlow,bdotv,b2,dum1,dum2,bxlow,bylow,bzlow CCTK_REAL :: bt,bx,by,bz,rhohstarW2,pstar logical, allocatable, dimension (:,:,:) :: force_spatial_second_order @@ -254,7 +254,7 @@ !$OMP dx_betax,dx_betay,dx_betaz,dy_betax,dy_betay,dy_betaz,& !$OMP dz_betax,dz_betay,dz_betaz,velxshift,velyshift,velzshift,& !$OMP vlowx,vlowy,vlowz,Bvecxlow,Bvecylow,Bveczlow, & - !$OMP bdotv,b2,dum,bxlow,bylow,bzlow,bt,bx,by,bz,rhohstarW2,pstar,& + !$OMP bdotv,b2,dum1,dum2,bxlow,bylow,bzlow,bt,bx,by,bz,rhohstarW2,pstar,& !$OMP t00,t0x,t0y,t0z,txx,txy,txz,tyy,tyz,tzz,& !$OMP dx_alp,dy_alp,dz_alp,& !$OMP tau_source,sx_source,sy_source,sz_source,& @@ -336,7 +336,7 @@ call calc_vlow_blow(localgxx,localgxy,localgxz,localgyy,localgyz,localgzz, & velx(i,j,k),vely(i,j,k),velz(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k), & vlowx,vlowy,vlowz,Bvecxlow,Bvecylow,Bveczlow, & - bdotv,b2,dum,dum,bxlow,bylow,bzlow) + bdotv,b2,dum1,dum2,bxlow,bylow,bzlow) !!$ These are the contravariant components bt = w_lorentz(i,j,k)/alp(i,j,k)*bdotv File [modified]: GRHydro_TmunuM.F90 Delta lines: +4 -3 =================================================================== --- trunk/src/GRHydro_TmunuM.F90 2012-11-09 01:53:05 UTC (rev 432) +++ trunk/src/GRHydro_TmunuM.F90 2012-11-09 01:53:08 UTC (rev 433) @@ -81,14 +81,14 @@ CCTK_REAL :: velxlow, velylow, velzlow CCTK_REAL :: betaxlow, betaylow, betazlow, beta2 CCTK_REAL :: Bvecxlow,Bvecylow,Bveczlow - CCTK_REAL :: bdotv,b2,bxlow,bylow,bzlow,btlow,dum + CCTK_REAL :: bdotv,b2,bxlow,bylow,bzlow,btlow,dum1,dum2 CCTK_REAL :: utlow,rhohstarw2,pstar CCTK_REAL :: bdotbeta,vdotbeta CCTK_INT :: i,j,k !$OMP PARALLEL DO PRIVATE(i,j,k,velxlow, velylow, velzlow,& - !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum,b2,bxlow,bylow,bzlow,& + !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum1,dum2,b2,bxlow,bylow,bzlow,& !$OMP betaxlow, betaylow, betazlow, beta2, bdotbeta,vdotbeta,utlow, btlow,& !$OMP rhohstarw2,pstar) @@ -96,11 +96,12 @@ do j = 1, cctk_lsh(2) do i = 1, cctk_lsh(1) + ! need separate dum1, dum2 b/c of Fortrans aliasing rules call calc_vlow_blow(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),& gyy(i,j,k),gyz(i,j,k),gzz(i,j,k), & velx(i,j,k),vely(i,j,k),velz(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k), & velxlow,velylow,velzlow,Bvecxlow,Bvecylow,Bveczlow, & - bdotv,b2,dum,dum,bxlow,bylow,bzlow) + bdotv,b2,dum1,dum2,bxlow,bylow,bzlow) !!$ Calculate lower components and square of shift vector.

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 432)
by rhaas＠tapir.caltech.edu 08 Nov '12

08 Nov '12

User: rhaas Date: 2012/11/08 07:53 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_Prim2ConM.F90 Log: GRHydro: add extensive information output to check_C2P_failed From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +55 -17 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2012-11-09 01:53:00 UTC (rev 431) +++ trunk/src/GRHydro_Con2Prim.F90 2012-11-09 01:53:05 UTC (rev 432) @@ -2231,18 +2231,8 @@ DECLARE_CCTK_ARGUMENTS integer :: i, j, k, nx, ny, nz - character(len=100) warnline + character(len=300) warnline - ! save memory when MP is not used - CCTK_INT :: GRHydro_UseGeneralCoordinates - CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup - - if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then - vup => lvel - else - vup => vel - end if - ! call CCTK_INFO("Checking the C2P failure mask.") nx = cctk_lsh(1) @@ -2273,16 +2263,64 @@ !$OMP CRITICAL call CCTK_WARN(1,'Con2Prim failed; stopping the code.') call CCTK_WARN(1,'Even with mesh refinement, this point is not restricted from a finer level, so this is really an error') - write(warnline,'(a28,i2)') 'on carpet reflevel: ',GRHydro_reflevel + write(warnline,'(a32,i2)') 'on carpet reflevel: ',GRHydro_reflevel call CCTK_WARN(1,warnline) - write(warnline,'(a20,3g16.7)') 'xyz location: ',& + write(warnline,'(a32,3g16.7)') 'xyz location: ',& x(i,j,k),y(i,j,k),z(i,j,k) call CCTK_WARN(1,warnline) - write(warnline,'(a20,g16.7)') 'radius: ',r(i,j,k) + write(warnline,'(a32,g16.7)') 'radius: ',r(i,j,k) call CCTK_WARN(1,warnline) - write(warnline,'(a20,3g16.7)') 'velocities: ',& - vup(i,j,k,1),vup(i,j,k,2),vup(i,j,k,3) - call CCTK_WARN(0,warnline) + write(warnline,'(a32,i3,1g16.7)') 'hydro_excision_mask, C2Pfailed: ', hydro_excision_mask(i,j,k), GRHydro_C2P_failed(i,j,k) + call CCTK_WARN(1,warnline) + write(warnline,'(a32,4g16.7)') 'rho eps press w_lorentz: ',& + rho(i,j,k),eps(i,j,k),press(i,j,k),w_lorentz(i,j,k) + call CCTK_WARN(1,warnline) + write(warnline,'(a32,3g16.7)') 'velocities: ',& + vel(i,j,k,1),vel(i,j,k,2),vel(i,j,k,3) + call CCTK_WARN(1,warnline) + write(warnline,'(a32,2g16.7)') 'dens tau: ',& + dens(i,j,k),tau(i,j,k) + call CCTK_WARN(1,warnline) + write(warnline,'(a32,3g16.7)') 'scon: ',& + scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3) + call CCTK_WARN(1,warnline) + if (evolve_mhd.ne.0) then + write(warnline,'(a32,3g16.7)') 'magnetic field: ',& + Bvec(i,j,k,1),Bvec(i,j,k,2),Bvec(i,j,k,3) + call CCTK_WARN(1,warnline) + write(warnline,'(a32,3g16.7)') 'cons magnetic field: ',& + Bcons(i,j,k,1),Bcons(i,j,k,2),Bcons(i,j,k,3) + call CCTK_WARN(1,warnline) + if (calculate_bcom.ne.0) then + write(warnline,'(a32,6g16.7)') 'com magnetic field, b0, b^2: ', bcom(i,j,k,1), bcom(i,j,k,2), bcom(i,j,k,3), bcom0(i,j,k), bcom_sq(i,j,k) + call CCTK_WARN(1,warnline) + end if + if (clean_divergence.ne.0) then + write(warnline,'(a32,g16.7)') 'psidc: ',psidc(i,j,k) + call CCTK_WARN(1,warnline) + end if + if (track_divB.ne.0) then + write(warnline,'(a32,g16.7)') 'divB: ',divB(i,j,k) + call CCTK_WARN(1,warnline) + end if + end if + if (evolve_temper.ne.0) then + write(warnline,'(a32,2g16.7)') 'temperature, entropy: ',& + temperature(i,j,k),entropy(i,j,k) + call CCTK_WARN(1,warnline) + end if + if (evolve_Y_e.ne.0) then + write(warnline,'(a32,2g16.7)') 'Y_e, Y_e_con: ',& + Y_e(i,j,k),Y_e_con(i,j,k) + call CCTK_WARN(1,warnline) + end if + write(warnline,'(a32,6g16.7)') '3-metric: ',& + gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k) + call CCTK_WARN(1,warnline) + write(warnline,'(a32,4g16.7)') 'lapse, shift: ',& + alp(i,j,k),betax(i,j,k),betay(i,j,k),betaz(i,j,k) + call CCTK_WARN(1,warnline) + call CCTK_WARN(0,"Aborting.") !$OMP END CRITICAL end if File [modified]: GRHydro_Prim2ConM.F90 Delta lines: +2 -0 =================================================================== --- trunk/src/GRHydro_Prim2ConM.F90 2012-11-09 01:53:00 UTC (rev 431) +++ trunk/src/GRHydro_Prim2ConM.F90 2012-11-09 01:53:05 UTC (rev 432) @@ -386,6 +386,8 @@ !!$ BEGIN: Check for NaN value if (w .ne. w) then !$OMP CRITICAL + write(NaN_WarnLine,'(a100,6g15.6)') 'NaN produced in sqrt(): (gxx,gxy,gxz,gyy,gyz,gzz)', gxx, gxy, gxz, gyy, gyz, gzz + call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine) write(NaN_WarnLine,'(a100,3g15.6)') 'NaN produced in sqrt(): (dvelx,dvely,dvelz)', dvelx, dvely, dvelz call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine) !$OMP END CRITICAL

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 431)
by rhaas＠tapir.caltech.edu 08 Nov '12

08 Nov '12

User: rhaas Date: 2012/11/08 07:53 PM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_EoSChangeGamma.F90 Log: GRHydro: remove usage of EOS_Module, poly_gamma_ini From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_EoSChangeGamma.F90 Delta lines: +8 -12 =================================================================== --- trunk/src/GRHydro_EoSChangeGamma.F90 2012-11-05 14:45:59 UTC (rev 430) +++ trunk/src/GRHydro_EoSChangeGamma.F90 2012-11-09 01:53:00 UTC (rev 431) @@ -39,8 +39,6 @@ subroutine GRHydro_EoSChangeGamma(CCTK_ARGUMENTS) - USE EOS_Omni_Module, only: press_gf, inv_rho_gf, poly_k_cgs, rho_gf - implicit none ! save memory when MP is not used @@ -95,8 +93,8 @@ end if local_Gamma = 1.0d0 + xpress/xeps - press = press_gf * poly_k_cgs * & - (rho * inv_rho_gf)**local_Gamma + press = poly_k * & + rho**local_Gamma eps = press / (rho * (local_Gamma - 1.d0)) !!$ Change the pressure and specific internal energy @@ -277,8 +275,6 @@ subroutine GRHydro_EoSChangeGammaK_Shibata(CCTK_ARGUMENTS) - USE EOS_Omni_Module, only: press_gf, inv_rho_gf, poly_k_cgs, rho_gf - implicit none ! save memory when MP is not used @@ -293,7 +289,7 @@ CCTK_INT :: i, j, k CCTK_REAL :: det - CCTK_REAL :: local_Gamma, local_k, eos_k_initial_cgs + CCTK_REAL :: local_Gamma, local_k, eos_k_initial CCTK_REAL, dimension(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)) :: Q @@ -340,13 +336,13 @@ local_Gamma = 1.0d0 + xpress/xeps local_K = xpress - eos_k_initial_cgs = initial_k * rho_gf**initial_Gamma / press_gf + eos_k_initial = initial_k - press = (local_Gamma - 1.d0) / (initial_Gamma - 1.0d0 ) * press_gf * eos_k_initial_cgs * & - (rho * rho_gf) ** initial_Gamma + press = (local_Gamma - 1.d0) / (initial_Gamma - 1.0d0 ) * eos_k_initial * & + rho ** initial_Gamma - eps = press_gf * eos_k_initial_cgs * & - (rho * inv_rho_gf) ** initial_Gamma / & + eps = eos_k_initial * & + rho ** initial_Gamma / & (rho * (initial_Gamma - 1.0d0)) do k = 1, cctk_lsh(3) Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +8 -0 =================================================================== --- trunk/param.ccl 2012-11-05 14:45:59 UTC (rev 430) +++ trunk/param.ccl 2012-11-09 01:53:00 UTC (rev 431) @@ -64,6 +64,14 @@ USES INT spatial_order +######################################### +### PARAMETRES SHARED FROM EOS_Omni ### +######################################### + +SHARES: EOS_Omni + +USES REAL poly_k + restricted: BOOLEAN GRHydro_enable_internal_excision "Set this to 'false' to disable the thorn-internal excision." STEERABLE=ALWAYS

1 0

[svn:einsteintoolkit] www/download/ (Rev. 913)
by knarf＠cct.lsu.edu 08 Nov '12

08 Nov '12

User: knarf Date: 2012/11/08 11:06 AM Modified: /download/ index.php Log: update for new release File Changes: Directory: /download/ ===================== File [modified]: index.php Delta lines: +3 -3 =================================================================== --- download/index.php 2012-11-08 17:02:00 UTC (rev 912) +++ download/index.php 2012-11-08 17:06:40 UTC (rev 913) @@ -15,10 +15,10 @@ commands listed below. This will create a directory called Cactus in which the components of the Einstein Toolkit are downloaded. -<h4>Current release: Lovelace (released on 28th May 2012)</h4> +<h4>Current release: Ãrsted (released on November 8th, 2012)</h4> This is the recommended version of the toolkit for most users. See -the <a href = "../about/releases/ET_2012_05_announcement.php">release +the <a href = "../about/releases/ET_2012_11_announcement.php">release notes</a> for more information. <pre> @@ -26,7 +26,7 @@ chmod a+x GetComponents -./GetComponents --parallel https://svn.einsteintoolkit.org/manifest/branches/ET_2012_05/einsteintoolki… +./GetComponents --parallel https://svn.einsteintoolkit.org/manifest/branches/ET_2012_11/einsteintoolki… </pre> <h4>Current development version</h4>

1 0

[svn:einsteintoolkit] www/ (Rev. 912)
by knarf＠cct.lsu.edu 08 Nov '12

08 Nov '12

User: knarf Date: 2012/11/08 11:02 AM Modified: / index.php Log: update main page File Changes: Directory: / ============ File [modified]: index.php Delta lines: +4 -4 =================================================================== --- index.php 2012-11-08 16:59:30 UTC (rev 911) +++ index.php 2012-11-08 17:02:00 UTC (rev 912) @@ -13,10 +13,10 @@ <a class = "downloadbutton" href = "download/">Download</a> -May 2012: We are pleased to -<a href="about/releases/ET_2012_05_announcement.php">announce the -fifth release</a> (code name -"<a href="http://en.wikipedia.org/wiki/Ada_Lovelace">Lovelace</a>") +November 2012: We are pleased to +<a href="about/releases/ET_2012_11_announcement.php">announce the +sixth release</a> (code name +"<a href="http://en.wikipedia.org/wiki/Hans_Christian_%C3%98rsted">Ãrsted</a>") of the Einstein Toolkit, an open, community developed software infrastructure for relativistic astrophysics.

1 0

[svn:einsteintoolkit] www/about/releases/ (Rev. 911)
by knarf＠cct.lsu.edu 08 Nov '12

08 Nov '12

User: knarf Date: 2012/11/08 10:59 AM Modified: /about/releases/ index.php Log: update for new release name File Changes: Directory: /about/releases/ =========================== File [modified]: index.php Delta lines: +4 -8 =================================================================== --- about/releases/index.php 2012-11-08 16:54:34 UTC (rev 910) +++ about/releases/index.php 2012-11-08 16:59:30 UTC (rev 911) @@ -2,8 +2,8 @@ include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> The most recent release of the Einstein Toolkit - is <tt>ET_2012_05</tt>, code named - "<a href="http://en.wikipedia.org/wiki/Ada_Lovelace">Lovelace</a>". + is <tt>ET_2012_11</tt>, code named + "<a href="http://en.wikipedia.org/wiki/Hans_Christian_%C3%98rsted">Ãrsted</a>". We expect this version to remain unchanged, save for correcting serious problems that may be detected. This version should be your first choice if you are new to the Einstein Toolkit, or if you want @@ -18,6 +18,8 @@ <hrule /> List of Einstein Toolkit releases: <ul> + <li><tt>ET_2012_11</tt> "Ãrsted", + <a href="ET_2012_11_announcement.php">released 2012-11-08</a></li> <li><tt>ET_2012_05</tt> "Lovelace", <a href="ET_2012_05_announcement.php">released 2012-05-28</a></li> <li><tt>ET_2011_10</tt> "Maxwell", @@ -29,11 +31,5 @@ <li><tt>ET_2010_06</tt> "Bohr", <a href="ET_2010_06_announcement.php">released 2010-06-17</a></li> </ul> - <?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?>

1 0

[svn:einsteintoolkit] www/about/releases/ (Rev. 910)
by knarf＠cct.lsu.edu 08 Nov '12

08 Nov '12

User: knarf Date: 2012/11/08 10:54 AM Modified: /about/releases/ ET_2012_11_announcement.php Log: boldface sub-headings File Changes: Directory: /about/releases/ =========================== File [modified]: ET_2012_11_announcement.php Delta lines: +5 -5 =================================================================== --- about/releases/ET_2012_11_announcement.php 2012-11-08 16:50:56 UTC (rev 909) +++ about/releases/ET_2012_11_announcement.php 2012-11-08 16:54:34 UTC (rev 910) @@ -50,7 +50,7 @@ <h4>The changes between this and the previous release include:</h4> -<ul><li>Compatibility</li></ul> +<ul><li>Compatibility</li></ul> The following changes lead to potential different numerical results or the need to modify optionlists: <ul> @@ -59,7 +59,7 @@ <li> McLachlan: Change the default value of apply_dissipation to "never". </ul> -<ul><li>Changes since last release</li></ul> +<ul><li>Changes since last release</li></ul> <ul> <li> GetComponents: <ul> @@ -152,13 +152,13 @@ <li> hwloc: new thorn for libhwloc </ul> -<ul><li>How to upgrade from Lovelace (ET_2012_05)</li></ul> +<ul><li>How to upgrade from Lovelace (ET_2012_05)</li></ul> To upgrade from the previous release, use GetComponents with the new thornlist to check out the new version. Due to the change of Carpet repository type, using GetComponents to upgrade a Lovelace checkout is unlikely to work. See the <a href="http://einsteintoolkit.org/download/">Download</a> page on the Einstein Toolkit website for download instructions. -<ul><li>Remaining issues with this release</li></ul> +<ul><li>Remaining issues with this release</li></ul> <ul> <li> Certain machines need to be configured specially in Simfactory because the remote directories cannot be determined automatically just from the username. See the Machine notes below. @@ -166,7 +166,7 @@ <li> Recovering with Carpet: The maximum number of timelevels that can be recovered is Carpet::prolongation_order_time+1. This is usually the case, but it is possible to write parameter files e.g. with prolongation_order_time=1 that use 3 timelevels. This bug manifests in an assert() failure when recovering from checkpoints. The work-around is to either increase Carpet::prolongation_order_time or decrease the number of timelevels for the grid function in question accordingly. </ul> -<ul><li>Machine notes</li></ul> +<ul><li>Machine notes</li></ul> <ul> <li>Kraken:

1 0

[svn:einsteintoolkit] www/about/releases/ (Rev. 909)
by knarf＠cct.lsu.edu 08 Nov '12

08 Nov '12

User: knarf Date: 2012/11/08 10:50 AM Added: /about/releases/ ET_2012_11_announcement.php, ET_2012_11_announcement.txt Log: release announcement files File Changes: Directory: /about/releases/ =========================== File [added]: ET_2012_11_announcement.php Delta lines: +200 -0 =================================================================== --- about/releases/ET_2012_11_announcement.php (rev 0) +++ about/releases/ET_2012_11_announcement.php 2012-11-08 16:50:56 UTC (rev 909) @@ -0,0 +1,200 @@ +<?php $title='Release Announcement'; +include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> + +We are pleased to announce the sixth release (code name +"<a href="http://en.wikipedia.org/wiki/Hans_Christian_%C3%98rsted">Ãrsted</a>") of the +Einstein Toolkit, an open, community developed software infrastructure for +relativistic astrophysics. This release includes various improvements to the +MHD handling within GRHydro, and improvements to the Cactus flesh. In addition, +bug fixes accumulated since the previous release in May 2012 have been +included. + +The Einstein Toolkit is a collection of software components and tools for +simulating and analyzing general relativistic astrophysical systems that builds +on numerous software efforts in the numerical relativity community including +CactusEinstein, the Carpet AMR infrastructure and the relativistic +hydrodynamics code GRHydro (an updated and extended version of the public +release of the Whisky code). The Cactus Framework is used as the underlying +computational infrastructure providing large-scale parallelization, general +computational components, and a model for collaborative, portable code +development. The toolkit includes modules to build complete codes for +simulating black hole spacetimes as well as systems governed by relativistic +hydrodynamics. + +The Einstein Toolkit uses a distributed software model and its different +modules are developed, distributed, and supported either by the core team of +Einstein Toolkit Maintainers, or by individual groups. Where modules are +provided by external groups, the Einstein Toolkit Maintainers provide quality +control for modules for inclusion in the toolkit and help coordinate support. +The Einstein Toolkit Maintainers currently involve postdocs and faculty from +six different institutions, and host weekly meetings that are open for anyone +to join in. + +Guiding principles for the design and implementation of the toolkit include: +open, community-driven software development; well thought out and stable +interfaces; separation of physics software from computational science +infrastructure; provision of complete working production code; training and +education for a new generation of researchers. + +For more information about using or contributing to the Einstein Toolkit, or to +join the Einstein Toolkit Consortium, please visit our web pages at +<<a href="http://einsteintoolkit.org">http://einsteintoolkit.org</a>>. + +The Einstein Toolkit is primarily supported by NSF +1212401/1212426/1212433/1212460, 0903973/0903782/0904015 (CIGR), and also by +0905046/0941653 (PetaCactus) and 0710874 (LONI Grid). + +The Einstein Toolkit contain over 170 regression test cases. On a large +portion of the tested machines, all of these testsuites pass, using both +MPI and OpenMP parallelization. + +<h4>The changes between this and the previous release include:</h4> + +<ul><li>Compatibility</li></ul> + +The following changes lead to potential different numerical results or the need to modify optionlists: +<ul> + <li> Cactus: Move MPI support from flesh to a thorn; use new MPI_DIR variable to control use of MPI (see below) + <li> McLachlan: Correct some errors in CCZ4, in particular the GammaShift parameter + <li> McLachlan: Change the default value of apply_dissipation to "never". +</ul> + +<ul><li>Changes since last release</li></ul> +<ul> + <li> GetComponents: + <ul> + <li> anonymous checkout is now used by default. Use --no-anonymous for authenticated checkouts + <li> no longer use "shallow git clones" by default + </ul> + <li> simfactory + <ul> + <li> updates of various machine configurations + <li> Replace large files by hard links during cleanup + <li> Basic configuration file for Ubuntu, Debian and Fedora operating systems. + </ul> + <li> utils: syntax highlight for vim and emacs + <li> Cactus flesh: + <ul> + <li> grid functions may be padded, i.e. the driver may allocate more grid points than requested. The allocated size of grid functions can be accessed using cctk_ash. Note however, that this can be different from cctk_lsh. The latter is that part (size) of cctk_ash that is used within a simulation, i.e. that was set as simulation size via parameters. This is necessary for vectorization where the size of arrays needs to be a multiple of some number, but the simulation domain might not be. There is no padding by default. + <li> correctly detect circular schedule items + <li> replace CCTK_RESTRICT by restrict, which is now defined by the flesh + <li> Change per thorn -DTHORN_IS_xxx to a per thorn -I bindings/include/xxx + <li> Correct error in determining dependencies + <li> have -reconfig target execute -cleandeps after reconfiguration + <li> add parameters to reverse un-enforced Cactus schedule ordering + <li> Move MPI support from flesh to a thorn + <li> Support various new Darwin minor versions + <li> increase default value of RELTOL to 1e-12 + <li> Enforce that arrangement names are legal C identifiers + <li> Update Fortran API for CCTK_LOOP macros + <li> Add autoconf check whether C99 is supported (now required) + <li> New option -P (--exit-after-param-check) to exit after param check + <li> allow ${parfile} in parameter files to allow ${parfile}/subdir + </ul> + <li> Carpet + <ul> + <li> Now lives in git again + <li> A parameter CarpetLib::interpolate_from_buffer_zones was introduced which, in order to get the correct behaviour (not to use the buffer zones for interpolations), should be set to "no". However, since this changes the way things had been done in the past and this change isn't well tested yet, the default is still "yes" (as like before the parameter was introduced). The default will very likely be changed for the next Einstein Toolkit release, after more testing. + <li> Added "balanced" domain decomposition + <li> Update setting CPU affinity + <li> change defaults of periodic_[xyz] to "yes" to match PUGH + <li> CarpetIOHDF5: add option checkpoint_every_divisor + <li> CarpetIOHDF5: teach the checkpointer to write index files + <li> CarpetIOHDF5: add support for index files to file reader + <li> LoopControl: Change printstats default to "no" + <li> CarpetInterp2: New ENO2 interpolator. + <li> Support more than one hierarchy of timers + <li> Add timers for refinement levels + <li> CarpetLib: Support higher-then 2nd order restriction operator for cell-centered meshes for Langrange prolongated grid functions. + </ul> + <li> EinsteinExact: New arrangement, to eventually replace the Exact thorn. The thorns in this arrangement are generated by Kranc + <li> McLachlan + <ul> + <li> Correct some errors in CCZ4, in particular the GammaShift parameter + <li> Change the default value of apply_dissipation to "never". + </ul> + <li> ExternalLibraries: MPI is now treated like any other library on which flesh and thorns may depend, such as e.g. the GSL. It is thus configured via thorn ExternalLibraries/MPI, and not via the flesh any more. In particular, the option "MPI = ..." does not work any more and must not be used; instead, use e.g. "MPI_DIR = <path>" to point to an existing installation. You may also say "MPI_DIR = BUILD" to have Cactus build OpenMPI automatically, which is a convenient way of using MPI on a system where it may not be installed. + <li> LocalInterp: The interpolation loop uses OpenMP + <li> MoL: add multirate capabilities and three new multirate rk schemes + <li> SphericalSurface: make maxnphi and maxntheta parameters steerable + <li> NaNChecker: Add options to use restriction mask to ignore unused points, output more information about the failed points, use OpenMP during NaN check + <li> Multipole: Implement Driscoll&Healy integration + <li> HydroBase: schedule hydrobase_InitExcisionmask global-early loop-local + <li> TmunuBase: support_old_CalcTmunu_mechanism is now deprecated and is likely to be removed soon + <li> EOS_Omni: make more parameters steerable, add option read_table_on_single_process to have a named processor read the full table and MPI_Bcast it to the other processors + <li> GRHydro + <ul> + <li> C2P: resort to bisection if Newton-Raphson fails. This is only activated if user requested this + <li> introduce parameters GRHydro_Y_e_min and GRHydro_Y_e_max to fix issues with EOS tables that don't extend sufficiently far + <li> Fixes to atmosphere mask and the way we sync + <li> various fixes to ePPM + <li> Option to use the slow sector of multirate RK methods + <li> fixed issues with hot EOS treatment + <li> Option to reconstruct W*vel instead of vel + <li> use an early atmosphere reset to ensure past timelevels are initialized correctly + <li> various MHD fixes + </ul> + <li> Exact: openmp fixes + <li> GRHydro_InitData: new poloidal magnetic field routine, Added option to smooth metric inside horizon, Alfven wave test + <li> SummationByParts: improve OpenMP parallelization + <li> PITTNull*: respect Cactus' IO::truncate parameter + <li> AEIThorns/Trigger now part of the Einstein Toolkit + +External libraries: + <li> MPI: new - use this instead of the flesh mechanism (see above) + <li> BLAS: update to 3.4.2 + <li> FFTW3: update to 3.3.2 + <li> HDF5: will likely not build fortran and C++ api for next release by default (but still does now) + <li> LAPACK: update to 3.4.2 + <li> LORENE: don't use OpenMP unless it's enabled within Cactus + <li> OpenSSL: improve search for existing library + <li> zlib: update to 1.2.7 + <li> hwloc: new thorn for libhwloc +</ul> + +<ul><li>How to upgrade from Lovelace (ET_2012_05)</li></ul> + +To upgrade from the previous release, use GetComponents with the new thornlist to check out the new version. Due to the change of Carpet repository type, using GetComponents to upgrade a Lovelace checkout is unlikely to work. + +See the <a href="http://einsteintoolkit.org/download/">Download</a> page on the Einstein Toolkit website for download instructions. + +<ul><li>Remaining issues with this release</li></ul> + +<ul> + <li> Certain machines need to be configured specially in Simfactory because the remote directories cannot be determined automatically just from the username. See the Machine notes below. + <li> CarpetIOHDF5: When the new parameter CarpetIOHDF5::output_symmetry_points is set to "yes", then symmetry points are not actually output if CarpetIOHDF5::output_buffer_points has been set to "no". + <li> Recovering with Carpet: The maximum number of timelevels that can be recovered is Carpet::prolongation_order_time+1. This is usually the case, but it is possible to write parameter files e.g. with prolongation_order_time=1 that use 3 timelevels. This bug manifests in an assert() failure when recovering from checkpoints. The work-around is to either increase Carpet::prolongation_order_time or decrease the number of timelevels for the grid function in question accordingly. +</ul> + +<ul><li>Machine notes</li></ul> + +<ul> + <li>Kraken: +defs.local.ini needs to have sourcebasedir = $HOME configured for this machine. You need to determine $HOME by logging in to the machine.</li> + <li>LoneStar and Ranger: +defs.local.ini needs to have sourcebasedir = $WORK and basedir = $SCRATCH/simulations configured for this machine. You need to determine $WORK and $SCRATCH by logging in to the machine.</li> + <li>Pandora: +The LSU system pandora has at the time of writing not only severe problems of compiling Cactus (can crash the filesystem), but in part due to this also outstanding are some failing testsuite results. Please contact Peter Diener at CCT (LSU) before trying to use this machine.</li> +</ul> + +All repositories participating in this release carry a branch ET_2012_11 +marking this release. These release branches will be updated if severe +errors are found. + +The "Ãrsted" Release Team on behalf of the Einstein Toolkit Consortium (2012-11-08) + + +   Eloisa Bentivegna +   Tanja Bode +   Peter Diener +   Roland Haas +   Ian Hinder +   Frank LÃ¶ffler +   Bruno Mundim +   Christian D. Ott +   Erik Schnetter + + +November 8, 2012 + +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?> File [added]: ET_2012_11_announcement.txt Delta lines: +183 -0 =================================================================== --- about/releases/ET_2012_11_announcement.txt (rev 0) +++ about/releases/ET_2012_11_announcement.txt 2012-11-08 16:50:56 UTC (rev 909) @@ -0,0 +1,183 @@ +We are pleased to announce the sixth release (code name "Ãrsted") of the +Einstein Toolkit, an open, community developed software infrastructure for +relativistic astrophysics. This release includes various improvements to the +MHD handling within GRHydro, and improvements to the Cactus flesh. In addition, +bug fixes accumulated since the previous release in May 2012 have been +included. + +The Einstein Toolkit is a collection of software components and tools for +simulating and analyzing general relativistic astrophysical systems that builds +on numerous software efforts in the numerical relativity community including +CactusEinstein, the Carpet AMR infrastructure and the relativistic +hydrodynamics code GRHydro (an updated and extended version of the public +release of the Whisky code). The Cactus Framework is used as the underlying +computational infrastructure providing large-scale parallelization, general +computational components, and a model for collaborative, portable code +development. The toolkit includes modules to build complete codes for +simulating black hole spacetimes as well as systems governed by relativistic +hydrodynamics. + +The Einstein Toolkit uses a distributed software model and its different +modules are developed, distributed, and supported either by the core team of +Einstein Toolkit Maintainers, or by individual groups. Where modules are +provided by external groups, the Einstein Toolkit Maintainers provide quality +control for modules for inclusion in the toolkit and help coordinate support. +The Einstein Toolkit Maintainers currently involve postdocs and faculty from +six different institutions, and host weekly meetings that are open for anyone +to join in. + +Guiding principles for the design and implementation of the toolkit include: +open, community-driven software development; well thought out and stable +interfaces; separation of physics software from computational science +infrastructure; provision of complete working production code; training and +education for a new generation of researchers. + +For more information about using or contributing to the Einstein Toolkit, or to +join the Einstein Toolkit Consortium, please visit our web pages at +<http://einsteintoolkit.org>. + +The Einstein Toolkit is primarily supported by NSF +1212401/1212426/1212433/1212460, 0903973/0903782/0904015 (CIGR), and also by +0905046/0941653 (PetaCactus) and 0710874 (LONI Grid). + +The Einstein Toolkit contain over 170 regression test cases. On a large +portion of the tested machines, all of these testsuites pass, using both +MPI and OpenMP parallelization. + +The changes between this and the previous release include: + +=== Compatibility === + +The following changes lead to potential different numerical results or the need to modify optionlists: + +* Cactus: Move MPI support from flesh to a thorn; use new MPI_DIR variable to control use of MPI (see below) +* McLachlan: Correct some errors in CCZ4, in particular the GammaShift parameter +* McLachlan: Change the default value of apply_dissipation to "never". + +=== Changes since last release === + +* GetComponents: +** anonymous checkout is now used by default. Use --no-anonymous for authenticated checkouts +** no longer use "shallow git clones" by default +* simfactory +** updates of various machine configurations +** Replace large files by hard links during cleanup +** Basic configuration file for Ubuntu, Debian and Fedora operating systems. +* utils: syntax highlight for vim and emacs +* Cactus flesh: +** grid functions may be padded, i.e. the driver may allocate more grid points than requested. The allocated size of grid functions can be accessed using cctk_ash. Note however, that this can be different from cctk_lsh. The latter is that part (size) of cctk_ash that is used within a simulation, i.e. that was set as simulation size via parameters. This is necessary for vectorization where the size of arrays needs to be a multiple of some number, but the simulation domain might not be. There is no padding by default. +** correctly detect circular schedule items +** replace CCTK_RESTRICT by restrict, which is now defined by the flesh +** Change per thorn -DTHORN_IS_xxx to a per thorn -I bindings/include/xxx +** Correct error in determining dependencies +** have -reconfig target execute -cleandeps after reconfiguration +** add parameters to reverse un-enforced Cactus schedule ordering +** Move MPI support from flesh to a thorn +** Support various new Darwin minor versions +** increase default value of RELTOL to 1e-12 +** Enforce that arrangement names are legal C identifiers +** Update Fortran API for CCTK_LOOP macros +** Add autoconf check whether C99 is supported (now required) +** New option -P (--exit-after-param-check) to exit after param check +** allow ${parfile} in parameter files to allow ${parfile}/subdir +* Carpet +** Now lives in git again +** A parameter CarpetLib::interpolate_from_buffer_zones was introduced which, in order to get the correct behaviour (not to use the buffer zones for interpolations), should be set to "no". However, since this changes the way things had been done in the past and this change isn't well tested yet, the default is still "yes" (as like before the parameter was introduced). The default will very likely be changed for the next Einstein Toolkit release, after more testing. +** Added "balanced" domain decomposition +** Update setting CPU affinity +** change defaults of periodic_[xyz] to "yes" to match PUGH +** CarpetIOHDF5: add option checkpoint_every_divisor +** CarpetIOHDF5: teach the checkpointer to write index files +** CarpetIOHDF5: add support for index files to file reader +** LoopControl: Change printstats default to "no" +** CarpetInterp2: New ENO2 interpolator. +** Support more than one hierarchy of timers +** Add timers for refinement levels +** CarpetLib: Support higher-then 2nd order restriction operator for cell-centered meshes for Langrange prolongated grid functions. +* EinsteinExact: New arrangement, to eventually replace the Exact thorn. The thorns in this arrangement are generated by Kranc +* McLachlan +** Correct some errors in CCZ4, in particular the GammaShift parameter +** Change the default value of apply_dissipation to "never". +* ExternalLibraries: MPI is now treated like any other library on which flesh and thorns may depend, such as e.g. the GSL. It is thus configured via thorn ExternalLibraries/MPI, and not via the flesh any more. In particular, the option "MPI = ..." does not work any more and must not be used; instead, use e.g. "MPI_DIR = <path>" to point to an existing installation. You may also say "MPI_DIR = BUILD" to have Cactus build OpenMPI automatically, which is a convenient way of using MPI on a system where it may not be installed. +* LocalInterp: The interpolation loop uses OpenMP +* MoL: add multirate capabilities and three new multirate rk schemes +* SphericalSurface: make maxnphi and maxntheta parameters steerable +* NaNChecker: Add options to use restriction mask to ignore unused points, output more information about the failed points, use OpenMP during NaN check +* Multipole: Implement Driscoll&Healy integration +* HydroBase: schedule hydrobase_InitExcisionmask global-early loop-local +* TmunuBase: support_old_CalcTmunu_mechanism is now deprecated and is likely to be removed soon +* EOS_Omni: make more parameters steerable, add option read_table_on_single_process to have a named processor read the full table and MPI_Bcast it to the other processors +* GRHydro +** C2P: resort to bisection if Newton-Raphson fails. This is only activated if user requested this +** introduce parameters GRHydro_Y_e_min and GRHydro_Y_e_max to fix issues with EOS tables that don't extend sufficiently far +** Fixes to atmosphere mask and the way we sync +** various fixes to ePPM +** Option to use the slow sector of multirate RK methods +** fixed issues with hot EOS treatment +** Option to reconstruct W*vel instead of vel +** use an early atmosphere reset to ensure past timelevels are initialized correctly +** various MHD fixes +* Exact: openmp fixes +* GRHydro_InitData: new poloidal magnetic field routine, Added option to smooth metric inside horizon, Alfven wave test +* SummationByParts: improve OpenMP parallelization +* PITTNull*: respect Cactus' IO::truncate parameter +* AEIThorns/Trigger now part of the Einstein Toolkit + +External libraries: +* MPI: new - use this instead of the flesh mechanism (see above) +* BLAS: update to 3.4.2 +* FFTW3: update to 3.3.2 +* HDF5: will likely not build fortran and C++ api for next release by default (but still does now) +* LAPACK: update to 3.4.2 +* LORENE: don't use OpenMP unless it's enabled within Cactus +* OpenSSL: improve search for existing library +* zlib: update to 1.2.7 +* hwloc: new thorn for libhwloc + +=== How to upgrade from Lovelace (ET_2012_05) === + +To upgrade from the previous release, use GetComponents with the new thornlist to check out the new version. Due to the change of Carpet repository type, using GetComponents to upgrade a Lovelace checkout is unlikely to work. + +See the [http://einsteintoolkit.org/download/ Download] page on the Einstein Toolkit website for download instructions. + +=== Remaining issues with this release === + +* Certain machines need to be configured specially in Simfactory because the remote directories cannot be determined automatically just from the username. See the Machine notes below. + +* CarpetIOHDF5: When the new parameter CarpetIOHDF5::output_symmetry_points is set to "yes", then symmetry points are not actually output if CarpetIOHDF5::output_buffer_points has been set to "no". + +* Recovering with Carpet: The maximum number of timelevels that can be recovered is Carpet::prolongation_order_time+1. This is usually the case, but it is possible to write parameter files e.g. with prolongation_order_time=1 that use 3 timelevels. This bug manifests in an assert() failure when recovering from checkpoints. The work-around is to either increase Carpet::prolongation_order_time or decrease the number of timelevels for the grid function in question accordingly. + +=== Machine notes === + +==== Kraken ==== +defs.local.ini needs to have sourcebasedir = $HOME configured for this machine. You need to determine $HOME by logging in to the machine. + +==== LoneStar and Ranger ==== +defs.local.ini needs to have sourcebasedir = $WORK and basedir = $SCRATCH/simulations configured for this machine. You need to determine $WORK and $SCRATCH by logging in to the machine. + +==== Pandora ==== +The LSU system pandora has at the time of writing not only severe problems of compiling Cactus (can crash the filesystem), but in part due to this also outstanding are some failing testsuite results. Please contact Peter Diener at CCT (LSU) before trying to use this machine. + +All repositories participating in this release carry a branch ET_2012_11 +marking this release. These release branches will be updated if severe +errors are found. + +For more detailed information about the "Ãrsted" release please read +the long release announcement on the Einstein Toolkit web pages: +http://einsteintoolkit.org/about/releases/ET_2012_11_announcement.php. + +The "Ãrsted" Release Team on behalf of the Einstein Toolkit Consortium (2012-11-08) + +Eloisa Bentivegna +Tanja Bode +Peter Diener +Roland Haas +Ian Hinder +Frank LÃ¶ffler +Bruno Mundim +Christian D. Ott +Erik Schnetter + +November 8, 2012 +

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 908)
by knarf＠cct.lsu.edu 08 Nov '12

08 Nov '12

User: knarf Date: 2012/11/08 10:09 AM Modified: /about/members/ members.txt Log: Bruno moved File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +1 -1 =================================================================== --- about/members/members.txt 2012-11-08 16:07:57 UTC (rev 907) +++ about/members/members.txt 2012-11-08 16:09:41 UTC (rev 908) @@ -7,6 +7,7 @@ Seth Hopper Michael Koppitz Gian Mario Manca + Bruno Mundim Francesco Pannarale David Radice Luciano Rezzolla @@ -97,7 +98,6 @@ Manuela Campanelli Joshua Faber Carlos Lousto - Bruno Mundim Scott Noble Marcelo Ponce Billy Vazquez

1 0

← Newer
1
2
3
4
5
6
7
8
9
Older →

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits November 2012