Commits September 2011

commits@lists.einsteintoolkit.org

6 participants
58 discussions

[svn:einsteintoolkit] manifest/trunk/ (Rev. 112)
by knarf＠cct.lsu.edu 20 Sep '11

20 Sep '11

User: knarf Date: 2011/09/20 12:29 PM Modified: /trunk/ einsteintoolkit.th Log: add FFTW3 File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.th Delta lines: +1 -0 =================================================================== --- trunk/einsteintoolkit.th 2011-09-19 17:11:22 UTC (rev 111) +++ trunk/einsteintoolkit.th 2011-09-20 17:29:26 UTC (rev 112) @@ -223,6 +223,7 @@ !CHECKOUT = ExternalLibraries/BLAS ExternalLibraries/GSL +ExternalLibraries/FFTW3 ExternalLibraries/HDF5 ExternalLibraries/LAPACK ExternalLibraries/libjpeg

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 277)
by roland.haas＠physics.gatech.edu 19 Sep '11

19 Sep '11

User: rhaas Date: 2011/09/19 06:16 PM Modified: /trunk/src/ GRHydro_Source.F90 Log: fix non-multipatch source terms File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Source.F90 Delta lines: +26 -11 =================================================================== --- trunk/src/GRHydro_Source.F90 2011-09-19 23:16:09 UTC (rev 276) +++ trunk/src/GRHydro_Source.F90 2011-09-19 23:16:38 UTC (rev 277) @@ -80,10 +80,13 @@ CCTK_INT :: GRHydro_UseGeneralCoordinates TARGET gaa, gab, gac, gbb, gbc, gcc TARGET gxx, gxy, gxz, gyy, gyz, gzz + TARGET kaa, kab, kac, kbb, kbc, kcc + TARGET kxx, kxy, kxz, kyy, kyz, kzz TARGET betaa, betab, betac TARGET betax, betay, betaz TARGET lvel, vel CCTK_REAL, DIMENSION(:,:,:), POINTER :: g11, g12, g13, g22, g23, g33 + CCTK_REAL, DIMENSION(:,:,:), POINTER :: k11, k12, k13, k22, k23, k33 CCTK_REAL, DIMENSION(:,:,:), POINTER :: beta1, beta2, beta3 CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup @@ -94,6 +97,12 @@ g22 => gbb g23 => gbc g33 => gcc + k11 => kaa + k12 => kab + k13 => kac + k22 => kbb + k23 => kbc + k33 => kcc beta1 => betaa beta2 => betab beta3 => betac @@ -105,6 +114,12 @@ g22 => gyy g23 => gyz g33 => gzz + k11 => kxx + k12 => kxy + k13 => kxz + k22 => kyy + k23 => kyz + k33 => kzz beta1 => betax beta2 => betay beta3 => betaz @@ -340,21 +355,21 @@ !!$ Contract the shift with the eatrinsic curvature - shiftshiftk = shifta*shifta*kaa(i,j,k) + & - shiftb*shiftb*kbb(i,j,k) + & - shiftc*shiftc*kcc(i,j,k) + & - two*(shifta*shiftb*kab(i,j,k) + & - shifta*shiftc*kac(i,j,k) + & - shiftb*shiftc*kbc(i,j,k)) + shiftshiftk = shifta*shifta*k11(i,j,k) + & + shiftb*shiftb*k22(i,j,k) + & + shiftc*shiftc*k33(i,j,k) + & + two*(shifta*shiftb*k12(i,j,k) + & + shifta*shiftc*k13(i,j,k) + & + shiftb*shiftc*k23(i,j,k)) - shiftka = shifta*kaa(i,j,k) + shiftb*kab(i,j,k) + shiftc*kac(i,j,k) - shiftkb = shifta*kab(i,j,k) + shiftb*kbb(i,j,k) + shiftc*kbc(i,j,k) - shiftkc = shifta*kac(i,j,k) + shiftb*kbc(i,j,k) + shiftc*kcc(i,j,k) + shiftka = shifta*k11(i,j,k) + shiftb*k12(i,j,k) + shiftc*k13(i,j,k) + shiftkb = shifta*k12(i,j,k) + shiftb*k22(i,j,k) + shiftc*k23(i,j,k) + shiftkc = shifta*k13(i,j,k) + shiftb*k23(i,j,k) + shiftc*k33(i,j,k) !!$ Contract the matter terms with the extrinsic curvature - sumTK = taa*kaa(i,j,k) + tbb*kbb(i,j,k) + tcc*kcc(i,j,k) & - + two*(tab*kab(i,j,k) + tac*kac(i,j,k) + tbc*kbc(i,j,k)) + sumTK = taa*k11(i,j,k) + tbb*k22(i,j,k) + tcc*k33(i,j,k) & + + two*(tab*k12(i,j,k) + tac*k13(i,j,k) + tbc*k23(i,j,k)) !!$ Update term for tau

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 276)
by roland.haas＠physics.gatech.edu 19 Sep '11

