Commits March 2012

commits@lists.einsteintoolkit.org

14 participants
128 discussions

[svn:einsteintoolkit] www/workshop2012/ (Rev. 783)
by tanja.bode＠physics.gatech.edu 08 Mar '12

08 Mar '12

User: tbode Date: 2012/03/08 10:29 AM Modified: /workshop2012/ .htaccess Log: Remove hanging slash in wiki site. Was pointed to an empty page as is. File Changes: Directory: /workshop2012/ ========================= File [modified]: .htaccess Delta lines: +1 -1 =================================================================== --- workshop2012/.htaccess 2012-02-25 09:02:05 UTC (rev 782) +++ workshop2012/.htaccess 2012-03-08 16:29:05 UTC (rev 783) @@ -1 +1 @@ -Redirect /workshop2012/ https://docs.einsteintoolkit.org/et-docs/ET_Workshop_Spring_2012/ +Redirect /workshop2012/ https://docs.einsteintoolkit.org/et-docs/ET_Workshop_Spring_2012

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 159)
by roland.haas＠physics.gatech.edu 08 Mar '12

08 Mar '12

User: rhaas Date: 2012/03/08 09:13 AM Modified: /trunk/ einsteintoolkit.bib Log: fix bibtex syntax in Christodoulou:1970wf File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: +1 -1 =================================================================== --- trunk/einsteintoolkit.bib 2012-03-07 17:22:49 UTC (rev 158) +++ trunk/einsteintoolkit.bib 2012-03-08 15:13:51 UTC (rev 159) @@ -905,7 +905,7 @@ @article{Christodoulou:1970wf, author = "Christodoulou, D.", - title = "{Reversible and Irreversible Transformations in Black-Hole Physics}" + title = "{Reversible and Irreversible Transformations in Black-Hole Physics}" , journal = "Phys.Rev.Lett.", volume = "25", pages = "1596-1597",

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/flyer/ (Rev. 3)
by knarf＠cct.lsu.edu 08 Mar '12

08 Mar '12

User: knarf Date: 2012/03/07 11:10 AM Added: /flyer/ etk2.pdf, etk2.svg, etk3.png Log: ETK flyer File Changes: Directory: /flyer/ ================== File [added]: etk2.pdf Delta lines: <Unavailable> File [added]: etk2.svg Delta lines: <Unavailable> File [added]: etk3.png Delta lines: <Unavailable>

3 2

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Mar '12

07 Mar '12

changeset: 3490:c81ef26b02ec tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Wed Mar 07 19:25:56 2012 -0500 summary: CarpetIOF5: Make generic reading work

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/ (Rev. 57)
by cott＠tapir.caltech.edu 07 Mar '12

07 Mar '12

User: cott Date: 2012/03/07 06:19 PM Modified: /trunk/ interface.ccl /trunk/src/ EOS_Omni_MultiVarCalls.F90 Log: * add aliased routine to compute deps/drho|p and deps/dp|rho File Changes: Directory: /trunk/src/ ====================== File [modified]: EOS_Omni_MultiVarCalls.F90 Delta lines: +101 -0 =================================================================== --- trunk/src/EOS_Omni_MultiVarCalls.F90 2012-02-06 02:50:28 UTC (rev 56) +++ trunk/src/EOS_Omni_MultiVarCalls.F90 2012-03-08 00:19:42 UTC (rev 57) @@ -202,3 +202,104 @@ end select end subroutine EOS_Omni_EOS_dpderho_dpdrhoe + + + subroutine EOS_Omni_EOS_DEpsByDRho_DEpsByDPress(eoskey,keytemp,rf_precision,npoints,& + rho,eps,temp,ye,depsdrho,depsdpress,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: depsdrho(npoints) + CCTK_REAL, intent(out) :: depsdpress(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + real*8 :: hybrid_local_gamma + real*8 :: hybrid_local_k_cgs + real*8 :: hybrid_dp_poly,hybrid_dp_th1,hybrid_dp_th2 + ! temporary vars for nuc_eos + real*8 :: xrho,xye,xtemp,xenr,xent + real*8 :: xprs,xmunu,xcs2 + real*8 :: xdedt,xdpderho,xdpdrhoe + + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + depsdpress(i) = 1.0d0/(poly_gamma - 1.0d0)/rho(i) + depsdrho(i) = depsdpress(i) * press_gf * poly_k_cgs * & + poly_gamma * inv_rho_gf * & + (rho(i)*inv_rho_gf) ** (poly_gamma - 1.d0) + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + depsdpress(i) = 1.0/( (gl_gamma - 1.0d0) * & + rho(i)) + depsdrho(i) = -eps(i)/rho(i) + enddo + case (3) + ! hybrid EOS + do i=1,npoints + if(rho(i).gt.hybrid_rho_nuc) then + hybrid_local_gamma = hybrid_gamma2 + hybrid_local_k_cgs = hybrid_k2_cgs + else + hybrid_local_gamma = hybrid_gamma1 + hybrid_local_k_cgs = hybrid_k1_cgs + endif + hybrid_dp_poly = hybrid_local_gamma * press_gf * & + hybrid_local_k_cgs * rho(i)**(hybrid_local_gamma - 1.0d0) * & + inv_rho_gf**hybrid_local_gamma + + hybrid_dp_th1 = - hybrid_local_gamma * press_gf * hybrid_local_k_cgs * & + (hybrid_gamma_th - 1.d0) / (hybrid_local_gamma - 1.d0) * & + rho(i)**(hybrid_local_gamma - 1.d0) * inv_rho_gf**hybrid_local_gamma + + hybrid_dp_th2 = (hybrid_gamma_th - 1.d0) * eps(i) & + - (hybrid_gamma_th - 1.d0) * (hybrid_local_gamma - hybrid_gamma1) / & + (hybrid_gamma1 - 1.d0) / (hybrid_gamma2 - 1.d0) * & + press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & + hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) + + xdpdrhoe = hybrid_dp_poly + hybrid_dp_th1 + hybrid_dp_th2 + xdpderho = (hybrid_gamma_th - 1.0d0) * rho(i) + depsdpress(i) = 1.0 / xdpderho + depsdrho(i) = - xdpdrhoe * depsdpress(i) + enddo + case (4) + write(warnstring,*) "depsdrho and depsdpress not implemented yet for hot nuclear EOS" + call CCTK_WARN(0,warnstring) + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + + end subroutine EOS_Omni_EOS_DEpsByDRho_DEpsByDPress Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +15 -0 =================================================================== --- trunk/interface.ccl 2012-02-06 02:50:28 UTC (rev 56) +++ trunk/interface.ccl 2012-03-08 00:19:42 UTC (rev 57) @@ -138,8 +138,23 @@ PROVIDES FUNCTION EOS_Omni_short WITH EOS_Omni_EOS_short LANGUAGE Fortran +################################################################################ +# the following routines are needed for MHD con2prim using the 2D Z-P scheme +void FUNCTION EOS_Omni_DEpsByDRho_DEpsByDPress(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY DEpsByDRho, \ + CCTK_REAL OUT ARRAY DEpsByDPress, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) +PROVIDES FUNCTION EOS_Omni_DEpsByDRho_DEpsByDPress WITH EOS_Omni_EOS_DEpsByDRho_DEpsByDPress LANGUAGE Fortran ################################################################################ # the following routines are needed for MHD con2prim based on the Harm Code

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 270)
by bcmsma＠astro.rit.edu 07 Mar '12

07 Mar '12

User: bmundim Date: 2012/03/07 03:24 PM Modified: / ET.tex Log: Fix grammar error in section name: "BH binary" --> "BH binary system". There are a few instances in the text where the adjective "binary" was used as a noun. I will fix this soon in the sections I am responsible for and later at night on the whole text. Please check your section if this error occurs. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +1 -1 =================================================================== --- ET.tex 2012-03-07 21:14:00 UTC (rev 269) +++ ET.tex 2012-03-07 21:24:22 UTC (rev 270) @@ -2449,7 +2449,7 @@ only second-order accurate overall, this is the expected result. The increase of about 0.22\% in the mass of the BH is caused solely by the increase in the irreducible mass. -\subsection{BH Binary} +\subsection{BH Binary System} \label{sec:bbh-example} To demonstrate the performance in the code for a current problem of wide scientific interest, we have evolved a non-spinning equal-mass

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 269)
by bcmsma＠astro.rit.edu 07 Mar '12

