Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 60)
by tanja.bode＠physics.gatech.edu 05 Apr '12

05 Apr '12

User: tbode Date: 2012/04/05 02:09 AM Added: /handson_sviz/ Makefile, handson_sviz.pdf, handson_sviz.tex, recompile /handson_sviz/extra_datafiles/ BH_diagnostics.ah1.gp, BH_diagnostics.ah2.gp /scientific_viz/ recompile, scientific_viz.pdf /scientific_viz/viz_images/ DiskMesh.png, ReductionHead.eps, ReductionHead.png, ReductionPlot.eps, ReductionPlot.png, heatmaps.png, sinc2d.png Modified: /scientific_viz/ Makefile, scientific_viz.tex Log: Getting there. File Changes: Directory: /handson_sviz/extra_datafiles/ ========================================= File [added]: BH_diagnostics.ah1.gp Delta lines: <Unavailable> File [added]: BH_diagnostics.ah2.gp Delta lines: <Unavailable> Directory: /handson_sviz/ ========================= File [added]: Makefile Delta lines: <Unavailable> File [added]: handson_sviz.pdf Delta lines: <Unavailable> File [added]: handson_sviz.tex Delta lines: <Unavailable> File [added]: recompile Delta lines: <Unavailable> Directory: /scientific_viz/ =========================== File [modified]: Makefile Delta lines: <Unavailable> File [added]: recompile Delta lines: <Unavailable> File [added]: scientific_viz.pdf Delta lines: <Unavailable> File [modified]: scientific_viz.tex Delta lines: <Unavailable> Directory: /scientific_viz/viz_images/ ====================================== File [added]: DiskMesh.png Delta lines: <Unavailable> File [added]: ReductionHead.eps Delta lines: <Unavailable> File [added]: ReductionHead.png Delta lines: <Unavailable> File [added]: ReductionPlot.eps Delta lines: <Unavailable> File [added]: ReductionPlot.png Delta lines: <Unavailable> File [added]: heatmaps.png Delta lines: <Unavailable> File [added]: sinc2d.png Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 59)
by roland.haas＠physics.gatech.edu 05 Apr '12

05 Apr '12

