Commits April 2013

commits@lists.einsteintoolkit.org

10 participants
99 discussions

[svn:einsteintoolkit] pyGWAnalysis/trunk/DataAnalysis/ (Rev. 17)
by reisswig＠tapir.caltech.edu 28 Apr '13

28 Apr '13

User: reisswig Date: 2013/04/28 11:56 AM Modified: /trunk/DataAnalysis/ DataAnalysis.py Log: Added Hubble constant. Added Luminosity distance caluclator. Added Big Bang Observer and DECIGO noise curves. File Changes: Directory: /trunk/DataAnalysis/ =============================== File [modified]: DataAnalysis.py Delta lines: +214 -1 =================================================================== --- trunk/DataAnalysis/DataAnalysis.py 2012-11-14 14:27:38 UTC (rev 16) +++ trunk/DataAnalysis/DataAnalysis.py 2013-04-28 16:56:56 UTC (rev 17) @@ -19,6 +19,7 @@ constant_c = 299792458 # speed of light in [m/sec] constant_G = 6.67428 * 1e-11 # gravitational constant in [m^3 kg^-1 s^-2] +constant_H0 = 70.0 # Hubble constant Msol = 1.9891*10**30 # solar mass[kg] ConversionFac = constant_G/constant_c**2 # [m/kg] (mass faktor G/c^2 to transform from geometrized units to SI units) @@ -379,7 +380,144 @@ return res +def LuminosityDistance(z, Omega_m=0.3175, Omega_lambda=0.6825): + """Computes the luminosity distance in Mpc for a given redshift z + Code from Ed Wright http://www.astro.ucla.edu/~wright/CosmoCalc.html""" + + H0 = constant_H0 + WM = Omega_m # Omega(matter) + WV = Omega_lambda # Omega(vacuum) or lambda + + WR = 0. # Omega(radiation) + WK = 0. # Omega curvaturve = 1-Omega(total) + c = 299792.458 # velocity of light in km/sec + Tyr = 977.8 # coefficent for converting 1/H into Gyr + DTT = 0.5 # time from z to now in units of 1/H0 + DTT_Gyr = 0.0 # value of DTT in Gyr + age = 0.5 # age of Universe in units of 1/H0 + age_Gyr = 0.0 # value of age in Gyr + zage = 0.1 # age of Universe at redshift z in units of 1/H0 + zage_Gyr = 0.0 # value of zage in Gyr + DCMR = 0.0 # comoving radial distance in units of c/H0 + DCMR_Mpc = 0.0 + DCMR_Gyr = 0.0 + DA = 0.0 # angular size distance + DA_Mpc = 0.0 + DA_Gyr = 0.0 + kpc_DA = 0.0 + DL = 0.0 # luminosity distance + DL_Mpc = 0.0 + DL_Gyr = 0.0 # DL in units of billions of light years + V_Gpc = 0.0 + a = 1.0 # 1/(1+z), the scale factor of the Universe + az = 0.5 # 1/(1+z(object)) + + h = H0/100. + WR = 4.165E-5/(h*h) # includes 3 massless neutrino species, T0 = 2.72528 + WK = 1-WM-WR-WV + az = 1.0/(1+1.0*z) + age = 0. + n=1000 # number of points in integrals + for i in range(n): + a = az*(i+0.5)/n + adot = sqrt(WK+(WM/a)+(WR/(a*a))+(WV*a*a)) + age = age + 1./adot + zage = az*age/n + zage_Gyr = (Tyr/H0)*zage + DTT = 0.0 + DCMR = 0.0 + + def DCMT(): + ratio = 1.00 + x = sqrt(abs(WK))*DCMR + if (x > 0.1): + if (WK > 0): ratio = 0.5*(exp(x)-exp(-x))/x + else: ratio = sin(x)/x; + y = ratio*DCMR + return y + y = x*x + # statement below fixed 13-Aug-03 to correct sign error in expansion + if (WK < 0): y = -y + ratio = 1 + y/6 + y*y/120 + y= ratio*DCMR + return y + + + def VCM(): + ratio = 1.00 + x = sqrt(abs(WK))*DCMR; + if (x > 0.1): + if (WK > 0): ratio = (0.125*(exp(2*x)-exp(-2*x))-x/2)/(x*x*x/3) + else: ratio = (x/2 - sin(2*x)/4)/(x*x*x/3) + y = ratio*DCMR*DCMR*DCMR/3 + return y + y = x*x + # statement below fixed 13-Aug-03 to correct sign error in expansion + if (WK < 0): y = -y + ratio = 1 + y/5 + (2/105)*y*y + y= ratio*DCMR*DCMR*DCMR/3 + return y + + + h = H0/100 + WR = 4.165E-5/(h*h) # includes 3 massless neutrino species, T0 = 2.72528 + WK = 1-WM-WR-WV + az = 1.0/(1+1.0*z) + age = 0; + for i in range (0, n+1): + a = az*(i+0.5)/n + adot = sqrt(WK+(WM/a)+(WR/(a*a))+(WV*a*a)) + age = age + 1/adot + zage = az*age/n + # correction for annihilations of particles not present now like e+/e- + # added 13-Aug-03 based on T_vs_t.f + lpz = log((1+1.0*z))/log(10.0) + dzage = 0 + if (lpz > 7.500): dzage = 0.002 * (lpz - 7.500) + if (lpz > 8.000): dzage = 0.014 * (lpz - 8.000) + 0.001 + if (lpz > 8.500): dzage = 0.040 * (lpz - 8.500) + 0.008 + if (lpz > 9.000): dzage = 0.020 * (lpz - 9.000) + 0.028 + if (lpz > 9.500): dzage = 0.019 * (lpz - 9.500) + 0.039 + if (lpz > 10.000): dzage = 0.048 + if (lpz > 10.775): dzage = 0.035 * (lpz - 10.775) + 0.048 + if (lpz > 11.851): dzage = 0.069 * (lpz - 11.851) + 0.086 + if (lpz > 12.258): dzage = 0.461 * (lpz - 12.258) + 0.114 + if (lpz > 12.382): dzage = 0.024 * (lpz - 12.382) + 0.171 + if (lpz > 13.055): dzage = 0.013 * (lpz - 13.055) + 0.188 + if (lpz > 14.081): dzage = 0.013 * (lpz - 14.081) + 0.201 + if (lpz > 15.107): dzage = 0.214 + zage = zage*10.0**dzage + + zage_Gyr = (Tyr/H0)*zage + DTT = 0.0 + DCMR = 0.0 + # do integral over a=1/(1+z) from az to 1 in n steps, midpoint rule + for i in range (0,n+1): + a = az+(1-az)*(i+0.5)/n + adot = sqrt(WK+(WM/a)+(WR/(a*a))+(WV*a*a)) + DTT = DTT + 1/adot + DCMR = DCMR + 1/(a*adot) + + DTT = (1-az)*DTT/n + DCMR = (1-az)*DCMR/n + age = DTT+zage + age_Gyr = age*(Tyr/H0) + DTT_Gyr = (Tyr/H0)*DTT + DCMR_Gyr = (Tyr/H0)*DCMR + DCMR_Mpc = (c/H0)*DCMR + DA = az*DCMT() + DA_Mpc = (c/H0)*DA + kpc_DA = DA_Mpc/206.264806 + DA_Gyr = (Tyr/H0)*DA + DL = DA/(az*az) + DL_Mpc = (c/H0)*DL + DL_Gyr = (Tyr/H0)*DL + V_Gpc = 4*pi*pow(0.001*c/H0,3)*VCM() + + return DL_Mpc + + def ShiftToZero(DF, trange): """Given a range [.,.] shift the wave according to that intervall to zero""" s = 0 @@ -427,7 +565,7 @@ res[-1].append([]) return res -def LoadModeArray(basename, lmax, col_time=1, col_re=2, col_im=3, filenameformat="cce_conversionscript"): +def LoadModeArray(basename, lmax, col_time=1, col_re=2, col_im=3, filenameformat="cce_conversion_script"): """Loads an array of modes up to lmax from disk. We assume that files begin with 'basename' and end with _l*m*.dat""" res = InitModeArray(lmax) @@ -551,7 +689,71 @@ return S +def SensitivityeLISA(f): + # Noise curve for eLISA-NGO. + # Function from eLISA yellow book, Sec. 4.3 + f0 = 150.0e-3 + a = 0.41 + T = sqrt(1.0+(f/(a*f0))**2) + + Da0 = 3.0 + fL = 0.1e-3 + fH = 8.0e-3 + da = Da0*1.0e-15 * sqrt(1.0+(f/fH)**4) * sqrt(1.0+fL/f) + + dx_DRS = 2.0*da / (2.0*pi*f)**2 + + Dx0 = 12.0 + f0 = 2.8e-3 + dx_IMS = Dx0*1e-12 * sqrt(1.0+(f0/f)**4) + + # Is this correct??? No... + S_IMS = dx_IMS + S_DRS = dx_DRS + + # arm length + L = 1e9 + + #h = sqrt(5.0)*2.0/sqrt(3.0)*T*sqrt(S_IMS + S_DRS) / L + # Use this in terms of dx_IMS, dx_DRS: + h = sqrt(5.0)*2.0/sqrt(3.0)*T*(dx_IMS + dx_DRS) / L + + S = h**2 + + return S + + +def SensitivityDECIGO(f): + """ From arXiv:1110.2865 + For this sensitivity curve, you don't need to select the + or x polarization component since the detector is a triangle. + The 60 degree abgle between detector arms is incoproated into the noise curve. + The antenna pattern is incorporated in the noise curve. + For SNRs, use an additional factor of sqrt(4) since BBO consists of 4 detectors!""" + + # White-dwarf confusion noise + SWD = 4.2e-47 * f**(-7.0/3.0) * exp(-2.0 * (f/5e-2)**2) + + fc = 7.69 + S = 3.30 * 1e-50*1.0/f**4 + 3.09*1e-47 * (1.0 + (f/fc)**2) + SWD + return S + + +def SensitivityBBO(f): + """ From arXiv:1110.2865 + For this sensitivity curve, you don't need to select the + or x polarization component since the detector is a triangle. + The 60 degree abgle between detector arms is incoproated into the noise curve. + The antenna pattern is incorporated in the noise curve. + For SNRs, use an additional factor of sqrt(4) since BBO consists of 4 detectors!""" + + # White-dwarf confusion noise + SWD = 4.2e-47 * f**(-7.0/3.0) * exp(-2.0 * (f/5e-2)**2) + + S = 6.15*1e-51*(1.0/f)**4 + 1.95*1e-48 + 1.20*1e-48*f**2 + SWD + return S + + + def SNR(fhSI, d, Sn, fLow, fHi): """Given a reconstructed h in SI units in Fourier domain and distance d (in Mpc) this computes the SNR for sensitivity curve Sn""" """Remember to select either h_+ or h_\cross! Otherise the SNR will be twice as large!""" @@ -561,6 +763,17 @@ return rho +def SNRavg(fhSIModes, lmax, d, Sn, fLow, fHi): + """Given a modes array in SI units in Fourier domain and distance d (in Mpc) this computes the angle averaged SNR for sensitivity curve Sn""" + S = map(Sn, fhSIModes[2][2].x) # create table of values for noise curve + + rho = 0.0 + for ll in range(2, lmax+1): + for mm in range(-ll, ll+1): + rho += 1.0/pi * WaveFunction(fhSIModes[ll][mm].x, abs(fhSIModes[ll][mm].f/(d*Mpc))**2 / S).Integrate(fLow, fHi, False).real + return sqrt(rho) + + def IntCharFreq(fhSI, d, Sn, fLow, fHi): """Given a reconstructed h in SI units in Fourier domain and distance d (in Mpc) this computes the integrated characteristic signal frequency for sensitivity curve Sn""" """Remember to select either h_+ or h_\cross! Otherise the SNR will be twice as large!"""

