Commits August 2011

commits@lists.einsteintoolkit.org

9 participants
99 discussions

[svn:einsteintoolkit] manifest/trunk/ (Rev. 107)
by knarf＠cct.lsu.edu 21 Aug '11

21 Aug '11

User: knarf Date: 2011/08/22 01:37 AM Added: /trunk/ CCE.th Log: CCE thornlist - distinct from einsteintoolkit.th for now (couldn't test it, some servers are down atm. However, this is using einsteintoolkit.th as include, so it's almost as good File Changes: Directory: /trunk/ ================== File [added]: CCE.th Delta lines: +22 -0 =================================================================== --- trunk/CCE.th (rev 0) +++ trunk/CCE.th 2011-08-22 06:37:42 UTC (rev 107) @@ -0,0 +1,22 @@ +!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 +

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 115)
by knarf＠cct.lsu.edu 21 Aug '11

21 Aug '11

User: knarf Date: 2011/08/22 12:47 AM Added: /examples/tov/ make_plots.sh, mode_labels.sh, mode_spectrum.pdf, mode_spectrum.plot, mode_spectrum.svg, modes.dat /examples/tov/tov_2/ hydrobase::rho.maximum.asc.bz2, hydrobase::rho.maximum.asc.spectrum.bz2 Modified: / ET.tex Log: add preliminary (still running) results for tov mode spectrum Directory Changes: Directory: /svn:executable/ =========================== + * Directory: /svn:mime-type/ ========================== + application/octet-stream File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +16 -0 =================================================================== --- ET.tex 2011-08-22 05:29:40 UTC (rev 114) +++ ET.tex 2011-08-22 05:47:08 UTC (rev 115) @@ -2446,8 +2446,24 @@ \end{figure} \subsection{TOV\pages{2 Frank}} +\todo{TOV bla} +In figure~\ref{fig:tov_mode_spectrum} we show the power spectrum of the central +density oscillations computed from a full 3D relativistic hydrodynamics +simulation using GRHydro compared to the cooresponding frequencies as obtained +with perturbative techniques\todo{thanks}. Clearly the 3D simulation can +correctly represent the fundamental mode frequency and its overtones of the +neutron star. +\begin{figure} + \label{fig:tov_mode_spectrum} + \includegraphics[width=0.8\textwidth]{examples/tov/mode_spectrum} + \caption{Eigenfrequency mode spectrum of TOV star. Shown is the spectrum + computed from a full 3D relativistic hydrodynamics simulation compared + to the values obtained by perturbation theory.} +\end{figure} + + \subsection{Collapse\pages{2 Christian}} Show TOV collapse and BH formation Directory: /examples/tov/ ========================= File [added]: make_plots.sh Delta lines: +2 -0 =================================================================== --- examples/tov/make_plots.sh (rev 0) +++ examples/tov/make_plots.sh 2011-08-22 05:47:08 UTC (rev 115) @@ -0,0 +1,2 @@ +gnuplot < mode_spectrum.plot +inkscape -A mode_spectrum.pdf mode_spectrum.svg Property changes on: examples/tov/make_plots.sh ___________________________________________________________________ File [added]: mode_labels.sh Delta lines: +1 -0 =================================================================== --- examples/tov/mode_labels.sh (rev 0) +++ examples/tov/mode_labels.sh 2011-08-22 05:47:08 UTC (rev 115) @@ -0,0 +1 @@ +perl -ane 'if ($F[2]) { print "set label \"$F[2]\" at ".($F[0]-0.15).",2.e-6\n"; }' modes.dat Property changes on: examples/tov/mode_labels.sh ___________________________________________________________________ File [added]: mode_spectrum.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: examples/tov/mode_spectrum.pdf ___________________________________________________________________ File [added]: mode_spectrum.plot Delta lines: +10 -0 =================================================================== --- examples/tov/mode_spectrum.plot (rev 0) +++ examples/tov/mode_spectrum.plot 2011-08-22 05:47:08 UTC (rev 115) @@ -0,0 +1,10 @@ +set size ratio 0.3 +set term svg enhanced size 600,250 +set output 'mode_spectrum.svg' +set logscale y +set xlabel 'f/[kHz]' +set xtics 2 +set mxtics 2 +unset ytics +load '<./mode_labels.sh' +plot [:10][1.e-10:1.e-5] '<bzcat tov_2/hydrobase::rho.maximum.asc.spectrum.bz2' u ($1/1000):2 w lp t '', 'modes.dat' w l t '' File [added]: mode_spectrum.svg Delta lines: +147 -0 =================================================================== --- examples/tov/mode_spectrum.svg (rev 0) +++ examples/tov/mode_spectrum.svg 2011-08-22 05:47:08 UTC (rev 115) @@ -0,0 +1,147 @@ +<?xml version="1.0" encoding="utf-8" standalone="no"?> +<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN" + "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd"> +<svg width="600" height="250" viewBox="0 0 600 250" + xmlns="http://www.w3.org/2000/svg" + xmlns:xlink="http://www.w3.org/1999/xlink"> +<desc>Produced by GNUPLOT 4.2 patchlevel 4 </desc> + +<defs> + + + <circle id='gpDot' r='1'/> + <path id='gpPt0' style='stroke-width:0.222' d='M-1,0 h2 M0,-1 v2'/> + <path id='gpPt1' style='stroke-width:0.222' d='M-1,-1 L1,1 M1,-1 L-1,1'/> + <path id='gpPt2' style='stroke-width:0.222' d='M-1,0 L1,0 M0,-1 L0,1 M-1,-1 L1,1 M-1,1 L1,-1'/> + <rect id='gpPt3' style='stroke-width:0.222' x='-1' y='-1' width='2' height='2'/> + <use xlink:href='#gpPt3' id='gpPt4' style='fill:currentColor; stroke:none'/> + <circle id='gpPt5' style='stroke-width:0.222' cx='0' cy='0' r='1'/> + <use xlink:href='#gpPt5' id='gpPt6' style='fill:currentColor; stroke:none'/> + <path id='gpPt7' style='stroke-width:0.222' d='M0,-1.33 L-1.33,0.67 L1.33,0.67 z'/> + <use xlink:href='#gpPt7' id='gpPt8' style='fill:currentColor; stroke:none'/> + <use xlink:href='#gpPt7' id='gpPt9' transform='rotate(180)'/> + <use xlink:href='#gpPt9' id='gpPt10' style='fill:currentColor; stroke:none'/> + <use xlink:href='#gpPt3' id='gpPt11' transform='rotate(45)'/> + <use xlink:href='#gpPt11' id='gpPt12' style='fill:currentColor; stroke:none'/> +</defs> +<g style="fill:none; color:white; stroke:currentColor; stroke-width:1.00; stroke-linecap:butt; stroke-linejoin:miter"> +</g> +<g style="fill:none; color:black; stroke:currentColor; stroke-width:1.00; stroke-linecap:butt; stroke-linejoin:miter"> + <path d='M25.6,188.1 L25.6,179.1 M25.6,20.9 L25.6,29.9 '></path> + <g transform="translate(25.6,210.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan> 0</tspan> + </text> + </g> + <path d='M81.3,188.1 L81.3,183.6 M81.3,20.9 L81.3,25.4 M137.0,188.1 L137.0,179.1 M137.0,20.9 L137.0,29.9 + '></path> + <g transform="translate(137.0,210.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan> 2</tspan> + </text> + </g> + <path d='M192.7,188.1 L192.7,183.6 M192.7,20.9 L192.7,25.4 M248.4,188.1 L248.4,179.1 M248.4,20.9 L248.4,29.9 + '></path> + <g transform="translate(248.4,210.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan> 4</tspan> + </text> + </g> + <path d='M304.2,188.1 L304.2,183.6 M304.2,20.9 L304.2,25.4 M359.9,188.1 L359.9,179.1 M359.9,20.9 L359.9,29.9 + '></path> + <g transform="translate(359.9,210.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan> 6</tspan> + </text> + </g> + <path d='M415.6,188.1 L415.6,183.6 M415.6,20.9 L415.6,25.4 M471.3,188.1 L471.3,179.1 M471.3,20.9 L471.3,29.9 + '></path> + <g transform="translate(471.3,210.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan> 8</tspan> + </text> + </g> + <path d='M527.0,188.1 L527.0,183.6 M527.0,20.9 L527.0,25.4 M582.7,188.1 L582.7,179.1 M582.7,20.9 L582.7,29.9 + '></path> + <g transform="translate(582.7,210.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan> 10</tspan> + </text> + </g> + <path d='M25.6,20.9 L25.6,188.1 L582.7,188.1 L582.7,20.9 L25.6,20.9 Z '></path> + <g transform="translate(304.1,237.6)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:middle"> + <text><tspan>f/[kHz]</tspan> + </text> + </g> + <g transform="translate(97.6,48.8)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:start"> + <text><tspan>F</tspan> + </text> + </g> + <g transform="translate(237.5,48.8)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:start"> + <text><tspan>H</tspan> + <tspan font-size="9.6pt" dy="3.60pt">1</tspan> + </text> + </g> + <g transform="translate(346.8,48.8)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:start"> + <text><tspan dy="-3.60pt">H</tspan> + <tspan font-size="9.6pt" dy="3.60pt">2</tspan> + </text> + </g> + <g transform="translate(450.3,48.8)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:start"> + <text><tspan dy="-3.60pt">H</tspan> + <tspan font-size="9.6pt" dy="3.60pt">3</tspan> + </text> + </g> + <g transform="translate(551.4,48.8)" style="stroke:none; fill:black; font-family:Arial; font-size:12.00pt; text-anchor:start"> + <text><tspan dy="-3.60pt">H</tspan> + <tspan font-size="9.6pt" dy="3.60pt">4</tspan> + </text> + </g> +</g> +<g style="fill:none; color:red; stroke:currentColor; stroke-width:1.00; stroke-linecap:butt; stroke-linejoin:miter"> + <path d='M37.5,20.9 L40.3,75.4 L55.0,94.0 L69.7,103.6 L84.5,98.6 L99.2,73.6 L113.9,89.2 L128.6,124.2 + L143.3,141.2 L158.0,161.8 L172.8,165.2 L187.5,170.0 L202.2,161.8 L216.9,155.7 L231.6,116.9 L246.3,66.2 + L261.1,126.2 L275.8,140.8 L290.5,139.8 L305.2,134.0 L319.9,124.7 L334.6,121.3 L349.3,87.1 L364.1,93.5 + L378.8,110.8 L393.5,117.3 L408.2,131.0 L422.9,119.6 L437.6,96.4 L452.4,100.2 L467.1,98.6 L481.8,107.8 + L496.5,116.0 L511.2,123.7 L525.9,115.1 L540.7,109.1 L555.4,105.7 L570.1,100.5 L582.7,112.4 '></path> + <use xlink:href='#gpPt0' transform='translate(40.3,75.4) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(55.0,94.0) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(69.7,103.6) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(84.5,98.6) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(99.2,73.6) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(113.9,89.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(128.6,124.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(143.3,141.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(158.0,161.8) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(172.8,165.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(187.5,170.0) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(202.2,161.8) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(216.9,155.7) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(231.6,116.9) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(246.3,66.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(261.1,126.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(275.8,140.8) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(290.5,139.8) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(305.2,134.0) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(319.9,124.7) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(334.6,121.3) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(349.3,87.1) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(364.1,93.5) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(378.8,110.8) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(393.5,117.3) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(408.2,131.0) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(422.9,119.6) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(437.6,96.4) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(452.4,100.2) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(467.1,98.6) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(481.8,107.8) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(496.5,116.0) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(511.2,123.7) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(525.9,115.1) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(540.7,109.1) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(555.4,105.7) scale(4.50)'/> + <use xlink:href='#gpPt0' transform='translate(570.1,100.5) scale(4.50)'/> +</g> +<g style="fill:none; color:green; stroke:currentColor; stroke-width:1.00; stroke-linecap:butt; stroke-linejoin:miter"> + <path d='M106.0,188.1 L106.0,54.3 M245.9,188.1 L245.9,54.3 M355.1,188.1 L355.1,54.3 M458.7,188.1 L458.7,54.3 + M559.7,188.1 L559.7,54.3 '></path> +</g> +<g style="fill:none; color:black; stroke:currentColor; stroke-width:1.00; stroke-linecap:butt; stroke-linejoin:miter"> + <path d='M25.6,20.9 L25.6,188.1 L582.7,188.1 L582.7,20.9 L25.6,20.9 Z '></path> +</g> +</svg> + File [added]: modes.dat Delta lines: +23 -0 =================================================================== --- examples/tov/modes.dat (rev 0) +++ examples/tov/modes.dat 2011-08-22 05:47:08 UTC (rev 115) @@ -0,0 +1,23 @@ +1.44252691528028 1.e-10 F +1.44252691528028 1.e-6 + +3.95408396916149 1.e-10 H_1 +3.95408396916149 1.e-6 + +5.91494894170517 1.e-10 H_2 +5.91494894170517 1.e-6 + +7.77382493405710 1.e-10 H_3 +7.77382493405710 1.e-6 + +9.58768293806276 1.e-10 H_4 +9.58768293806276 1.e-6 + +11.3772129983097 1.e-10 H_5 +11.3772129983097 1.e-6 + +13.1520241905666 1.e-10 H_6 +13.1520241905666 1.e-6 + +14.9172321735655 1.e-10 H_7 +14.9172321735655 1.e-6 Directory: /examples/tov/tov_2/ =============================== File [added]: hydrobase::rho.maximum.asc.bz2 Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: examples/tov/tov_2/hydrobase::rho.maximum.asc.bz2 ___________________________________________________________________ File [added]: hydrobase::rho.maximum.asc.spectrum.bz2 Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: examples/tov/tov_2/hydrobase::rho.maximum.asc.spectrum.bz2 ___________________________________________________________________

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 114)
by bcmsma＠astro.rit.edu 21 Aug '11

