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}