Commits March 2011

commits@lists.einsteintoolkit.org

6 participants
48 discussions

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 37)
by cott＠tapir.caltech.edu 26 Mar '11

26 Mar '11

User: cott Date: 2011/03/26 10:14 AM Added: / widetext.sty /local_bibtex/ ott_references.bib Modified: / ET.tex Log: * bunch of changes; copied in some text from another paper that needs massaging File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +282 -4 =================================================================== --- ET.tex 2011-03-22 15:00:14 UTC (rev 36) +++ ET.tex 2011-03-26 15:14:54 UTC (rev 37) @@ -16,6 +16,7 @@ \usepackage[hyphens]{url} \usepackage{wrapfig} \usepackage{multicol} +\usepackage{widetext} \usepackage[bookmarks, bookmarksopen, bookmarksnumbered]{hyperref} \usepackage[all]{hypcap} @@ -388,6 +389,7 @@ supercomputers to run the same simulation there. \subsection{Kranc\pages{1 Ian}} +\label{sec:kranc} Kranc\cite{Husa:2004ip,kranc04,krancweb} is a Mathematica application which converts a high-level continuum description of a PDE into a highly optimised module for @@ -542,13 +544,289 @@ the \codename{TwoPunctures} module can also be used to construct neutron star black hole binary initial data, when being coupled with \codename{TOVSolver}. -\subsection{Evolution Methods\pages{4 Christian}} +\subsection{Spacetime Curvature and Hydrodynamics Evolution} +\todo{Christian in charge} -\paragraph{Spacetime} +In the following, we assume that the reader is familiar with the +basics of numerical relativity and GR hydrodynamics. Detailed +introductions to numerical relativity have recently been given by +Alcubierre~\cite{alcubierre:08} and Baumgarte \& Shapiro~\cite{baumgarte:10}. +GR hydrodynamics has been reviewed by Font~\cite{font:08}. -\paragraph{Relativistic Matter} +\subsubsection{Spacetime Curvature Evolution} The Einstein Toolkit +curvature evolution code \codename{McLachlan}~\cite{brown:09} is +auto-generated from tensor equations via \codename{Kranc} +(Sec.~\ref{sec:kranc}). It implements the Einstein equations in a +$3+1$ split as a Cauchy initial boundary value problem +\cite{york:79}. For this, the Baumgarte-Shapiro-Shibata-Nakamura +(BSSN) conformal-tracefree reformulation +\cite{shibata:95,baumgarte:95,alcubierre:00} of the original +Arnowitt-Deser-Misner formalism~\cite{adm:62} is employed. +This leads to the following set of evolved variables: +\todo{continue editing here} + +\begin{eqnarray} + \phi & := & \log \left[ \frac{1}{12} \det \gamma_{ij} \right] + \\ + \tilde\gamma_{ij} & := & e^{-4\phi}\; \gamma_{ij} + \\ + K & := & g^{ij} K_{ij} + \\ + \tilde A_{ij} & := & e^{-4\phi} \left[ K_{ij} - \frac{1}{3} g_{ij} K + \right] + \\ + \tilde\Gamma^i & := & \tilde\gamma^{jk} \tilde\Gamma^i_{jk} . +\end{eqnarray} +Our exact evolution equations are as described by Eqs.~(3) to (10) of +\cite{ES-Brown2007b}, which we list here for completeness: +\begin{widetext} +\begin{eqnarray} + \partial_0 \alpha & = & -\alpha^2 f(\alpha, \phi, x^\mu) (K - + K_0(x^\mu)) + \\ + \partial_0 K & = & -e^{-4\phi} \left[ \tilde{D}^i \tilde{D}_i \alpha + + 2 \partial_i \phi \cdot \tilde{D}^i \alpha \right] + \alpha + \left( \tilde{A}^{ij} \tilde{A}_{ij} + \frac{1}{3} K^2 \right) - \alpha S + \\ + \partial_0 \beta^i & = & \alpha^2 G(\alpha,\phi,x^\mu) B^i + \\ + \partial_0 B^i & = & e^{-4\phi} H(\alpha,\phi,x^\mu) + \partial_0\tilde{\Gamma}^i - \eta^i(B^i,\alpha,x^\mu) + \\ + \partial_0 \phi & = & -\frac{\alpha}{6}\, K + + \frac{1}{6}\partial_k\beta^k + \\ + \partial_0 \tilde{\gamma}_{ij} & = & -2\alpha\tilde{A}_{ij} + + 2\tilde{\gamma}_{k(i}\partial_{j)}\beta^k + - \frac{2}{3}\tilde{\gamma}_{ij}\partial_k\beta^k + \\ + \partial_0 \tilde{A}_{ij} & = & e^{-4\phi}\left[ + \alpha\tilde{R}_{ij} + \alpha R^\phi_{ij} - \tilde{D}_i\tilde{D}_j\alpha + + 4\partial_{(i}\phi\cdot\tilde{D}_{j)}\alpha\right]^{TF} + \nonumber\\ + & & {} + \alpha K\tilde{A}_{ij} - 2\alpha\tilde{A}_{ik}\tilde{A}^k_{\; j} + + 2\tilde{A}_{k(i}\partial_{j)}\beta^k + - \frac{2}{3}\tilde{A}_{ij}\partial_k\beta^k + - \alpha e^{-4\phi} \hat{S}_{ij} + \\ + \partial_0\tilde{\Gamma}^i & = & + \tilde{\gamma}^{kl}\partial_k\partial_l\beta^i + + \frac{1}{3} \tilde{\gamma}^{ij}\partial_j\partial_k\beta^k + + \partial_k\tilde{\gamma}^{kj} \cdot \partial_j\beta^i + - \frac{2}{3}\partial_k\tilde{\gamma}^{ki} \cdot \partial_j\beta^j + \nonumber\\ + & & {} - 2\tilde{A}^{ij}\partial_j\alpha + + 2\alpha\left[ (m-1)\partial_k\tilde{A}^{ki} - \frac{2m}{3}\tilde{D}^i K + + m(\tilde{\Gamma}^i_{\; kl}\tilde{A}^{kl} + 6\tilde{A}^{ij}\partial_j\phi) + \right] - S^i +\end{eqnarray} +\end{widetext} +with the momentum constraint damping constant set to $m=1$. The stress +energy tensor $T_{\mu\nu}$ is incorporated via the projections +\begin{eqnarray} + \rho & := & \frac{1}{\alpha^2} \left( T_{00} - 2 \beta^i T_{0i} + + \beta^i \beta^j T^{ij} \right) + \\ + S & := & g^{ij} T_{ij} + \\ + S_i & := & - \frac{1}{\alpha} \left( T_{0i} - \beta^j T_{ij} \right) . +\end{eqnarray} +We have introduced the notation $\partial_0 = \partial_t - +\beta^j\partial_j$. All quantities with a tilde $\tilde{~}$ refer to +the conformal 3-metric $\tilde{\gamma}_{ij}$, which is used to +raise and lower indices. In particular, $\tilde{D}_i$ and +$\tilde{\Gamma}^k_{ij}$ refer to the covariant derivative and the +Christoffel symbols with respect to $\tilde{\gamma}_{ij}$. The +expression $[ \cdots ]^{TF}$ denotes the trace-free part of the +expression inside the parentheses, and we define the Ricci tensor +contributions +\begin{widetext} +\begin{eqnarray} +\tilde{R}_{ij} + &=& -\frac{1}{2} \tilde{\gamma}^{kl}\partial_k\partial_l\tilde{\gamma}_{ij} + + \tilde{\gamma}_{k(i}\partial_{j)}\tilde{\Gamma}^k + - \tilde{\Gamma}_{(ij)k}\partial_l\tilde{\gamma}^{lk} + + \tilde{\gamma}^{ls}\left( 2\tilde{\Gamma}^k_{\; l(i}\tilde{\Gamma}_{j)ks} + + \tilde{\Gamma}^k_{\; is}\tilde{\Gamma}_{klj} \right) +\\ +R^\phi_{ij} &=& -2\tilde{D}_i\tilde{D}_j\phi + - 2\tilde{\gamma}_{ij} \tilde{D}^k\tilde{D}_k\phi + + 4\tilde{D}_i\phi\, \tilde{D}_j\phi + - 4\tilde{\gamma}_{ij}\tilde{D}^k\phi\, \tilde{D}_k\phi . +\end{eqnarray} +\end{widetext} +This is a so-called $\phi$-variant of BSSN. The evolved gauge +variables are lapse $\alpha$, shift $\beta^i$, and a quantity $B^i$ +related to the time derivative of the shift. The gauge parameters $f$, +$G$, $H$, and $\eta$ are determined by our choice of a $1+\log$ +slicing: +\begin{eqnarray} + f(\alpha,\phi,x^\mu) & := & 2/\alpha + \\ + K_0(x^\mu) & := & 0 +\end{eqnarray} +and $\Gamma$-driver shift condition: +\begin{eqnarray} + G(\alpha,\phi,x^\mu) & := & (3/4)\, \alpha^{-2} + \\ + H(\alpha,\phi,x^\mu) & := & \exp\{4\phi\} + \\\label{eq:eta} + \eta(B^i,\alpha,x^\mu) & := & (1/2)\, B^i q(r) . +\end{eqnarray} +The expression $q(r)$ attenuates the $\Gamma$-driver depending on the +radius as described below. + +The $\Gamma$-driver shift condition is symmetry-seeking, +driving the shift $\beta^i$ to a state that renders the conformal +connection functions $\tilde\Gamma^i$ +stationary. Of course, such a stationary state cannot be achieved +while the metric is evolving, but in a stationary spacetime the time +evolution of the shift +$\beta^i$ and thus that of the spatial coordinates $x^i$ will be exponentially +damped. This damping time scale is set by the gauge parameter $\eta$ +(see Eq.~\ref{eq:eta}) which has dimension $1/T$ (inverse time). +As described, e.g., in \cite{Muller:2009jx, ES-Schnetter2010a}, this +time scale may need to be adapted in different regions of the domain +to avoid spurious high-frequency behavior in regions that otherwise +evolve only very slowly, e.g., far away from the source. + +Here we use the simple damping mechanism described in Eq.~(12) of +\cite{ES-Schnetter2010a}, which is defined as +\begin{eqnarray} + \label{eq:varying-simple} + q(r) & := & \left\{ + \begin{array}{llll} + 1 & \mathrm{for} & r \le R & \textrm{(near the origin)} + \\ + R/r & \mathrm{for} & r \ge R & \textrm{(far away)} + \end{array} + \right. +\end{eqnarray} +with a constant $R$ defining the transition radius between the +interior, where $q\approx1$, and the exterior, where $q$ falls off as +$1/r$. Eq.~\ref{eq:eta} describes how $q$ appears in the gauge +parameters. In this paper we use $R=250\,M_\odot$ ($R = +369.2\,\mathrm{km}$). + +We implement the above BSSN equations and gauge conditions in the +{\tt McLachlan} code \cite{ES-Brown2007b, ES-mclachlanweb} which is +freely available as part of the {\tt EinsteinToolkit}. {\tt McLachlan} is +auto-generated from the definition of the variables and equations in the +{\tt Mathematica} format by the {\tt Kranc} code generator \cite{kranc04, + Husa:2004ip, krancweb}. {\tt Kranc} is a suite of {\tt Mathematica} packages +comprising a computer algebra toolbox for numerical relativists. {\tt Kranc} +can be used as a ``rapid prototyping'' system for physicists or +mathematicians handling complex systems of partial differential +equations, and through integration into the {\tt Cactus} framework one can +also produce efficient production codes. + +We use fourth-order accurate finite differencing for the spacetime +variables and add a fifth-order Kreiss-Oliger dissipation term to +remove high frequency noise. We use a fourth-order Runge-Kutta time +integrator for all evolved variables. + +\subsubsection{Initial Conditions} + +We set up our initial condition from the ADM variables $g_{ij}$, +$K_{ij}$, lapse $\alpha$, and shift $\beta^i$, as provided by the +initial data discussed in Sec.~\ref{sec:initial_models}. From these we +calculate the BSSN quantities via their definition, setting $B^i=0$, +and using cubic extrapolation for $\tilde\Gamma^i$ at the outer +boundary. This extrapolation is necessary since the $\tilde\Gamma^i$ are +calculated from derivatives of the metric, and one cannot use centered +finite differencing stencils near the outer boundary. We assume that +one could instead also use one-sided derivatives to calculate +$\tilde\Gamma^i$ on the boundary. + +The extrapolation stencils distinguish between points on the faces, +edges, and corners of the grid. Points on the faces are extrapolated +via stencils perpendicular to that face, while points on the edges and +corners are extrapolated with stencils aligned with the average of the +normals of the adjoining faces. For example, points on the $(+x,+y)$ +edge are extrapolated in the $(1,1,0)$ direction, while points in the +$(+x,+y+z)$ corner are extrapolated in the $(1,1,1)$ direction. Since +several layers of boundary points have to be filled for higher order +schemes (e.g., three layers for a fourth order scheme), we proceed +outwards starting from the innermost layer. Each subsequent layer is +then defined via the points in the interior and the previously +calculated layers. + +\subsubsection{Boundary Conditions} + +During time evolution, we apply a Sommerfeld-type radiative boundary +condition to all components of the evolved BSSN variables as described +in \cite{Alcubierre2000}. The main feature of this boundary condition +is that it assumes approximate spherical symmetry of the solution, +while applying the actual boundary condition on the boundary of a +cubic grid where the face normals are not aligned with the radial +direction. This boundary condition defines the right hand side +of the BSSN state vector on the outer boundary, which is then +integrated in time as well, so that the boundary and interior are +calculated with the same order of accuracy. + +The main part of the boundary condition assumes that we have an +outgoing radial wave with some speed $v_0$: +\begin{eqnarray} + X & = & X_0 + \frac{u(r - v_0 t)}{r} +\end{eqnarray} +where $X$ is any of the tensor components of evolved variables, $X_0$ +the value at infinity, and $u$ a spherically symmetric perturbation. +Both $X_0$ and $v_0$ depend on the particular variable and have to be +specified. This implies the following differential equation: +\begin{eqnarray} + \partial_t X & = & - v^i \partial_i X - v_0\, \frac{X - X_0}{r}\,, +\end{eqnarray} +where $v^i = v_0\, x^i/r$. The spatial derivatives $\partial_i$ are +evaluated using centered finite differencing where possible, and +one-sided finite differencing elsewhere. We use second order stencils +in our implementation. + +In addition to this main part, we also account for those parts of the +solution that do not behave as a pure wave, e.g., Coulomb type terms +caused by infall of the coordinate lines. We assume that these parts +decay with a certain power $p$ of the radius. We implement this by +considering the radial derivative of the source term above, and +extrapolating according to this power-law decay. + +Given a source term $(\partial_t X)$, we define the corrected source +term $(\partial_t X)^*$ via +\begin{eqnarray} + (\partial_t X)^* & = & (\partial_t X) + \left( \frac{r}{r - n^i + \partial_i r} \right)^p\; n^i \partial_i (\partial_t X)\,, +\end{eqnarray} +where $n^i$ is the normal vector of the corresponding boundary face. +The spatial derivatives $\partial_i$ are evaluated by comparing +neighbouring grid points, corresponding to a second-order stencil +evaluated in the middle between the two neighbouring grid points. We +assume a second-order decay, i.e., we choose $p=2$. + +As with the initial conditions above, this boundary condition is +evaluated on several layers of grid points, starting from the +innermost layer. Both the extrapolation and radiative boundary +condition algorithms are implemented in the publicly available +\texttt{NewRad} component of the Einstein Toolkit. + +This boundary condition is only a coarse approximation of the actual +decay behavior of the BSSN state vector, and it does not capture the +correct behavior of the evolved variables. However, we observe that +this boundary condition leads to stable evolutions if applied +sufficiently far from the source. Errors introduced at the boundary +(both errors in the geometry and constraint violations) propagate +inwards with the speed of light \cite{ES-Brown2007b}. Gauge changes +introduced by the boundary condition, which are physically not +observable, propagate faster, with a speed up to $\sqrt{2}$ for our +gauge conditions. + + + + + +\subsubsection{Hydrodynamics: \codename{GRHydro}} + + + \subsection{Analysis\pages{4 Tanja}} Any simulated system may have many applicable quantities for analysis which @@ -696,6 +974,6 @@ \bibliographystyle{amsplain-url} -\bibliography{manifest/einsteintoolkit} +\bibliography{manifest/einsteintoolkit,local_bibtex/ott_references} \end{document} File [added]: widetext.sty Delta lines: +86 -0 =================================================================== --- widetext.sty (rev 0) +++ widetext.sty 2011-03-26 15:14:54 UTC (rev 37) @@ -0,0 +1,86 @@ +\NeedsTeXFormat{LaTeX2e} +\ProvidesPackage{widetext} + +%% Mimics the widetext environment of revtex4 for any other class package +%% Eg: article.cls +%% +%% Compiled by: Anjishnu Sarkar +%% +%% Advantages: +%% *) Supports float (eg: figure) in two column format (Advantage over +%% multicol package) +%% *) One and twocolumn exist on the same page +%% *) Flow of text shown via rule +%% *) Equal height of text when in two column format +%% +%% Acknowledgment(s): +%% 1. Instead of re-inventing the wheel, two packages (flushend, cuted) of +%% the sttools bundle are used. The sttools bundle is available from CTAN. +%% Lisence of these packages rests with their corresponding author. +%% Any bug/problem with flushend and cuted should be forwarded to their +%% corresponding package authors. +%% 2. The idea of the rule came from the following latex community website +%% http://www.latex-community.org/forum/viewtopic.php?f=5&t=2770 +%% +%% This package just defines the widetext environment and the rules. +%% +%% Usage: +%% \documentclass[a4paper,12pt,twocolumn]{article} +%% \usepackage{widetext} +%% +%% \begin{document} +%% +%% Some text in twocolumn +%% +%% \begin{widetext} +%% Text in onecolumn format. +%% \end{widetext} +%% +%% Some more text in twocolumn +%% +%% \end{document} +%%%%%%%%%%%%%%%%%%%% + +%% Package required for equal height while in 2 columns format +\IfFileExists{flushend.sty} + {\RequirePackage{flushend}} + {\typeout{} + \typeout{Package widetext error: Install the flushend package which is + a part of sttools bundle. Available from CTAN.} + \typeout{} + \stop + } + +%% Package required for onecolumn and twocolumn to exist on the same page. +%% and also required for widetext environment. +\IfFileExists{cuted.sty} + {\RequirePackage{cuted}} + {\typeout{} + \typeout{Package widetext error: Install the cuted package which is + a part of sttools bundle. Available from CTAN.} + \typeout{} + \stop + } + + +\newlength\@parindent +\setlength\@parindent{\parindent} + +\if@twocolumn + \newenvironment{widetext} + {% + \begin{strip} + \rule{\dimexpr(0.5\textwidth-0.5\columnsep-0.4pt)}{0.4pt}% + \rule{0.4pt}{6pt} + \par %\vspace{6pt} + \parindent \@parindent + }% + {% + \par + \hfill\rule[-6pt]{0.4pt}{6.4pt}% + \rule{\dimexpr(0.5\textwidth-0.5\columnsep-1pt)}{0.4pt} + \end{strip} + } +\else + \newenvironment{widetext}{}{} +\fi Directory: /local_bibtex/ ========================= File [added]: ott_references.bib Delta lines: +104 -0 =================================================================== --- local_bibtex/ott_references.bib (rev 0) +++ local_bibtex/ott_references.bib 2011-03-26 15:14:54 UTC (rev 37) @@ -0,0 +1,104 @@ +@incollection{york:79, + author = "York, James W.", + title = "Kinematics and Dynamics of General Relativity", + pages = "83--126", + editor = "Larry L. Smarr", + booktitle = "Sources of Gravitational Radiation", + publisher = "Cambridge University Press", + address = "Cambridge, UK", + year = 1979, + snote = "Proceedings of the Battelle Seattle Workshop, + 24 July -- 4 August, 1978", +} + +@Book{adm:62, +author = {R.~Arnowitt and S.~Deser and C.~Misner}, +title = {Dynamics of General Relativity, in "Gravitation: An Introduction to Current Research"}, +publisher = {Wiley, New York}, +year = {1962}, +} + +@article{shibata:95, + key = {Shibata95}, + author = {Masaru Shibata and Takashi Nakamura}, + journal = {Phys. Rev. D}, + title = {Evolution of three-dimensional gravitational waves: + {H}armonic slicing case}, + volume = {52}, + year = {1995}, + pages = {5428} +} + +@article{baumgarte:95, + key = {Baumgarte99}, + author = {Thomas W. Baumgarte and Stuart L. Shapiro}, + title = {On the numerical integration of {E}instein's field + equations}, + journal = {Phys. Rev. D}, + year = 1999, + volume = 59, + pages = 024007, +} + +@article{alcubierre:00, + author = {M. Alcubierre and B. Br\"{u}gmann and T. Dramlitsch and J. A. Font + and Philippos Papadopoulos and E. Seidel and N. Stergioulas and R. Takahashi}, + title = {Towards a stable numerical evolution of strongly gravitating + systems in general relativity: The conformal treatments}, + journal = {Phys. Rev. D}, + volume = {62}, + year = 2000, + pages = 044034, +} + + +@article{alcubierre:03a, + author = {Miguel Alcubierre and Bernd Br\"ugmann and Peter Diener and + Michael Koppitz and Denis Pollney and Edward Seidel and Ryoji Takahashi}, + title = "Gauge conditions for long-term numerical black hole + evolutions without excision", + journal = {Phys. Rev. D}, + volume = {67}, + pages = {084023}, + year = 2003, +} + +@BOOK{alcubierre:08, + author = {{Alcubierre}, M.}, + title = "{Introduction to 3+1 Numerical Relativity}", +booktitle = {Introduction to 3+1 Numerical Relativity}, + year = 2008, + editor = "{Alcubierre, M.}", +publisher = {Oxford University Press}, +} + +@BOOK{baumgarte:10, + author = {{Baumgarte}, T.~W. and {Shapiro}, S.~L.}, + title = "{Numerical Relativity: Solving Einstein's Equations on the Computer}", + year = 2010, + publisher = "Cambridge University Press", +} + + +@ARTICLE{font:08, + author = {{Font}, J.~A.}, + title = "{Numerical Hydrodynamics and Magnetohydrodynamics in General Relativity}", + journal = {Liv. Rev. Rel.}, + keywords = {Hydrodynamics, Relativistic hydrodynamics, Numerical relativity, Magnetohydrodynamics}, + year = 2008, + month = sep, + volume = 11, + pages = {7}, +} + +@Article{brown:09, + author = {David Brown and Peter Diener and Olivier Sarbach and + Erik Schnetter and Manuel Tiglio}, + title = {Turduckening black holes: an analytical and + computational study}, + journal = {Phys. Rev. D}, + year = 2009, + volume = 79, + pages = 044023, +} +

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_LR/ (Rev. 1)
by knarf＠cct.lsu.edu 23 Mar '11