19 Sep '11

User: rhaas Date: 2011/09/19 06:16 PM Modified: /trunk/src/ GRHydro_ENOReconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_TVDReconstruct_drv.F90 Log: fix non-multipatch support for MHD File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: +7 -4 =================================================================== --- trunk/src/GRHydro_ENOReconstruct_drv.F90 2011-09-19 23:15:08 UTC (rev 275) +++ trunk/src/GRHydro_ENOReconstruct_drv.F90 2011-09-19 23:16:09 UTC (rev 276) @@ -20,9 +20,9 @@ #define sx(i,j,k) scon(i,j,k,1) #define sy(i,j,k) scon(i,j,k,2) #define sz(i,j,k) scon(i,j,k,3) -#define Bvecx(i,j,k) lBvec(i,j,k,1) -#define Bvecy(i,j,k) lBvec(i,j,k,2) -#define Bvecz(i,j,k) lBvec(i,j,k,3) +#define Bvecx(i,j,k) Bprim(i,j,k,1) +#define Bvecy(i,j,k) Bprim(i,j,k,2) +#define Bvecz(i,j,k) Bprim(i,j,k,3) #define Bconsx(i,j,k) Bcons(i,j,k,1) #define Bconsy(i,j,k) Bcons(i,j,k,2) #define Bconsz(i,j,k) Bcons(i,j,k,3) @@ -67,12 +67,15 @@ ! save memory when MP is not used CCTK_INT :: GRHydro_UseGeneralCoordinates TARGET vel, lvel - CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + TARGET Bvec, lBvec + CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup, Bprim if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then vup => lvel + Bprim => lBvec else vup => vel + Bprim => Bvec end if allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: +7 -4 =================================================================== --- trunk/src/GRHydro_PPMMReconstruct_drv.F90 2011-09-19 23:15:08 UTC (rev 275) +++ trunk/src/GRHydro_PPMMReconstruct_drv.F90 2011-09-19 23:16:09 UTC (rev 276) @@ -20,9 +20,9 @@ #define sx(i,j,k) scon(i,j,k,1) #define sy(i,j,k) scon(i,j,k,2) #define sz(i,j,k) scon(i,j,k,3) -#define Bvecx(i,j,k) lBvec(i,j,k,1) -#define Bvecy(i,j,k) lBvec(i,j,k,2) -#define Bvecz(i,j,k) lBvec(i,j,k,3) +#define Bvecx(i,j,k) Bprim(i,j,k,1) +#define Bvecy(i,j,k) Bprim(i,j,k,2) +#define Bvecz(i,j,k) Bprim(i,j,k,3) #define Bconsx(i,j,k) Bcons(i,j,k,1) #define Bconsy(i,j,k) Bcons(i,j,k,2) #define Bconsz(i,j,k) Bcons(i,j,k,3) @@ -71,9 +71,10 @@ TARGET betaa, betab, betac TARGET betax, betay, betaz TARGET lvel, vel + TARGET lBvec, Bvec CCTK_REAL, DIMENSION(:,:,:), POINTER :: g11, g12, g13, g22, g23, g33 CCTK_REAL, DIMENSION(:,:,:), POINTER :: beta1, beta2, beta3 - CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup, Bprim if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then g11 => gaa @@ -86,6 +87,7 @@ beta2 => betab beta3 => betac vup => lvel + Bprim => lBvec else g11 => gxx g12 => gxy @@ -97,6 +99,7 @@ beta2 => betay beta3 => betaz vup => vel + Bprim => Bvec end if allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: +10 -7 =================================================================== --- trunk/src/GRHydro_TVDReconstruct_drv.F90 2011-09-19 23:15:08 UTC (rev 275) +++ trunk/src/GRHydro_TVDReconstruct_drv.F90 2011-09-19 23:16:09 UTC (rev 276) @@ -20,9 +20,9 @@ #define sx(i,j,k) scon(i,j,k,1) #define sy(i,j,k) scon(i,j,k,2) #define sz(i,j,k) scon(i,j,k,3) -#define Bvecx(i,j,k) lBvec(i,j,k,1) -#define Bvecy(i,j,k) lBvec(i,j,k,2) -#define Bvecz(i,j,k) lBvec(i,j,k,3) +#define Bvecx(i,j,k) Bprim(i,j,k,1) +#define Bvecy(i,j,k) Bprim(i,j,k,2) +#define Bvecz(i,j,k) Bprim(i,j,k,3) #define Bconsx(i,j,k) Bcons(i,j,k,1) #define Bconsy(i,j,k) Bcons(i,j,k,2) #define Bconsz(i,j,k) Bcons(i,j,k,3) @@ -65,12 +65,15 @@ ! save memory when MP is not used CCTK_INT :: GRHydro_UseGeneralCoordinates TARGET vel, lvel - CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + TARGET Bvec, lBvec + CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup, Bprim if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then vup => lvel + Bprim => lBvec else vup => vel + Bprim => Bvec end if allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) @@ -175,11 +178,11 @@ eps, epsplus, epsminus, trivial_rp, hydro_excision_mask) if(evolve_mhd.ne.0) then call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & - lBvec(:,:,:,1), Bvecxplus, Bvecxminus, trivial_rp, hydro_excision_mask) + Bprim(:,:,:,1), Bvecxplus, Bvecxminus, trivial_rp, hydro_excision_mask) call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & - lBvec(:,:,:,2), Bvecyplus, Bvecyminus, trivial_rp, hydro_excision_mask) + Bprim(:,:,:,2), Bvecyplus, Bvecyminus, trivial_rp, hydro_excision_mask) call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & - lBvec(:,:,:,3), Bveczplus, Bveczminus, trivial_rp, hydro_excision_mask) + Bprim(:,:,:,3), Bveczplus, Bveczminus, trivial_rp, hydro_excision_mask) endif else if (CCTK_EQUALS(recon_vars,"conservative")) then

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 275)
by roland.haas＠physics.gatech.edu 19 Sep '11