1 0

[Carpet] branch, master, updated. carpet-6-init-228-g8c63348
by Carpet git repositories 27 Apr '13

27 Apr '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "Carpet — Adaptive Mesh Refinement for the Cactus Framework <http://carpetcode.org/>". The branch, master has been updated via 8c63348a8570b61f802fd366812c74b6d5a6153e (commit) from 9b4e5841e2f30bded4cfed65e37c0ff2393250e8 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 8c63348a8570b61f802fd366812c74b6d5a6153e Author: Roland Haas <rhaas(a)tapir.caltech.edu> Date: Sat Apr 27 15:52:13 2013 -0700 Requirements: add missing include file ----------------------------------------------------------------------- Summary of changes: Carpet/Requirements/src/Requirements.cc | 1 + Carpet/Requirements/src/all_state.cc | 2 ++ Carpet/Requirements/src/gridpoint.cc | 2 ++ 3 files changed, 5 insertions(+), 0 deletions(-) hooks/post-receive -- Carpet — Adaptive Mesh Refinement for the Cactus Framework <http://carpetcode.org/>

1 0

[Carpet] branch, master, updated. carpet-6-init-227-g9b4e584
by Carpet git repositories 26 Apr '13

26 Apr '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "Carpet — Adaptive Mesh Refinement for the Cactus Framework <http://carpetcode.org/>". The branch, master has been updated via 9b4e5841e2f30bded4cfed65e37c0ff2393250e8 (commit) from fa53c27c1de5d1afd5bc9d904e8a1858b06bebc4 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 9b4e5841e2f30bded4cfed65e37c0ff2393250e8 Author: Erik Schnetter <schnetter(a)gmail.com> Date: Sat Apr 27 00:16:44 2013 -0400 CarpetLib: Correct alignment mechanism when allocating memory Handle the case where new returns memory with an alignment less than the type's size. ----------------------------------------------------------------------- Summary of changes: Carpet/CarpetLib/src/mem.cc | 4 +++- 1 files changed, 3 insertions(+), 1 deletions(-) hooks/post-receive -- Carpet — Adaptive Mesh Refinement for the Cactus Framework <http://carpetcode.org/>

1 0

[Carpet] branch, master, updated. carpet-6-init-226-gfa53c27
by Carpet git repositories 25 Apr '13

25 Apr '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "Carpet — Adaptive Mesh Refinement for the Cactus Framework <http://carpetcode.org/>". The branch, master has been updated via fa53c27c1de5d1afd5bc9d904e8a1858b06bebc4 (commit) from c35f19d61efd8d93b4d58882cf1ad7c8481c2f3a (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit fa53c27c1de5d1afd5bc9d904e8a1858b06bebc4 Author: Erik Schnetter <schnetter(a)gmail.com> Date: Thu Apr 25 14:12:14 2013 -0400 CarpetReduce: Offset reduction weights when periodic boundaries are used Offset reduction weights by 1/2 grid point when periodic boundaries are used with vertex centering. This means that all interior points have then weight 1, and all boundary points have then weight 0. ----------------------------------------------------------------------- Summary of changes: Carpet/CarpetReduce/src/mask_coords.c | 38 ++++++++++++++++++++++++++++---- 1 files changed, 33 insertions(+), 5 deletions(-) hooks/post-receive -- Carpet — Adaptive Mesh Refinement for the Cactus Framework <http://carpetcode.org/>

1 0

[McLachlan]McLachlan branch, master, updated. 966f1fdcffb6c327d4083a285edf895e3cd51da3
by Carpet git repositories 25 Apr '13