User: rhaas Date: 2012/04/05 12:36 AM Added: /handson_viz/ Makefile, add_isosurface.png, add_pseudocolor.png, add_reflection.png, annotation_menu.png, change_colour_scheme.png, disable_x_axis.png, final_plot_snapshot.png, handson_viz.tex, initial_plot_snapshot.png, isocontour_options.png, movie_frame_options.png, open_file.png, pseudocolour_options.png, reflection_options.png, render_frames_menu.png, rendering_options.png Log: add first version of visualization hands on this is really just a collection of images of where to click on not quite sure how else to guide them through VisIt and still be done on time. Directory Changes: Directory: /svn:mime-type/ ========================== + application/octet-stream File Changes: Directory: /handson_viz/ ======================== File [added]: Makefile Delta lines: +10 -0 =================================================================== --- handson_viz/Makefile (rev 0) +++ handson_viz/Makefile 2012-04-05 05:36:03 UTC (rev 59) @@ -0,0 +1,10 @@ +all: handson_viz.pdf + +%.pdf: %.tex + pdflatex $* + +clean: + /bin/bash -exec "rm -f handson_viz.{aux,log,nav,out,snm,pdf,toc} *~" + +.PHONY: all clean + File [added]: add_isosurface.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/add_isosurface.png ___________________________________________________________________ File [added]: add_pseudocolor.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/add_pseudocolor.png ___________________________________________________________________ File [added]: add_reflection.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/add_reflection.png ___________________________________________________________________ File [added]: annotation_menu.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/annotation_menu.png ___________________________________________________________________ File [added]: change_colour_scheme.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/change_colour_scheme.png ___________________________________________________________________ File [added]: disable_x_axis.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/disable_x_axis.png ___________________________________________________________________ File [added]: final_plot_snapshot.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/final_plot_snapshot.png ___________________________________________________________________ File [added]: handson_viz.tex Delta lines: +245 -0 =================================================================== --- handson_viz/handson_viz.tex (rev 0) +++ handson_viz/handson_viz.tex 2012-04-05 05:36:03 UTC (rev 59) @@ -0,0 +1,245 @@ +\input{../preamble.tex} +\usepackage{ragged2e} +\usepackage{amsmath} +\usepackage{listings} +\newcommand{\Heaviside}{\operatorname{\Theta}} + +\subtitle[Visualization hands-on]{{\large Einstein Toolkit +Workshop}\\*[0.3em]Hands-on exercise: visualizing density oscillations and +gravitational waves} +\author{Roland Haas} +\date{Apr 05 2012} + +\usecolortheme[RGB={200,200,200}]{structure} + +\begin{document} + +\begin{frame} +\titlepage +\end{frame} + +\begin{frame} +\frametitle{Installing software} +To turn the rendered frames into a movie we will use \texttt{ffmpeg} which +unfortunately did not make it into the virtual machine. To install it, enter +on a command prompt: +\begin{alltt} + sudo apt-get install ffmpeg +\end{alltt} +Enter \texttt{etk-user} (ie. your password) when asked for the password. +\end{frame} + +\begin{frame} +\frametitle{TOV files to visualize} +You can either use the output files from the \texttt{fast\_tov.par} run you +started yesterday, or download a (slightly different) dataset from +\url{http://www.physics.gatech.edu/~rhaas3/Workshop\_Spring\_2012/rho.rl4.h5.xz} +(about 4MB). Copy the file to a folder of your choice and decompress it via +\begin{alltt} + xz -d -f --keep path-to-folder/rho.rl4.h5.xz +\end{alltt} + +The former file will work but will be a bit boring. The later one I tweaked a +bit to get nice large oscillations. +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} + \item type \texttt{\$HOME/VisIt/bin/visit} at a command prompt to start + visit + \item You can try enabling 3D acceleration, but at least on Intel graphics + chips this tends to make VisIt fail at startup. + \item open the output file of your choice +\end{itemize} +\begin{center} +\includegraphics[height=48mm]{open_file} +\end{center} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item add a pseudocolor plot from the "Plots" drop-down button +\item we will display a number of isocontour lines in the star to show its oscillations +\item we will use the \texttt{IsoSurface} operator to compute the contour surfaces +\item \texttt{VisIt} also offers a \texttt{Contour} plot which however is limited in the transparency options compared to the \texttt{IsoSurface}, \texttt{PseudoColor} combination we use +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{add_pseudocolor} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item add \texttt{IsoSurface} operator from \texttt{Operators/Slicing/Isosurface} +\item add \texttt{Reflect} operator from \texttt{Operators/Transforms/Reflect}, to populate the ``missing'' octants +\end{itemize} +\begin{columns} +\column{60mm} +\includegraphics[width=60mm]{add_reflection} +\column{60mm} +\includegraphics[width=60mm]{add_isosurface} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item double-click on the \texttt{Pseudocolor} entry in the list of plots to display its options page. Your might have to expand the list using the triangle shaped button on the left first +\item add transparency by reducing the opacity to about 25\% and set limits of $[10^{-8},10^{-3}]$ for the plot range. We will use the same range when creating isocontours in the next step +\item you have to click on \texttt{Apply} then on \texttt{Dismiss} to apply the new settings +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{pseudocolour_options} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item double click on the \texttt{Reflect} operator, change to \texttt{3D} and activate all octants except the original ($+X+Y+Z$) one and the $+X-Y+Z$ one, change the reflection point to the origin. +\end{itemize} +\begin{center} +\includegraphics[height=60mm]{reflection_options} +\end{center} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item double click on the \texttt{IsoSurface} operator and select the range used for the \texttt{Pseudocolor} plot, use $6$ levels and leave the scaling at \texttt{linear} +\end{itemize} +\begin{center} +\includegraphics[width=60mm]{isocontour_options} +\end{center} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item click on \texttt{Draw} to see how your plot looks like +\item rotate to an orientation that looks nice, by grabbing the plot with the mouse +\end{itemize} +\begin{center} +\includegraphics[width=60mm]{initial_plot_snapshot} +\end{center} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item next we will change \texttt{VisIt}'s rendering options to add shadows and change the background to something more attractive +\item unfortunately this means we will eventually loose the transparancy (at least in the virtual machine) +\item first though open the \texttt{Annotation} control page +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{annotation_menu} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item in the \texttt{X-Axis} subtab of the \texttt{3D} tab, disable the axis title +\item repeat for the other two axes +\item in the \texttt{Colors} tab, change the background to black and the foreground to white +\item add a gradient to the background, with a dark gray as the second colour +\end{itemize} +\begin{center} +\includegraphics[width=58mm]{disable_x_axis} +\hfill +\includegraphics[width=58mm]{change_colour_scheme} +\end{center} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item open \texttt{VisIt}'s rendering options from the \texttt{Options/Rendering...} menu +\item on the \texttt{Advanced} tab, enable \texttt{scalable rendering}, \texttt{Shadows} and \texttt{Depth Cueing} +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{rendering_options} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item click on \texttt{Draw} to see how your plot looks like +\end{itemize} +\begin{center} +\includegraphics[width=60mm]{final_plot_snapshot} +\end{center} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item in this penultimate step we will let \texttt{VisIt} rendering frames of a small movie showing the oscillations +\item to start the movie wizard, select \texttt{Save movie...} from the \texttt{File} menu +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{render_frames_menu} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item in this penultimate step we will let \texttt{VisIt} rendering frames of a small movie showing the oscillations +\item to start the movie wizard, select \texttt{Save movie...} from the \texttt{File} menu +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{render_frames_menu} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{columns} +\column{60mm} +\begin{itemize} +\item Finally, inside the wizard, change the frame size to approximately $640\times640$ pixels (multiples of $16$ are best) +\item choose \texttt{png} output +\item add the movie format to the list on the right using the \texttt{->} button +\item click on next and accept the defaults for all other steps, maybe choosing a different output directory for the frames +\end{itemize} +\column{60mm} +\includegraphics[width=60mm]{movie_frame_options} +\end{columns} +\end{frame} + +\begin{frame} +\frametitle{Visualizing density oscillations} +\begin{itemize} +\item once \texttt{VisIt} is done, which will take some minutes, we can render the frames into an actual movie +\item for this, \texttt{cd} into the folder containing the frames and enter +\begin{alltt} +ffmpeg -loop\_input -vframes 200 -qscale 1 -b 1000 -i movie\%04d.png movie.avi +\end{alltt} +\begin{itemize} +\item \texttt{-loop\_input} makes \texttt{ffmpeg} wrap around when looking for frames +\item \texttt{-vframes 200} is the number of frames to include in the movie, we set this to twice as many as we have actual frames, ie. the movie will loop twice +\item \texttt{-qscale 1} sets the quality scale for the encoder. $1$ is best quality. +\item \texttt{-b 1000} asks for a bitrate of $1000$ kilobits per second +\item \texttt{-i movie\%04d.png} tells \texttt{ffmpeg} the names of the input files, it accepts \texttt{printf} format specifiers to construct a file name based on the frame number. \texttt{ffmpeg} stops encoding once it generates a name for which no or no valid file is found. +\end{itemize} +\item play the movie using \texttt{ffplay} +\end{itemize} +\end{frame} +\end{document} File [added]: initial_plot_snapshot.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/initial_plot_snapshot.png ___________________________________________________________________ File [added]: isocontour_options.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/isocontour_options.png ___________________________________________________________________ File [added]: movie_frame_options.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/movie_frame_options.png ___________________________________________________________________ File [added]: open_file.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/open_file.png ___________________________________________________________________ File [added]: pseudocolour_options.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/pseudocolour_options.png ___________________________________________________________________ File [added]: reflection_options.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/reflection_options.png ___________________________________________________________________ File [added]: render_frames_menu.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/render_frames_menu.png ___________________________________________________________________ File [added]: rendering_options.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: handson_viz/rendering_options.png ___________________________________________________________________

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/example_wavetoy/ (Rev. 58)
by roland.haas＠physics.gatech.edu 04 Apr '12

04 Apr '12

User: rhaas Date: 2012/04/04 03:45 PM Modified: /example_wavetoy/ example_wavetoy.pdf, example_wavetoy.tex Log: fix drawings to reflect what the code actually does wrt boundaries gray out very old timelevels highlight only initial data on current level File Changes: Directory: /example_wavetoy/ ============================ File [modified]: example_wavetoy.pdf Delta lines: +0 -0 =================================================================== (Binary files differ) File [modified]: example_wavetoy.tex Delta lines: +11 -9 =================================================================== --- example_wavetoy/example_wavetoy.tex 2012-04-04 18:43:17 UTC (rev 57) +++ example_wavetoy/example_wavetoy.tex 2012-04-04 20:45:32 UTC (rev 58) @@ -260,8 +260,8 @@ % Step 7: Layers rotation \uncover<7->{ \scope[very thick,color=green,xstep=0.6cm,ystep=1cm] - \draw[xshift=0cm,yshift=0cm] (0cm,-1pt) grid (3.0cm,1pt); - \draw[xshift=2.45cm,yshift=0.05cm] (0cm,-1pt) grid (3.0cm,1pt); + \draw[xshift=0cm,yshift=0cm,color=gray] (0cm,-1pt) grid (3.0cm,1pt); + \draw[xshift=2.45cm,yshift=0.05cm,color=gray] (0cm,-1pt) grid (3.0cm,1pt); \draw[xshift=0cm,yshift=0.4cm] (0cm,-1pt) grid (3.0cm,1pt); \draw[xshift=2.45cm,yshift=0.45cm] (0cm,-1pt) grid (3.0cm,1pt); \draw[xshift=0cm,yshift=0.8cm] (0cm,-1pt) grid (3.0cm,1pt); @@ -325,16 +325,16 @@ \scope[very thick,color=green,xstep=0.6cm,ystep=1cm] \draw[xshift=0.6cm,yshift=1.2cm] (0cm,-1pt) grid (2.4cm,1pt); \draw[xshift=2.45cm,yshift=1.25cm] (0cm,-1pt) grid (2.4cm,1pt); - \draw[->,color=red] (0.0,0.8) -- (0.0,1.2); - \draw[->,color=red] (5.45,0.85) -- (5.45,1.25); + \draw[->,color=red] (-0.4,1.2) -- (0.0,1.2); + \draw[->,color=red] (5.85,1.25) -- (5.45,1.25); \endscope }% % Step 11: Output data \uncover<11->{ \scope[very thick,color=green,xstep=0.6cm,ystep=1cm] - \draw[xshift=0cm,yshift=0.8cm] (0cm,-1pt) grid (3.0cm,1pt); - \draw[xshift=2.45cm,yshift=0.85cm] (0cm,-1pt) grid (3.0cm,1pt); + \draw[xshift=0cm,yshift=1.2cm] (0cm,-1pt) grid (3.0cm,1pt); + \draw[xshift=2.45cm,yshift=1.25cm] (0cm,-1pt) grid (3.0cm,1pt); \endscope } @@ -346,6 +346,8 @@ \draw[xshift=2.45cm,yshift=1.65cm,color=black] (0cm,-1pt) grid (3.0cm,1pt); \draw[xshift=0cm,yshift=1.2cm] (0cm,-1pt) grid (3.0cm,1pt); \draw[xshift=2.45cm,yshift=1.25cm] (0cm,-1pt) grid (3.0cm,1pt); + \draw[xshift=0cm,yshift=0.4cm,color=gray] (0cm,-1pt) grid (3.0cm,1pt); + \draw[xshift=2.45cm,yshift=0.45cm,color=gray] (0cm,-1pt) grid (3.0cm,1pt); \endscope \scope[color=blue] \draw (-0.5,0.8) node[anchor=east] {\small phi\_p\_p}; @@ -596,15 +598,15 @@ CCTK\_REAL Z = z[vindex]; CCTK\_REAL R = sqrt(X*X + Y*Y + Z*Z); - phi[vindex] = amplitude*exp( - SQ( (R - radius) / sigma ) ); + \textcolor{blue}{phi[vindex] = amplitude*exp( - SQ( (R - radius) / sigma ) );} if (R == 0.0) \{ phi\_p[vindex] = amplitude*(1.0 - 2.0*dt*dt/sigma)*exp(-dt*dt/sigma); \} else \{ - \textcolor{blue}{phi\_p[vindex] = amplitude/2.0*(R-dt)/R* + phi\_p[vindex] = amplitude/2.0*(R-dt)/R* exp( - SQ( (R - radius - dt)/ sigma ) ) + amplitude/2.0*(R+dt)/R* - exp( - SQ( (R - radius + dt)/ sigma ) );} + exp( - SQ( (R - radius + dt)/ sigma ) ); \textcolor{Green}{\} \} \}}

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 04 Apr '12