19 Sep '11

User: rhaas Date: 2011/09/19 06:15 PM Modified: /trunk/src/ GRHydro_ENOReconstruct_drv.F90 Log: fix non-multipatch velocities in ENO reconstruction File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: +14 -3 =================================================================== --- trunk/src/GRHydro_ENOReconstruct_drv.F90 2011-09-15 16:51:30 UTC (rev 274) +++ trunk/src/GRHydro_ENOReconstruct_drv.F90 2011-09-19 23:15:08 UTC (rev 275) @@ -14,9 +14,9 @@ #include "SpaceMask.h" -#define velx(i,j,k) lvel(i,j,k,1) -#define vely(i,j,k) lvel(i,j,k,2) -#define velz(i,j,k) lvel(i,j,k,3) +#define velx(i,j,k) vup(i,j,k,1) +#define vely(i,j,k) vup(i,j,k,2) +#define velz(i,j,k) vup(i,j,k,3) #define sx(i,j,k) scon(i,j,k,1) #define sy(i,j,k) scon(i,j,k,2) #define sz(i,j,k) scon(i,j,k,3) @@ -64,6 +64,17 @@ CCTK_INT :: ierr + ! save memory when MP is not used + CCTK_INT :: GRHydro_UseGeneralCoordinates + TARGET vel, lvel + CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + + if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then + vup => lvel + else + vup => vel + end if + allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp")

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 111)
by knarf＠cct.lsu.edu 19 Sep '11

19 Sep '11