25 Apr '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "McLachlan". The branch, master has been updated via 966f1fdcffb6c327d4083a285edf895e3cd51da3 (commit) from d28dc270949233910c696d3860fbe3d0d6fe67de (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 966f1fdcffb6c327d4083a285edf895e3cd51da3 Author: Frank Löffler <knarf(a)cct.lsu.edu> Date: Thu Apr 25 10:05:23 2013 -0500 remove outdated thornlists (closes #1156) ----------------------------------------------------------------------- Summary of changes: doc/mclachlan-public.th | 136 ------------------------------------------- doc/qc0-mclachlan-public.th | 133 ------------------------------------------ 2 files changed, 0 insertions(+), 269 deletions(-) delete mode 100644 doc/mclachlan-public.th delete mode 100644 doc/qc0-mclachlan-public.th hooks/post-receive -- McLachlan

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 971)
by eloisa.bentivegna＠aei.mpg.de 25 Apr '13

25 Apr '13

User: bentivegna Date: 2013/04/25 08:55 AM Added: /release-info/ supermuc__1_8.log, supermuc__2_8.log Log: First test results for supermuc. File Changes: Directory: /release-info/ ========================= File [added]: supermuc__1_8.log Delta lines: <Unavailable> File [added]: supermuc__2_8.log Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 232)
by knarf＠cct.lsu.edu 23 Apr '13

23 Apr '13

User: knarf Date: 2013/04/23 08:27 PM Modified: /trunk/ einsteintoolkit.bib Log: add Moesta:2013dna File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: +14 -0 =================================================================== --- trunk/einsteintoolkit.bib 2013-04-08 22:47:47 UTC (rev 231) +++ trunk/einsteintoolkit.bib 2013-04-24 01:27:37 UTC (rev 232) @@ -2396,6 +2396,20 @@ primaryClass = "astro-ph.HE", SLACcitation = "%%CITATION = ARXIV:1001.2832;%%", } + +@article{Moesta:2013dna, + author = "Moesta, Philipp and Mundim, Bruno C. and Faber, Joshua A. + and Haas, Roland and Noble, Scott C. and Bode, Tanja and + L{\"o}ffler, Frank and Ott, Christian D. and Reisswig, Christian + and Schnetter, Erik", + title = "{{GRHydro}: {A} new open source general-relativistic + magnetohydrodynamics code for the {E}instein {T}oolkit}", + year = "2013", + eprint = "1304.5544", + archivePrefix = "arXiv", + primaryClass = "gr-qc", +} + @Article{Moncrief:1974am, author = "Moncrief, V.", title = "{Gravitational perturbations of spherically

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 970)
by knarf＠cct.lsu.edu 23 Apr '13

23 Apr '13

User: knarf Date: 2013/04/23 03:36 PM Modified: /about/members/ members.txt Log: remove Josh File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +0 -1 =================================================================== --- about/members/members.txt 2013-04-23 20:35:32 UTC (rev 969) +++ about/members/members.txt 2013-04-23 20:36:14 UTC (rev 970) @@ -62,7 +62,6 @@ http://www.cct.lsu.edu/~sbrandt/ Steven Brandt Dennis Castleberry Peter Diener - Joshua Harris Hari Krishnan Frank LÃ¶ffler Jian Tao

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 969)
by knarf＠cct.lsu.edu 23 Apr '13

23 Apr '13

User: knarf Date: 2013/04/23 03:35 PM Modified: /about/members/ members.txt Log: correct name File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +1 -1 =================================================================== --- about/members/members.txt 2013-04-23 20:32:41 UTC (rev 968) +++ about/members/members.txt 2013-04-23 20:35:32 UTC (rev 969) @@ -113,7 +113,7 @@ Scott Noble Marcelo Ponce Billy Vazquez - Miguel Zilhao + Miguel ZilhÃ£o Yosef Zlochower http://www.rutgers.edu/ Rutgers University

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 968)
by knarf＠cct.lsu.edu 23 Apr '13

23 Apr '13

User: knarf Date: 2013/04/23 03:32 PM Modified: /about/members/ members.txt Log: new member File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +3 -0 =================================================================== --- about/members/members.txt 2013-04-23 03:15:06 UTC (rev 967) +++ about/members/members.txt 2013-04-23 20:32:41 UTC (rev 968) @@ -75,6 +75,9 @@ Bernard Kelly Jennifer Seiler +https://www.camk.edu.pl/en/ Nicolaus Copernicus Astronomical Center (NCAC) + Bhupendra Prakash Mishra + http://www.ncsu.edu/ North Carolina State University Cody Simmons David Brown

1 0

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

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits April 2013