21 Aug '11

User: bmundim Date: 2011/08/22 12:29 AM Modified: / ET.tex Log: A few typos fixed. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +6 -5 =================================================================== --- ET.tex 2011-08-22 04:21:20 UTC (rev 113) +++ ET.tex 2011-08-22 05:29:40 UTC (rev 114) @@ -349,7 +349,7 @@ In addition, the Einstein Toolkit provides computer scientists an ideal platform to perform state-of-the-art research, which directly -benefits research in other areas of science and providines an +benefits research in other areas of science and provides an immediate application of their research. One of the most prominent examples within the Einstein Toolkit is the Cactus Computational Toolkit, a framework developed by computer scientists and now used by researchers @@ -360,7 +360,7 @@ The mechanisms for the development and support of the Einstein Toolkit are designed to be open, transparent and community-driven. The complete source code, documentation and tools included in the Einstein Toolkit are distributed -under open-source licenses. The Einstein Toolkit maintains an version +under open-source licenses. The Einstein Toolkit maintains a version control system ({\tt svn.einsteintoolkit.org}) with open access that contains software supported by the Einstein Toolkit, the toolkit web pages, and documentation. An open wiki for documentation @@ -414,7 +414,7 @@ develop new application level tools for performance and correctness. The structure of the Cactus framework is completely modular, with -only a very small core proving the interfaces between modules, +only a very small core providing the interfaces between modules, both at compile- and run-time. The Cactus modules, called ``thorns'', may, and typically do, specify inter-module dependencies, e.g., to share or extend configuration information, common variables or runtime parameters. @@ -663,11 +663,12 @@ The Einstein Toolkit provides modules with the sole purpose to declare commonly used variables, and to define their meaning and units. The latter is not strictly enforced, and instead documented as part of the -module documentation. Two of the these base modules are described in the +module documentation. Two of these base modules are described in the following in more detail, because especially these quantities might be defined differently in various simulation codes. \subsubsection{ADMBase} +%BCM: the following sentence seems to be an inadvertently truncated one: Relativistic spacetime evolution methods used within the Cactus framework use different formulations, but essentially all are based on the $3+1$ ADM construction~\cite{Arnowitt:1962hi}, in which @@ -704,7 +705,7 @@ \end{itemize} The type of initial data chosen for a simulation is specified by the -parameters{\tt initial\_data} (3-metric and extrinsic curvature), +parameters {\tt initial\_data} (3-metric and extrinsic curvature), {\tt initial\_lapse}, {\tt initial\_shift}, and by parameters for their first time derivatives, {\tt initial\_dtlapse} and {\tt initial\_dtshift}, respectively.

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/tov/ (Rev. 113)
by knarf＠cct.lsu.edu 21 Aug '11