23 Mar '11

User: knarf Date: 2011/03/23 01:00 PM Added: / test Log: test commit messages File Changes: Directory: / ============ File [added]: test Delta lines: +0 -0 Added: test

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_LR/ (Rev. 2)
by knarf＠cct.lsu.edu 23 Mar '11

23 Mar '11

User: knarf Date: 2011/03/23 12:22 PM Removed: / test Log: test commit messages File Changes: Directory: / ============ File [removed]: test Delta lines: +0 -0 Deleted: test

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/figures/ (Rev. 36)
by knarf＠cct.lsu.edu 22 Mar '11

22 Mar '11

User: knarf Date: 2011/03/22 10:00 AM Added: /figures/ Log: add figures directory

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 35)
by knarf＠cct.lsu.edu 22 Mar '11

22 Mar '11

User: knarf Date: 2011/03/22 09:55 AM Modified: / ET.tex Log: add comments with names and length to sections File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +37 -54 =================================================================== --- ET.tex 2011-03-14 15:08:29 UTC (rev 34) +++ ET.tex 2011-03-22 14:55:47 UTC (rev 35) @@ -29,7 +29,9 @@ \def\gt#1{Georgia Tech#1 (GaTech#1)\gdef\gt{GaTech}} \newcommand{\codename}[1]{\texttt{#1}} -\newcommand{\todo}[1]{{\color{blue}$\blacksquare$~\textsf{[TODO: #1]}}} +\newcommand{\todo}[1]{{\color{red}$\blacksquare$~\textsf{[TODO: #1]}}} +\newcommand{\comment}[1]{{\color{blue}$\blacksquare$~\textsf{[Comment: #1]}}} +\newcommand{\pages}[1]{{\color{blue}$\blacksquare$~\textsf{[#1]}}} % Don't use tt font for urls \urlstyle{rm} @@ -69,19 +71,16 @@ -\thanks{{}$^{(1)}$ Department of Computer Science, - Louisiana State University, Baton Rouge, LA 70803, USA} +%\thanks{{}$^{(1)}$ Department of Computer Science, +% Louisiana State University, Baton Rouge, LA 70803, USA} +%\thanks{{}$^{(2)}$ Department of Physics \& Astronomy, +% Louisiana State University, Baton Rouge, LA 70803, USA} +%\thanks{{}$^{(3)}$ Center for Computation \& Technology, +% Louisiana State University, Baton Rouge, LA 70803, USA} +%\thanks{{}$^{(4)}$ Center for Computational Relativity \& Gravitation, +% School of Mathematical Sciences, +% Rochester Institute of Technology, Rochester, New York 14623, USA} -\thanks{{}$^{(2)}$ Department of Physics \& Astronomy, - Louisiana State University, Baton Rouge, LA 70803, USA} - -\thanks{{}$^{(3)}$ Center for Computation \& Technology, - Louisiana State University, Baton Rouge, LA 70803, USA} - -\thanks{{}$^{(4)}$ Center for Computational Relativity \& Gravitation, - School of Mathematical Sciences, - Rochester Institute of Technology, Rochester, New York 14623, USA} - % TODO: Already add a few affiliations which we know are going % to be used @@ -91,8 +90,11 @@ \maketitle +% We will remove this at the end, but we leave it in for ourselfes now +\tableofcontents -\section{Introduction} +\section{Introduction\pages{1 Frank}} +\comment{About 1 page} Scientific progress in the field of numerical relativity has always been closely tied with the availability and ease-of-use of enabling @@ -168,7 +170,7 @@ -\section{Requirements} +\section{Requirements\pages{2 Frank}} \subsection{Scientific} @@ -246,7 +248,7 @@ other areas of science and hereby providing an immediate application of their research. -\section{Design and Strategy for the Einstein Toolkit} +\section{Design and Strategy\pages{1 Frank}} The mechanisms for the development and support of the Einstein Toolkit are designed to be open, transparent and community driven. All the source code, @@ -274,10 +276,9 @@ Einstein Toolkit Maintainer, and is generally open for discussion on the users mailing list. -\section{Core Technologies} +\section{Core Technologies\pages{0.5 Frank}} -\subsection{Cactus Framework} -\todo{1/2 page Gab} +\subsection{Cactus Framework\pages{1 Gab}} The \textbf{Cactus Framework}~\cite{CS_cactus_web,CS_Goodale02a,CS_cactususersguide} is @@ -311,7 +312,7 @@ The Einstein Toolkit now collects the widely used parts of CactusEinstein, combined with contributions from the community. -\subsection{Adaptive Mesh Refinement} +\subsection{Adaptive Mesh Refinement\pages{1 Erik}} \todo{ES: add/check citations} \todo{ES: remove names} @@ -354,7 +355,7 @@ % E. Schnetter \todo{ES: do we want names here?} and is funded via several % NSF awards \todo{ES: list them somewhere}. -\subsection{Simulation Factory} +\subsection{Simulation Factory\pages{1 Erik/Michael}} It is unfortunate that today's supercomputers differ significantly in their hardware configuration, available software, directory structure, @@ -386,8 +387,7 @@ The exact same set of commands can be used on all other supported supercomputers to run the same simulation there. -\subsection{Kranc} -\todo{1/3 page Ian} +\subsection{Kranc\pages{1 Ian}} Kranc\cite{Husa:2004ip,kranc04,krancweb} is a Mathematica application which converts a high-level continuum description of a PDE into a highly optimised module for @@ -424,9 +424,9 @@ as well as the wave extraction thorn WeylScal4, are both generated using Kranc, and hence support all the above features. -\section{Components} +\section{Components\pages{0.5 Frank}} -\subsection{Base Modules} +\subsection{Base Modules\pages{4 Frank}} \label{sec:base_modules} Modular designs have proven to be essential for distributed development of complex software systems, and require the use of well defined interfaces. @@ -485,8 +485,7 @@ to the Faraday tensor and the unit normal of the foliation of spacetime $n^\mu$. -\subsection{Initial Data} -\todo{1/2 page Josh, Bruno} +\subsection{Initial Data\pages{4 Josh/Bruno}} The Einstein Toolkit contains many modules used to generate initial data for general relativistic simulations, including both vacuum and hydrodynamical configurations. @@ -543,16 +542,14 @@ the \codename{TwoPunctures} module can also be used to construct neutron star black hole binary initial data, when being coupled with \codename{TOVSolver}. -\subsection{Evolution Methods} -\todo{1/2 page Christian} +\subsection{Evolution Methods\pages{4 Christian}} \paragraph{Spacetime} \paragraph{Relativistic Matter} -\subsection{Analysis} -%\todo{1/2 page Tanja} +\subsection{Analysis\pages{4 Tanja}} Any simulated system may have many applicable quantities for analysis which are most efficiently calculated during evolution. The Einstein Toolkit @@ -589,8 +586,7 @@ Weyl scalars, as well as any other quantity, can be projected onto a basis of spin-weighted spherical harmonics using the \codename{Multipole} module. -\subsection{Relativity Tools} -\todo{1/2 page Peter} +\subsection{Relativity Tools\pages{3 Peter}} \paragraph{Black Hole Excision} @@ -621,8 +617,7 @@ -%\subsection{Infrastructure and Numerical Methods} -\subsection{Simulation Domain, Symmetries, Boundaries} +\subsection{Simulation Domain, Symmetries, Boundaries\pages{3 Roland?}} \paragraph{Domains and Coordinates.} The simulation domain is specified via the parameter file \todo{ES: @@ -671,36 +666,24 @@ -\section{Examples} -\todo{2 pages in total, 1 plot per paragraph, Frank coordinates} +\section{Examples\pages{0.5 Frank}} -\paragraph{Kerr-Schild} -\todo{1/2 page, Bruno} -Show stable Kerr-evolution. What is commonly shown for this nowadays, the -spectrum of the ringdown of a perturbation? +\subsection{Vacuum\pages{2 Bruno}} -\paragraph{BH-binary} -\todo{1/2 page, Peter} -Show binary merger and plot of waveform in either time or frequency domain. +\subsection{TOV\pages{2 Frank}} +Show stable TOV evolution -\paragraph{TOV+mEOS} -\todo{1/2 page, Frank} -Show stable TOV evolution, once with Polytropes and once with mEOS, -both including oscillation spectra - -\paragraph{Collapse} -\todo{1/2 page, Frank, Christian} +\subsection{Collapse\pages{2 Christian}} Show TOV collapse and BH formation -\section{Future Work} -\todo{Frank} +\section{Future Work\pages{1 Frank}} - MHD - massively parallel AMR - radiation transport - advanced wave extraction -\section*{Acknowledgments} +\section*{Acknowledgments\pages{0.5 All}} \todo{Frank, ADDME} The Einstein Toolkit is directly supported by the National Science Foundation under the grant numbers 0903973/0903782/0904015 (CIGR).

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 38)
by cott＠tapir.caltech.edu 20 Mar '11