04 Apr '12

changeset: 3508:ded3e3f657fe user: Erik Schnetter <schnetter(a)gmail.com> date: Wed Apr 04 15:55:18 2012 -0400 summary: Carpet: Time OutputGH operations changeset: 3509:752bc38fb8a5 user: Erik Schnetter <schnetter(a)gmail.com> date: Wed Apr 04 15:55:57 2012 -0400 summary: CarpetLib: Improve comment of dh data structures changeset: 3510:6e7df7a98136 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Wed Apr 04 16:09:23 2012 -0400 summary: CarpetRegrid2: Ensure that the refinement hierarchy is symmetry about the origin, if desired

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 04 Apr '12

04 Apr '12

changeset: 3507:e970c8c56a63 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Wed Apr 04 15:54:40 2012 -0400 summary: Carpet: Do not count ghost zones when calculating speed

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/supercomputing/ (Rev. 57)
by jfaber＠einsteintoolkit.org 04 Apr '12

04 Apr '12

User: jfaber Date: 2012/04/04 01:43 PM Added: /supercomputing/ supercomputing.aux, supercomputing.log, supercomputing.nav, supercomputing.out, supercomputing.pdf, supercomputing.snm, supercomputing.tex, supercomputing.toc Log: Contributions to the notes on supercomputing/Xsede resources File Changes: Directory: /supercomputing/ =========================== File [added]: supercomputing.aux Delta lines: <Unavailable> File [added]: supercomputing.log Delta lines: <Unavailable> File [added]: supercomputing.nav Delta lines: <Unavailable> File [added]: supercomputing.out Delta lines: <Unavailable> File [added]: supercomputing.pdf Delta lines: <Unavailable> File [added]: supercomputing.snm Delta lines: <Unavailable> File [added]: supercomputing.tex Delta lines: <Unavailable> File [added]: supercomputing.toc Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/example_wavetoy/Cactus/arrangements/Introduction/WaveToy/ (Rev. 56)
by roland.haas＠physics.gatech.edu 04 Apr '12