21 Aug '11

User: knarf Date: 2011/08/21 11:21 PM Modified: /examples/tov/ tov.rpar Log: update to current rpar File Changes: Directory: /examples/tov/ ========================= File [modified]: tov.rpar Delta lines: +29 -2 =================================================================== --- examples/tov/tov.rpar 2011-08-21 20:30:39 UTC (rev 112) +++ examples/tov/tov.rpar 2011-08-22 04:21:20 UTC (rev 113) @@ -178,8 +178,35 @@ ML_BSSN::ML_shift_bound = "none" ML_BSSN::ML_dtshift_bound = "none" +ActiveThorns = "SphericalSurface Dissipation" +SphericalSurface::nsurfaces = 2 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +SphericalSurface::ntheta [1] = 39 +SphericalSurface::nphi [1] = 76 +SphericalSurface::nghoststheta[1] = 2 +SphericalSurface::nghostsphi [1] = 2 + +Dissipation::verbose = "no" +Dissipation::epsdis = 0.01 +Dissipation::vars = " + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_curv + ML_BSSN::ML_trace_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift +" + + # init parameters InitBase::initial_data_setup_method = "init_some_levels" @@ -191,8 +218,8 @@ ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" -#TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::TOV_Mass[0] = 1.4 +TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +#TOVSolver::TOV_Mass[0] = 1.4 TOVSolver::TOV_Gamma[0] = 2.0 TOVSolver::TOV_K[0] = $K_TOV $K_change

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 112)
by diener＠cct.lsu.edu 21 Aug '11