User: knarf Date: 2011/09/19 12:11 PM Removed: /trunk/ CCE.th Modified: /trunk/ einsteintoolkit.th Log: add CCE to thornlist File Changes: Directory: /trunk/ ================== File [removed]: CCE.th Delta lines: +0 -22 =================================================================== --- trunk/CCE.th 2011-09-15 17:19:41 UTC (rev 110) +++ trunk/CCE.th 2011-09-19 17:11:22 UTC (rev 111) @@ -1,22 +0,0 @@ -!CRL_VERSION = 1.0 - -!INCLUDE = http://svn.einsteintoolkit.org/manifest/trunk/einsteintoolkit.th - -# CCE -!TARGET = $ARR -!TYPE = svn -!AUTH_URL = https://svn.einsteintoolkit.org/cactus/$1/$2/trunk -!URL = http://svn.einsteintoolkit.org/cactus/$1/$2/trunk -!CHECKOUT = -PITTNullCode/NullConstr -PITTNullCode/NullDecomp -PITTNullCode/NullEvolve -PITTNullCode/NullExact -PITTNullCode/NullGrid -PITTNullCode/NullInterp -PITTNullCode/NullNews -PITTNullCode/NullPsiInt -PITTNullCode/NullSHRExtract -PITTNullCode/NullVars -PITTNullCode/SphericalHarmonicRecon - File [modified]: einsteintoolkit.th Delta lines: +20 -0 =================================================================== --- trunk/einsteintoolkit.th 2011-09-15 17:19:41 UTC (rev 110) +++ trunk/einsteintoolkit.th 2011-09-19 17:11:22 UTC (rev 111) @@ -317,6 +317,26 @@ +# CCE +!TARGET = $ARR +!TYPE = svn +!AUTH_URL = https://svn.einsteintoolkit.org/cactus/$1/$2/trunk +!URL = http://svn.einsteintoolkit.org/cactus/$1/$2/trunk +!CHECKOUT = +PITTNullCode/NullConstr +PITTNullCode/NullDecomp +PITTNullCode/NullEvolve +PITTNullCode/NullExact +PITTNullCode/NullGrid +PITTNullCode/NullInterp +PITTNullCode/NullNews +PITTNullCode/NullPsiInt +PITTNullCode/NullSHRExtract +PITTNullCode/NullVars +PITTNullCode/SphericalHarmonicRecon + + + # Private thorns !TARGET = $ARR !TYPE = ignore

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/bbh/ (Rev. 134)
by bcmsma＠astro.rit.edu 16 Sep '11

16 Sep '11

User: bmundim Date: 2011/09/16 05:43 PM Added: /examples/bbh/BBHHighRes/ amp_phs_frq_mp_psi4_l2_m2_r30.00.asc, mp_psi4_l2_m2_r30.00.asc /examples/bbh/BBHHigherRes/ amp_phs_frq_mp_psi4_l2_m2_r30.00.asc, amp_phs_frq_mp_psi4_l2_m2_r60.00.asc, hi4th_amp_phs_mp_psi4_l2_m2_r60.00.asc, hi6th_amp_phs_mp_psi4_l2_m2_r60.00.asc, hi8th_amp_phs_mp_psi4_l2_m2_r60.00.asc, ml_admconstraints::ml_ham.x.asc, mp_psi4_l2_m2_r30.00.asc, mp_psi4_l2_m2_r60.00.asc, puncturetracker::pt_loc..asc /examples/bbh/BBHLowRes/ amp_phs_frq_mp_psi4_l2_m2_r30.00.asc, mp_psi4_l2_m2_r30.00.asc /examples/bbh/BBHMedHRes/ amp_phs_frq_mp_psi4_l2_m2_r30.00.asc, mp_psi4_l2_m2_r30.00.asc /examples/bbh/BBHMedRes/ amp_phs_frq_mp_psi4_l2_m2_r30.00.asc, mp_psi4_l2_m2_r30.00.asc /examples/bbh/figs/ amp_convergence_mehhihir_4th.eps, amp_convergence_mehhihir_6th.eps, amp_convergence_mehhihir_8th.eps, phase_convergence_mehhihir_4th.eps, phase_convergence_mehhihir_6th.eps, phase_convergence_mehhihir_8th.eps /examples/bbh/sm_scripts/ sm_conv_amp_mehhihir_4th, sm_conv_amp_mehhihir_6th, sm_conv_amp_mehhihir_8th, sm_conv_phs_mehhihir_4th, sm_conv_phs_mehhihir_6th, sm_conv_phs_mehhihir_8th Modified: /examples/bbh/sm_scripts/ .smhist Log: Data, scripts and plots for bbh example taking into account a higher resolution run. Comments: 1) No improvement at all on determining a convergence rate close to 4th order. Indeed data gets noisier. 2) Add waveform data at r=30M. Not that helpful either. 3) I will focus on showing convergence and not convergence rate. File Changes: Directory: /examples/bbh/BBHHighRes/ ==================================== File [added]: amp_phs_frq_mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHHigherRes/ ====================================== File [added]: amp_phs_frq_mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> File [added]: amp_phs_frq_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: hi4th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: hi6th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: hi8th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: ml_admconstraints::ml_ham.x.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: puncturetracker::pt_loc..asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHLowRes/ =================================== File [added]: amp_phs_frq_mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHMedHRes/ ==================================== File [added]: amp_phs_frq_mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHMedRes/ =================================== File [added]: amp_phs_frq_mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m2_r30.00.asc Delta lines: <Unavailable> Directory: /examples/bbh/figs/ ============================== File [added]: amp_convergence_mehhihir_4th.eps Delta lines: <Unavailable> File [added]: amp_convergence_mehhihir_6th.eps Delta lines: <Unavailable> File [added]: amp_convergence_mehhihir_8th.eps Delta lines: <Unavailable> File [added]: phase_convergence_mehhihir_4th.eps Delta lines: <Unavailable> File [added]: phase_convergence_mehhihir_6th.eps Delta lines: <Unavailable> File [added]: phase_convergence_mehhihir_8th.eps Delta lines: <Unavailable> Directory: /examples/bbh/sm_scripts/ ==================================== File [modified]: .smhist Delta lines: <Unavailable> File [added]: sm_conv_amp_mehhihir_4th Delta lines: <Unavailable> File [added]: sm_conv_amp_mehhihir_6th Delta lines: <Unavailable> File [added]: sm_conv_amp_mehhihir_8th Delta lines: <Unavailable> File [added]: sm_conv_phs_mehhihir_4th Delta lines: <Unavailable> File [added]: sm_conv_phs_mehhihir_6th Delta lines: <Unavailable> File [added]: sm_conv_phs_mehhihir_8th Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] HydroBase/trunk/doc/ (Rev. 42)
by bcmsma＠astro.rit.edu 15 Sep '11