04 Apr '12

User: rhaas Date: 2012/04/04 01:22 PM Modified: /example_wavetoy/Cactus/arrangements/Introduction/WaveToy/ interface.ccl, schedule.ccl /example_wavetoy/Cactus/arrangements/Introduction/WaveToy/src/ WaveToy.c Log: add routines to compute wave energy to wavetoy example File Changes: Directory: /example_wavetoy/Cactus/arrangements/Introduction/WaveToy/ ===================================================================== File [modified]: interface.ccl Delta lines: +2 -0 =================================================================== --- example_wavetoy/Cactus/arrangements/Introduction/WaveToy/interface.ccl 2012-04-04 18:04:36 UTC (rev 55) +++ example_wavetoy/Cactus/arrangements/Introduction/WaveToy/interface.ccl 2012-04-04 18:22:04 UTC (rev 56) @@ -8,7 +8,9 @@ phi } "The evolved scalar field" +CCTK_REAL energy TYPE=gf TIMELEVELS=3 "Energy in wave" + CCTK_INT FUNCTION Boundary_SelectVarForBC(CCTK_POINTER_TO_CONST IN GH, \ CCTK_INT IN faces, CCTK_INT IN boundary_width, CCTK_INT IN table_handle, \ CCTK_STRING IN var_name, CCTK_STRING IN bc_name) File [modified]: schedule.ccl Delta lines: +6 -1 =================================================================== --- example_wavetoy/Cactus/arrangements/Introduction/WaveToy/schedule.ccl 2012-04-04 18:04:36 UTC (rev 55) +++ example_wavetoy/Cactus/arrangements/Introduction/WaveToy/schedule.ccl 2012-04-04 18:22:04 UTC (rev 56) @@ -1,4 +1,4 @@ -STORAGE: scalarevolve[3] +STORAGE: scalarevolve[3] energy[3] SCHEDULE WaveToy_InitialData AT CCTK_INITIAL { @@ -11,6 +11,11 @@ SYNC: scalarevolve } "Evolution of 3D wave equation" +SCHEDULE WaveToy_Energy AT CCTK_EVOL AFTER WaveToy_Evolution +{ + LANG: C +} "Compute wave energy" + SCHEDULE WaveToy_Boundaries AT CCTK_EVOL AFTER WaveToy_Evolution { LANG: C Directory: /example_wavetoy/Cactus/arrangements/Introduction/WaveToy/src/ ========================================================================= File [modified]: WaveToy.c Delta lines: +19 -0 =================================================================== --- example_wavetoy/Cactus/arrangements/Introduction/WaveToy/src/WaveToy.c 2012-04-04 18:04:36 UTC (rev 55) +++ example_wavetoy/Cactus/arrangements/Introduction/WaveToy/src/WaveToy.c 2012-04-04 18:22:04 UTC (rev 56) @@ -70,7 +70,26 @@ } } +/* + * Routine to compute energy + */ +void WaveToy_Energy(CCTK_ARGUMENTS) { + DECLARE_CCTK_ARGUMENTS; + /* Energy density is just square of wave */ + for (int k=1; k<cctk_lsh[2]-1; k++) { + for (int j=1; j<cctk_lsh[1]-1; j++) { + for (int i=1; i<cctk_lsh[0]-1; i++) { + energy[CCTK_GFINDEX3D(cctkGH,i,j,k)] = + SQ(phi[CCTK_GFINDEX3D(cctkGH,i,j,k)]) * + CCTK_DELTA_SPACE(0) * CCTK_DELTA_SPACE(1) * + CCTK_DELTA_SPACE(2); + } + } + } +} + + /* * Select boundary condition */

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/example_tov/ (Rev. 55)
by knarf＠cct.lsu.edu 04 Apr '12