21 Aug '11

User: diener Date: 2011/08/21 03:30 PM Modified: / ET.tex Log: Add the first couple of plots for the distorted rotating black hole example. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +20 -1 =================================================================== --- ET.tex 2011-08-21 20:29:34 UTC (rev 111) +++ ET.tex 2011-08-21 20:30:39 UTC (rev 112) @@ -2252,8 +2252,27 @@ resolution on the finest grid and the second number in parenthesis is the resolution on the coarsest grid. The runs where performed using the tapering evolution scheme in \codename{Carpet} in order to avoid interpolation in -time during prolongation. +time during prolongation. The initial data corresponds to a rotating black +hole perturbed by an $\ell =2, m=0$ Brill wave and as such has a non-zero +gravitational wave content. We evolved using 4th order finite differencing from +$T=0M$ until the black hole had settled down to a stationary state at $T=120M$. +\begin{figure} + \includegraphics[width=0.33\textwidth,angle=-90]{examples/kerr/figs/waves} + \includegraphics[width=0.33\textwidth,angle=-90]{examples/kerr/figs/waves_conv} + \caption{The right plot shows the extracted $\ell =2, m=0$ mode of $\Psi_4$ + as function of time from the high resultion run. The extraction was + done at $R=30M$. Shown is both the real (solid red line) and the + imaginary (dashed green line) part of the waveform. The left plot + shows for the real part of the $\ell =2, m=0$ waveforms the + difference between the medium and low resolution runs (solid red + line), the difference between the high and medium resolution runs + (dashed green line) as well as the scaled (for 4th order + convergence) difference between the medium and low resolution runs + (dotted blue line).} + \label{fig:kerr_waves} +\end{figure} + \subsection{BHB\pages{2 Bruno}} \label{sec:bbh-example} %BCM: DONE:ID parameters, TP number of collocation points, grid structure,

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/kerr/ (Rev. 111)
by diener＠cct.lsu.edu 21 Aug '11