20 Mar '11

User: cott Date: 2011/03/20 02:46 PM Modified: /trunk/src/nuc_eos/ nuc_eos.F90 Log: * improve error handling in nuc_eos a bit. File Changes: Directory: /trunk/src/nuc_eos/ ============================== File [modified]: nuc_eos.F90 Delta lines: +10 -5 =================================================================== --- trunk/src/nuc_eos/nuc_eos.F90 2011-02-07 05:36:36 UTC (rev 37) +++ trunk/src/nuc_eos/nuc_eos.F90 2011-03-20 19:46:09 UTC (rev 38) @@ -48,11 +48,13 @@ keyerrr = 0 if(xrho.gt.eos_rhomax) then - call CCTK_WARN(0, "nuc_eos: rho > rhomax") + keyerr = 103 + return endif if(xrho.lt.eos_rhomin) then - call CCTK_WARN (0, "nuc_eos: rho < rhomin") + keyerr = 104 + return endif if(xye.gt.eos_yemax) then @@ -228,7 +230,8 @@ keyerrr = 0 if(xrho.gt.eos_rhomax) then - call CCTK_WARN (0, "nuc_eos: rho > rhomax") + keyerr = 103 + return endif if(xrho.lt.eos_rhomin*1.126d0) then @@ -351,7 +354,8 @@ keyerrt = 0 if(xrho.gt.eos_rhomax) then - call CCTK_WARN (0, "nuc_eos: rho > rhomax") + keyerr = 103 + return endif if(xrho.lt.eos_rhomin*1.126d0) then @@ -430,7 +434,8 @@ keyerrt = 0 if(xrho.gt.eos_rhomax) then - call CCTK_WARN (0, "nuc_eos: rho > rhomax") + keyerr = 103 + return endif if(xrho.lt.eos_rhomin*1.126d0) then

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 34)
by schnetter＠cct.lsu.edu 14 Mar '11

14 Mar '11

User: eschnett Date: 2011/03/14 10:08 AM Modified: / ET.tex Log: Remove Schnetter's name from code description File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +10 -4 =================================================================== --- ET.tex 2011-03-07 15:53:14 UTC (rev 33) +++ ET.tex 2011-03-14 15:08:29 UTC (rev 34) @@ -289,7 +289,7 @@ formats. Generic interfaces are used, making it possible to use external packages and improved modules which are immediately available to its users. Cactus is involved in the NSF Blue Waters consortium -for petascale computing, has funding from NSF SDCI (PI Schnetter) to +for petascale computing, has funding from NSF SDCI to develop new application level tools for performance and correctness, as well as the current NSF XiRel award for scalable AMR\@. \S\ref{broader} describes the growing impact of Cactus in other @@ -313,6 +313,11 @@ \subsection{Adaptive Mesh Refinement} +\todo{ES: add/check citations} +\todo{ES: remove names} +\todo{ES: describe grants in acknowledgements, also in introduction, + but not in individual sections} + In Cactus, infrastructure capabilities such as memory management, parallelisation, time evolution, mesh refinement, and I/O are delegated to a set of special \emph{driver} components. This helps @@ -344,9 +349,10 @@ hours) were used via Carpet by more than a dozen research groups world-wide. To date, more than 90 peer-reviewed publication and more than 15 student theses are based on Carpet \todo{ES: update these numbers} -\cite{CarpetCode:web}. Carpet's continued development is overseen by -E. Schnetter \todo{ES: do we want names here?} and is funded via several -NSF awards \todo{ES: list them somewhere}. +\cite{CarpetCode:web}. +% Carpet's continued development is overseen by +% E. Schnetter \todo{ES: do we want names here?} and is funded via several +% NSF awards \todo{ES: list them somewhere}. \subsection{Simulation Factory}

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 224)
by cott＠tapir.caltech.edu 12 Mar '11

12 Mar '11

User: cott Date: 2011/03/12 09:05 AM Modified: /trunk/src/ GRHydro_Prim2Con.F90 Log: * optimize hot P2C routine a bit. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Prim2Con.F90 Delta lines: +14 -13 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2011-03-12 14:25:11 UTC (rev 223) +++ trunk/src/GRHydro_Prim2Con.F90 2011-03-12 15:05:33 UTC (rev 224) @@ -227,6 +227,7 @@ deps(1), dpress(1), w, vlowx, vlowy, vlowz CCTK_REAL :: temp(1),ye(1), x, y, z CCTK_INT :: handle, GRHydro_reflevel, ii, jj, kk + CCTK_REAL :: sdet, h character(len=512) warnline ! begin EOS Omni vars @@ -253,7 +254,7 @@ drho,deps,temp,ye,dpress,keyerr,anyerr) keytemp=0 if(anyerr.ne.0) then - !OMP CRITICAL + !$OMP CRITICAL(p2c1) call CCTK_WARN(1,"EOS error in prim2con_hot: lev 2") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(1,warnline) @@ -263,14 +264,13 @@ call CCTK_WARN(1,warnline) write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(1,warnline) - !OMP END CRITICAL + !$OMP END CRITICAL(p2c1) endif else ! This is a way of recovering even on finer refinement levels: ! Use the average temperature at the interface instead of the ! reconstructed specific internal energy. -!! !OMP CRITICAL -#if 0 + !$OMP CRITICAL(p2c2) call CCTK_WARN(1,"EOS error in prim2con_hot: NOW using averaged temp!") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(1,warnline) @@ -280,15 +280,14 @@ call CCTK_WARN(1,warnline) write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(1,warnline) -#endif -!!! !OMP END CRITICAL + !$OMP END CRITICAL(p2c2) keytemp=1 temp = temp0 call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n,& drho,deps,temp,ye,dpress,keyerr,anyerr) keytemp=0 if(anyerr.ne.0) then - !OMP CRITICAL + !$OMP CRITICAL(p2c3) call CCTK_WARN(1,"EOS error in prim2con_hot") write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z call CCTK_WARN(1,warnline) @@ -299,7 +298,7 @@ write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel call CCTK_WARN(1,warnline) call CCTK_WARN(0,"Aborting!!!") - !OMP END CRITICAL + !$OMP END CRITICAL(p2c3) endif endif endif @@ -309,11 +308,13 @@ vlowy = gxy*dvelx + gyy*dvely + gyz*dvelz vlowz = gxz*dvelx + gyz*dvely + gzz*dvelz - ddens = sqrt(det) * drho(1) * w - dsx = sqrt(det) * (drho(1)*(1+deps(1))+dpress(1))*w*w * vlowx - dsy = sqrt(det) * (drho(1)*(1+deps(1))+dpress(1))*w*w * vlowy - dsz = sqrt(det) * (drho(1)*(1+deps(1))+dpress(1))*w*w * vlowz - dtau = sqrt(det) * ((drho(1)*(1+deps(1))+dpress(1))*w*w - dpress(1)) - ddens + sdet = sqrt(det) + h = drho(1)*(1.0d0+deps(1))+dpress(1) + ddens = sdet * drho(1) * w + dsx = sdet * h*w*w * vlowx + dsy = sdet * h*w*w * vlowy + dsz = sdet * h*w*w * vlowz + dtau = sdet * (h*w*w - dpress(1)) - ddens end subroutine prim2con_hot

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 223)
by cott＠tapir.caltech.edu 12 Mar '11

12 Mar '11

User: cott Date: 2011/03/12 08:25 AM Modified: /trunk/ interface.ccl, param.ccl, schedule.ccl /trunk/src/ GRHydro_Boundaries.F90, GRHydro_BoundariesM.F90, make.code.defn Log: * remove GZ patch system code * test suites pass File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Boundaries.F90 Delta lines: +0 -156 =================================================================== --- trunk/src/GRHydro_Boundaries.F90 2011-03-12 04:06:44 UTC (rev 222) +++ trunk/src/GRHydro_Boundaries.F90 2011-03-12 14:25:11 UTC (rev 223) @@ -257,159 +257,3 @@ end subroutine GRHydro_Boundaries - /*@@ - @routine GRHydro_OutflowBoundaries - @date Tue May 17 16:58:06 2005 - @author Ian Hawke - @desc - Set outflow boundaries over only part of the domain. - This is designed to be used with GZPatchSystem. - @enddesc - @calls - @calledby - @history - - @endhistory - -@@*/ - -subroutine GRHydro_OutflowBoundaries(CCTK_ARGUMENTS) - - implicit none - - DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - DECLARE_CCTK_FUNCTIONS - - integer, dimension(3) :: sw - integer :: ierr - integer :: i,j,k - - CCTK_REAL, dimension(3) :: posn - - CCTK_REAL :: det,psi4pt - - sw = GRHydro_stencil - - if (r(1,1,1) < r(1,1,cctk_lsh(3))) then - - if (cctk_bbox(6) .ne. 0) then - - do k = cctk_lsh(3) - sw(3), cctk_lsh(3) - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) - - posn(1) = x(i,j,k) - posn(2) = y(i,j,k) - posn(3) = z(i,j,k) - - if (dot_product(outflowboundary_normal, posn) > 0.d0) then - - rho(i,j,k) = rho(i,j,k-1) - velx(i,j,k) = velx(i,j,k-1) - vely(i,j,k) = vely(i,j,k-1) - velz(i,j,k) = velz(i,j,k-1) - eps(i,j,k) = eps(i,j,k-1) - press(i,j,k) = press(i,j,k-1) - w_lorentz(i,j,k) = w_lorentz(i,j,k-1) - if(evolve_mhd.ne.0) then - Bvecx(i,j,k) = Bvecx(i,j,k-1) - Bvecy(i,j,k) = Bvecy(i,j,k-1) - Bvecz(i,j,k) = Bvecz(i,j,k-1) - if(clean_divergence.ne.0) then - psidc(i,j,k) = psidc(i,j,k-1) - endif - endif - - psi4pt = 1.0d0 - det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\ - gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) - - if(evolve_mhd.ne.0) then - call prim2conM(GRHydro_eos_handle,psi4pt*gxx(i,j,k),& - psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),& - psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),& - det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),& - tau(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),& - rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),& - eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)) - else - call prim2con(GRHydro_eos_handle,psi4pt*gxx(i,j,k),& - psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),& - psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),& - det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),& - tau(i,j,k),rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),& - eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)) - endif - end if - - end do - end do - end do - - end if - - else - - if (cctk_bbox(5) .ne. 0) then - - do k = sw(3), 1, -1 - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) - - posn(1) = x(i,j,k) - posn(2) = y(i,j,k) - posn(3) = z(i,j,k) - - if (dot_product(outflowboundary_normal, posn) > 0.d0) then - - rho(i,j,k) = rho(i,j,k+1) - velx(i,j,k) = velx(i,j,k+1) - vely(i,j,k) = vely(i,j,k+1) - velz(i,j,k) = velz(i,j,k+1) - eps(i,j,k) = eps(i,j,k+1) - press(i,j,k) = press(i,j,k+1) - w_lorentz(i,j,k) = w_lorentz(i,j,k+1) - - if(evolve_mhd.ne.0) then - Bvecx(i,j,k) = Bvecx(i,j,k+1) - Bvecy(i,j,k) = Bvecy(i,j,k+1) - Bvecz(i,j,k) = Bvecz(i,j,k+1) - if(clean_divergence.ne.0) then - psidc(i,j,k) = psidc(i,j,k+1) - endif - endif - - psi4pt = 1.0d0 - det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\ - gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) - if(evolve_mhd.ne.0) then - call prim2conM(GRHydro_eos_handle,psi4pt*gxx(i,j,k),& - psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),& - psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),& - det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),& - tau(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),& - rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),& - eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)) - else - call prim2con(GRHydro_eos_handle,psi4pt*gxx(i,j,k),& - psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),& - psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),& - det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),& - tau(i,j,k),rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),& - eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)) - endif - end if - - end do - end do - end do - - end if - - end if - -end subroutine GRHydro_OutflowBoundaries - - - File [modified]: GRHydro_BoundariesM.F90 Delta lines: +0 -131 =================================================================== --- trunk/src/GRHydro_BoundariesM.F90 2011-03-12 04:06:44 UTC (rev 222) +++ trunk/src/GRHydro_BoundariesM.F90 2011-03-12 14:25:11 UTC (rev 223) @@ -233,137 +233,6 @@ end subroutine GRHydro_BoundariesM - /*@@ - @routine GRHydro_OutflowBoundariesM - @date Aug 30, 2010 - @author Joshua Faber, Ian Hawke - @desc - Set outflow boundaries over only part of the domain. - This is designed to be used with GZPatchSystem. - @enddesc - @calls - @calledby - @history - - @endhistory -@@*/ -subroutine GRHydro_OutflowBoundariesM(CCTK_ARGUMENTS) - - implicit none - - DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - DECLARE_CCTK_FUNCTIONS - - integer, dimension(3) :: sw - integer :: ierr - integer :: i,j,k - - CCTK_REAL, dimension(3) :: posn - CCTK_REAL :: det,psi4pt - - sw = GRHydro_stencil - - if (r(1,1,1) < r(1,1,cctk_lsh(3))) then - - if (cctk_bbox(6) .ne. 0) then - - do k = cctk_lsh(3) - sw(3), cctk_lsh(3) - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) - - posn(1) = x(i,j,k) - posn(2) = y(i,j,k) - posn(3) = z(i,j,k) - - if (dot_product(outflowboundary_normal, posn) > 0.d0) then - - rho(i,j,k) = rho(i,j,k-1) - velx(i,j,k) = velx(i,j,k-1) - vely(i,j,k) = vely(i,j,k-1) - velz(i,j,k) = velz(i,j,k-1) - eps(i,j,k) = eps(i,j,k-1) - press(i,j,k) = press(i,j,k-1) - w_lorentz(i,j,k) = w_lorentz(i,j,k-1) - Bvecx(i,j,k) = Bvecx(i,j,k-1) - Bvecy(i,j,k) = Bvecy(i,j,k-1) - Bvecz(i,j,k) = Bvecz(i,j,k-1) - if(clean_divergence.ne.0) then - psidc(i,j,k) = psidc(i,j,k-1) - endif - - psi4pt = 1.d0 - det=SPATIAL_DETERMINANT(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 prim2conM(GRHydro_eos_handle,psi4pt*gxx(i,j,k),& - psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),& - psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),& - det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),& - tau(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),& - rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),& - eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)) - - end if - - end do - end do - end do - - end if - - else - - if (cctk_bbox(5) .ne. 0) then - - do k = sw(3), 1, -1 - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) - - posn(1) = x(i,j,k) - posn(2) = y(i,j,k) - posn(3) = z(i,j,k) - - if (dot_product(outflowboundary_normal, posn) > 0.d0) then - - rho(i,j,k) = rho(i,j,k+1) - velx(i,j,k) = velx(i,j,k+1) - vely(i,j,k) = vely(i,j,k+1) - velz(i,j,k) = velz(i,j,k+1) - eps(i,j,k) = eps(i,j,k+1) - press(i,j,k) = press(i,j,k+1) - w_lorentz(i,j,k) = w_lorentz(i,j,k+1) - Bvecx(i,j,k) = Bvecx(i,j,k+1) - Bvecy(i,j,k) = Bvecy(i,j,k+1) - Bvecz(i,j,k) = Bvecz(i,j,k+1) - if(clean_divergence.ne.0) then - psidc(i,j,k) = psidc(i,j,k+1) - endif - - psi4pt = 1.d0 - det=SPATIAL_DETERMINANT(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 prim2conM(GRHydro_eos_handle,psi4pt*gxx(i,j,k),& - psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),& - psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),& - det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),& - tau(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),& - rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),& - eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)) - - end if - - end do - end do - end do - - end if - - end if - -end subroutine GRHydro_OutflowBoundariesM - - - File [modified]: make.code.defn Delta lines: +0 -1 =================================================================== --- trunk/src/make.code.defn 2011-03-12 04:06:44 UTC (rev 222) +++ trunk/src/make.code.defn 2011-03-12 14:25:11 UTC (rev 223) @@ -28,7 +28,6 @@ GRHydro_Reconstruct.F90 \ GRHydro_ReconstructPoly.F90 \ GRHydro_RegisterMask.c \ - GRHydro_RegisterGZ.cc \ GRHydro_RegisterVars.cc \ GRHydro_RiemannSolve.F90 \ GRHydro_RoeAverage.F90 \ Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +0 -11 =================================================================== --- trunk/interface.ccl 2011-03-12 04:06:44 UTC (rev 222) +++ trunk/interface.ccl 2011-03-12 14:25:11 UTC (rev 223) @@ -206,18 +206,7 @@ CCTK_REAL IN z ) USES FUNCTION SetDriftCorrectPosition -############################################################## -### Functions provided by GZPatchSystem for registration ### -############################################################## -CCTK_INT FUNCTION GZPatchSystem_register_sync(CCTK_STRING IN cactus_group_name) -USES FUNCTION GZPatchSystem_register_sync -CCTK_INT FUNCTION \ - GZPatchSystem_register_cxform(CCTK_INT IN time_level, \ - CCTK_STRING IN cactus_group_name) -USES FUNCTION GZPatchSystem_register_cxform - - ####################################################### ### Functions provided by the Omni EOS interface ### ####################################################### File [modified]: param.ccl Delta lines: +0 -11 =================================================================== --- trunk/param.ccl 2011-03-12 04:06:44 UTC (rev 222) +++ trunk/param.ccl 2011-03-12 14:25:11 UTC (rev 223) @@ -515,17 +515,6 @@ ".*" :: "Should be in the full Thorn::GF format" } "" -#Setting outflow boundaries over only some of the domain for multipatch runs - -BOOLEAN outflow_boundaries "Use the outflow boundary conditions?" -{ -} "no" - -REAL outflowboundary_normal[3] "The normal to the plane; to the right of this outflow conditions will be set" -{ - *:* :: "Anything" -} 0.0 - # NaN detection level INT GRHydro_NaN_verbose "The warning level for NaNs occuring within GRHydro" File [modified]: schedule.ccl Delta lines: +0 -21 =================================================================== --- trunk/schedule.ccl 2011-03-12 04:06:44 UTC (rev 222) +++ trunk/schedule.ccl 2011-03-12 14:25:11 UTC (rev 223) @@ -212,20 +212,7 @@ # LANG: C # } "Set up the weighting terms for the fluxes" -################################### -### Register with GZPatchSystem ### -################################### -SCHEDULE GROUP GZPatchSystem_register -{ -} "Tell Cactus that this group exists, but is not scheduled from here" - -SCHEDULE GRHydro_register_GZPatchSystem IN GZPatchSystem_register -{ - LANG: C - OPTIONS: meta -} "register to-be-interpatch-synchronized variables with GZPatchSystem" - ################################# ### Set the handle for the EOS ### ################################## @@ -1002,14 +989,6 @@ STORAGE:GRHydro_tracer_rhs } -if (outflow_boundaries) -{ - schedule GRHydro_OutflowBoundaries IN HydroBase_PostStep BEFORE GRHydro_Bound - { - LANG: Fortran - } "Outflow boundaries over only some of the domain" -} - # This should not be used anymore and should be removed after some time schedule group Do_GRHydro_Boundaries IN HydroBase_Boundaries {

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 222)
by cott＠tapir.caltech.edu 11 Mar '11

11 Mar '11

User: cott Date: 2011/03/11 10:06 PM Modified: /trunk/test/ GRHydro_test_shock.par, GRHydro_test_shock_ppm.par, GRHydro_test_tov_ppm.par, GRHydro_test_tov_ppm_ML.par, GRHydro_test_tov_ppm_ML_disable_internal_excision.par, GRHydro_test_tov_ppm_disable_internal_excision.par, GRHydro_test_tov_ppm_no_trp.par, GRHydro_test_tov_ppm_no_trp_ML.par Log: * modify test suite parameter files: * those using carpet now activate loopcontrol * change GRHydro::w_lorentz to HydroBase::w_lorentz File Changes: Directory: /trunk/test/ ======================= File [modified]: GRHydro_test_shock.par Delta lines: +1 -1 =================================================================== --- trunk/test/GRHydro_test_shock.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_shock.par 2011-03-12 04:06:44 UTC (rev 222) @@ -28,7 +28,7 @@ IOBasic::outInfo_vars = "HydroBase::rho" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon" ActiveThorns = "nanchecker mol" File [modified]: GRHydro_test_shock_ppm.par Delta lines: +1 -1 =================================================================== --- trunk/test/GRHydro_test_shock_ppm.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_shock_ppm.par 2011-03-12 04:06:44 UTC (rev 222) @@ -28,7 +28,7 @@ IOBasic::outInfo_vars = "HydroBase::rho" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon" ActiveThorns = "nanchecker mol" File [modified]: GRHydro_test_tov_ppm.par Delta lines: +2 -2 =================================================================== --- trunk/test/GRHydro_test_tov_ppm.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_tov_ppm.par 2011-03-12 04:06:44 UTC (rev 222) @@ -74,8 +74,8 @@ IO::out_dir = "GRHydro_test_tov_ppm" IOBasic::outScalar_every = 1 -IOBasic::outScalar_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" +IOBasic::outScalar_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" IOASCII::out_format = ".10e" File [modified]: GRHydro_test_tov_ppm_ML.par Delta lines: +3 -3 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_ML.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_tov_ppm_ML.par 2011-03-12 04:06:44 UTC (rev 222) @@ -1,6 +1,6 @@ !DESC "TOV from TOVSolverC, ICN, Marquina, PPM, Polytrope" -ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce InitBase" +ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce InitBase LoopControl" TmunuBase::stress_energy_storage = yes TmunuBase::stress_energy_at_RHS = yes @@ -89,8 +89,8 @@ IO::out_dir = $parfile IOBasic::outScalar_every = 1 -IOBasic::outScalar_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" +IOBasic::outScalar_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" IOASCII::out_format = ".10e" File [modified]: GRHydro_test_tov_ppm_ML_disable_internal_excision.par Delta lines: +3 -3 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_ML_disable_internal_excision.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_tov_ppm_ML_disable_internal_excision.par 2011-03-12 04:06:44 UTC (rev 222) @@ -1,6 +1,6 @@ !DESC "TOV from TOVSolverC, ICN, Marquina, PPM, Polytrope" -ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce InitBase" +ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce InitBase LoopControl" TmunuBase::stress_energy_storage = yes TmunuBase::stress_energy_at_RHS = yes @@ -91,8 +91,8 @@ IO::out_dir = $parfile IOBasic::outScalar_every = 1 -IOBasic::outScalar_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" +IOBasic::outScalar_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" IOASCII::out_format = ".10e" File [modified]: GRHydro_test_tov_ppm_disable_internal_excision.par Delta lines: +2 -2 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_disable_internal_excision.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_tov_ppm_disable_internal_excision.par 2011-03-12 04:06:44 UTC (rev 222) @@ -76,8 +76,8 @@ IO::out_dir = $parfile IOBasic::outScalar_every = 1 -IOBasic::outScalar_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" +IOBasic::outScalar_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" IOASCII::out_format = ".10e" File [modified]: GRHydro_test_tov_ppm_no_trp.par Delta lines: +2 -2 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_no_trp.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_tov_ppm_no_trp.par 2011-03-12 04:06:44 UTC (rev 222) @@ -75,8 +75,8 @@ IO::out_fileinfo="none" IOBasic::outScalar_every = 1 -IOBasic::outScalar_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" +IOBasic::outScalar_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" IOASCII::out_format = ".10e" File [modified]: GRHydro_test_tov_ppm_no_trp_ML.par Delta lines: +3 -3 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_no_trp_ML.par 2011-03-12 03:16:20 UTC (rev 221) +++ trunk/test/GRHydro_test_tov_ppm_no_trp_ML.par 2011-03-12 04:06:44 UTC (rev 222) @@ -1,6 +1,6 @@ !DESC "TOV from TOVSolverC, ICN, Marquina, PPM, Polytrope" -ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce Tmunubase InitBase" +ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce Tmunubase InitBase LoopControl" TmunuBase::stress_energy_storage = yes TmunuBase::stress_energy_at_RHS = yes @@ -90,8 +90,8 @@ IO::out_fileinfo="none" IOBasic::outScalar_every = 1 -IOBasic::outScalar_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" +IOBasic::outScalar_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum" IOASCII::out1D_every = 1 -IOASCII::out1D_vars = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" +IOASCII::out1D_vars = "hydrobase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum" IOASCII::out_format = ".10e"

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits March 2011