07 Mar '12

User: bmundim Date: 2012/03/07 03:14 PM Modified: / ET.tex Log: Fixes related to the ID table. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +6 -5 =================================================================== --- ET.tex 2012-03-07 18:11:55 UTC (rev 268) +++ ET.tex 2012-03-07 21:14:00 UTC (rev 269) @@ -1023,7 +1023,7 @@ {\tt Lorene} is often used to generate initial data involving neutron stars, where the initial data problem is sufficiently hard that finding a numerical solution requires a significant amount of computer time as well as some human -vversight in order to guide the routine to the correct solution. +oversight in order to guide the routine to the correct solution. Therefore, recommended practice is to let Lorene output data into files, and then read those into the Einstein Toolkit at the beginning of a run. @@ -2478,10 +2478,11 @@ \hline QC3 & 3.0 & -3.0 & 0.0 & 0.13808 & 0.47656 & 0.984618 \\ \end{tabular}\\} -\caption{Initial data parameters for a non-spinning equal mass -BH binary. The punctures are located on the -$x$-axis at positions $x_1$ and $x_2$, with puncture bare -mass parameters $m_1 = m_2 = m$, and momenta $\pm\vec p$. +\caption{Initial data parameters and initial ADM mass +for a non-spinning equal mass BH binary system. The punctures are located +on the $x$-axis at positions $x_1$ and $x_2$, with puncture bare +mass parameters $m_1 = m_2 = m$, and momenta %$\pm\vec p$. +$\vec p_1 = -\vec p_2 = \vec p$. } \end{table}

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 268)
by diener＠cct.lsu.edu 07 Mar '12