21 Aug '11

User: diener Date: 2011/08/21 03:29 PM Added: /examples/kerr/by-a0.7-hi/ mp_psi4_l2_m0_r30.00.asc, mp_psi4_l2_m0_r40.00.asc, mp_psi4_l2_m0_r50.00.asc, mp_psi4_l2_m0_r60.00.asc /examples/kerr/by-a0.7-med/ mp_psi4_l2_m0_r30.00.asc, mp_psi4_l2_m0_r40.00.asc, mp_psi4_l2_m0_r50.00.asc, mp_psi4_l2_m0_r60.00.asc /examples/kerr/by-a0.7/ mp_psi4_l2_m0_r30.00.asc, mp_psi4_l2_m0_r40.00.asc, mp_psi4_l2_m0_r50.00.asc, mp_psi4_l2_m0_r60.00.asc /examples/kerr/figs/ waves.eps, waves.gp, waves.pdf, waves_conv.eps, waves_conv.gp, waves_conv.pdf, waves_conv_imag.eps, waves_conv_imag.gp, waves_conv_imag.pdf Log: Add waveform data and gnuplot script for generating some figures. File Changes: Directory: /examples/kerr/by-a0.7/ ================================== File [added]: mp_psi4_l2_m0_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r40.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r50.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r60.00.asc Delta lines: <Unavailable> Directory: /examples/kerr/by-a0.7-hi/ ===================================== File [added]: mp_psi4_l2_m0_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r40.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r50.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r60.00.asc Delta lines: <Unavailable> Directory: /examples/kerr/by-a0.7-med/ ====================================== File [added]: mp_psi4_l2_m0_r30.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r40.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r50.00.asc Delta lines: <Unavailable> File [added]: mp_psi4_l2_m0_r60.00.asc Delta lines: <Unavailable> Directory: /examples/kerr/figs/ =============================== File [added]: waves.eps Delta lines: <Unavailable> File [added]: waves.gp Delta lines: <Unavailable> File [added]: waves.pdf Delta lines: <Unavailable> File [added]: waves_conv.eps Delta lines: <Unavailable> File [added]: waves_conv.gp Delta lines: <Unavailable> File [added]: waves_conv.pdf Delta lines: <Unavailable> File [added]: waves_conv_imag.eps Delta lines: <Unavailable> File [added]: waves_conv_imag.gp Delta lines: <Unavailable> File [added]: waves_conv_imag.pdf Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 106)
by bcmsma＠astro.rit.edu 21 Aug '11

21 Aug '11

User: bmundim Date: 2011/08/21 12:51 PM Modified: /trunk/ einsteintoolkit.bib Log: Update .bib file File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 110)
by bcmsma＠astro.rit.edu 21 Aug '11

21 Aug '11