15 Sep '11

User: bmundim Date: 2011/09/15 04:54 PM Modified: /trunk/doc/ documentation.tex Log: Better form. File Changes: Directory: /trunk/doc/ ====================== File [modified]: documentation.tex Delta lines: +7 -5 =================================================================== --- trunk/doc/documentation.tex 2011-09-15 20:48:00 UTC (rev 41) +++ trunk/doc/documentation.tex 2011-09-15 21:54:40 UTC (rev 42) @@ -215,11 +215,13 @@ Inserting the SI units into the above unit correspondences, we find the following conversion factors: -\begin{equation} - M_\mathrm{sun} \equiv 1.477\, \mathrm{km} - \equiv 4.92673\, \mathrm{\mu s} - \equiv 8.344956\cdot 10^{15}\, \mathrm{T}, -\end{equation} +\begin{eqnarray} + [L] & = & 1\, M_\mathrm{sun} \equiv 1.477\, \mathrm{km} + \\{} + [T] & = & 1\, M_\mathrm{sun} \equiv 4.92673\, \mathrm{\mu s} + \\{} + [B] & = & 1\, \sqrt{\mu_0} /M_\mathrm{sun} \equiv 8.344956\cdot 10^{15}\, \mathrm{T}, +\end{eqnarray} where T (Tesla) is the magnetic field unit in SI, $1\,\mathrm{T}=1\,\mathrm{N/(A\cdot m)}$, or in its cgs equivalent, $1\,\mathrm{Tesla} = 10^4\,\mathrm{Gauss}$.

1 0

[svn:einsteintoolkit] HydroBase/trunk/doc/ (Rev. 41)
by bcmsma＠astro.rit.edu 15 Sep '11

15 Sep '11

User: bmundim Date: 2011/09/15 03:48 PM Modified: /trunk/doc/ documentation.tex Log: Add a few lines about units adopted, in particular the conversion factors for your convenience. Please review these numbers though. File Changes: Directory: /trunk/doc/ ====================== File [modified]: documentation.tex Delta lines: +15 -1 =================================================================== --- trunk/doc/documentation.tex 2011-04-21 04:46:32 UTC (rev 40) +++ trunk/doc/documentation.tex 2011-09-15 20:48:00 UTC (rev 41) @@ -132,7 +132,7 @@ a common set of variables by different hydrodynamics codes creates the possibility to share parts of the code, e.g.\ initial data solvers or analysis routines. -Currently the defined primitive variables are (see~\cite{livrevgrfd} for +Currently the defined primitive variables are (see~\cite{livrevgrrfd} for details): \begin{itemize} \item \verb|rho|: rest mass density $\varrho$ @@ -185,6 +185,8 @@ M_\mathrm{sun} = 1 &;& c = G = 1 \end{eqnarray} which are commonly used in astrophysics and in relativity. +The former sets the mass scale to the solar one and the latter +adopts the same units for time, length and mass. We assume the following definitions and constants of nature: \begin{eqnarray} @@ -211,6 +213,18 @@ [B] & = & 1/[L]\; / \sqrt{\epsilon_0 G / c^2} \end{eqnarray} +Inserting the SI units into the above unit correspondences, +we find the following conversion factors: +\begin{equation} + M_\mathrm{sun} \equiv 1.477\, \mathrm{km} + \equiv 4.92673\, \mathrm{\mu s} + \equiv 8.344956\cdot 10^{15}\, \mathrm{T}, +\end{equation} +where T (Tesla) is the magnetic field unit in SI, +$1\,\mathrm{T}=1\,\mathrm{N/(A\cdot m)}$, or in its cgs equivalent, +$1\,\mathrm{Tesla} = 10^4\,\mathrm{Gauss}$. + + \section{Acknowledgments} This thorn was produced by Tanja Bode, Roland Haas, Frank L\"offler, and Erik Schnetter.

1 0

[svn:einsteintoolkit] manifest/branches/ET_2011_05/ (Rev. 110)
by knarf＠cct.lsu.edu 15 Sep '11

15 Sep '11

User: knarf Date: 2011/09/15 12:19 PM Modified: /branches/ET_2011_05/ einsteintoolkit.th Log: reflect move of two Cactus thorns File Changes: Directory: /branches/ET_2011_05/ ================================ File [modified]: einsteintoolkit.th Delta lines: +2 -2 =================================================================== --- branches/ET_2011_05/einsteintoolkit.th 2011-09-06 21:01:51 UTC (rev 109) +++ branches/ET_2011_05/einsteintoolkit.th 2011-09-15 17:19:41 UTC (rev 110) @@ -73,8 +73,6 @@ CactusBase/IOASCII CactusBase/IOBasic CactusBase/IOUtil -CactusBase/LocalInterp -CactusBase/LocalReduce CactusBase/SymBase CactusBase/Time @@ -91,6 +89,8 @@ CactusNumerical/Cartoon2D CactusNumerical/Dissipation CactusNumerical/InterpToArray +CactusNumerical/LocalInterp +CactusNumerical/LocalReduce CactusNumerical/MoL CactusNumerical/Noise CactusNumerical/Norms

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 274)
by roland.haas＠physics.gatech.edu 15 Sep '11

15 Sep '11

User: rhaas Date: 2011/09/15 11:51 AM Added: /trunk/src/ GRHydro_TransformTensorBasis.c Removed: /trunk/src/ GRHydro_TransformTensorBasis.F90 Modified: /trunk/ interface.ccl, schedule.ccl /trunk/src/ GRHydro_Boundaries.F90, GRHydro_Con2Prim.F90, GRHydro_EoSChangeGamma.F90, GRHydro_HLLE.F90, GRHydro_Jacobian_state.c, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_ParamCheck.F90, GRHydro_Prim2Con.F90, GRHydro_RegisterVars.cc, GRHydro_Source.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_UpdateMask.F90, make.code.defn Log: remove explicit dependence on Multipatch thorns GRHydro now compiles and runs without MP thorns present and behaves identical to the old version File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Boundaries.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [modified]: GRHydro_EoSChangeGamma.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLE.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Jacobian_state.c Delta lines: <Unavailable> File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPMReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ParamCheck.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2Con.F90 Delta lines: <Unavailable> File [modified]: GRHydro_RegisterVars.cc Delta lines: <Unavailable> File [modified]: GRHydro_Source.F90 Delta lines: <Unavailable> File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: <Unavailable> File [removed]: GRHydro_TransformTensorBasis.F90 Delta lines: <Unavailable> File [added]: GRHydro_TransformTensorBasis.c Delta lines: <Unavailable> File [modified]: GRHydro_UpdateMask.F90 Delta lines: <Unavailable> File [modified]: make.code.defn Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: <Unavailable> File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits September 2011