04 Apr '12

User: knarf Date: 2012/04/04 01:04 PM Added: /example_tov/ fast_tov.par Log: fast tov parameter file File Changes: Directory: /example_tov/ ======================== File [added]: fast_tov.par Delta lines: +221 -0 =================================================================== --- example_tov/fast_tov.par (rev 0) +++ example_tov/fast_tov.par 2012-04-04 18:04:36 UTC (rev 55) @@ -0,0 +1,221 @@ +## mclachlan tov_static +ActiveThorns = "Time MoL" +ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal" +ActiveThorns = "SymBase ADMBase TmunuBase HydroBase InitBase ADMCoupling ADMMacros" +ActiveThorns = "IOUtil Formaline" +ActiveThorns = "SpaceMask CoordGauge Constants LocalReduce aeilocalinterp LoopControl" +ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetRegrid2 CarpetInterp" +ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic" + +ActiveThorns = "ADMConstraints NaNChecker" + +# grid parameters +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" +CoordBase::xmin = 0.0 +CoordBase::ymin = 0.0 +CoordBase::zmin = 0.0 +CoordBase::xmax = 15.0 +CoordBase::ymax = 15.0 +CoordBase::zmax = 15.0 +CoordBase::dx = 1 +CoordBase::dy = 1 +CoordBase::dz = 1 + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + + + + +ActiveThorns = "ReflectionSymmetry" + +ReflectionSymmetry::reflection_x = "yes" +ReflectionSymmetry::reflection_y = "yes" +ReflectionSymmetry::reflection_z = "yes" +ReflectionSymmetry::avoid_origin_x = "no" +ReflectionSymmetry::avoid_origin_y = "no" +ReflectionSymmetry::avoid_origin_z = "no" + + + +TmunuBase::stress_energy_storage = yes +TmunuBase::stress_energy_at_RHS = yes +TmunuBase::timelevels = 1 +TmunuBase::prolongation_type = none + +HydroBase::timelevels = 3 + +ADMMacros::spatial_order = 4 + +ADMBase::metric_type = "physical" + +ADMConstraints::bound = "static" +ADMConstraints::constraints_timelevels = 3 +ADMConstraints::constraints_persist = yes + +SpaceMask::use_mask = "yes" + +Cactus::terminate = "time" +Cactus::cctk_final_time = 1000 + +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = no +Carpet::use_buffer_zones = "yes" + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" +Carpet::poison_value = 113 + +Carpet::init_3_timelevels = no +Carpet::init_fill_timelevels = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 + +# system specific Carpet paramters +Carpet::max_refinement_levels = 10 +driver::ghost_size = 3 +Carpet::prolongation_order_space = 3 +Carpet::prolongation_order_time = 2 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::num_centres = 1 + +CarpetRegrid2::num_levels_1 = 1 +CarpetRegrid2::radius_1[1] = 15.0 + + +time::dtfac = 0.25 +MoL::ODE_Method = "rk4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +# check all physical variables for NaNs +NaNChecker::check_every = 1 +NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" + +## Lapse Condition: \partial_t alpha = - alpha K +## Shift Condition: \partial_t beta^i = 0 + +# Hydro paramters +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro" + +HydroBase::evolution_method = "GRHydro" + +GRHydro::riemann_solver = "Marquina" +GRHydro::GRHydro_eos_type = "Polytype" +GRHydro::GRHydro_eos_table = "2D_Polytrope" +GRHydro::recon_method = "ppm" +GRHydro::GRHydro_stencil = 3 +GRHydro::bound = "none" +GRHydro::rho_abs_min = 1.e-10 +#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau +#GRHydro::GRHydro = 10 # gij(6), alpha, beta(3) + + +ActiveThorns = "GenericFD NewRad" +ActiveThorns = "ML_BSSN ML_BSSN_Helper" +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::shift_evolution_method = "ML_BSSN" +ADMBase::dtlapse_evolution_method= "ML_BSSN" +ADMBase::dtshift_evolution_method= "ML_BSSN" + +TmunuBase::support_old_CalcTmunu_mechanism = "no" + +ML_BSSN::timelevels = 3 + +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 1.0 # 1+log +ML_BSSN::LapseACoeff = 1.0 +ML_BSSN::ShiftBCoeff = 1.0 +ML_BSSN::ShiftGammaCoeff = 0.0 +ML_BSSN::AlphaDriver = 0.0 +ML_BSSN::BetaDriver = 0.0 +ML_BSSN::LapseAdvectionCoeff = 0.0 +ML_BSSN::ShiftAdvectionCoeff = 0.0 + +ML_BSSN::MinimumLapse = 1.0e-8 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" + +ML_BSSN::ML_log_confac_bound = "none" +ML_BSSN::ML_metric_bound = "none" +ML_BSSN::ML_Gamma_bound = "none" +ML_BSSN::ML_trace_curv_bound = "none" +ML_BSSN::ML_curv_bound = "none" +ML_BSSN::ML_lapse_bound = "none" +ML_BSSN::ML_dtlapse_bound = "none" +ML_BSSN::ML_shift_bound = "none" +ML_BSSN::ML_dtshift_bound = "none" + + + +# init parameters +InitBase::initial_data_setup_method = "init_some_levels" + +ActiveThorns = "TOVSolver" + +ADMBase::initial_data = "tov" +ADMBase::initial_lapse = "tov" +ADMBase::initial_shift = "tov" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +TOVSolver::TOV_Gamma[0] = 2.0 +TOVSolver::TOV_K[0] = 100.0 + + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = "HydroBase::rho ADMBase::lapse" + +IO::out_dir = $parfile + +IOScalar::outScalar_every = 32 +IOScalar::one_file_per_group = yes +IOScalar::outScalar_vars = " + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ADMConstraints::ham + ADMConstraints::momentum + +" + +IOASCII::out1D_every = 256 +IOASCII::one_file_per_group = yes +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no +IOASCII::out3D_outer_ghosts = no +IOASCII::out1D_vars = " + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ADMConstraints::ham + ADMConstraints::momentum + +"

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 54)
by knarf＠cct.lsu.edu 04 Apr '12

04 Apr '12

User: knarf Date: 2012/04/04 12:59 PM Added: /example_tov/ 2_rho.maximum.asc.xz, 4_rho.maximum.asc.xz Log: some precomputed tov data Directory Changes: Directory: /svn:mime-type/ ========================== + application/octet-stream File Changes: Directory: /example_tov/ ======================== File [added]: 2_rho.maximum.asc.xz Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: example_tov/2_rho.maximum.asc.xz ___________________________________________________________________ File [added]: 4_rho.maximum.asc.xz Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: example_tov/4_rho.maximum.asc.xz ___________________________________________________________________

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 53)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 12:51 PM Modified: /numerical_relativity/ numerical_relativity.tex Log: Outline. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +9 -0 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 17:35:10 UTC (rev 52) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 17:51:11 UTC (rev 53) @@ -23,6 +23,15 @@ \item Finite differencing. \item The method of lines. \item McLachlan. + \item GRMHD - Valencia Formulation. + \item Conservation Laws. + \item Finite Volume Methods. + \item Riemann Problem. + \item Example: Burgers Equation. + \item Approximate Riemann Solvers. + \item Slope Limiters + \item Algorithm Overview. + \item Shocktube Test. \end{itemize} }

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits