[Commits] [svn:einsteintoolkit] Workshop_Spring_2012/handson_analysis/ (Rev. 26)

[roland.haas at physics.gatech.edu]

User: rhaas
Date: 2012/04/02 04:39 PM

Modified:
 /handson_analysis/
  Makefile, handson_analysis.tex

Log:
 add slides describing coretracker code units

File Changes:

Directory: /handson_analysis/
=============================

File [modified]: Makefile
Delta lines: +1 -1
===================================================================
--- handson_analysis/Makefile	2012-04-02 15:24:59 UTC (rev 25)
+++ handson_analysis/Makefile	2012-04-02 21:39:51 UTC (rev 26)
@@ -1,4 +1,4 @@
-all: handson_aanalysis.pdf
+all: handson_analysis.pdf
 
 %.pdf: %.tex
 	pdflatex $*

File [modified]: handson_analysis.tex
Delta lines: +58 -7
===================================================================
--- handson_analysis/handson_analysis.tex	2012-04-02 15:24:59 UTC (rev 25)
+++ handson_analysis/handson_analysis.tex	2012-04-02 21:39:51 UTC (rev 26)
@@ -1,6 +1,7 @@
 \input{../preamble.tex}
 \usepackage{ragged2e}
 \usepackage{amsmath}
+\usepackage{listings}
 \newcommand{\Heaviside}{\operatorname{\Theta}}
 
 \subtitle[Analysis hands-on]{{\large Einstein Toolkit
@@ -12,9 +13,12 @@
 
 \begin{document}
 
-\frame{\titlepage}
+\begin{frame}
+\titlepage
+\end{frame}
 
-\frame{\frametitle{Problem description}
+\begin{frame}
+\frametitle{Problem description}
 \includegraphics[width=\linewidth]{inspiral_image}
 {\\\scriptsize image courtesy of NASA}
 \begin{block}{Motivation}
@@ -24,9 +28,10 @@
         \item need to track location of star for moving boxed
     \end{itemize}
 \end{block}
-}
+\end{frame}
 
-\frame{\frametitle{Problem description}
+\begin{frame}
+\frametitle{Problem description}
 \begin{block}{Requirements}
     \begin{itemize}
         \item be able to pick out a single neutron star, among several
@@ -44,9 +49,10 @@
             influence of outflowing tails
     \end{itemize}
 \end{block}
-}
+\end{frame}
 
-\frame[containsverbatim]{\frametitle{Recipe in formulae}
+\begin{frame}[containsverbatim]
+\frametitle{Recipe in formulae}
 \begin{block}{What to compute}
     \begin{align*}
         M  &= \int \Heaviside(\rho - \rho_{\text{min}})
@@ -65,6 +71,51 @@
         \item $\Heaviside$ $\rightarrow$ \verb|if(...)|
     \end{itemize}
 \end{block}
-}
+\end{frame}
 
+\begin{frame}
+\frametitle{Entities to define}
+\begin{itemize}
+    \item<1-> \texttt{param.ccl}
+        \begin{alltt}\color{blue}
+            CCTK\_REAL density\_threshold\\
+            CCTK\_REAL core\_radius\\
+            CCTK\_REAL inital\_core\_\{x,y,z\}
+        \end{alltt}
+    \item<2-> \texttt{interface.ccl}
+        \begin{alltt}\color{blue}
+            CCTK\_REAL core\_density\_gf type=gf\\
+            \hspace{1cm}core\_density\_\{x,y,z\}\\
+            \hspace{1cm}core\_density\\
+            CCTK\_REAL core\_location type=scalar\\
+            \hspace{1cm}core\_[xyz]
+        \end{alltt}
+\end{itemize}
+\end{frame}
+
+\begin{frame}[fragile]
+\frametitle{Algorithm outline}
+\begin{itemize}
+    \item<1-> \texttt{CCTK\_INITIAL} copy initial location into evolution variables
+        \begin{alltt}\color{blue}
+*core\_x = initial\_core\_x\
+        \end{alltt}
+\item<2-> \texttt{MoL\_PseudoEvolution} evaluate integrand on current time-slice
+\begin{alltt}\color{blue}
+if(distance2 < SQR(core\_radius)
+  rho[idx] >= threshold\_density) \{
+    core\_density\_x[idx] = x[idx] * rho[idx];
+    core\_density[idx]   = rho[idx];
+\} else \{
+    core\_density\_x[idx] = 0.;
+    core\_density[idx]   = 0.;
+\}
+\end{alltt}
+\item<3-> \texttt{CCTK\_ANALYSIS} compute \texttt{sum} reduction on integrand
+\begin{alltt}\color{blue}
+CCTK\_Reduce(cctkGH, -1, handle\_sum, 1, CCTK\_VARIABLE\_REAL,
+             core\_x, 1, varindex);
+\end{alltt}
+\end{itemize}
+\end{frame}
 \end{document}