07 Mar '12

User: diener Date: 2012/03/07 12:11 PM Modified: / ET.tex Log: Add a statement that Multipole can handle arbitrary spin-weights. And that s=-2 should be used for Psi_4. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +3 -0 =================================================================== --- ET.tex 2012-03-07 17:53:09 UTC (rev 267) +++ ET.tex 2012-03-07 18:11:55 UTC (rev 268) @@ -2039,6 +2039,9 @@ r^2 d\Omega \end{equation} of a projection onto spin-weighted spherical harmonics ${}_s Y_{\ell m}$. +Multipole provides spin-weighted spherical harmonics of arbitrary spin-weight +and can thus be used to decompose fields with different spin-weights. For +$\Psi_4$ the appropriate spin-weight is $s=-2$. \subsubsection{Object tracking} \label{sec:object-tracking}

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 267)
by diener＠cct.lsu.edu 07 Mar '12

07 Mar '12

User: diener Date: 2012/03/07 11:53 AM Modified: / ET.tex Log: Add a citation to Zerilli and explicitly name the Regge-Wheeler and Zerilli equations. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +5 -3 =================================================================== --- ET.tex 2012-03-07 17:25:31 UTC (rev 266) +++ ET.tex 2012-03-07 17:53:09 UTC (rev 267) @@ -1950,8 +1950,9 @@ on a Schwarzschild background or the calculation of the Weyl scalar $\Psi_4$. The module \codename{Extract} uses the Moncrief formalism~\cite{ -Moncrief:1974am} to extract gauge-invariant wave functions $Q_{\ell m}^\times$ and $Q_{\ell -m}^+$ given spherical surfaces of constant coordinate +Moncrief:1974am} to extract gauge-invariant wave functions $Q_{\ell m}^\times$ +(see~\cite{Regge:1957td}) and $Q_{\ell +m}^+$ (see~\cite{Zerilli:1970se}) given spherical surfaces of constant coordinate radius. The spatial metric is expressed as a perturbation on Schwarzschild and expanded into a tensor basis of the Regge-Wheeler harmonics~\cite{Regge:1957td} described by six standard @@ -1969,7 +1970,8 @@ + \Lambda r K^{+\ell m} \Bigg) \end{eqnarray} are calculated, where $S=1-2M/r$ and $\Lambda=(\ell-1)(\ell+2)+6M/r$. -These functions then satisfy the wave equations: +These functions then satisfy the Regge-Wheeler ($Q_{\ell m}^\times$) and +Zerilli ($Q_{\ell m}^+$) wave equations: \begin{eqnarray} (\partial_t^2-\partial_{r^*}^2)Q_{\ell m}^\times &=& - S \left[ \frac{\ell(\ell+1)}{r^2}-\frac{6M}{r^3} \right]

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 266)
by diener＠cct.lsu.edu 07 Mar '12

07 Mar '12

User: diener Date: 2012/03/07 11:25 AM Modified: / ET.tex Log: a citation to Steve's thesis where he shows in an appendix that a axi-symmetric spacetime cannot radiate angular momentum. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +1 -1 =================================================================== --- ET.tex 2012-03-07 09:25:16 UTC (rev 265) +++ ET.tex 2012-03-07 17:25:31 UTC (rev 266) @@ -2433,7 +2433,7 @@ \end{figure} In both cases the black (solid) curve is for high, blue (dashed) for medium and red (dotted) for low resolution. Since the spacetime is -axisymmetric the gravitational waves cannot radiate angular momentum. Thus +axisymmetric the gravitational waves cannot radiate angular momentum~\cite{Brandt:1996th}. Thus any change in the spin must be due to numerical error and $\Delta S$ should converge to zero with increasing resolution. This is clearly shown in the bottom plot of Figure~\ref{fig:ah_mass_spin}; the green

1 0

← Newer
1
...
7
8
9
10
11
12
13
Older →

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits March 2012