User: bmundim Date: 2011/08/21 12:40 PM Added: /examples/bbh/ BBHHigherRes.par, BBHMedHRes.par /examples/bbh/BBHHighRes/ meh4th_amp_phs_mp_psi4_l2_m2_r60.00.asc, meh6th_amp_phs_mp_psi4_l2_m2_r60.00.asc, meh8th_amp_phs_mp_psi4_l2_m2_r60.00.asc /examples/bbh/BBHMedHRes/ amp_phs_frq_mp_psi4_l2_m2_r60.00.asc, me4th_amp_phs_mp_psi4_l2_m2_r60.00.asc, me6th_amp_phs_mp_psi4_l2_m2_r60.00.asc, me8th_amp_phs_mp_psi4_l2_m2_r60.00.asc, ml_admconstraints::ml_ham.x.asc, mp_psi4_l2_m2_r60.00.asc, puncturetracker::pt_loc..asc /examples/bbh/figs/ amp_convergence_lomehi_4th.eps, amp_convergence_lomehi_6th.eps, amp_convergence_lomehi_8th.eps, amp_convergence_memehhi_4th.eps, amp_convergence_memehhi_6th.eps, amp_convergence_memehhi_8th.eps, amp_convergence_memehhi_8th.pdf, phase_convergence_lomehi_4th.eps, phase_convergence_lomehi_6th.eps, phase_convergence_lomehi_8th.eps, phase_convergence_memehhi_4th.eps, phase_convergence_memehhi_6th.eps, phase_convergence_memehhi_8th.eps, phase_convergence_memehhi_8th.pdf /examples/bbh/sm_scripts/ sm_conv_amp_lomehi_4th, sm_conv_amp_lomehi_6th, sm_conv_amp_lomehi_8th, sm_conv_amp_memehhi_4th, sm_conv_amp_memehhi_6th, sm_conv_amp_memehhi_8th, sm_conv_phs_lomehi_4th, sm_conv_phs_lomehi_6th, sm_conv_phs_lomehi_8th, sm_conv_phs_memehhi_4th, sm_conv_phs_memehhi_6th, sm_conv_phs_memehhi_8th Removed: /examples/bbh/figs/ amp_convergence_4th.eps, amp_convergence_6th.eps, amp_convergence_8th.eps, amp_convergence_8th.pdf, phase_convergence_4th.eps, phase_convergence_6th.eps, phase_convergence_8th.eps, phase_convergence_8th.pdf /examples/bbh/sm_scripts/ sm_conv_amp_4th, sm_conv_amp_6th, sm_conv_amp_8th, sm_conv_phs_4th, sm_conv_phs_6th, sm_conv_phs_8th Modified: / ET.tex /examples/bbh/BBHHighRes/ ml_admconstraints::ml_ham.x.asc /examples/bbh/BBHLowRes/ ml_admconstraints::ml_ham.x.asc, puncturetracker::pt_loc..asc /examples/bbh/BBHMedRes/ ml_admconstraints::ml_ham.x.asc, puncturetracker::pt_loc..asc /examples/bbh/sm_scripts/ .smhist Log: Update data, scripts and plots. First draft of the section. Text still needs some polishing. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: <Unavailable> Directory: /examples/bbh/BBHHighRes/ ==================================== File [added]: meh4th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: meh6th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: meh8th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [modified]: ml_admconstraints::ml_ham.x.asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHLowRes/ =================================== File [modified]: ml_admconstraints::ml_ham.x.asc Delta lines: <Unavailable> File [modified]: puncturetracker::pt_loc..asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHMedHRes/ ==================================== File [added]: amp_phs_frq_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: me4th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: me6th_amp_phs_mp_psi4_l2_m2_r60.00.asc Delta lines: <Unavailable> File [added]: me8th_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_r60.00.asc Delta lines: <Unavailable> File [added]: puncturetracker::pt_loc..asc Delta lines: <Unavailable> Directory: /examples/bbh/BBHMedRes/ =================================== File [modified]: ml_admconstraints::ml_ham.x.asc Delta lines: <Unavailable> File [modified]: puncturetracker::pt_loc..asc Delta lines: <Unavailable> Directory: /examples/bbh/ ========================= File [added]: BBHHigherRes.par Delta lines: <Unavailable> File [added]: BBHMedHRes.par Delta lines: <Unavailable> Directory: /examples/bbh/figs/ ============================== File [removed]: amp_convergence_4th.eps Delta lines: <Unavailable> File [removed]: amp_convergence_6th.eps Delta lines: <Unavailable> File [removed]: amp_convergence_8th.eps Delta lines: <Unavailable> File [removed]: amp_convergence_8th.pdf Delta lines: <Unavailable> File [added]: amp_convergence_lomehi_4th.eps Delta lines: <Unavailable> File [added]: amp_convergence_lomehi_6th.eps Delta lines: <Unavailable> File [added]: amp_convergence_lomehi_8th.eps Delta lines: <Unavailable> File [added]: amp_convergence_memehhi_4th.eps Delta lines: <Unavailable> File [added]: amp_convergence_memehhi_6th.eps Delta lines: <Unavailable> File [added]: amp_convergence_memehhi_8th.eps Delta lines: <Unavailable> File [added]: amp_convergence_memehhi_8th.pdf Delta lines: <Unavailable> File [removed]: phase_convergence_4th.eps Delta lines: <Unavailable> File [removed]: phase_convergence_6th.eps Delta lines: <Unavailable> File [removed]: phase_convergence_8th.eps Delta lines: <Unavailable> File [removed]: phase_convergence_8th.pdf Delta lines: <Unavailable> File [added]: phase_convergence_lomehi_4th.eps Delta lines: <Unavailable> File [added]: phase_convergence_lomehi_6th.eps Delta lines: <Unavailable> File [added]: phase_convergence_lomehi_8th.eps Delta lines: <Unavailable> File [added]: phase_convergence_memehhi_4th.eps Delta lines: <Unavailable> File [added]: phase_convergence_memehhi_6th.eps Delta lines: <Unavailable> File [added]: phase_convergence_memehhi_8th.eps Delta lines: <Unavailable> File [added]: phase_convergence_memehhi_8th.pdf Delta lines: <Unavailable> Directory: /examples/bbh/sm_scripts/ ==================================== File [modified]: .smhist Delta lines: <Unavailable> File [removed]: sm_conv_amp_4th Delta lines: <Unavailable> File [removed]: sm_conv_amp_6th Delta lines: <Unavailable> File [removed]: sm_conv_amp_8th Delta lines: <Unavailable> File [added]: sm_conv_amp_lomehi_4th Delta lines: <Unavailable> File [added]: sm_conv_amp_lomehi_6th Delta lines: <Unavailable> File [added]: sm_conv_amp_lomehi_8th Delta lines: <Unavailable> File [added]: sm_conv_amp_memehhi_4th Delta lines: <Unavailable> File [added]: sm_conv_amp_memehhi_6th Delta lines: <Unavailable> File [added]: sm_conv_amp_memehhi_8th Delta lines: <Unavailable> File [removed]: sm_conv_phs_4th Delta lines: <Unavailable> File [removed]: sm_conv_phs_6th Delta lines: <Unavailable> File [removed]: sm_conv_phs_8th Delta lines: <Unavailable> File [added]: sm_conv_phs_lomehi_4th Delta lines: <Unavailable> File [added]: sm_conv_phs_lomehi_6th Delta lines: <Unavailable> File [added]: sm_conv_phs_lomehi_8th Delta lines: <Unavailable> File [added]: sm_conv_phs_memehhi_4th Delta lines: <Unavailable> File [added]: sm_conv_phs_memehhi_6th Delta lines: <Unavailable> File [added]: sm_conv_phs_memehhi_8th Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 109)
by knarf＠cct.lsu.edu 20 Aug '11

20 Aug '11

User: knarf Date: 2011/08/20 10:36 PM Modified: / ET.tex Log: more on future work and acknowl. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +30 -15 =================================================================== --- ET.tex 2011-08-15 15:18:25 UTC (rev 108) +++ ET.tex 2011-08-21 03:36:17 UTC (rev 109) @@ -2330,21 +2330,36 @@ \end{figure} -\section{Future Work\pages{1 Frank}} +\section{Future Work\pagesdone{1 Frank}} In this paper we illustrated the current state of the ``Einstein Toolkit'', a collection of freely available and easy to use computational codes for numerical relativity. However, there is room for improvement on both the underlying infrastructure and the included physics. One of the desirable additions of physics is a proper treatment of -radiation, in particular neutrinos, but possibly even some approximation -of emmission of electro-magnetic waves. Radiation transport is, even compared -to the full GRMHD problem, computationally very expensive, especially in -three dimensions. \todo{Erik: mention PetaCactus} +radiation, in particular neutrinos. Another desirable addition would be some +approximation of emmission of electro-magnetic waves. Radiation transport is, +even compared to the full GRMHD problem, computationally very expensive, +especially in three dimensions. Several members of the CIGR team are involved +in the NSF project PetaCactus, which aims to explore these possibilities. -\todo{CCE people: mention advanced wave extraction} +Besides new additions, existing techniques require to be improved. One such +example is the improvement of current wave extraction techniques within the +Einstein Toolkit to include Cauchy characteristic wave extraction techniques, +as recently studied in~\cite{todo}. The authors of these research codes +recently agreed to make their work available to the whole community within +the Einstein toolkit. This is expected to happen for the second release of +the toolkit in 2011. -Another important goal is to increase the scalability of the Carpet AMR +A second example for an improvement of an existing method is extending the +current implementation of ideal magneto hydrodynamics to the case of non-zero +resistivity. Such implementations exist +and are described in the literature~\cite{todo}, but none are provided freely +within the Einstein Toolkit yet. Resistive MHD is important for a number of +astrophysical scenarios, one of which is the merger of two magnetized neutron +stars, a candidate for short Gamma ray bursts. + +Yet another important goal is to increase the scalability of the Carpet AMR infrastructure. It has been shown that good scaling is limited to less than a few thousend proccesses, for some of the most used simulation scenarios. In this case, work is already in progress to eliminate this bottle-neck. @@ -2357,17 +2372,17 @@ for the forseeable future. \section*{Acknowledgments\pages{0.5 All}} -\todo{Frank, ADDME} +\todo{Frank, ADD} The Einstein Toolkit is directly supported by the National Science Foundation under the grant numbers 0903973/0903782/0904015 (CIGR). -\todo{NSF OCI 0721915 ``Alpaca'', NSF OCI 0725070 ``Blue Waters'', NSF - OCI 0905046 ``PetaCactus'', NSF OCI 0941653 ``PRAC''.} -Results presented +Related projects contributed directly and indirectly to the success of CIGR, +including the NSF OCI 0721915 project ``Alpaca'', the NSF OCI 0725070 project +``Blue Waters'', the NSF OCI project 0905046 ``PetaCactus'' and the +NSF OCI 0941653 ``PRAC'' award. Results presented in this article were obtained through computations on the Louisiana Optical -Network Initiative under allocation loni\_cactus05, as well as the NSF Teragrid -under allocation TG-MCA02N014 \todo{DOE repository m152, - HLRB, LRZ, Sharcnet project - project cfz-411-aa}. +Network Initiative under allocation loni\_cactus05 and loni\_numrel06, as well +as the NSF Teragrid under allocation TG-MCA02N014, DOE repository m152, +HLRB, at the LRZ, and using the Sharcnet project cfz-411-aa. \bibliographystyle{amsplain-url} % TODO: all references should be in manifest/einsteintoolkit

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 265)
by schnetter＠cct.lsu.edu 19 Aug '11

19 Aug '11

User: eschnett Date: 2011/08/19 06:58 PM Modified: /trunk/ schedule.ccl Log: Schedule certain items before MoL_PostStepModify Schedule all those items which were scheduled "before MoL_PostStep" also "before MoL_PostStepModify". File Changes: Directory: /trunk/ ================== File [modified]: schedule.ccl Delta lines: +4 -4 =================================================================== --- trunk/schedule.ccl 2011-08-05 18:07:18 UTC (rev 264) +++ trunk/schedule.ccl 2011-08-19 23:58:27 UTC (rev 265) @@ -271,11 +271,11 @@ if (CCTK_IsThornActive("PUGH" )) { - schedule GRHydro_Rho_Minima_Setup_Final_PUGH AT CCTK_PostInitial BEFORE MoL_PostStep + schedule GRHydro_Rho_Minima_Setup_Final_PUGH AT CCTK_PostInitial BEFORE (MoL_PostStepModify MoL_PostStep) { LANG: C } "Set the value of the rest-mass density of the atmosphere which will be used during the evolution (PUGH)" - schedule GRHydro_InitialAtmosphereResetM AT CCTK_PostInitial BEFORE MoL_PostStep AFTER GRHydro_Rho_Minima_Setup_Final_PUGH + schedule GRHydro_InitialAtmosphereResetM AT CCTK_PostInitial BEFORE (MoL_PostStepModify MoL_PostStep) AFTER GRHydro_Rho_Minima_Setup_Final_PUGH { LANG: Fortran } "Use mask to enforce atmosphere at initial time" @@ -296,11 +296,11 @@ if (CCTK_IsThornActive("PUGH" )) { - schedule GRHydro_Rho_Minima_Setup_Final_PUGH AT CCTK_PostInitial BEFORE MoL_PostStep + schedule GRHydro_Rho_Minima_Setup_Final_PUGH AT CCTK_PostInitial BEFORE (MoL_PostStepModify MoL_PostStep) { LANG: C } "Set the value of the rest-mass density of the atmosphere which will be used during the evolution (PUGH)" - schedule GRHydro_InitialAtmosphereReset AT CCTK_PostInitial BEFORE MoL_PostStep AFTER GRHydro_Rho_Minima_Setup_Final_PUGH + schedule GRHydro_InitialAtmosphereReset AT CCTK_PostInitial BEFORE (MoL_PostStepModify MoL_PostStep) AFTER GRHydro_Rho_Minima_Setup_Final_PUGH { LANG: Fortran } "Use mask to enforce atmosphere at initial time"

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 August 2011