Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] incoming/ (Rev. 37)
by reisswig＠tapir.caltech.edu 13 Jun '11

13 Jun '11

User: reisswig Date: 2011/06/13 04:44 PM Added: /PITTNullCode/WorldTube/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/WorldTube/doc/ documentation.tex /PITTNullCode/WorldTube/src/ areal_sphere.cc, decompose.cc, make.code.defn, paramcheck.cc, register.cc, register.h Removed: /WorldTube/ Log: Place thorn for metric decomposition in PITTNullCode repository. File Changes: Directory: /PITTNullCode/WorldTube/ =================================== File [added]: README Delta lines: +12 -0 =================================================================== --- PITTNullCode/WorldTube/README (rev 0) +++ PITTNullCode/WorldTube/README 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,12 @@ +CVS info : $Header:$ + +Cactus Code Thorn WorldTube +Thorn Author(s) : Christian Reisswig <reisswig(a)aei.mpg.de> +Thorn Maintainer(s) : Christian Reisswig <reisswig(a)aei.mpg.de> +-------------------------------------------------------------------------- + +Purpose of the thorn: + +This thorn is ment as an interface between Cauchy and characteristic codes. + +It can also just output metric components and find spheres of constant areal radius. \ No newline at end of file File [added]: configuration.ccl Delta lines: +1 -0 =================================================================== --- PITTNullCode/WorldTube/configuration.ccl (rev 0) +++ PITTNullCode/WorldTube/configuration.ccl 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1 @@ +REQUIRES THORNS: ADMBase ADMDerivatives File [added]: interface.ccl Delta lines: +65 -0 =================================================================== --- PITTNullCode/WorldTube/interface.ccl (rev 0) +++ PITTNullCode/WorldTube/interface.ccl 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,65 @@ +# Interface definition for thorn WorldTube +# $Header:$ + +implements: WorldTube +inherits: SphericalSlice ADMBase + + +uses include: slices.hh + + + +# Register a Cactus variable for being sliced and return registration number +CCTK_INT FUNCTION \ + SphericalSlice_RegisterVariable \ + (CCTK_POINTER_TO_CONST IN varname, \ + CCTK_INT IN sn, \ + CCTK_INT IN timelevels, \ + CCTK_POINTER_TO_CONST IN distrib_method) +REQUIRES FUNCTION SphericalSlice_RegisterVariable + + +# Return an array of all values of scalar product with all spin-weighted spherical harmonics sYlm up to lmax +CCTK_INT FUNCTION \ + SphericalSlice_ContractVariableWithAllsYlm \ + (CCTK_INT IN varno, \ + CCTK_INT IN timelevel, \ + CCTK_INT IN s, \ + CCTK_INT IN lmin, \ + CCTK_INT IN lmax, \ + CCTK_COMPLEX OUT ARRAY coeffs) +REQUIRES FUNCTION SphericalSlice_ContractVariableWithAllsYlm \ + + +# Synchronzie slice, i.e. interpolate a whole set of variables on the same slice +CCTK_INT FUNCTION \ + SphericalSlice_CollectiveSyncVariables \ + (CCTK_POINTER_TO_CONST IN cctkGH_, \ + CCTK_POINTER_TO_CONST IN varno, \ + CCTK_INT IN number_of_vars) +REQUIRES FUNCTION SphericalSlice_CollectiveSyncVariables \ + + + +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) +CCTK_INT FUNCTION Boundary_SelectGroupForBC(CCTK_POINTER_TO_CONST IN GH, \ + CCTK_INT IN faces, CCTK_INT IN boundary_width, CCTK_INT IN table_handle, \ + CCTK_STRING IN group_name, CCTK_STRING IN bc_name) + +REQUIRES FUNCTION Boundary_SelectVarForBC +REQUIRES FUNCTION Boundary_SelectGroupForBC + + +PUBLIC: + +COMPLEX extracted_vars[50] TYPE=ARRAY TIMELEVELS=1 DIM=2 SIZE=((lmax+1)*(lmax+1)),ntubes DISTRIB=CONSTANT +{ + extraction_vars +} "Stores extracted lm-modes for each metric component" + + + + + File [added]: param.ccl Delta lines: +67 -0 =================================================================== --- PITTNullCode/WorldTube/param.ccl (rev 0) +++ PITTNullCode/WorldTube/param.ccl 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,67 @@ +# Parameter definitions for thorn WorldTube +# $Header:$ + + +KEYWORD boundary_behavior "The type of boundary" +{ + "none" :: "don't extract/inject anything" + "extract metric" :: "extract metric components, lapse and shift" + "CCE" :: "extract metric components, lapse, shift and derivatives for CCE" + "CCE_cart" :: "extract metric components, lapse, shift and (cartesian) derivatives for CCE" + "CCM" :: "to be implemented" +} "CCE" + + +INT lmax "Maximum extracted l-mode" +{ + 0:* :: "any positive number" +} 8 + + +CCTK_INT ntubes "The number of worldtubes to be placed" +{ + 0:42 :: "Any number between 0 and 42" +} 1 + + + +CCTK_INT which_slice_to_take[42] "Which SphericalSlice to use" +{ + 0:* :: "Any positive number" +} 0 + + +BOOLEAN seek_const_areal_radius[42] "Set slice radius to a surface of constant areal radius before extracting/injecting" +{ +} no + + +REAL tolerance "Tolerance for finding areal radius" +{ + 0:* :: "Any positive number" +} 1e-6 + + +INT max_it "Maximum number of iterations for finding areal radius" +{ + 1:* :: "any positive number greater than 1" +} 10 + + +CCTK_INT radial_deriv_order "Order of radial derivatives" +{ + -1 :: "Use whatever is specified by SummationByParts" + 2:8 :: "Use this order" +} -1 + + + +SHARES: SphericalSlice +USES KEYWORD type[42] +USES REAL radius[42] +USES INT nslices AS ss_nslices +USES INT ntheta +USES INT nphi +USES INT nghostzones + + File [added]: schedule.ccl Delta lines: +45 -0 =================================================================== --- PITTNullCode/WorldTube/schedule.ccl (rev 0) +++ PITTNullCode/WorldTube/schedule.ccl 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,45 @@ +# Schedule definitions for thorn WorldTube +# $Header:$ + + +STORAGE: extracted_vars + + +SCHEDULE WorldTube_ParamCheck AT PARAMCHECK +{ + LANG: C + OPTIONS: GLOBAL +} "Check parameters" + + +SCHEDULE WorldTube_RegisterSlices AT CCTK_INITIAL AFTER SphericalSlice_HasBeenSet +{ + LANG: C + OPTIONS: GLOBAL +} "Register variables for extraction" + + +SCHEDULE WorldTube_RegisterSlices AT CCTK_POST_RECOVER_VARIABLES AFTER SphericalSlice_HasBeenSet +{ + LANG: C + OPTIONS: GLOBAL +} "Register variables for extraction" + + + +SCHEDULE WorldTube_ArealSphere AT CCTK_POSTSTEP AFTER SphericalSlice_HasBeenSet +{ + LANG: C + OPTIONS: GLOBAL +} "Calculate surfaces of constant areal radius if selected" + + + +SCHEDULE WorldTube_Decompose AT CCTK_ANALYSIS AFTER (SphericalSlice_HasBeenSet, WorldTube_ArealSphere) +{ + LANG: C + OPTIONS: GLOBAL + TRIGGERS: extracted_vars +} "Contract sliced variables with sYlm's" + + Directory: /PITTNullCode/WorldTube/doc/ ======================================= File [added]: documentation.tex Delta lines: +144 -0 =================================================================== --- PITTNullCode/WorldTube/doc/documentation.tex (rev 0) +++ PITTNullCode/WorldTube/doc/documentation.tex 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,144 @@ +% *======================================================================* +% Cactus Thorn template for ThornGuide documentation +% Author: Ian Kelley +% Date: Sun Jun 02, 2002 +% $Header: /cactusdevcvs/Cactus/doc/ThornGuide/template.tex,v 1.12 2004/01/07 20:12:39 rideout Exp $ +% +% Thorn documentation in the latex file doc/documentation.tex +% will be included in ThornGuides built with the Cactus make system. +% The scripts employed by the make system automatically include +% pages about variables, parameters and scheduling parsed from the +% relevant thorn CCL files. +% +% This template contains guidelines which help to assure that your +% documentation will be correctly added to ThornGuides. More +% information is available in the Cactus UsersGuide. +% +% Guidelines: +% - Do not change anything before the line +% % START CACTUS THORNGUIDE", +% except for filling in the title, author, date, etc. fields. +% - Each of these fields should only be on ONE line. +% - Author names should be separated with a \\ or a comma. +% - You can define your own macros, but they must appear after +% the START CACTUS THORNGUIDE line, and must not redefine standard +% latex commands. +% - To avoid name clashes with other thorns, 'labels', 'citations', +% 'references', and 'image' names should conform to the following +% convention: +% ARRANGEMENT_THORN_LABEL +% For example, an image wave.eps in the arrangement CactusWave and +% thorn WaveToyC should be renamed to CactusWave_WaveToyC_wave.eps +% - Graphics should only be included using the graphicx package. +% More specifically, with the "\includegraphics" command. Do +% not specify any graphic file extensions in your .tex file. This +% will allow us to create a PDF version of the ThornGuide +% via pdflatex. +% - References should be included with the latex "\bibitem" command. +% - Use \begin{abstract}...\end{abstract} instead of \abstract{...} +% - Do not use \appendix, instead include any appendices you need as +% standard sections. +% - For the benefit of our Perl scripts, and for future extensions, +% please use simple latex. +% +% *======================================================================* +% +% Example of including a graphic image: +% \begin{figure}[ht] +% \begin{center} +% \includegraphics[width=6cm]{MyArrangement_MyThorn_MyFigure} +% \end{center} +% \caption{Illustration of this and that} +% \label{MyArrangement_MyThorn_MyLabel} +% \end{figure} +% +% Example of using a label: +% \label{MyArrangement_MyThorn_MyLabel} +% +% Example of a citation: +% \cite{MyArrangement_MyThorn_Author99} +% +% Example of including a reference +% \bibitem{MyArrangement_MyThorn_Author99} +% {J. Author, {\em The Title of the Book, Journal, or periodical}, 1 (1999), +% 1--16. {\tt http://www.nowhere.com/}} +% +% *======================================================================* + +% If you are using CVS use this line to give version information +% $Header: /cactusdevcvs/Cactus/doc/ThornGuide/template.tex,v 1.12 2004/01/07 20:12:39 rideout Exp $ + +\documentclass{article} + +% Use the Cactus ThornGuide style file +% (Automatically used from Cactus distribution, if you have a +% thorn without the Cactus Flesh download this from the Cactus +% homepage at www.cactuscode.org) +\usepackage{../../../../doc/latex/cactus} + +\begin{document} + +% The author of the documentation +\author{Christian Reisswig \textless reisswig(a)aei.mpg.de\textgreater} + +% The title of the document (not necessarily the name of the Thorn) +\title{WorldTube} + +% the date your document was last changed, if your document is in CVS, +% please use: +% \date{$ $Date: 2004/01/07 20:12:39 $ $} +\date{January 30 2009} + +\maketitle + +% Do not delete next line +% START CACTUS THORNGUIDE + +% Add all definitions used in this documentation here +% \def\mydef etc + +% Add an abstract for this thorn's documentation +\begin{abstract} + +\end{abstract} + +% The following sections are suggestive only. +% Remove them or add your own. + +\section{Introduction} + +\section{Physical System} + +\section{Numerical Implementation} + +\section{Using This Thorn} + +\subsection{Obtaining This Thorn} + +\subsection{Basic Usage} + +\subsection{Special Behaviour} + +\subsection{Interaction With Other Thorns} + +\subsection{Examples} + +\subsection{Support and Feedback} + +\section{History} + +\subsection{Thorn Source Code} + +\subsection{Thorn Documentation} + +\subsection{Acknowledgements} + + +\begin{thebibliography}{9} + +\end{thebibliography} + +% Do not delete next line +% END CACTUS THORNGUIDE + +\end{document} Directory: /PITTNullCode/WorldTube/src/ ======================================= File [added]: areal_sphere.cc Delta lines: +226 -0 =================================================================== --- PITTNullCode/WorldTube/src/areal_sphere.cc (rev 0) +++ PITTNullCode/WorldTube/src/areal_sphere.cc 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,226 @@ +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" +#include "cctk_Functions.h" +#include <assert.h> +#include <iostream> + +#include "register.h" +#include "slices.hh" + + + +using namespace WorldTube; +using namespace SPS; + +double norm_inf(double const * const f, int const ni, int const nj); +int get_1patch_metric_component(double g_sph[4][4], double const r, double const theta, double const phi, + const double lapse, const double shiftx, const double shifty, const double shiftz, + const double gxx, const double gxy, const double gxz, + const double gyy, const double gyz, + const double gzz); +double calc_R(double const target_radius, + double const rp, double const tp, double const pp, const double g_sph[4][4], + double const diff_theta, double const diff_phi); + +extern "C" void WorldTube_ArealSphere(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + int mode_dim = (lmax+1)*(lmax+1); + + + for (int ii=0; ii < ntubes; ++ii) + { + const int si = which_slice_to_take[ii]; + + if (seek_const_areal_radius[ii]) + { + int it = 0; + double res_Linf = 2.0*tolerance; + while ((res_Linf > tolerance) && (it < max_it)) + { + // first of all we do a collective interpolation for the current iterated radius onto the sphere + SphericalSlice_CollectiveSyncVariables(cctkGH, &sid[si].front(), 10); + + + if (CCTK_Equals(type[si], "1patch")) + { + vector<CCTK_REAL> R(slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[0] * slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[1], 0); + vector<CCTK_REAL> res(slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[0] * slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[1], 0); + + for (spheredata_1patch<CCTK_REAL>::const_iter i=slices_1patch(INDEX1P(ss_radius_id[si]), 0).begin(); !i.done(); ++i) + { + // dont't integrate in ghostzones + if (i.ghostzone()) + continue; + + int const idx = i.idx().ij; + + // get metric components at current point in polar-spherical basis + CCTK_REAL g_sph[4][4]; + const CCTK_REAL lapse = slices_1patch(INDEX1P(sid[si][0]), 0)(i); + const CCTK_REAL shiftx = slices_1patch(INDEX1P(sid[si][1]), 0)(i); + const CCTK_REAL shifty = slices_1patch(INDEX1P(sid[si][2]), 0)(i); + const CCTK_REAL shiftz = slices_1patch(INDEX1P(sid[si][3]), 0)(i); + const CCTK_REAL gxx = slices_1patch(INDEX1P(sid[si][4]), 0)(i); + const CCTK_REAL gxy = slices_1patch(INDEX1P(sid[si][5]), 0)(i); + const CCTK_REAL gxz = slices_1patch(INDEX1P(sid[si][6]), 0)(i); + const CCTK_REAL gyy = slices_1patch(INDEX1P(sid[si][7]), 0)(i); + const CCTK_REAL gyz = slices_1patch(INDEX1P(sid[si][8]), 0)(i); + const CCTK_REAL gzz = slices_1patch(INDEX1P(sid[si][9]), 0)(i); + get_1patch_metric_component(g_sph, *i, i.idx().theta, i.idx().phi, + lapse, + shiftx, shifty, shiftz, + gxx, gxy, gxz, + gyy, gyz, + gzz); + + // get angular derivatives of radius + CCTK_REAL diff_theta = slices_1patch(INDEX1P(ss_radius_id[si]), 0).dx(i); + CCTK_REAL diff_phi = slices_1patch(INDEX1P(ss_radius_id[si]), 0).dy(i); + + R[idx] = calc_R(radius[si], *i, i.idx().theta, i.idx().phi, g_sph, diff_theta, diff_phi); + res[idx] = fabs(R[idx] - r[idx]); + } + + res_Linf = norm_inf(&res.front(), slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[0], + slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[1]); + + for (spheredata_1patch<CCTK_REAL>::iter i=slices_1patch(INDEX1P(ss_radius_id[si]), 0).begin(); !i.done(); ++i) + *i = R[i.idx().ij]; + + ++it; + } + if (CCTK_Equals(type[si], "6patch")) + CCTK_WARN(0, "Sorry: Finding surface of constant areal radius not implemented for 6-patches yet :("); + + + } + } + + } +} + + + +/* + * Infinity norm over the slice. + */ +double +norm_inf(double const * const f, int const ni, int const nj) +{ + double norm = f[0]; + for (int idx=0; idx<ni*nj; ++idx) + if (fabs(f[idx]) > norm) + norm = fabs(f[idx]); + double global_norm; + MPI_Allreduce (&norm, &global_norm, 1, MPI_DOUBLE, MPI_MAX, MPI_COMM_WORLD ); + return global_norm; +} + + +/* + * Get the metric at a point on the sphere (1-patch coordinate basis). + * + * The metric components should be given in cartesian components to + * fill the g_cart[] array. This is transformed to spherical polar + * coordinates g_sph[] using the jacobian J[]. + * + * In practice, the jacobian J[] will also need to recognise the possibility + * of alternate patch systems that are in use, eg. thornburg04 style 6-patch + * systems. + * + * The metric components should be obtained by interpolation from the + * numerically generated spacetime. + * + */ +int +get_1patch_metric_component(double g_sph[4][4], double const r, double const theta, double const phi, + const double lapse, const double shiftx, const double shifty, const double shiftz, + const double gxx, const double gxy, const double gxz, + const double gyy, const double gyz, + const double gzz) +{ + double const J[4][4] = { + { 1, 0, 0, 0 }, + { 0, sin(theta)*cos(phi), r*cos(theta)*cos(phi), -r*sin(theta)*sin(phi) }, + { 0, sin(theta)*sin(phi), r*cos(theta)*sin(phi), r*sin(theta)*cos(phi) }, + { 0, cos(theta), -r*sin(theta), 0 } + }; + + double const dt2 = -lapse*lapse + + gxx*shiftx*shiftx + + 2*gxy*shiftx*shifty + + 2*gxz*shiftx*shiftz + + gyy*shifty*shifty + + 2*gyz*shifty*shiftz + + gzz*shiftz*shiftz; + + double const dtdx = gxx*shiftx+gxy*shifty+gxz*shiftz; + double const dtdy = gxy*shiftx+gyy*shifty+gyz*shiftz; + double const dtdz = gxz*shiftx+gyz*shifty+gzz*shiftz; + + double const g_cart[4][4] = { + { dt2, dtdx, dtdy, dtdz }, + { dtdx, gxx, gxy, gxz }, + { dtdy, gxy, gyy, gyz }, + { dtdz, gxz, gyz, gzz } + }; + + for (int a=0; a<4; ++a) + for (int b=0; b<4; ++b) + { + g_sph[a][b] = 0.0; + for (int ix=1; ix<4; ++ix) + for (int jx=1; jx<4; ++jx) + g_sph[a][b] += g_cart[ix][jx]*J[ix][a]*J[jx][b]; + } + + + return 0; +} + + + +/* + * Calculate the updated radius, R, at a point given an initial guess + * for the radius, r, and the metric on the sphere. See Nigel's notes + * for details on the computation. + */ +double +calc_R(double const target_radius, + double const rp, double const tp, double const pp, const double g_sph[4][4], + double const diff_theta, double const diff_phi) +{ + double const grr = g_sph[1][1]; + double const hr2 = g_sph[1][2] / rp; + double const hr3 = g_sph[1][3] / rp; + double const h22 = g_sph[2][2] / (rp*rp); + double const h23 = g_sph[2][3] / (rp*rp); + double const h33 = g_sph[3][3] / (rp*rp); + + double const log_r_theta = diff_theta / rp; + double const log_r_phi = diff_phi / rp; + + double const cg = sin(tp)*sin(tp) / + (h22*h33 - h23 * h23 + + 2.0 * log_r_theta * hr2 * log_r_phi * hr3 + + 2.0 * log_r_theta * hr2 * h33 + + h22 * log_r_phi * log_r_phi * grr + + log_r_theta * log_r_theta * grr * h33 + - 2.0 * log_r_theta * log_r_phi * grr * h23 + - log_r_phi * log_r_phi * hr2 * hr2 + + 2.0 * h22 * log_r_phi * hr3 + - log_r_theta * log_r_theta * hr3 * hr3 + - 2.0 * log_r_theta * hr3 * h23 + - 2.0 * log_r_phi * hr2 * h23); + + return target_radius * pow(cg, 0.25); +} + + + + + File [added]: decompose.cc Delta lines: +45 -0 =================================================================== --- PITTNullCode/WorldTube/src/decompose.cc (rev 0) +++ PITTNullCode/WorldTube/src/decompose.cc 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,45 @@ +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" +#include "cctk_Functions.h" +#include <assert.h> +#include <iostream> + +#include "register.h" + + + + + +namespace WorldTube { + + +extern "C" void WorldTube_Decompose(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + int mode_dim = (lmax+1)*(lmax+1); + + // set all coeffs to zero initially + for (int i=0; i < 50*ntubes*(lmax+1)*(lmax+1); ++i) + extraction_vars[i] = CCTK_Cmplx(0,0); + + for (int si=0; si < ntubes; ++si) + { + + // first of all we do a collective interpolation onto the sphere + SphericalSlice_CollectiveSyncVariables(cctkGH, &sid[si].front(), extr_vars.size()); + + for (int j=0; j < extr_vars.size(); ++j) + { + SphericalSlice_ContractVariableWithAllsYlm(sid[si][j], 0, 0, 0, lmax, &extraction_vars[mode_dim*(si + ntubes*j)]); + } + } +} + + +} + + + File [added]: make.code.defn Delta lines: +8 -0 =================================================================== --- PITTNullCode/WorldTube/src/make.code.defn (rev 0) +++ PITTNullCode/WorldTube/src/make.code.defn 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,8 @@ +# Main make.code.defn file for thorn WorldTube +# $Header:$ + +# Source files in this directory +SRCS = register.cc decompose.cc areal_sphere.cc paramcheck.cc + +# Subdirectories containing source files +SUBDIRS = File [added]: paramcheck.cc Delta lines: +73 -0 =================================================================== --- PITTNullCode/WorldTube/src/paramcheck.cc (rev 0) +++ PITTNullCode/WorldTube/src/paramcheck.cc 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,73 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + + + +extern "C" void WorldTube_ParamCheck(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + // do this only on processor 0 + if (CCTK_MyProc(cctkGH) != 0) + return; + + for (int i=0; i < ntubes; ++i) + { + if (which_slice_to_take[i] >= ss_nslices) + CCTK_PARAMWARN("which_slice_to_take[i]: Slice number is greater than available slices!"); + + for (int j=0; j < ntubes; ++j) + if (i != j && which_slice_to_take[i] == which_slice_to_take[j]) + CCTK_PARAMWARN("which_slice_to_take[i]: You have selected the same slice twice!"); + + // check, if the number of points required for the (lmax, mmax)-mode is met + // Eg: l=2, m=2 has 4 maxima -> we need at least twice as many points in order to resolve that + int required_points = lmax*4; + int recommended_points = lmax*8; + + if (CCTK_Equals(type[which_slice_to_take[i]], "1patch")) + { + if (nphi[which_slice_to_take[i]] < required_points || 2*ntheta[which_slice_to_take[i]] < required_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d you need at least nphi=2*ntheta=%d to resolve it.", lmax, required_points); + if (nphi[which_slice_to_take[i]] < recommended_points || 2*ntheta[which_slice_to_take[i]] < recommended_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d I recommend nphi=2*ntheta=%d points.", lmax, recommended_points); + } + if (CCTK_Equals(type[which_slice_to_take[i]], "6patch")) + { + if (6*nphi[which_slice_to_take[i]] < required_points || 6*ntheta[which_slice_to_take[i]] < required_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d you need at least 6*nphi=6*ntheta=%d to resolve it.", lmax, required_points); + if (6*nphi[which_slice_to_take[i]] < recommended_points || 6*ntheta[which_slice_to_take[i]] < recommended_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d I recommend 6*nphi=6*ntheta=%d points.", lmax, recommended_points); + } + + if (seek_const_areal_radius[i]) + if (nghostzones[which_slice_to_take[i]] != 2) + CCTK_PARAMWARN("Ghostsize must be >= 2."); + } + + if (CCTK_Equals(boundary_behavior, "CCE")) + { + int calc_dr = *((int*) CCTK_ParameterGet("store_radial_derivatives", "ADMDerivatives", NULL)); + int calc_dt = *((int*) CCTK_ParameterGet("store_time_derivatives", "ADMDerivatives", NULL)); + + if (calc_dr == 0 || calc_dt == 0) + CCTK_WARN(0, "ADMDerivatives has to be set to calculate radial and time derivatives!"); + } + + if (CCTK_Equals(boundary_behavior, "CCE_cart")) + { + int calc_dx = *((int*) CCTK_ParameterGet("store_cartesian_derivatives", "ADMDerivatives", NULL)); + int calc_dt = *((int*) CCTK_ParameterGet("store_time_derivatives", "ADMDerivatives", NULL)); + + if (calc_dx == 0 || calc_dt == 0) + CCTK_WARN(0, "ADMDerivatives has to be set to calculate cartesian and time derivatives!"); + } + +} File [added]: register.cc Delta lines: +168 -0 =================================================================== --- PITTNullCode/WorldTube/src/register.cc (rev 0) +++ PITTNullCode/WorldTube/src/register.cc 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,168 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +#include "register.h" + +namespace WorldTube { + + +using namespace std; + + +// the slice ids for each slice and variable +vector<vector<int> > sid; +// all Cauchy extracted variables +vector<string> extr_vars; + + + +extern "C" void WorldTube_RegisterSlices(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + if (CCTK_Equals(boundary_behavior, "extract metric")) + { + extr_vars = vector<string>(10); + extr_vars[0] = "ADMBase::alp"; + extr_vars[1] = "ADMBase::betax"; + extr_vars[2] = "ADMBase::betay"; + extr_vars[3] = "ADMBase::betaz"; + extr_vars[4] = "ADMBase::gxx"; + extr_vars[5] = "ADMBase::gxy"; + extr_vars[6] = "ADMBase::gxz"; + extr_vars[7] = "ADMBase::gyy"; + extr_vars[8] = "ADMBase::gyz"; + extr_vars[9] = "ADMBase::gzz"; + } + if (CCTK_Equals(boundary_behavior, "CCE")) + { + extr_vars = vector<string>(30); + extr_vars[0] = "ADMBase::alp"; + extr_vars[1] = "ADMBase::betax"; + extr_vars[2] = "ADMBase::betay"; + extr_vars[3] = "ADMBase::betaz"; + extr_vars[4] = "ADMBase::gxx"; + extr_vars[5] = "ADMBase::gxy"; + extr_vars[6] = "ADMBase::gxz"; + extr_vars[7] = "ADMBase::gyy"; + extr_vars[8] = "ADMBase::gyz"; + extr_vars[9] = "ADMBase::gzz"; + + // radial derivatives + + extr_vars[10] = "ADMDerivatives::alp_dr"; + extr_vars[11] = "ADMDerivatives::betax_dr"; + extr_vars[12] = "ADMDerivatives::betay_dr"; + extr_vars[13] = "ADMDerivatives::betaz_dr"; + extr_vars[14] = "ADMDerivatives::gxx_dr"; + extr_vars[15] = "ADMDerivatives::gxy_dr"; + extr_vars[16] = "ADMDerivatives::gxz_dr"; + extr_vars[17] = "ADMDerivatives::gyy_dr"; + extr_vars[18] = "ADMDerivatives::gyz_dr"; + extr_vars[19] = "ADMDerivatives::gzz_dr"; + + // time derivatives + + // this can only be checked here (and not in parametercheck) because this is initialized in basegrid! + if (*dtlapse_state == 0 || *dtshift_state == 0) + CCTK_WARN(0, "You need to activate storage for ADMBase::dtalp and ADMBase::dtshift for CCE!"); + + extr_vars[20] = "ADMBase::dtalp"; + extr_vars[21] = "ADMBase::dtbetax"; + extr_vars[22] = "ADMBase::dtbetay"; + extr_vars[23] = "ADMBase::dtbetaz"; + extr_vars[24] = "ADMDerivatives::gxx_dt"; + extr_vars[25] = "ADMDerivatives::gxy_dt"; + extr_vars[26] = "ADMDerivatives::gxz_dt"; + extr_vars[27] = "ADMDerivatives::gyy_dt"; + extr_vars[28] = "ADMDerivatives::gyz_dt"; + extr_vars[29] = "ADMDerivatives::gzz_dt"; + } + + if (CCTK_Equals(boundary_behavior, "CCE_cart")) + { + extr_vars = vector<string>(50); + extr_vars[0] = "ADMBase::alp"; + extr_vars[1] = "ADMBase::betax"; + extr_vars[2] = "ADMBase::betay"; + extr_vars[3] = "ADMBase::betaz"; + extr_vars[4] = "ADMBase::gxx"; + extr_vars[5] = "ADMBase::gxy"; + extr_vars[6] = "ADMBase::gxz"; + extr_vars[7] = "ADMBase::gyy"; + extr_vars[8] = "ADMBase::gyz"; + extr_vars[9] = "ADMBase::gzz"; + + // cartesian derivatives + + extr_vars[10] = "ADMDerivatives::alp_dx"; + extr_vars[11] = "ADMDerivatives::betax_dx"; + extr_vars[12] = "ADMDerivatives::betay_dx"; + extr_vars[13] = "ADMDerivatives::betaz_dx"; + extr_vars[14] = "ADMDerivatives::gxx_dx"; + extr_vars[15] = "ADMDerivatives::gxy_dx"; + extr_vars[16] = "ADMDerivatives::gxz_dx"; + extr_vars[17] = "ADMDerivatives::gyy_dx"; + extr_vars[18] = "ADMDerivatives::gyz_dx"; + extr_vars[19] = "ADMDerivatives::gzz_dx"; + + extr_vars[20] = "ADMDerivatives::alp_dy"; + extr_vars[21] = "ADMDerivatives::betax_dy"; + extr_vars[22] = "ADMDerivatives::betay_dy"; + extr_vars[23] = "ADMDerivatives::betaz_dy"; + extr_vars[24] = "ADMDerivatives::gxx_dy"; + extr_vars[25] = "ADMDerivatives::gxy_dy"; + extr_vars[26] = "ADMDerivatives::gxz_dy"; + extr_vars[27] = "ADMDerivatives::gyy_dy"; + extr_vars[28] = "ADMDerivatives::gyz_dy"; + extr_vars[29] = "ADMDerivatives::gzz_dy"; + + extr_vars[30] = "ADMDerivatives::alp_dz"; + extr_vars[31] = "ADMDerivatives::betax_dz"; + extr_vars[32] = "ADMDerivatives::betay_dz"; + extr_vars[33] = "ADMDerivatives::betaz_dz"; + extr_vars[34] = "ADMDerivatives::gxx_dz"; + extr_vars[35] = "ADMDerivatives::gxy_dz"; + extr_vars[36] = "ADMDerivatives::gxz_dz"; + extr_vars[37] = "ADMDerivatives::gyy_dz"; + extr_vars[38] = "ADMDerivatives::gyz_dz"; + extr_vars[39] = "ADMDerivatives::gzz_dz"; + + // time derivatives + + // this can only be checked here (and not in parametercheck) because this is initialized in basegrid! + if (*dtlapse_state == 0 || *dtshift_state == 0) + CCTK_WARN(0, "You need to activate storage for ADMBase::dtalp and ADMBase::dtshift for CCE!"); + + extr_vars[40] = "ADMBase::dtalp"; + extr_vars[41] = "ADMBase::dtbetax"; + extr_vars[42] = "ADMBase::dtbetay"; + extr_vars[43] = "ADMBase::dtbetaz"; + extr_vars[44] = "ADMDerivatives::gxx_dt"; + extr_vars[45] = "ADMDerivatives::gxy_dt"; + extr_vars[46] = "ADMDerivatives::gxz_dt"; + extr_vars[47] = "ADMDerivatives::gyy_dt"; + extr_vars[48] = "ADMDerivatives::gyz_dt"; + extr_vars[49] = "ADMDerivatives::gzz_dt"; + } + + + sid = vector<vector<int> >(ntubes); + + for (int i=0; i < ntubes; ++i) + { + sid[i] = vector<int>(extr_vars.size(), 0); + for (int j=0; j < extr_vars.size(); ++j) + { + sid[i][j] = SphericalSlice_RegisterVariable(extr_vars[j].c_str(), which_slice_to_take[i], 1, "split"); + } + } +} + + + +} + File [added]: register.h Delta lines: +22 -0 =================================================================== --- PITTNullCode/WorldTube/src/register.h (rev 0) +++ PITTNullCode/WorldTube/src/register.h 2011-06-13 21:44:24 UTC (rev 37) @@ -0,0 +1,22 @@ +#ifndef _WORLDTUBE_ +#define _WORLDTUBE_ + +#include <vector> +#include <iostream> + +namespace WorldTube { + +using namespace std; + +// the slice ids for each slice and variable +extern vector<vector<int> > sid; + +// all Cauchy extracted variables +extern vector<string> extr_vars; + + + +} + + +#endif

1 0

[svn:einsteintoolkit] incoming/WorldTube/ (Rev. 36)
by reisswig＠tapir.caltech.edu 13 Jun '11

13 Jun '11

User: reisswig Date: 2011/06/13 01:37 PM Added: /WorldTube/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /WorldTube/doc/ documentation.tex /WorldTube/src/ areal_sphere.cc, decompose.cc, make.code.defn, paramcheck.cc, register.cc, register.h Log: A thorn capable of decomposing metric data on spheres. File Changes: Directory: /WorldTube/ ====================== File [added]: README Delta lines: +12 -0 =================================================================== --- WorldTube/README (rev 0) +++ WorldTube/README 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,12 @@ +CVS info : $Header:$ + +Cactus Code Thorn WorldTube +Thorn Author(s) : Christian Reisswig <reisswig(a)aei.mpg.de> +Thorn Maintainer(s) : Christian Reisswig <reisswig(a)aei.mpg.de> +-------------------------------------------------------------------------- + +Purpose of the thorn: + +This thorn is ment as an interface between Cauchy and characteristic codes. + +It can also just output metric components and find spheres of constant areal radius. \ No newline at end of file File [added]: configuration.ccl Delta lines: +1 -0 =================================================================== --- WorldTube/configuration.ccl (rev 0) +++ WorldTube/configuration.ccl 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1 @@ +REQUIRES THORNS: ADMBase ADMDerivatives File [added]: interface.ccl Delta lines: +65 -0 =================================================================== --- WorldTube/interface.ccl (rev 0) +++ WorldTube/interface.ccl 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,65 @@ +# Interface definition for thorn WorldTube +# $Header:$ + +implements: WorldTube +inherits: SphericalSlice ADMBase + + +uses include: slices.hh + + + +# Register a Cactus variable for being sliced and return registration number +CCTK_INT FUNCTION \ + SphericalSlice_RegisterVariable \ + (CCTK_POINTER_TO_CONST IN varname, \ + CCTK_INT IN sn, \ + CCTK_INT IN timelevels, \ + CCTK_POINTER_TO_CONST IN distrib_method) +REQUIRES FUNCTION SphericalSlice_RegisterVariable + + +# Return an array of all values of scalar product with all spin-weighted spherical harmonics sYlm up to lmax +CCTK_INT FUNCTION \ + SphericalSlice_ContractVariableWithAllsYlm \ + (CCTK_INT IN varno, \ + CCTK_INT IN timelevel, \ + CCTK_INT IN s, \ + CCTK_INT IN lmin, \ + CCTK_INT IN lmax, \ + CCTK_COMPLEX OUT ARRAY coeffs) +REQUIRES FUNCTION SphericalSlice_ContractVariableWithAllsYlm \ + + +# Synchronzie slice, i.e. interpolate a whole set of variables on the same slice +CCTK_INT FUNCTION \ + SphericalSlice_CollectiveSyncVariables \ + (CCTK_POINTER_TO_CONST IN cctkGH_, \ + CCTK_POINTER_TO_CONST IN varno, \ + CCTK_INT IN number_of_vars) +REQUIRES FUNCTION SphericalSlice_CollectiveSyncVariables \ + + + +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) +CCTK_INT FUNCTION Boundary_SelectGroupForBC(CCTK_POINTER_TO_CONST IN GH, \ + CCTK_INT IN faces, CCTK_INT IN boundary_width, CCTK_INT IN table_handle, \ + CCTK_STRING IN group_name, CCTK_STRING IN bc_name) + +REQUIRES FUNCTION Boundary_SelectVarForBC +REQUIRES FUNCTION Boundary_SelectGroupForBC + + +PUBLIC: + +COMPLEX extracted_vars[50] TYPE=ARRAY TIMELEVELS=1 DIM=2 SIZE=((lmax+1)*(lmax+1)),ntubes DISTRIB=CONSTANT +{ + extraction_vars +} "Stores extracted lm-modes for each metric component" + + + + + File [added]: param.ccl Delta lines: +67 -0 =================================================================== --- WorldTube/param.ccl (rev 0) +++ WorldTube/param.ccl 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,67 @@ +# Parameter definitions for thorn WorldTube +# $Header:$ + + +KEYWORD boundary_behavior "The type of boundary" +{ + "none" :: "don't extract/inject anything" + "extract metric" :: "extract metric components, lapse and shift" + "CCE" :: "extract metric components, lapse, shift and derivatives for CCE" + "CCE_cart" :: "extract metric components, lapse, shift and (cartesian) derivatives for CCE" + "CCM" :: "to be implemented" +} "CCE" + + +INT lmax "Maximum extracted l-mode" +{ + 0:* :: "any positive number" +} 8 + + +CCTK_INT ntubes "The number of worldtubes to be placed" +{ + 0:42 :: "Any number between 0 and 42" +} 1 + + + +CCTK_INT which_slice_to_take[42] "Which SphericalSlice to use" +{ + 0:* :: "Any positive number" +} 0 + + +BOOLEAN seek_const_areal_radius[42] "Set slice radius to a surface of constant areal radius before extracting/injecting" +{ +} no + + +REAL tolerance "Tolerance for finding areal radius" +{ + 0:* :: "Any positive number" +} 1e-6 + + +INT max_it "Maximum number of iterations for finding areal radius" +{ + 1:* :: "any positive number greater than 1" +} 10 + + +CCTK_INT radial_deriv_order "Order of radial derivatives" +{ + -1 :: "Use whatever is specified by SummationByParts" + 2:8 :: "Use this order" +} -1 + + + +SHARES: SphericalSlice +USES KEYWORD type[42] +USES REAL radius[42] +USES INT nslices AS ss_nslices +USES INT ntheta +USES INT nphi +USES INT nghostzones + + File [added]: schedule.ccl Delta lines: +45 -0 =================================================================== --- WorldTube/schedule.ccl (rev 0) +++ WorldTube/schedule.ccl 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,45 @@ +# Schedule definitions for thorn WorldTube +# $Header:$ + + +STORAGE: extracted_vars + + +SCHEDULE WorldTube_ParamCheck AT PARAMCHECK +{ + LANG: C + OPTIONS: GLOBAL +} "Check parameters" + + +SCHEDULE WorldTube_RegisterSlices AT CCTK_INITIAL AFTER SphericalSlice_HasBeenSet +{ + LANG: C + OPTIONS: GLOBAL +} "Register variables for extraction" + + +SCHEDULE WorldTube_RegisterSlices AT CCTK_POST_RECOVER_VARIABLES AFTER SphericalSlice_HasBeenSet +{ + LANG: C + OPTIONS: GLOBAL +} "Register variables for extraction" + + + +SCHEDULE WorldTube_ArealSphere AT CCTK_POSTSTEP AFTER SphericalSlice_HasBeenSet +{ + LANG: C + OPTIONS: GLOBAL +} "Calculate surfaces of constant areal radius if selected" + + + +SCHEDULE WorldTube_Decompose AT CCTK_ANALYSIS AFTER (SphericalSlice_HasBeenSet, WorldTube_ArealSphere) +{ + LANG: C + OPTIONS: GLOBAL + TRIGGERS: extracted_vars +} "Contract sliced variables with sYlm's" + + Directory: /WorldTube/doc/ ========================== File [added]: documentation.tex Delta lines: +144 -0 =================================================================== --- WorldTube/doc/documentation.tex (rev 0) +++ WorldTube/doc/documentation.tex 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,144 @@ +% *======================================================================* +% Cactus Thorn template for ThornGuide documentation +% Author: Ian Kelley +% Date: Sun Jun 02, 2002 +% $Header: /cactusdevcvs/Cactus/doc/ThornGuide/template.tex,v 1.12 2004/01/07 20:12:39 rideout Exp $ +% +% Thorn documentation in the latex file doc/documentation.tex +% will be included in ThornGuides built with the Cactus make system. +% The scripts employed by the make system automatically include +% pages about variables, parameters and scheduling parsed from the +% relevant thorn CCL files. +% +% This template contains guidelines which help to assure that your +% documentation will be correctly added to ThornGuides. More +% information is available in the Cactus UsersGuide. +% +% Guidelines: +% - Do not change anything before the line +% % START CACTUS THORNGUIDE", +% except for filling in the title, author, date, etc. fields. +% - Each of these fields should only be on ONE line. +% - Author names should be separated with a \\ or a comma. +% - You can define your own macros, but they must appear after +% the START CACTUS THORNGUIDE line, and must not redefine standard +% latex commands. +% - To avoid name clashes with other thorns, 'labels', 'citations', +% 'references', and 'image' names should conform to the following +% convention: +% ARRANGEMENT_THORN_LABEL +% For example, an image wave.eps in the arrangement CactusWave and +% thorn WaveToyC should be renamed to CactusWave_WaveToyC_wave.eps +% - Graphics should only be included using the graphicx package. +% More specifically, with the "\includegraphics" command. Do +% not specify any graphic file extensions in your .tex file. This +% will allow us to create a PDF version of the ThornGuide +% via pdflatex. +% - References should be included with the latex "\bibitem" command. +% - Use \begin{abstract}...\end{abstract} instead of \abstract{...} +% - Do not use \appendix, instead include any appendices you need as +% standard sections. +% - For the benefit of our Perl scripts, and for future extensions, +% please use simple latex. +% +% *======================================================================* +% +% Example of including a graphic image: +% \begin{figure}[ht] +% \begin{center} +% \includegraphics[width=6cm]{MyArrangement_MyThorn_MyFigure} +% \end{center} +% \caption{Illustration of this and that} +% \label{MyArrangement_MyThorn_MyLabel} +% \end{figure} +% +% Example of using a label: +% \label{MyArrangement_MyThorn_MyLabel} +% +% Example of a citation: +% \cite{MyArrangement_MyThorn_Author99} +% +% Example of including a reference +% \bibitem{MyArrangement_MyThorn_Author99} +% {J. Author, {\em The Title of the Book, Journal, or periodical}, 1 (1999), +% 1--16. {\tt http://www.nowhere.com/}} +% +% *======================================================================* + +% If you are using CVS use this line to give version information +% $Header: /cactusdevcvs/Cactus/doc/ThornGuide/template.tex,v 1.12 2004/01/07 20:12:39 rideout Exp $ + +\documentclass{article} + +% Use the Cactus ThornGuide style file +% (Automatically used from Cactus distribution, if you have a +% thorn without the Cactus Flesh download this from the Cactus +% homepage at www.cactuscode.org) +\usepackage{../../../../doc/latex/cactus} + +\begin{document} + +% The author of the documentation +\author{Christian Reisswig \textless reisswig(a)aei.mpg.de\textgreater} + +% The title of the document (not necessarily the name of the Thorn) +\title{WorldTube} + +% the date your document was last changed, if your document is in CVS, +% please use: +% \date{$ $Date: 2004/01/07 20:12:39 $ $} +\date{January 30 2009} + +\maketitle + +% Do not delete next line +% START CACTUS THORNGUIDE + +% Add all definitions used in this documentation here +% \def\mydef etc + +% Add an abstract for this thorn's documentation +\begin{abstract} + +\end{abstract} + +% The following sections are suggestive only. +% Remove them or add your own. + +\section{Introduction} + +\section{Physical System} + +\section{Numerical Implementation} + +\section{Using This Thorn} + +\subsection{Obtaining This Thorn} + +\subsection{Basic Usage} + +\subsection{Special Behaviour} + +\subsection{Interaction With Other Thorns} + +\subsection{Examples} + +\subsection{Support and Feedback} + +\section{History} + +\subsection{Thorn Source Code} + +\subsection{Thorn Documentation} + +\subsection{Acknowledgements} + + +\begin{thebibliography}{9} + +\end{thebibliography} + +% Do not delete next line +% END CACTUS THORNGUIDE + +\end{document} Directory: /WorldTube/src/ ========================== File [added]: areal_sphere.cc Delta lines: +226 -0 =================================================================== --- WorldTube/src/areal_sphere.cc (rev 0) +++ WorldTube/src/areal_sphere.cc 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,226 @@ +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" +#include "cctk_Functions.h" +#include <assert.h> +#include <iostream> + +#include "register.h" +#include "slices.hh" + + + +using namespace WorldTube; +using namespace SPS; + +double norm_inf(double const * const f, int const ni, int const nj); +int get_1patch_metric_component(double g_sph[4][4], double const r, double const theta, double const phi, + const double lapse, const double shiftx, const double shifty, const double shiftz, + const double gxx, const double gxy, const double gxz, + const double gyy, const double gyz, + const double gzz); +double calc_R(double const target_radius, + double const rp, double const tp, double const pp, const double g_sph[4][4], + double const diff_theta, double const diff_phi); + +extern "C" void WorldTube_ArealSphere(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + int mode_dim = (lmax+1)*(lmax+1); + + + for (int ii=0; ii < ntubes; ++ii) + { + const int si = which_slice_to_take[ii]; + + if (seek_const_areal_radius[ii]) + { + int it = 0; + double res_Linf = 2.0*tolerance; + while ((res_Linf > tolerance) && (it < max_it)) + { + // first of all we do a collective interpolation for the current iterated radius onto the sphere + SphericalSlice_CollectiveSyncVariables(cctkGH, &sid[si].front(), 10); + + + if (CCTK_Equals(type[si], "1patch")) + { + vector<CCTK_REAL> R(slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[0] * slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[1], 0); + vector<CCTK_REAL> res(slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[0] * slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[1], 0); + + for (spheredata_1patch<CCTK_REAL>::const_iter i=slices_1patch(INDEX1P(ss_radius_id[si]), 0).begin(); !i.done(); ++i) + { + // dont't integrate in ghostzones + if (i.ghostzone()) + continue; + + int const idx = i.idx().ij; + + // get metric components at current point in polar-spherical basis + CCTK_REAL g_sph[4][4]; + const CCTK_REAL lapse = slices_1patch(INDEX1P(sid[si][0]), 0)(i); + const CCTK_REAL shiftx = slices_1patch(INDEX1P(sid[si][1]), 0)(i); + const CCTK_REAL shifty = slices_1patch(INDEX1P(sid[si][2]), 0)(i); + const CCTK_REAL shiftz = slices_1patch(INDEX1P(sid[si][3]), 0)(i); + const CCTK_REAL gxx = slices_1patch(INDEX1P(sid[si][4]), 0)(i); + const CCTK_REAL gxy = slices_1patch(INDEX1P(sid[si][5]), 0)(i); + const CCTK_REAL gxz = slices_1patch(INDEX1P(sid[si][6]), 0)(i); + const CCTK_REAL gyy = slices_1patch(INDEX1P(sid[si][7]), 0)(i); + const CCTK_REAL gyz = slices_1patch(INDEX1P(sid[si][8]), 0)(i); + const CCTK_REAL gzz = slices_1patch(INDEX1P(sid[si][9]), 0)(i); + get_1patch_metric_component(g_sph, *i, i.idx().theta, i.idx().phi, + lapse, + shiftx, shifty, shiftz, + gxx, gxy, gxz, + gyy, gyz, + gzz); + + // get angular derivatives of radius + CCTK_REAL diff_theta = slices_1patch(INDEX1P(ss_radius_id[si]), 0).dx(i); + CCTK_REAL diff_phi = slices_1patch(INDEX1P(ss_radius_id[si]), 0).dy(i); + + R[idx] = calc_R(radius[si], *i, i.idx().theta, i.idx().phi, g_sph, diff_theta, diff_phi); + res[idx] = fabs(R[idx] - r[idx]); + } + + res_Linf = norm_inf(&res.front(), slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[0], + slices_1patch(INDEX1P(ss_radius_id[si]), 0).lsh(0)[1]); + + for (spheredata_1patch<CCTK_REAL>::iter i=slices_1patch(INDEX1P(ss_radius_id[si]), 0).begin(); !i.done(); ++i) + *i = R[i.idx().ij]; + + ++it; + } + if (CCTK_Equals(type[si], "6patch")) + CCTK_WARN(0, "Sorry: Finding surface of constant areal radius not implemented for 6-patches yet :("); + + + } + } + + } +} + + + +/* + * Infinity norm over the slice. + */ +double +norm_inf(double const * const f, int const ni, int const nj) +{ + double norm = f[0]; + for (int idx=0; idx<ni*nj; ++idx) + if (fabs(f[idx]) > norm) + norm = fabs(f[idx]); + double global_norm; + MPI_Allreduce (&norm, &global_norm, 1, MPI_DOUBLE, MPI_MAX, MPI_COMM_WORLD ); + return global_norm; +} + + +/* + * Get the metric at a point on the sphere (1-patch coordinate basis). + * + * The metric components should be given in cartesian components to + * fill the g_cart[] array. This is transformed to spherical polar + * coordinates g_sph[] using the jacobian J[]. + * + * In practice, the jacobian J[] will also need to recognise the possibility + * of alternate patch systems that are in use, eg. thornburg04 style 6-patch + * systems. + * + * The metric components should be obtained by interpolation from the + * numerically generated spacetime. + * + */ +int +get_1patch_metric_component(double g_sph[4][4], double const r, double const theta, double const phi, + const double lapse, const double shiftx, const double shifty, const double shiftz, + const double gxx, const double gxy, const double gxz, + const double gyy, const double gyz, + const double gzz) +{ + double const J[4][4] = { + { 1, 0, 0, 0 }, + { 0, sin(theta)*cos(phi), r*cos(theta)*cos(phi), -r*sin(theta)*sin(phi) }, + { 0, sin(theta)*sin(phi), r*cos(theta)*sin(phi), r*sin(theta)*cos(phi) }, + { 0, cos(theta), -r*sin(theta), 0 } + }; + + double const dt2 = -lapse*lapse + + gxx*shiftx*shiftx + + 2*gxy*shiftx*shifty + + 2*gxz*shiftx*shiftz + + gyy*shifty*shifty + + 2*gyz*shifty*shiftz + + gzz*shiftz*shiftz; + + double const dtdx = gxx*shiftx+gxy*shifty+gxz*shiftz; + double const dtdy = gxy*shiftx+gyy*shifty+gyz*shiftz; + double const dtdz = gxz*shiftx+gyz*shifty+gzz*shiftz; + + double const g_cart[4][4] = { + { dt2, dtdx, dtdy, dtdz }, + { dtdx, gxx, gxy, gxz }, + { dtdy, gxy, gyy, gyz }, + { dtdz, gxz, gyz, gzz } + }; + + for (int a=0; a<4; ++a) + for (int b=0; b<4; ++b) + { + g_sph[a][b] = 0.0; + for (int ix=1; ix<4; ++ix) + for (int jx=1; jx<4; ++jx) + g_sph[a][b] += g_cart[ix][jx]*J[ix][a]*J[jx][b]; + } + + + return 0; +} + + + +/* + * Calculate the updated radius, R, at a point given an initial guess + * for the radius, r, and the metric on the sphere. See Nigel's notes + * for details on the computation. + */ +double +calc_R(double const target_radius, + double const rp, double const tp, double const pp, const double g_sph[4][4], + double const diff_theta, double const diff_phi) +{ + double const grr = g_sph[1][1]; + double const hr2 = g_sph[1][2] / rp; + double const hr3 = g_sph[1][3] / rp; + double const h22 = g_sph[2][2] / (rp*rp); + double const h23 = g_sph[2][3] / (rp*rp); + double const h33 = g_sph[3][3] / (rp*rp); + + double const log_r_theta = diff_theta / rp; + double const log_r_phi = diff_phi / rp; + + double const cg = sin(tp)*sin(tp) / + (h22*h33 - h23 * h23 + + 2.0 * log_r_theta * hr2 * log_r_phi * hr3 + + 2.0 * log_r_theta * hr2 * h33 + + h22 * log_r_phi * log_r_phi * grr + + log_r_theta * log_r_theta * grr * h33 + - 2.0 * log_r_theta * log_r_phi * grr * h23 + - log_r_phi * log_r_phi * hr2 * hr2 + + 2.0 * h22 * log_r_phi * hr3 + - log_r_theta * log_r_theta * hr3 * hr3 + - 2.0 * log_r_theta * hr3 * h23 + - 2.0 * log_r_phi * hr2 * h23); + + return target_radius * pow(cg, 0.25); +} + + + + + File [added]: decompose.cc Delta lines: +45 -0 =================================================================== --- WorldTube/src/decompose.cc (rev 0) +++ WorldTube/src/decompose.cc 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,45 @@ +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" +#include "cctk_Functions.h" +#include <assert.h> +#include <iostream> + +#include "register.h" + + + + + +namespace WorldTube { + + +extern "C" void WorldTube_Decompose(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + int mode_dim = (lmax+1)*(lmax+1); + + // set all coeffs to zero initially + for (int i=0; i < 50*ntubes*(lmax+1)*(lmax+1); ++i) + extraction_vars[i] = CCTK_Cmplx(0,0); + + for (int si=0; si < ntubes; ++si) + { + + // first of all we do a collective interpolation onto the sphere + SphericalSlice_CollectiveSyncVariables(cctkGH, &sid[si].front(), extr_vars.size()); + + for (int j=0; j < extr_vars.size(); ++j) + { + SphericalSlice_ContractVariableWithAllsYlm(sid[si][j], 0, 0, 0, lmax, &extraction_vars[mode_dim*(si + ntubes*j)]); + } + } +} + + +} + + + File [added]: make.code.defn Delta lines: +8 -0 =================================================================== --- WorldTube/src/make.code.defn (rev 0) +++ WorldTube/src/make.code.defn 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,8 @@ +# Main make.code.defn file for thorn WorldTube +# $Header:$ + +# Source files in this directory +SRCS = register.cc decompose.cc areal_sphere.cc paramcheck.cc + +# Subdirectories containing source files +SUBDIRS = File [added]: paramcheck.cc Delta lines: +73 -0 =================================================================== --- WorldTube/src/paramcheck.cc (rev 0) +++ WorldTube/src/paramcheck.cc 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,73 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + + + +extern "C" void WorldTube_ParamCheck(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + // do this only on processor 0 + if (CCTK_MyProc(cctkGH) != 0) + return; + + for (int i=0; i < ntubes; ++i) + { + if (which_slice_to_take[i] >= ss_nslices) + CCTK_PARAMWARN("which_slice_to_take[i]: Slice number is greater than available slices!"); + + for (int j=0; j < ntubes; ++j) + if (i != j && which_slice_to_take[i] == which_slice_to_take[j]) + CCTK_PARAMWARN("which_slice_to_take[i]: You have selected the same slice twice!"); + + // check, if the number of points required for the (lmax, mmax)-mode is met + // Eg: l=2, m=2 has 4 maxima -> we need at least twice as many points in order to resolve that + int required_points = lmax*4; + int recommended_points = lmax*8; + + if (CCTK_Equals(type[which_slice_to_take[i]], "1patch")) + { + if (nphi[which_slice_to_take[i]] < required_points || 2*ntheta[which_slice_to_take[i]] < required_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d you need at least nphi=2*ntheta=%d to resolve it.", lmax, required_points); + if (nphi[which_slice_to_take[i]] < recommended_points || 2*ntheta[which_slice_to_take[i]] < recommended_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d I recommend nphi=2*ntheta=%d points.", lmax, recommended_points); + } + if (CCTK_Equals(type[which_slice_to_take[i]], "6patch")) + { + if (6*nphi[which_slice_to_take[i]] < required_points || 6*ntheta[which_slice_to_take[i]] < required_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d you need at least 6*nphi=6*ntheta=%d to resolve it.", lmax, required_points); + if (6*nphi[which_slice_to_take[i]] < recommended_points || 6*ntheta[which_slice_to_take[i]] < recommended_points) + CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, + "For mode lmax=%d I recommend 6*nphi=6*ntheta=%d points.", lmax, recommended_points); + } + + if (seek_const_areal_radius[i]) + if (nghostzones[which_slice_to_take[i]] != 2) + CCTK_PARAMWARN("Ghostsize must be >= 2."); + } + + if (CCTK_Equals(boundary_behavior, "CCE")) + { + int calc_dr = *((int*) CCTK_ParameterGet("store_radial_derivatives", "ADMDerivatives", NULL)); + int calc_dt = *((int*) CCTK_ParameterGet("store_time_derivatives", "ADMDerivatives", NULL)); + + if (calc_dr == 0 || calc_dt == 0) + CCTK_WARN(0, "ADMDerivatives has to be set to calculate radial and time derivatives!"); + } + + if (CCTK_Equals(boundary_behavior, "CCE_cart")) + { + int calc_dx = *((int*) CCTK_ParameterGet("store_cartesian_derivatives", "ADMDerivatives", NULL)); + int calc_dt = *((int*) CCTK_ParameterGet("store_time_derivatives", "ADMDerivatives", NULL)); + + if (calc_dx == 0 || calc_dt == 0) + CCTK_WARN(0, "ADMDerivatives has to be set to calculate cartesian and time derivatives!"); + } + +} File [added]: register.cc Delta lines: +168 -0 =================================================================== --- WorldTube/src/register.cc (rev 0) +++ WorldTube/src/register.cc 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,168 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +#include "register.h" + +namespace WorldTube { + + +using namespace std; + + +// the slice ids for each slice and variable +vector<vector<int> > sid; +// all Cauchy extracted variables +vector<string> extr_vars; + + + +extern "C" void WorldTube_RegisterSlices(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + if (CCTK_Equals(boundary_behavior, "extract metric")) + { + extr_vars = vector<string>(10); + extr_vars[0] = "ADMBase::alp"; + extr_vars[1] = "ADMBase::betax"; + extr_vars[2] = "ADMBase::betay"; + extr_vars[3] = "ADMBase::betaz"; + extr_vars[4] = "ADMBase::gxx"; + extr_vars[5] = "ADMBase::gxy"; + extr_vars[6] = "ADMBase::gxz"; + extr_vars[7] = "ADMBase::gyy"; + extr_vars[8] = "ADMBase::gyz"; + extr_vars[9] = "ADMBase::gzz"; + } + if (CCTK_Equals(boundary_behavior, "CCE")) + { + extr_vars = vector<string>(30); + extr_vars[0] = "ADMBase::alp"; + extr_vars[1] = "ADMBase::betax"; + extr_vars[2] = "ADMBase::betay"; + extr_vars[3] = "ADMBase::betaz"; + extr_vars[4] = "ADMBase::gxx"; + extr_vars[5] = "ADMBase::gxy"; + extr_vars[6] = "ADMBase::gxz"; + extr_vars[7] = "ADMBase::gyy"; + extr_vars[8] = "ADMBase::gyz"; + extr_vars[9] = "ADMBase::gzz"; + + // radial derivatives + + extr_vars[10] = "ADMDerivatives::alp_dr"; + extr_vars[11] = "ADMDerivatives::betax_dr"; + extr_vars[12] = "ADMDerivatives::betay_dr"; + extr_vars[13] = "ADMDerivatives::betaz_dr"; + extr_vars[14] = "ADMDerivatives::gxx_dr"; + extr_vars[15] = "ADMDerivatives::gxy_dr"; + extr_vars[16] = "ADMDerivatives::gxz_dr"; + extr_vars[17] = "ADMDerivatives::gyy_dr"; + extr_vars[18] = "ADMDerivatives::gyz_dr"; + extr_vars[19] = "ADMDerivatives::gzz_dr"; + + // time derivatives + + // this can only be checked here (and not in parametercheck) because this is initialized in basegrid! + if (*dtlapse_state == 0 || *dtshift_state == 0) + CCTK_WARN(0, "You need to activate storage for ADMBase::dtalp and ADMBase::dtshift for CCE!"); + + extr_vars[20] = "ADMBase::dtalp"; + extr_vars[21] = "ADMBase::dtbetax"; + extr_vars[22] = "ADMBase::dtbetay"; + extr_vars[23] = "ADMBase::dtbetaz"; + extr_vars[24] = "ADMDerivatives::gxx_dt"; + extr_vars[25] = "ADMDerivatives::gxy_dt"; + extr_vars[26] = "ADMDerivatives::gxz_dt"; + extr_vars[27] = "ADMDerivatives::gyy_dt"; + extr_vars[28] = "ADMDerivatives::gyz_dt"; + extr_vars[29] = "ADMDerivatives::gzz_dt"; + } + + if (CCTK_Equals(boundary_behavior, "CCE_cart")) + { + extr_vars = vector<string>(50); + extr_vars[0] = "ADMBase::alp"; + extr_vars[1] = "ADMBase::betax"; + extr_vars[2] = "ADMBase::betay"; + extr_vars[3] = "ADMBase::betaz"; + extr_vars[4] = "ADMBase::gxx"; + extr_vars[5] = "ADMBase::gxy"; + extr_vars[6] = "ADMBase::gxz"; + extr_vars[7] = "ADMBase::gyy"; + extr_vars[8] = "ADMBase::gyz"; + extr_vars[9] = "ADMBase::gzz"; + + // cartesian derivatives + + extr_vars[10] = "ADMDerivatives::alp_dx"; + extr_vars[11] = "ADMDerivatives::betax_dx"; + extr_vars[12] = "ADMDerivatives::betay_dx"; + extr_vars[13] = "ADMDerivatives::betaz_dx"; + extr_vars[14] = "ADMDerivatives::gxx_dx"; + extr_vars[15] = "ADMDerivatives::gxy_dx"; + extr_vars[16] = "ADMDerivatives::gxz_dx"; + extr_vars[17] = "ADMDerivatives::gyy_dx"; + extr_vars[18] = "ADMDerivatives::gyz_dx"; + extr_vars[19] = "ADMDerivatives::gzz_dx"; + + extr_vars[20] = "ADMDerivatives::alp_dy"; + extr_vars[21] = "ADMDerivatives::betax_dy"; + extr_vars[22] = "ADMDerivatives::betay_dy"; + extr_vars[23] = "ADMDerivatives::betaz_dy"; + extr_vars[24] = "ADMDerivatives::gxx_dy"; + extr_vars[25] = "ADMDerivatives::gxy_dy"; + extr_vars[26] = "ADMDerivatives::gxz_dy"; + extr_vars[27] = "ADMDerivatives::gyy_dy"; + extr_vars[28] = "ADMDerivatives::gyz_dy"; + extr_vars[29] = "ADMDerivatives::gzz_dy"; + + extr_vars[30] = "ADMDerivatives::alp_dz"; + extr_vars[31] = "ADMDerivatives::betax_dz"; + extr_vars[32] = "ADMDerivatives::betay_dz"; + extr_vars[33] = "ADMDerivatives::betaz_dz"; + extr_vars[34] = "ADMDerivatives::gxx_dz"; + extr_vars[35] = "ADMDerivatives::gxy_dz"; + extr_vars[36] = "ADMDerivatives::gxz_dz"; + extr_vars[37] = "ADMDerivatives::gyy_dz"; + extr_vars[38] = "ADMDerivatives::gyz_dz"; + extr_vars[39] = "ADMDerivatives::gzz_dz"; + + // time derivatives + + // this can only be checked here (and not in parametercheck) because this is initialized in basegrid! + if (*dtlapse_state == 0 || *dtshift_state == 0) + CCTK_WARN(0, "You need to activate storage for ADMBase::dtalp and ADMBase::dtshift for CCE!"); + + extr_vars[40] = "ADMBase::dtalp"; + extr_vars[41] = "ADMBase::dtbetax"; + extr_vars[42] = "ADMBase::dtbetay"; + extr_vars[43] = "ADMBase::dtbetaz"; + extr_vars[44] = "ADMDerivatives::gxx_dt"; + extr_vars[45] = "ADMDerivatives::gxy_dt"; + extr_vars[46] = "ADMDerivatives::gxz_dt"; + extr_vars[47] = "ADMDerivatives::gyy_dt"; + extr_vars[48] = "ADMDerivatives::gyz_dt"; + extr_vars[49] = "ADMDerivatives::gzz_dt"; + } + + + sid = vector<vector<int> >(ntubes); + + for (int i=0; i < ntubes; ++i) + { + sid[i] = vector<int>(extr_vars.size(), 0); + for (int j=0; j < extr_vars.size(); ++j) + { + sid[i][j] = SphericalSlice_RegisterVariable(extr_vars[j].c_str(), which_slice_to_take[i], 1, "split"); + } + } +} + + + +} + File [added]: register.h Delta lines: +22 -0 =================================================================== --- WorldTube/src/register.h (rev 0) +++ WorldTube/src/register.h 2011-06-13 18:37:44 UTC (rev 36) @@ -0,0 +1,22 @@ +#ifndef _WORLDTUBE_ +#define _WORLDTUBE_ + +#include <vector> +#include <iostream> + +namespace WorldTube { + +using namespace std; + +// the slice ids for each slice and variable +extern vector<vector<int> > sid; + +// all Cauchy extracted variables +extern vector<string> extr_vars; + + + +} + + +#endif

1 0

[svn:einsteintoolkit] incoming/ADMDerivatives/ (Rev. 35)
by reisswig＠tapir.caltech.edu 13 Jun '11

13 Jun '11

User: reisswig Date: 2011/06/13 01:36 PM Added: /ADMDerivatives/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /ADMDerivatives/doc/ documentation.tex /ADMDerivatives/src/ calc_derivs.cc, derivative_operators.h, derivatives.h, make.code.defn, sync.cc Log: Store and compute the derivatives (time, radial, Cartesian) of the metric. File Changes: Directory: /ADMDerivatives/ =========================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /ADMDerivatives/doc/ =============================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /ADMDerivatives/src/ =============================== File [added]: calc_derivs.cc Delta lines: <Unavailable> File [added]: derivative_operators.h Delta lines: <Unavailable> File [added]: derivatives.h Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: sync.cc Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] incoming/SphericalSlice/ (Rev. 34)
by reisswig＠tapir.caltech.edu 13 Jun '11

13 Jun '11

User: reisswig Date: 2011/06/13 01:35 PM Added: /SphericalSlice/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /SphericalSlice/doc/ Doxyfile, documentation.tex, refman.pdf /SphericalSlice/mma/ 1patch.nb, 6patch.nb /SphericalSlice/src/ check.cc, commstack.cc, commstack.hh, contract.cc, int_coeffs.hh, integrate.cc, make.code.defn, paramcheck.cc, register.cc, setup.cc, setup.hh, slices.cc, slices.hh, spheredata.cc, spheredata.hh, spheredata_1patch.cc, spheredata_1patch.hh, spheredata_2patch.cc, spheredata_2patch.hh, spheredata_6patch.cc, spheredata_6patch.hh, sync.cc /SphericalSlice/src/IO/ attributes.h, io.h, io_ascii.h, io_hdf5.h, make.code.defn /SphericalSlice/src/sYlm/ make.code.defn, sYlm.cc, sYlm.h Log: SphericalSlice: A sophisticated version of SphericalSurface. File Changes: Directory: /SphericalSlice/src/IO/ ================================== File [added]: attributes.h Delta lines: <Unavailable> File [added]: io.h Delta lines: <Unavailable> File [added]: io_ascii.h Delta lines: <Unavailable> File [added]: io_hdf5.h Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /SphericalSlice/ =========================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /SphericalSlice/doc/ =============================== File [added]: Doxyfile Delta lines: <Unavailable> File [added]: documentation.tex Delta lines: <Unavailable> File [added]: refman.pdf Delta lines: <Unavailable> Directory: /SphericalSlice/mma/ =============================== File [added]: 1patch.nb Delta lines: <Unavailable> File [added]: 6patch.nb Delta lines: <Unavailable> Directory: /SphericalSlice/src/sYlm/ ==================================== File [added]: make.code.defn Delta lines: <Unavailable> File [added]: sYlm.cc Delta lines: <Unavailable> File [added]: sYlm.h Delta lines: <Unavailable> Directory: /SphericalSlice/src/ =============================== File [added]: check.cc Delta lines: <Unavailable> File [added]: commstack.cc Delta lines: <Unavailable> File [added]: commstack.hh Delta lines: <Unavailable> File [added]: contract.cc Delta lines: <Unavailable> File [added]: int_coeffs.hh Delta lines: <Unavailable> File [added]: integrate.cc Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: paramcheck.cc Delta lines: <Unavailable> File [added]: register.cc Delta lines: <Unavailable> File [added]: setup.cc Delta lines: <Unavailable> File [added]: setup.hh Delta lines: <Unavailable> File [added]: slices.cc Delta lines: <Unavailable> File [added]: slices.hh Delta lines: <Unavailable> File [added]: spheredata.cc Delta lines: <Unavailable> File [added]: spheredata.hh Delta lines: <Unavailable> File [added]: spheredata_1patch.cc Delta lines: <Unavailable> File [added]: spheredata_1patch.hh Delta lines: <Unavailable> File [added]: spheredata_2patch.cc Delta lines: <Unavailable> File [added]: spheredata_2patch.hh Delta lines: <Unavailable> File [added]: spheredata_6patch.cc Delta lines: <Unavailable> File [added]: spheredata_6patch.hh Delta lines: <Unavailable> File [added]: sync.cc Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] incoming/PITTNullCode/ (Rev. 33)
by reisswig＠tapir.caltech.edu 13 Jun '11

13 Jun '11

User: reisswig Date: 2011/06/13 11:55 AM Added: /PITTNullCode/ README /PITTNullCode/DebugBondiR/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/DebugBondiR/doc/ documentation.tex /PITTNullCode/DebugBondiR/src/ ComputeBondiR.F90, make.code.defn /PITTNullCode/HDF5Extra/ HDF5Extra.sh, README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/HDF5Extra/src/ make.code.defn /PITTNullCode/NonCactusVer/ interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NonCactusVer/src/ make.code.defn /PITTNullCode/NonCactusVer/src/PittNullCode/ ForwardBackward-HOWTO, Makefile, StereoIntegration.f90, TS_et_d.f, affine.f90, ascii_io.f90, axiio.f90, boundary.f90, boundary_robust.f90, calling.tree, dpsort.f, estimateTX.f90, evolve2D.f90, geomview.f90, gridtranslator.f90, hdriver.f90, hio.f90, horizon.f90, horizon_eth.f90, makemake, model.f90, model.in, news2.f90, null.f90, null.in, null_analytic.f90, null_cfl_test.f90, null_code.f90, null_eth.f90, null_evol.f90, null_grid.f90, null_hyper_beta.f90, null_hyper_u.f90, null_hyper_w.f90, null_interp.f90, null_interp2.f90, null_io.f90, null_mask.f90, null_params.f90, null_vars.f90, numservicef90.f90, odedriver.f90, particle.f90, particle.in, rksuite.f, rksuite_90.f90, set_precision.f90, spheroid_initial.f90 /PITTNullCode/Null2DWave/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/Null2DWave/doc/ documentation.tex /PITTNullCode/Null2DWave/par/ test.par, testl1m1_r65_NI.par, testl2m2_r65ND.par /PITTNullCode/Null2DWave/src/ Null2DWave_Analysis.F90, Null2DWave_Bdry.F90, Null2DWave_DissipMask.F90, Null2DWave_Eth2.F90, Null2DWave_Evol.F90, Null2DWave_IMask.F90, Null2DWave_Init.F90, Null2DWave_MoLRegister.c, Null2DWave_Mod.F90, Null2DWave_Sync.F90, make.code.defn, my_YlmYst.mw /PITTNullCode/NullConstr/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullConstr/doc/ documentation.tex /PITTNullCode/NullConstr/par/ constr.par, constraints.par /PITTNullCode/NullConstr/src/ NullConstr_AuxConstr.F90, NullConstr_Driver.F90, NullConstr_R00.F90, NullConstr_R01.F90, NullConstr_R0A.F90, NullConstr_Util.F90, make.code.defn /PITTNullCode/NullConstr/src/maple-files/ READ.me, constraint.map, constraint.out, gamma.map, gamma.out, k.map /PITTNullCode/NullDecomp/ README, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullDecomp/doc/ documentation.tex /PITTNullCode/NullDecomp/src/ NullDecomp_IO.F90, NullDecomp_IO.F90.bak, NullDecomp_SpinDecomp.F90, NullDecomp_Vars.F90, NullDecomp_sYlm.F90, Startup.F90, make.code.defn, make.code.deps /PITTNullCode/NullEvolve/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullEvolve/doc/ documentation.tex /PITTNullCode/NullEvolve/par/ test_Jdissip.par, test_circular.par /PITTNullCode/NullEvolve/src/ NullEvol_Bdry.F90, NullEvol_Diag.F90, NullEvol_DiagMod.F90, NullEvol_DissipMask.F90, NullEvol_Evol.F90, NullEvol_InitialData.F90, NullEvol_InitialSlice.F90, NullEvol_Mask.F90, NullEvol_Pulse.F90, NullEvol_Step.F90, NullEvol_TimeStep.F90, NullEvol_cfl_test.F90, NullEvol_hyper_beta.F90, NullEvol_hyper_q.F90, NullEvol_hyper_u.F90, NullEvol_hyper_w.F90, NullEvol_sYlm.F90, gauge.h, make.code.defn /PITTNullCode/NullEvolve/src/evolQU/ evolB.f90, evolQU.f90, redo.sh /PITTNullCode/NullExact/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullExact/doc/ documentation.tex /PITTNullCode/NullExact/par/ random.par /PITTNullCode/NullExact/src/ NullExact_AnaFields.F90, NullExact_Analytic.F90, NullExact_Boundary.F90, NullExact_Error.F90, NullExact_InitMoL.F90, NullExact_Initial.F90, NullExact_MoLRegister.c, NullExact_PreCalcYlm.F90, NullExact_dotNewsMoL.F90, make.code.defn /PITTNullCode/NullGrid/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullGrid/par/ DemoInterp.par /PITTNullCode/NullGrid/src/ NullGrid_InitCoord.F90, NullGrid_RegisterRanges.c, NullGrid_Setup.F90, NullGrid_Vars.F90, make.code.defn /PITTNullCode/NullGrid/src/ret/ InitCoord.F90, InitCoord.c, InterpCoord.F90_eth1, InterpCoord.F90_eth2_failed, InterpCoord.F90_work, getset.f90 /PITTNullCode/NullInterp/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullInterp/src/ NullInterp.F90, NullInterp_Deriv.F90, NullInterp_Eth.F90, NullInterp_Interp.F90, NullInterp_InterpUtil.F90, NullInterp_MaskInit.F90, NullInterp_ParamCheck.F90, NullInterp_TempInit.F90, ethdef.maple, make.code.defn /PITTNullCode/NullInterp/test/ news_test_c.par, null_test_c.par, scalar_test_c.par /PITTNullCode/NullInterpTest/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullInterpTest/doc/ documentation.tex /PITTNullCode/NullInterpTest/par/ conv_check_Cn.agr, e-1_conv.agr, e-2_conv.agr, e1_conv.agr, e2_conv.agr, err_cbcn.agr, eth_conv.agr, nullvars_conv.agr, t.par, test.par /PITTNullCode/NullInterpTest/src/ NullInterpTest_Test.F90, NullInterpTest_sYlm.F90, make.code.defn /PITTNullCode/NullInterpTest/src/scripts/ makemake /PITTNullCode/NullNews/ configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullNews/src/ NullNews_Bondi.F90, NullNews_CalcPsi4.F90, NullNews_GetNews.F90, NullNews_Init.F90, NullNews_Omega.F90, NullNews_ScriUtil.F90, NullNews_ScriVals.F90, NullNews_SphHarmDecomp.F90, NullNews_YlmScri.F90, NullNews_interp_to_constant_uBondi.cc, make.code.defn, set_pointers.h, set_pointers.h.old /PITTNullCode/NullNews/src/maple-files/ J_uPsi4.map, d2Udl2.map, jeff-weyl.f90, jeff-weyl.map, jeff-weyl.out, ruleFort.map /PITTNullCode/NullNewsAdv/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullNewsAdv/par/ t1d2.par, t1ev.par /PITTNullCode/NullNewsAdv/src/ NullNewsAdv_Init.F90, NullNewsAdv_Omega.F90, make.code.defn /PITTNullCode/NullNewsIO/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullNewsIO/src/ NullNewsIO_DumpNewsToFlexIOFile.c, NullNewsIO_Observers.F90, NullNewsIO_Test.F90, NullNewsIO_WriteToFile.c, NullNewsIO_sYlm.F90, make.code.defn, make.code.deps, make.configuration.defn, make.configuration.deps /PITTNullCode/NullNewsIO/src/util/ Decomp.f90, driver.f90, ieeeread.f90, interpnews.f90, newsread.f90, sYlm.f90 /PITTNullCode/NullNewsPostProcess/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullNewsPostProcess/doc/ documentation.tex /PITTNullCode/NullNewsPostProcess/par/ postprocess.par, postprocess2.par /PITTNullCode/NullNewsPostProcess/scripts/ GenAnalysis.sh, GenParamFile.sh /PITTNullCode/NullNewsPostProcess/src/ Cleanup.F90, NullNewsPostProcess_CalcCoeffs.F90, NullNewsPostProcess_Dump.F90, NullNewsPostProcess_SpinDecomp.F90, NullNewsPostProcess_Vars.F90, NullNewsPostProcess_ieee_read.F90, NullNewsPostProcess_interpnews.F90, NullNewsPostProcess_newsread.F90, NullNewsPostProcess_sYlm.F90, Startup.F90, String.c, Terminate.c, comm.c, make.code.defn, make.code.deps, make.configuration.defn, make.configuration.deps /PITTNullCode/NullNewsPostProcess/src/util/ GetGridSize.f90, ieee_merge.c /PITTNullCode/NullNewsTest/ README, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullNewsTest/doc/ documentation.tex /PITTNullCode/NullNewsTest/math/ distort.m, distort_u_q_p.nb /PITTNullCode/NullNewsTest/par/ distort_u_q_p.par, nullchk2.par /PITTNullCode/NullNewsTest/src/ NullNewsTest_Err.F90, NullNewsTest_Init.F90, NullNewsTest_Setup.c, NullNewsTest_Step.F90, NullNewsTest_rtmetric.F90, make.code.defn, make.code.deps /PITTNullCode/NullPsi4IO/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullPsi4IO/src/ NullPsi4IO_DumpPsi4ToFlexIOFile.c, NullPsi4IO_Observers.F90, NullPsi4IO_Test.F90, NullPsi4IO_WriteToFile.c, NullPsi4IO_sYlm.F90, make.code.defn, make.code.deps, make.configuration.defn, make.configuration.deps /PITTNullCode/NullPsi4IO/src/util/ Decomp.f90, Psi4read.f90, driver.f90, ieeeread.f90, interpPsi4.f90, sYlm.f90 /PITTNullCode/NullPsiInt/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullPsiInt/doc/ documentation.tex /PITTNullCode/NullPsiInt/src/ NullPsiInt_InitMoL.F90, NullPsiInt_MoLRegister.c, NullPsiInt_dotNewsPsi.F90, make.code.defn /PITTNullCode/NullPsiIntIO/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullPsiIntIO/src/ NullPsiIntIO_DumpPsiIntToFlexIOFile.c, NullPsiIntIO_Observers.F90, NullPsiIntIO_Test.F90, NullPsiIntIO_WriteToFile.c, NullPsiIntIO_sYlm.F90, make.code.defn, make.code.deps, make.configuration.defn, make.configuration.deps /PITTNullCode/NullPsiIntIO/src/util/ Decomp.f90, PsiIntread.f90, driver.f90, ieeeread.f90, interpPsiInt.f90, sYlm.f90 /PITTNullCode/NullSHRExtract/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullSHRExtract/doc/ documentation.tex /PITTNullCode/NullSHRExtract/src/ NullSHRE_AnaSchMetricRecon.F90, NullSHRE_ErrorSch.F90, NullSHRE_ExtractWT.F90, NullSHRE_FullMetricRecon.F90, NullSHRE_InitialWT.F90, NullSHRE_MetricReconPast.F90, NullSHRE_Pointing.F90, NullSHRE_PreCalcNSYlm.F90, NullSHRE_SchMetricRecon.F90, NullSHRE_SchPreCalc.F90, NullSHRE_SchReadClm.F90, NullSHRE_WTSphHarm.F90, NullSHRE_mod4Metric.F90, NullSHRE_modAnaCoord.F90, NullSHRE_modBoundary.F90, NullSHRE_modEta.F90, NullSHRE_modGFdef.F90, NullSHRE_modInverse.F90, NullSHRE_modNullDir.F90, NullSHRE_modRadius.F90, NullSHRE_modSchClm.F90, NullSHRE_modSchYlm.F90, NullSHRE_modVars.F90, NullSHRE_modsinSchClm.F90, NullSHRE_modvibSchClm.F90, make.code.defn, wt_SphHarmDecomp.F90 /PITTNullCode/NullSHRExtract/src/Maple/ B10_3.f90, B11_3.f90, B1m1_3.f90, C00_3.f90, C20_3.f90, C21_3.f90, C22_3.f90, C2m1_3.f90, C2m2_3.f90, Cpq_gN.f90, Cxyz_g.f90, ExtractionSchwarzIEF.mw, ExtractionSchwarzIEF_MAS.bak, JacobianCtoW.mw, Schwarz3Ylm_NS.mw, SchwarzYlm_NS.mw, SchwarzYlm_NS_Deriv.mw, SchwarzYlm_NS_MAS.bak, WT_dRg3N.f90, WT_dRg3S.f90, WT_g3N.f90, WT_g3S.f90, WT_gN.f90, WT_gS.f90, Y00.f90, Y10.f90, Y11.f90, Y1m1.f90, Y20.f90, Y21.f90, Y22.f90, Y2m1.f90, Y2m2.f90, dRB10_3.f90, dRB11_3.f90, dRB1m1_3.f90, dRC00_3.f90, dRC20_3.f90, dRC21_3.f90, dRC22_3.f90, dRC2m1_3.f90, dRC2m2_3.f90, g_N.f90 /PITTNullCode/NullSHRExtract/src/draft/ NullSHRE_AnaSchWT.h, NullSHRE_BrokenAnalysis.F90, NullSHRE_FullMetricReconBK.F90, NullSHRE_JInitial.F90, NullSHRE_OldSchMetricRecon.F90, NullSHRE_SchMetricReconBK.F90, NullSHRE_SchPreCalc.F90, NullSHRE_SchReadClm.F90, NullSHRE_modBrokenAnalytic.F90, NullSHRE_modSchClm.F90, NullSHRE_modSchYlm.F90 /PITTNullCode/NullSHRExtract/src/include/ NullSHRE_PointingLib.h, NullSHRE_PointingNull.h /PITTNullCode/NullSHRExtract/src/include/Maple/ 1YlmDeriv.mw, B10_3.f90, B11_3.f90, B1m1_3.f90, C00_3.f90, C20_3.f90, C21_3.f90, C22_3.f90, C2m1_3.f90, C2m2_3.f90, Cpq_gN.f90, Cxyz_dxg.f90, Cxyz_dzg.f90, Cxyz_g.f90, DxSchwarz3Ylm_NS.mw, JacobianCtoW.mw, JacobianCtoW_MAS.bak, MinkYlm_NS.mw, Schwarz3Ylm_NS.mw, Schwarz3Ylm_NS_MAS.bak, SchwarzYlm_0.mw, SchwarzYlm_NS.mw, SchwarzYlm_NS_Deriv.mw, SchwarzYlm_NS_Deriv_MAS.bak, SchwarzYlm_NS_MAS.bak, WT_dRg3N.f90, WT_dRg3S.f90, WT_dxgS.f90, WT_g3N.f90, WT_g3S.f90, WT_gN.f90, WT_gS.f90, Y00.f90, Y10.f90, Y11.f90, Y1m1.f90, Y20.f90, Y21.f90, Y22.f90, Y2m1.f90, Y2m2.f90, dRB10_3.f90, dRB11_3.f90, dRB1m1_3.f90, dRC00_3.f90, dRC20_3.f90, dRC21_3.f90, dRC22_3.f90, dRC2m1_3.f90, dRC2m2_3.f90, dxC00.f90, dxC11.f90, dxC20.f90, dxC21.f90, dxC22.f90, dxC2m1.f90, dxC2m2.f90, dxC30.f90, dxC31.f90, dxC32.f90, dxC33.f90, dxC3m1.f90, dxC3m2.f90, dxC3m3.f90, dxCpq_gN.f90, dxY10.f90, g_N.f90 /PITTNullCode/NullScalar/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullScalar/doc/ documentation.tex /PITTNullCode/NullScalar/par/ test.par /PITTNullCode/NullScalar/src/ NullScalar_Analysis.F90, NullScalar_Bdry.F90, NullScalar_Evol.F90, NullScalar_Init.F90, NullScalar_Loop.F90, NullScalar_Step.F90, NullScalar_Sync.F90, NullScalar_Testbed.F90, make.code.defn /PITTNullCode/NullScalar/src/MAPLE/ SphWave.mw /PITTNullCode/NullScalar/test/ t1.par, t2.par /PITTNullCode/NullScalar/test/t1/ /PITTNullCode/NullVars/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/NullVars/doc/ documentation.tex /PITTNullCode/NullVars/src/ make.code.defn /PITTNullCode/SphericalHarmonicDecomp/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/SphericalHarmonicDecomp/doc/ documentation.tex /PITTNullCode/SphericalHarmonicDecomp/src/ Chebyshev.cc, Chebyshev.hh, ExtractMetric.cc, Legendre.cc, Legendre.hh, SphericalHarmonicDecomp.h, alias.cc, decomp.cc, decomp.hh, h5read.cc, make.code.defn, matrix.cc, matrix.hh, myassert.h, sYlm.cc, sYlm.hh /PITTNullCode/SphericalHarmonicDecomp/src/TestLeastSquare/ Chebyshev.cc, Chebyshev.hh, Legendre.cc, Legendre.hh, Makefile, TestLeastSquare.h, decomp.cc, decomp.hh, driver.cc, matrix.cc, matrix.hh, myassert.h, sYlm.cc, sYlm.hh, testmatrix.cc /PITTNullCode/SphericalHarmonicDecomp/src/util/ Legendre.cc, Legendre.hh, Makefile, h5read.cc, run.sh, sYlm.cc, sYlm.hh /PITTNullCode/SphericalHarmonicRecon/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/SphericalHarmonicRecon/codes/ Makefile, ascii_output.c, dampfilter.c, fftwfilter.c, filter.c, findlast.c, printtime.c, readmeta.c, setmeta.c, tderiv.c, testfilter.c, util.c, util.h /PITTNullCode/SphericalHarmonicRecon/doc/ documentation.tex /PITTNullCode/SphericalHarmonicRecon/par/ /PITTNullCode/SphericalHarmonicRecon/src/ Chebyshev.cc, Chebyshev.hh, Recon.cc, make.code.defn /PITTNullCode/SphericalHarmonicRecon/test/ /PITTNullCode/SphericalHarmonicReconASCII/ configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/SphericalHarmonicReconASCII/src/ get.cc, make.code.defn, read.cc, sph_database.cc, sph_database.hh, startup.cc, vars.hh Log: The PITTNullCode and thorns capable of doing Cauchy-characteristic extraction (CCE). File Changes: Directory: /PITTNullCode/DebugBondiR/ ===================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/HDF5Extra/ =================================== File [added]: HDF5Extra.sh Delta lines: <Unavailable> File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullScalar/src/MAPLE/ ============================================== File [added]: SphWave.mw Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/src/Maple/ ================================================== File [added]: B10_3.f90 Delta lines: <Unavailable> File [added]: B11_3.f90 Delta lines: <Unavailable> File [added]: B1m1_3.f90 Delta lines: <Unavailable> File [added]: C00_3.f90 Delta lines: <Unavailable> File [added]: C20_3.f90 Delta lines: <Unavailable> File [added]: C21_3.f90 Delta lines: <Unavailable> File [added]: C22_3.f90 Delta lines: <Unavailable> File [added]: C2m1_3.f90 Delta lines: <Unavailable> File [added]: C2m2_3.f90 Delta lines: <Unavailable> File [added]: Cpq_gN.f90 Delta lines: <Unavailable> File [added]: Cxyz_g.f90 Delta lines: <Unavailable> File [added]: ExtractionSchwarzIEF.mw Delta lines: <Unavailable> File [added]: ExtractionSchwarzIEF_MAS.bak Delta lines: <Unavailable> File [added]: JacobianCtoW.mw Delta lines: <Unavailable> File [added]: Schwarz3Ylm_NS.mw Delta lines: <Unavailable> File [added]: SchwarzYlm_NS.mw Delta lines: <Unavailable> File [added]: SchwarzYlm_NS_Deriv.mw Delta lines: <Unavailable> File [added]: SchwarzYlm_NS_MAS.bak Delta lines: <Unavailable> File [added]: WT_dRg3N.f90 Delta lines: <Unavailable> File [added]: WT_dRg3S.f90 Delta lines: <Unavailable> File [added]: WT_g3N.f90 Delta lines: <Unavailable> File [added]: WT_g3S.f90 Delta lines: <Unavailable> File [added]: WT_gN.f90 Delta lines: <Unavailable> File [added]: WT_gS.f90 Delta lines: <Unavailable> File [added]: Y00.f90 Delta lines: <Unavailable> File [added]: Y10.f90 Delta lines: <Unavailable> File [added]: Y11.f90 Delta lines: <Unavailable> File [added]: Y1m1.f90 Delta lines: <Unavailable> File [added]: Y20.f90 Delta lines: <Unavailable> File [added]: Y21.f90 Delta lines: <Unavailable> File [added]: Y22.f90 Delta lines: <Unavailable> File [added]: Y2m1.f90 Delta lines: <Unavailable> File [added]: Y2m2.f90 Delta lines: <Unavailable> File [added]: dRB10_3.f90 Delta lines: <Unavailable> File [added]: dRB11_3.f90 Delta lines: <Unavailable> File [added]: dRB1m1_3.f90 Delta lines: <Unavailable> File [added]: dRC00_3.f90 Delta lines: <Unavailable> File [added]: dRC20_3.f90 Delta lines: <Unavailable> File [added]: dRC21_3.f90 Delta lines: <Unavailable> File [added]: dRC22_3.f90 Delta lines: <Unavailable> File [added]: dRC2m1_3.f90 Delta lines: <Unavailable> File [added]: dRC2m2_3.f90 Delta lines: <Unavailable> File [added]: g_N.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/src/include/Maple/ ========================================================== File [added]: 1YlmDeriv.mw Delta lines: <Unavailable> File [added]: B10_3.f90 Delta lines: <Unavailable> File [added]: B11_3.f90 Delta lines: <Unavailable> File [added]: B1m1_3.f90 Delta lines: <Unavailable> File [added]: C00_3.f90 Delta lines: <Unavailable> File [added]: C20_3.f90 Delta lines: <Unavailable> File [added]: C21_3.f90 Delta lines: <Unavailable> File [added]: C22_3.f90 Delta lines: <Unavailable> File [added]: C2m1_3.f90 Delta lines: <Unavailable> File [added]: C2m2_3.f90 Delta lines: <Unavailable> File [added]: Cpq_gN.f90 Delta lines: <Unavailable> File [added]: Cxyz_dxg.f90 Delta lines: <Unavailable> File [added]: Cxyz_dzg.f90 Delta lines: <Unavailable> File [added]: Cxyz_g.f90 Delta lines: <Unavailable> File [added]: DxSchwarz3Ylm_NS.mw Delta lines: <Unavailable> File [added]: JacobianCtoW.mw Delta lines: <Unavailable> File [added]: JacobianCtoW_MAS.bak Delta lines: <Unavailable> File [added]: MinkYlm_NS.mw Delta lines: <Unavailable> File [added]: Schwarz3Ylm_NS.mw Delta lines: <Unavailable> File [added]: Schwarz3Ylm_NS_MAS.bak Delta lines: <Unavailable> File [added]: SchwarzYlm_0.mw Delta lines: <Unavailable> File [added]: SchwarzYlm_NS.mw Delta lines: <Unavailable> File [added]: SchwarzYlm_NS_Deriv.mw Delta lines: <Unavailable> File [added]: SchwarzYlm_NS_Deriv_MAS.bak Delta lines: <Unavailable> File [added]: SchwarzYlm_NS_MAS.bak Delta lines: <Unavailable> File [added]: WT_dRg3N.f90 Delta lines: <Unavailable> File [added]: WT_dRg3S.f90 Delta lines: <Unavailable> File [added]: WT_dxgS.f90 Delta lines: <Unavailable> File [added]: WT_g3N.f90 Delta lines: <Unavailable> File [added]: WT_g3S.f90 Delta lines: <Unavailable> File [added]: WT_gN.f90 Delta lines: <Unavailable> File [added]: WT_gS.f90 Delta lines: <Unavailable> File [added]: Y00.f90 Delta lines: <Unavailable> File [added]: Y10.f90 Delta lines: <Unavailable> File [added]: Y11.f90 Delta lines: <Unavailable> File [added]: Y1m1.f90 Delta lines: <Unavailable> File [added]: Y20.f90 Delta lines: <Unavailable> File [added]: Y21.f90 Delta lines: <Unavailable> File [added]: Y22.f90 Delta lines: <Unavailable> File [added]: Y2m1.f90 Delta lines: <Unavailable> File [added]: Y2m2.f90 Delta lines: <Unavailable> File [added]: dRB10_3.f90 Delta lines: <Unavailable> File [added]: dRB11_3.f90 Delta lines: <Unavailable> File [added]: dRB1m1_3.f90 Delta lines: <Unavailable> File [added]: dRC00_3.f90 Delta lines: <Unavailable> File [added]: dRC20_3.f90 Delta lines: <Unavailable> File [added]: dRC21_3.f90 Delta lines: <Unavailable> File [added]: dRC22_3.f90 Delta lines: <Unavailable> File [added]: dRC2m1_3.f90 Delta lines: <Unavailable> File [added]: dRC2m2_3.f90 Delta lines: <Unavailable> File [added]: dxC00.f90 Delta lines: <Unavailable> File [added]: dxC11.f90 Delta lines: <Unavailable> File [added]: dxC20.f90 Delta lines: <Unavailable> File [added]: dxC21.f90 Delta lines: <Unavailable> File [added]: dxC22.f90 Delta lines: <Unavailable> File [added]: dxC2m1.f90 Delta lines: <Unavailable> File [added]: dxC2m2.f90 Delta lines: <Unavailable> File [added]: dxC30.f90 Delta lines: <Unavailable> File [added]: dxC31.f90 Delta lines: <Unavailable> File [added]: dxC32.f90 Delta lines: <Unavailable> File [added]: dxC33.f90 Delta lines: <Unavailable> File [added]: dxC3m1.f90 Delta lines: <Unavailable> File [added]: dxC3m2.f90 Delta lines: <Unavailable> File [added]: dxC3m3.f90 Delta lines: <Unavailable> File [added]: dxCpq_gN.f90 Delta lines: <Unavailable> File [added]: dxY10.f90 Delta lines: <Unavailable> File [added]: g_N.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/NonCactusVer/ ====================================== File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/Null2DWave/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullConstr/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullDecomp/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullEvolve/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullExact/ =================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullGrid/ ================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterp/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterpTest/ ======================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullNews/ ================================== File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsAdv/ ===================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsIO/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsPostProcess/ ============================================= File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsTest/ ====================================== File [added]: README Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsi4IO/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsiInt/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsiIntIO/ ====================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/ ======================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullScalar/ ==================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/NullVars/ ================================== File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/ ========================= File [added]: README Delta lines: <Unavailable> Directory: /PITTNullCode/NonCactusVer/src/PittNullCode/ ======================================================= File [added]: ForwardBackward-HOWTO Delta lines: <Unavailable> File [added]: Makefile Delta lines: <Unavailable> File [added]: StereoIntegration.f90 Delta lines: <Unavailable> File [added]: TS_et_d.f Delta lines: <Unavailable> File [added]: affine.f90 Delta lines: <Unavailable> File [added]: ascii_io.f90 Delta lines: <Unavailable> File [added]: axiio.f90 Delta lines: <Unavailable> File [added]: boundary.f90 Delta lines: <Unavailable> File [added]: boundary_robust.f90 Delta lines: <Unavailable> File [added]: calling.tree Delta lines: <Unavailable> File [added]: dpsort.f Delta lines: <Unavailable> File [added]: estimateTX.f90 Delta lines: <Unavailable> File [added]: evolve2D.f90 Delta lines: <Unavailable> File [added]: geomview.f90 Delta lines: <Unavailable> File [added]: gridtranslator.f90 Delta lines: <Unavailable> File [added]: hdriver.f90 Delta lines: <Unavailable> File [added]: hio.f90 Delta lines: <Unavailable> File [added]: horizon.f90 Delta lines: <Unavailable> File [added]: horizon_eth.f90 Delta lines: <Unavailable> File [added]: makemake Delta lines: <Unavailable> File [added]: model.f90 Delta lines: <Unavailable> File [added]: model.in Delta lines: <Unavailable> File [added]: news2.f90 Delta lines: <Unavailable> File [added]: null.f90 Delta lines: <Unavailable> File [added]: null.in Delta lines: <Unavailable> File [added]: null_analytic.f90 Delta lines: <Unavailable> File [added]: null_cfl_test.f90 Delta lines: <Unavailable> File [added]: null_code.f90 Delta lines: <Unavailable> File [added]: null_eth.f90 Delta lines: <Unavailable> File [added]: null_evol.f90 Delta lines: <Unavailable> File [added]: null_grid.f90 Delta lines: <Unavailable> File [added]: null_hyper_beta.f90 Delta lines: <Unavailable> File [added]: null_hyper_u.f90 Delta lines: <Unavailable> File [added]: null_hyper_w.f90 Delta lines: <Unavailable> File [added]: null_interp.f90 Delta lines: <Unavailable> File [added]: null_interp2.f90 Delta lines: <Unavailable> File [added]: null_io.f90 Delta lines: <Unavailable> File [added]: null_mask.f90 Delta lines: <Unavailable> File [added]: null_params.f90 Delta lines: <Unavailable> File [added]: null_vars.f90 Delta lines: <Unavailable> File [added]: numservicef90.f90 Delta lines: <Unavailable> File [added]: odedriver.f90 Delta lines: <Unavailable> File [added]: particle.f90 Delta lines: <Unavailable> File [added]: particle.in Delta lines: <Unavailable> File [added]: rksuite.f Delta lines: <Unavailable> File [added]: rksuite_90.f90 Delta lines: <Unavailable> File [added]: set_precision.f90 Delta lines: <Unavailable> File [added]: spheroid_initial.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicDecomp/ ================================================= File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/ ================================================ File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicReconASCII/ ===================================================== File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicDecomp/src/TestLeastSquare/ ===================================================================== File [added]: Chebyshev.cc Delta lines: <Unavailable> File [added]: Chebyshev.hh Delta lines: <Unavailable> File [added]: Legendre.cc Delta lines: <Unavailable> File [added]: Legendre.hh Delta lines: <Unavailable> File [added]: Makefile Delta lines: <Unavailable> File [added]: TestLeastSquare.h Delta lines: <Unavailable> File [added]: decomp.cc Delta lines: <Unavailable> File [added]: decomp.hh Delta lines: <Unavailable> File [added]: driver.cc Delta lines: <Unavailable> File [added]: matrix.cc Delta lines: <Unavailable> File [added]: matrix.hh Delta lines: <Unavailable> File [added]: myassert.h Delta lines: <Unavailable> File [added]: sYlm.cc Delta lines: <Unavailable> File [added]: sYlm.hh Delta lines: <Unavailable> File [added]: testmatrix.cc Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/codes/ ====================================================== File [added]: Makefile Delta lines: <Unavailable> File [added]: ascii_output.c Delta lines: <Unavailable> File [added]: dampfilter.c Delta lines: <Unavailable> File [added]: fftwfilter.c Delta lines: <Unavailable> File [added]: filter.c Delta lines: <Unavailable> File [added]: findlast.c Delta lines: <Unavailable> File [added]: printtime.c Delta lines: <Unavailable> File [added]: readmeta.c Delta lines: <Unavailable> File [added]: setmeta.c Delta lines: <Unavailable> File [added]: tderiv.c Delta lines: <Unavailable> File [added]: testfilter.c Delta lines: <Unavailable> File [added]: util.c Delta lines: <Unavailable> File [added]: util.h Delta lines: <Unavailable> Directory: /PITTNullCode/DebugBondiR/doc/ ========================================= File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/Null2DWave/doc/ ======================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullConstr/doc/ ======================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullDecomp/doc/ ======================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullEvolve/doc/ ======================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullExact/doc/ ======================================= File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterpTest/doc/ ============================================ File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsPostProcess/doc/ ================================================= File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsTest/doc/ ========================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsiInt/doc/ ======================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/doc/ ============================================ File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullScalar/doc/ ======================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullVars/doc/ ====================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicDecomp/doc/ ===================================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/doc/ ==================================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/src/draft/ ================================================== File [added]: NullSHRE_AnaSchWT.h Delta lines: <Unavailable> File [added]: NullSHRE_BrokenAnalysis.F90 Delta lines: <Unavailable> File [added]: NullSHRE_FullMetricReconBK.F90 Delta lines: <Unavailable> File [added]: NullSHRE_JInitial.F90 Delta lines: <Unavailable> File [added]: NullSHRE_OldSchMetricRecon.F90 Delta lines: <Unavailable> File [added]: NullSHRE_SchMetricReconBK.F90 Delta lines: <Unavailable> File [added]: NullSHRE_SchPreCalc.F90 Delta lines: <Unavailable> File [added]: NullSHRE_SchReadClm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modBrokenAnalytic.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modSchClm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modSchYlm.F90 Delta lines: <Unavailable> Directory: /PITTNullCode/NullEvolve/src/evolQU/ =============================================== File [added]: evolB.f90 Delta lines: <Unavailable> File [added]: evolQU.f90 Delta lines: <Unavailable> File [added]: redo.sh Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/src/include/ ==================================================== File [added]: NullSHRE_PointingLib.h Delta lines: <Unavailable> File [added]: NullSHRE_PointingNull.h Delta lines: <Unavailable> Directory: /PITTNullCode/NullConstr/src/maple-files/ ==================================================== File [added]: READ.me Delta lines: <Unavailable> File [added]: constraint.map Delta lines: <Unavailable> File [added]: constraint.out Delta lines: <Unavailable> File [added]: gamma.map Delta lines: <Unavailable> File [added]: gamma.out Delta lines: <Unavailable> File [added]: k.map Delta lines: <Unavailable> Directory: /PITTNullCode/NullNews/src/maple-files/ ================================================== File [added]: J_uPsi4.map Delta lines: <Unavailable> File [added]: d2Udl2.map Delta lines: <Unavailable> File [added]: jeff-weyl.f90 Delta lines: <Unavailable> File [added]: jeff-weyl.map Delta lines: <Unavailable> File [added]: jeff-weyl.out Delta lines: <Unavailable> File [added]: ruleFort.map Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsTest/math/ =========================================== File [added]: distort.m Delta lines: <Unavailable> File [added]: distort_u_q_p.nb Delta lines: <Unavailable> Directory: /PITTNullCode/Null2DWave/par/ ======================================== File [added]: test.par Delta lines: <Unavailable> File [added]: testl1m1_r65_NI.par Delta lines: <Unavailable> File [added]: testl2m2_r65ND.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullConstr/par/ ======================================== File [added]: constr.par Delta lines: <Unavailable> File [added]: constraints.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullEvolve/par/ ======================================== File [added]: test_Jdissip.par Delta lines: <Unavailable> File [added]: test_circular.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullExact/par/ ======================================= File [added]: random.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullGrid/par/ ====================================== File [added]: DemoInterp.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterpTest/par/ ============================================ File [added]: conv_check_Cn.agr Delta lines: <Unavailable> File [added]: e-1_conv.agr Delta lines: <Unavailable> File [added]: e-2_conv.agr Delta lines: <Unavailable> File [added]: e1_conv.agr Delta lines: <Unavailable> File [added]: e2_conv.agr Delta lines: <Unavailable> File [added]: err_cbcn.agr Delta lines: <Unavailable> File [added]: eth_conv.agr Delta lines: <Unavailable> File [added]: nullvars_conv.agr Delta lines: <Unavailable> File [added]: t.par Delta lines: <Unavailable> File [added]: test.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsAdv/par/ ========================================= File [added]: t1d2.par Delta lines: <Unavailable> File [added]: t1ev.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsPostProcess/par/ ================================================= File [added]: postprocess.par Delta lines: <Unavailable> File [added]: postprocess2.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsTest/par/ ========================================== File [added]: distort_u_q_p.par Delta lines: <Unavailable> File [added]: nullchk2.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullScalar/par/ ======================================== File [added]: test.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullGrid/src/ret/ ========================================== File [added]: InitCoord.F90 Delta lines: <Unavailable> File [added]: InitCoord.c Delta lines: <Unavailable> File [added]: InterpCoord.F90_eth1 Delta lines: <Unavailable> File [added]: InterpCoord.F90_eth2_failed Delta lines: <Unavailable> File [added]: InterpCoord.F90_work Delta lines: <Unavailable> File [added]: getset.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterpTest/src/scripts/ ==================================================== File [added]: makemake Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsPostProcess/scripts/ ===================================================== File [added]: GenAnalysis.sh Delta lines: <Unavailable> File [added]: GenParamFile.sh Delta lines: <Unavailable> Directory: /PITTNullCode/DebugBondiR/src/ ========================================= File [added]: ComputeBondiR.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/HDF5Extra/src/ ======================================= File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NonCactusVer/src/ ========================================== File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/Null2DWave/src/ ======================================== File [added]: Null2DWave_Analysis.F90 Delta lines: <Unavailable> File [added]: Null2DWave_Bdry.F90 Delta lines: <Unavailable> File [added]: Null2DWave_DissipMask.F90 Delta lines: <Unavailable> File [added]: Null2DWave_Eth2.F90 Delta lines: <Unavailable> File [added]: Null2DWave_Evol.F90 Delta lines: <Unavailable> File [added]: Null2DWave_IMask.F90 Delta lines: <Unavailable> File [added]: Null2DWave_Init.F90 Delta lines: <Unavailable> File [added]: Null2DWave_MoLRegister.c Delta lines: <Unavailable> File [added]: Null2DWave_Mod.F90 Delta lines: <Unavailable> File [added]: Null2DWave_Sync.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: my_YlmYst.mw Delta lines: <Unavailable> Directory: /PITTNullCode/NullConstr/src/ ======================================== File [added]: NullConstr_AuxConstr.F90 Delta lines: <Unavailable> File [added]: NullConstr_Driver.F90 Delta lines: <Unavailable> File [added]: NullConstr_R00.F90 Delta lines: <Unavailable> File [added]: NullConstr_R01.F90 Delta lines: <Unavailable> File [added]: NullConstr_R0A.F90 Delta lines: <Unavailable> File [added]: NullConstr_Util.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullDecomp/src/ ======================================== File [added]: NullDecomp_IO.F90 Delta lines: <Unavailable> File [added]: NullDecomp_IO.F90.bak Delta lines: <Unavailable> File [added]: NullDecomp_SpinDecomp.F90 Delta lines: <Unavailable> File [added]: NullDecomp_Vars.F90 Delta lines: <Unavailable> File [added]: NullDecomp_sYlm.F90 Delta lines: <Unavailable> File [added]: Startup.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: make.code.deps Delta lines: <Unavailable> Directory: /PITTNullCode/NullEvolve/src/ ======================================== File [added]: NullEvol_Bdry.F90 Delta lines: <Unavailable> File [added]: NullEvol_Diag.F90 Delta lines: <Unavailable> File [added]: NullEvol_DiagMod.F90 Delta lines: <Unavailable> File [added]: NullEvol_DissipMask.F90 Delta lines: <Unavailable> File [added]: NullEvol_Evol.F90 Delta lines: <Unavailable> File [added]: NullEvol_InitialData.F90 Delta lines: <Unavailable> File [added]: NullEvol_InitialSlice.F90 Delta lines: <Unavailable> File [added]: NullEvol_Mask.F90 Delta lines: <Unavailable> File [added]: NullEvol_Pulse.F90 Delta lines: <Unavailable> File [added]: NullEvol_Step.F90 Delta lines: <Unavailable> File [added]: NullEvol_TimeStep.F90 Delta lines: <Unavailable> File [added]: NullEvol_cfl_test.F90 Delta lines: <Unavailable> File [added]: NullEvol_hyper_beta.F90 Delta lines: <Unavailable> File [added]: NullEvol_hyper_q.F90 Delta lines: <Unavailable> File [added]: NullEvol_hyper_u.F90 Delta lines: <Unavailable> File [added]: NullEvol_hyper_w.F90 Delta lines: <Unavailable> File [added]: NullEvol_sYlm.F90 Delta lines: <Unavailable> File [added]: gauge.h Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullExact/src/ ======================================= File [added]: NullExact_AnaFields.F90 Delta lines: <Unavailable> File [added]: NullExact_Analytic.F90 Delta lines: <Unavailable> File [added]: NullExact_Boundary.F90 Delta lines: <Unavailable> File [added]: NullExact_Error.F90 Delta lines: <Unavailable> File [added]: NullExact_InitMoL.F90 Delta lines: <Unavailable> File [added]: NullExact_Initial.F90 Delta lines: <Unavailable> File [added]: NullExact_MoLRegister.c Delta lines: <Unavailable> File [added]: NullExact_PreCalcYlm.F90 Delta lines: <Unavailable> File [added]: NullExact_dotNewsMoL.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullGrid/src/ ====================================== File [added]: NullGrid_InitCoord.F90 Delta lines: <Unavailable> File [added]: NullGrid_RegisterRanges.c Delta lines: <Unavailable> File [added]: NullGrid_Setup.F90 Delta lines: <Unavailable> File [added]: NullGrid_Vars.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterp/src/ ======================================== File [added]: NullInterp.F90 Delta lines: <Unavailable> File [added]: NullInterp_Deriv.F90 Delta lines: <Unavailable> File [added]: NullInterp_Eth.F90 Delta lines: <Unavailable> File [added]: NullInterp_Interp.F90 Delta lines: <Unavailable> File [added]: NullInterp_InterpUtil.F90 Delta lines: <Unavailable> File [added]: NullInterp_MaskInit.F90 Delta lines: <Unavailable> File [added]: NullInterp_ParamCheck.F90 Delta lines: <Unavailable> File [added]: NullInterp_TempInit.F90 Delta lines: <Unavailable> File [added]: ethdef.maple Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterpTest/src/ ============================================ File [added]: NullInterpTest_Test.F90 Delta lines: <Unavailable> File [added]: NullInterpTest_sYlm.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullNews/src/ ====================================== File [added]: NullNews_Bondi.F90 Delta lines: <Unavailable> File [added]: NullNews_CalcPsi4.F90 Delta lines: <Unavailable> File [added]: NullNews_GetNews.F90 Delta lines: <Unavailable> File [added]: NullNews_Init.F90 Delta lines: <Unavailable> File [added]: NullNews_Omega.F90 Delta lines: <Unavailable> File [added]: NullNews_ScriUtil.F90 Delta lines: <Unavailable> File [added]: NullNews_ScriVals.F90 Delta lines: <Unavailable> File [added]: NullNews_SphHarmDecomp.F90 Delta lines: <Unavailable> File [added]: NullNews_YlmScri.F90 Delta lines: <Unavailable> File [added]: NullNews_interp_to_constant_uBondi.cc Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: set_pointers.h Delta lines: <Unavailable> File [added]: set_pointers.h.old Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsAdv/src/ ========================================= File [added]: NullNewsAdv_Init.F90 Delta lines: <Unavailable> File [added]: NullNewsAdv_Omega.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsIO/src/ ======================================== File [added]: NullNewsIO_DumpNewsToFlexIOFile.c Delta lines: <Unavailable> File [added]: NullNewsIO_Observers.F90 Delta lines: <Unavailable> File [added]: NullNewsIO_Test.F90 Delta lines: <Unavailable> File [added]: NullNewsIO_WriteToFile.c Delta lines: <Unavailable> File [added]: NullNewsIO_sYlm.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: make.code.deps Delta lines: <Unavailable> File [added]: make.configuration.defn Delta lines: <Unavailable> File [added]: make.configuration.deps Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsPostProcess/src/ ================================================= File [added]: Cleanup.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_CalcCoeffs.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_Dump.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_SpinDecomp.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_Vars.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_ieee_read.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_interpnews.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_newsread.F90 Delta lines: <Unavailable> File [added]: NullNewsPostProcess_sYlm.F90 Delta lines: <Unavailable> File [added]: Startup.F90 Delta lines: <Unavailable> File [added]: String.c Delta lines: <Unavailable> File [added]: Terminate.c Delta lines: <Unavailable> File [added]: comm.c Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: make.code.deps Delta lines: <Unavailable> File [added]: make.configuration.defn Delta lines: <Unavailable> File [added]: make.configuration.deps Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsTest/src/ ========================================== File [added]: NullNewsTest_Err.F90 Delta lines: <Unavailable> File [added]: NullNewsTest_Init.F90 Delta lines: <Unavailable> File [added]: NullNewsTest_Setup.c Delta lines: <Unavailable> File [added]: NullNewsTest_Step.F90 Delta lines: <Unavailable> File [added]: NullNewsTest_rtmetric.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: make.code.deps Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsi4IO/src/ ======================================== File [added]: NullPsi4IO_DumpPsi4ToFlexIOFile.c Delta lines: <Unavailable> File [added]: NullPsi4IO_Observers.F90 Delta lines: <Unavailable> File [added]: NullPsi4IO_Test.F90 Delta lines: <Unavailable> File [added]: NullPsi4IO_WriteToFile.c Delta lines: <Unavailable> File [added]: NullPsi4IO_sYlm.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: make.code.deps Delta lines: <Unavailable> File [added]: make.configuration.defn Delta lines: <Unavailable> File [added]: make.configuration.deps Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsiInt/src/ ======================================== File [added]: NullPsiInt_InitMoL.F90 Delta lines: <Unavailable> File [added]: NullPsiInt_MoLRegister.c Delta lines: <Unavailable> File [added]: NullPsiInt_dotNewsPsi.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsiIntIO/src/ ========================================== File [added]: NullPsiIntIO_DumpPsiIntToFlexIOFile.c Delta lines: <Unavailable> File [added]: NullPsiIntIO_Observers.F90 Delta lines: <Unavailable> File [added]: NullPsiIntIO_Test.F90 Delta lines: <Unavailable> File [added]: NullPsiIntIO_WriteToFile.c Delta lines: <Unavailable> File [added]: NullPsiIntIO_sYlm.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: make.code.deps Delta lines: <Unavailable> File [added]: make.configuration.defn Delta lines: <Unavailable> File [added]: make.configuration.deps Delta lines: <Unavailable> Directory: /PITTNullCode/NullSHRExtract/src/ ============================================ File [added]: NullSHRE_AnaSchMetricRecon.F90 Delta lines: <Unavailable> File [added]: NullSHRE_ErrorSch.F90 Delta lines: <Unavailable> File [added]: NullSHRE_ExtractWT.F90 Delta lines: <Unavailable> File [added]: NullSHRE_FullMetricRecon.F90 Delta lines: <Unavailable> File [added]: NullSHRE_InitialWT.F90 Delta lines: <Unavailable> File [added]: NullSHRE_MetricReconPast.F90 Delta lines: <Unavailable> File [added]: NullSHRE_Pointing.F90 Delta lines: <Unavailable> File [added]: NullSHRE_PreCalcNSYlm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_SchMetricRecon.F90 Delta lines: <Unavailable> File [added]: NullSHRE_SchPreCalc.F90 Delta lines: <Unavailable> File [added]: NullSHRE_SchReadClm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_WTSphHarm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_mod4Metric.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modAnaCoord.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modBoundary.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modEta.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modGFdef.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modInverse.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modNullDir.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modRadius.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modSchClm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modSchYlm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modVars.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modsinSchClm.F90 Delta lines: <Unavailable> File [added]: NullSHRE_modvibSchClm.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: wt_SphHarmDecomp.F90 Delta lines: <Unavailable> Directory: /PITTNullCode/NullScalar/src/ ======================================== File [added]: NullScalar_Analysis.F90 Delta lines: <Unavailable> File [added]: NullScalar_Bdry.F90 Delta lines: <Unavailable> File [added]: NullScalar_Evol.F90 Delta lines: <Unavailable> File [added]: NullScalar_Init.F90 Delta lines: <Unavailable> File [added]: NullScalar_Loop.F90 Delta lines: <Unavailable> File [added]: NullScalar_Step.F90 Delta lines: <Unavailable> File [added]: NullScalar_Sync.F90 Delta lines: <Unavailable> File [added]: NullScalar_Testbed.F90 Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/NullVars/src/ ====================================== File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicDecomp/src/ ===================================================== File [added]: Chebyshev.cc Delta lines: <Unavailable> File [added]: Chebyshev.hh Delta lines: <Unavailable> File [added]: ExtractMetric.cc Delta lines: <Unavailable> File [added]: Legendre.cc Delta lines: <Unavailable> File [added]: Legendre.hh Delta lines: <Unavailable> File [added]: SphericalHarmonicDecomp.h Delta lines: <Unavailable> File [added]: alias.cc Delta lines: <Unavailable> File [added]: decomp.cc Delta lines: <Unavailable> File [added]: decomp.hh Delta lines: <Unavailable> File [added]: h5read.cc Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: matrix.cc Delta lines: <Unavailable> File [added]: matrix.hh Delta lines: <Unavailable> File [added]: myassert.h Delta lines: <Unavailable> File [added]: sYlm.cc Delta lines: <Unavailable> File [added]: sYlm.hh Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/src/ ==================================================== File [added]: Chebyshev.cc Delta lines: <Unavailable> File [added]: Chebyshev.hh Delta lines: <Unavailable> File [added]: Recon.cc Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicReconASCII/src/ ========================================================= File [added]: get.cc Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable> File [added]: read.cc Delta lines: <Unavailable> File [added]: sph_database.cc Delta lines: <Unavailable> File [added]: sph_database.hh Delta lines: <Unavailable> File [added]: startup.cc Delta lines: <Unavailable> File [added]: vars.hh Delta lines: <Unavailable> Directory: /PITTNullCode/NullInterp/test/ ========================================= File [added]: news_test_c.par Delta lines: <Unavailable> File [added]: null_test_c.par Delta lines: <Unavailable> File [added]: scalar_test_c.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullScalar/test/ ========================================= File [added]: t1.par Delta lines: <Unavailable> File [added]: t2.par Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsIO/src/util/ ============================================= File [added]: Decomp.f90 Delta lines: <Unavailable> File [added]: driver.f90 Delta lines: <Unavailable> File [added]: ieeeread.f90 Delta lines: <Unavailable> File [added]: interpnews.f90 Delta lines: <Unavailable> File [added]: newsread.f90 Delta lines: <Unavailable> File [added]: sYlm.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/NullNewsPostProcess/src/util/ ====================================================== File [added]: GetGridSize.f90 Delta lines: <Unavailable> File [added]: ieee_merge.c Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsi4IO/src/util/ ============================================= File [added]: Decomp.f90 Delta lines: <Unavailable> File [added]: Psi4read.f90 Delta lines: <Unavailable> File [added]: driver.f90 Delta lines: <Unavailable> File [added]: ieeeread.f90 Delta lines: <Unavailable> File [added]: interpPsi4.f90 Delta lines: <Unavailable> File [added]: sYlm.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/NullPsiIntIO/src/util/ =============================================== File [added]: Decomp.f90 Delta lines: <Unavailable> File [added]: PsiIntread.f90 Delta lines: <Unavailable> File [added]: driver.f90 Delta lines: <Unavailable> File [added]: ieeeread.f90 Delta lines: <Unavailable> File [added]: interpPsiInt.f90 Delta lines: <Unavailable> File [added]: sYlm.f90 Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicDecomp/src/util/ ========================================================== File [added]: Legendre.cc Delta lines: <Unavailable> File [added]: Legendre.hh Delta lines: <Unavailable> File [added]: Makefile Delta lines: <Unavailable> File [added]: h5read.cc Delta lines: <Unavailable> File [added]: run.sh Delta lines: <Unavailable> File [added]: sYlm.cc Delta lines: <Unavailable> File [added]: sYlm.hh Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 251)
by schnetter＠cct.lsu.edu 13 Jun '11

13 Jun '11

User: eschnett Date: 2011/06/13 10:40 AM Modified: /trunk/src/ GRHydro_CalcUpdate.F90, GRHydro_Con2Prim.F90, GRHydro_Con2PrimM.F90, GRHydro_HLLE.F90, GRHydro_Marquina.F90, GRHydro_Prim2Con.F90, GRHydro_Prim2ConM.F90, GRHydro_Reconstruct.F90, GRHydro_ReconstructM.F90, GRHydro_ReconstructPoly.F90, GRHydro_ReconstructPolyM.F90, GRHydro_Source.F90, GRHydro_SourceM.F90, GRHydro_TVDReconstruct.F90, GRHydro_Tmunu.F90, GRHydro_TmunuM.F90, GRHydro_UpdateMask.F90 Log: Declare all OpenMP private variables explicitly as private In Fortran OpenMP, loop variables of parallel loops are implicitly private and do not need to be declared as private. To improve code clarity, all such variables are now explicitly declared to be private. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_CalcUpdate.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2PrimM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLE.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Marquina.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2Con.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2ConM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Reconstruct.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructPoly.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructPolyM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Source.F90 Delta lines: <Unavailable> File [modified]: GRHydro_SourceM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_TVDReconstruct.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Tmunu.F90 Delta lines: <Unavailable> File [modified]: GRHydro_TmunuM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_UpdateMask.F90 Delta lines: <Unavailable>

1 0

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

11 Jun '11

User: cott Date: 2011/06/11 04:54 PM Modified: /trunk/ param.ccl Log: * add new parameters for hot EOS atmosphere handling File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +10 -0 =================================================================== --- trunk/param.ccl 2011-06-11 21:51:15 UTC (rev 249) +++ trunk/param.ccl 2011-06-11 21:54:46 UTC (rev 250) @@ -268,6 +268,16 @@ { } "no" +REAL GRHydro_hot_atmo_temp "Temperature of the hot atmosphere in MeV" STEERABLE=ALWAYS +{ + 0:* :: "Larger than 0 MeV" +} 0.5e0 + +REAL GRHydro_hot_atmo_Y_e "Y_e of the hot atmosphere" STEERABLE=ALWAYS +{ + 0:* :: "Larger than 0" +} 0.5e0 + ####################### Other Parameters ############################## #Parameters controlling conservative <-> primitive change.

1 0

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

11 Jun '11

User: cott Date: 2011/06/11 04:51 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_Prim2Con.F90, GRHydro_UpdateMask.F90 Log: * improve atmosphere handling so that it is now possible to use hot EOS with an atmosphere File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +26 -11 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2011-05-23 22:31:18 UTC (rev 248) +++ trunk/src/GRHydro_Con2Prim.F90 2011-06-11 21:51:15 UTC (rev 249) @@ -88,7 +88,7 @@ GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) !$OMP PARALLEL DO PRIVATE(i,j,itracer,& - !$OMP uxx, uxy, uxz, uyy, uyz, uzz, det, epsnegative) + !$OMP uxx, uxy, uxz, uyy, uyz, uzz, det, epsnegative, anyerr, keyerr, keytemp) do k = 1, nz do j = 1, ny do i = 1, nx @@ -128,7 +128,6 @@ Y_e(i,j,k) = Y_e_con(i,j,k) / dens(i,j,k) endif - if ( dens(i,j,k) .le. sqrt(det)*GRHydro_rho_min*(1.d0+GRHydro_atmo_tolerance) ) then dens(i,j,k) = sqrt(det)*GRHydro_rho_min !/(1.d0+GRHydro_atmo_tolerance) @@ -137,15 +136,26 @@ vel(i,j,k,:) = 0.d0 w_lorentz(i,j,k) = 1.d0 - call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& - rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) + if(evolve_temper.ne.0) then + ! set hot atmosphere values + temperature(i,j,k) = grhydro_hot_atmo_temp + y_e(i,j,k) = grhydro_hot_atmo_Y_e + y_e_con(i,j,k) = y_e(i,j,k) * dens(i,j,k) + keytemp = 1 + call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& + rho(i,j,k),eps(i,j,k),temperature(i,j,k),y_e(i,j,k),& + press(i,j,k),keyerr,anyerr) + else + ! use polytropic EOS + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& + rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) - call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& - rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),eps(i,j,k),keyerr,anyerr) + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& + rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),eps(i,j,k),keyerr,anyerr) + endif - ! w_lorentz=1, so the expression for tau reduces to: tau(i,j,k) = sqrt(det) * (rho(i,j,k)+rho(i,j,k)*eps(i,j,k)) - dens(i,j,k) - + cycle end if @@ -196,6 +206,7 @@ endif endif + if (epsnegative) then #if 0 @@ -848,10 +859,14 @@ rho = GRHydro_rho_min udens = rho dens = sqrt(det) * rho -! epsilon = (sqrt( (utau + pnew + udens)**2) - pnew - udens)/udens - epsilon = epsminl + temp = GRHydro_hot_atmo_temp + ye = GRHydro_hot_atmo_Y_e + keytemp=1 + call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n, & + rho,epsilon,temp,ye,xpress,keyerr,anyerr) + keytemp=0 ! w_lorentz=1, so the expression for utau reduces to: - utau = rho + rho*epsminl - udens + utau = rho + rho*epsilon - udens sx = 0.d0 sy = 0.d0 sz = 0.d0 File [modified]: GRHydro_Prim2Con.F90 Delta lines: +0 -1 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2011-05-23 22:31:18 UTC (rev 248) +++ trunk/src/GRHydro_Prim2Con.F90 2011-06-11 21:51:15 UTC (rev 249) @@ -303,7 +303,6 @@ endif endif - vlowx = gxx*dvelx + gxy*dvely + gxz*dvelz vlowy = gxy*dvelx + gyy*dvely + gyz*dvelz vlowz = gxz*dvelx + gyz*dvely + gzz*dvelz File [modified]: GRHydro_UpdateMask.F90 Delta lines: +50 -23 =================================================================== --- trunk/src/GRHydro_UpdateMask.F90 2011-05-23 22:31:18 UTC (rev 248) +++ trunk/src/GRHydro_UpdateMask.F90 2011-06-11 21:51:15 UTC (rev 249) @@ -219,13 +219,19 @@ CCTK_REAL :: det, psi4pt CCTK_INT :: type_bits, atmosphere, not_atmosphere +! begin EOS Omni vars + integer :: n,keytemp,anyerr,keyerr(1) + real*8 :: xpress,xeps,xtemp,xye + n=1;keytemp=0;anyerr=0;keyerr(1)=0 + xpress=0.0d0;xye=0.0d0;xeps=0.0d0;xtemp=0.0d0 +! end EOS Omni vars call SpaceMask_GetTypeBits(type_bits, "Hydro_Atmosphere") call SpaceMask_GetStateBits(atmosphere, "Hydro_Atmosphere",& "in_atmosphere") call SpaceMask_GetStateBits(not_atmosphere, "Hydro_Atmosphere",& "not_in_atmosphere") - +!$OMP PARALLEL DO PRIVATE(det,keytemp,i,j,k,anyerr,keyerr) do k = 1, cctk_lsh(3) do j = 1, cctk_lsh(2) do i = 1, cctk_lsh(1) @@ -235,37 +241,52 @@ atmosphere)) & &) then -!!$ write(*,*) 'Resetting at ',i,j,k, atmosphere_mask(i, j, k), & -!!$ & (SpaceMask_CheckStateBitsF90(space_mask,i, j, k, type_bits,\ -!!$ atmosphere)) - rho(i,j,k) = GRHydro_rho_min velx(i,j,k) = 0.0d0 vely(i,j,k) = 0.0d0 velz(i,j,k) = 0.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)) - call prim2conpolytype(GRHydro_polytrope_handle, & - gxx(i,j,k), gxy(i,j,k), gxz(i,j,k), & - gyy(i,j,k), gyz(i,j,k), gzz(i,j,k), det, & - dens(i,j,k), scon(i,j,k,1), scon(i,j,k,2), scon(i,j,k,3), & - 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)) - if (wk_atmosphere .eq. 0) then - atmosphere_mask(i, j, k) = 0 - SpaceMask_SetStateBitsF90(space_mask, i, j, k, type_bits,\ - not_atmosphere) - end if - end if + if(evolve_temper.ne.0) then +! ! set the temperature to be relatively low + temperature(i,j,k) = grhydro_hot_atmo_temp + y_e(i,j,k) = grhydro_hot_atmo_Y_e + keytemp = 1 + call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& + rho(i,j,k),eps(i,j,k),temperature(i,j,k),y_e(i,j,k),& + press(i,j,k),keyerr,anyerr) + call prim2con_hot(GRHydro_eos_handle, GRHydro_reflevel,& + i,j,k,x(i,j,k),y(i,j,k),z(i,j,k),gxx(i,j,k),gxy(i,j,k),& + gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k), & + det,dens(i,j,k),scon(i,j,k,1), scon(i,j,k,2), scon(i,j,k,3), & + 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),& + temperature(i,j,k),y_e(i,j,k)) + y_e_con(i,j,k) = dens(i,j,k) * y_e(i,j,k) + else + call prim2conpolytype(GRHydro_polytrope_handle, & + gxx(i,j,k), gxy(i,j,k), gxz(i,j,k), & + gyy(i,j,k), gyz(i,j,k), gzz(i,j,k), det, & + dens(i,j,k), scon(i,j,k,1), scon(i,j,k,2), scon(i,j,k,3), & + 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)) + if (wk_atmosphere .eq. 0) then + atmosphere_mask(i, j, k) = 0 + SpaceMask_SetStateBitsF90(space_mask, i, j, k, type_bits,\ + not_atmosphere) + end if + endif + + end if + end do end do end do +!$OMP END PARALLEL DO -!!$ call GRHydro_Boundaries(CCTK_PASS_FTOF) - end subroutine GRHydro_AtmosphereReset subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) @@ -316,8 +337,9 @@ velz(i,j,k) = 0.0d0 if(evolve_temper.ne.0) then - ! set the temperature to be relatively low - temperature(i,j,k) = 0.1d0 +! ! set the temperature to be relatively low + temperature(i,j,k) = grhydro_hot_atmo_temp + y_e(i,j,k) = grhydro_hot_atmo_Y_e keytemp = 1 call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& rho(i,j,k),eps(i,j,k),temperature(i,j,k),y_e(i,j,k),& @@ -330,6 +352,7 @@ 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),& temperature(i,j,k),y_e(i,j,k)) + y_e_con(i,j,k) = dens(i,j,k) * y_e(i,j,k) else call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) @@ -355,7 +378,8 @@ if(evolve_temper.ne.0) then ! set the temperature to be relatively low - temperature_p(i,j,k) = 0.1d0 + temperature_p(i,j,k) = grhydro_hot_atmo_temp + y_e_p(i,j,k) = grhydro_hot_atmo_Y_e keytemp = 1 call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& rho_p(i,j,k),eps_p(i,j,k),temperature_p(i,j,k),y_e_p(i,j,k),& @@ -368,6 +392,7 @@ tau_p(i,j,k), rho_p(i,j,k), velx_p(i,j,k), vely_p(i,j,k), & velz_p(i,j,k), eps_p(i,j,k), press_p(i,j,k), w_lorentz_p(i,j,k),& temperature_p(i,j,k),y_e_p(i,j,k)) + y_e_con_p(i,j,k) = dens_p(i,j,k) * y_e_p(i,j,k) else call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho_p(i,j,k),eps_p(i,j,k),xtemp,xye,press_p(i,j,k),keyerr,anyerr) @@ -394,7 +419,8 @@ if(evolve_temper.ne.0) then ! set the temperature to be relatively low - temperature_p_p(i,j,k) = 0.1d0 + temperature_p_p(i,j,k) = grhydro_hot_atmo_temp + y_e_p_p(i,j,k) = grhydro_hot_atmo_Y_e keytemp = 1 call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& rho_p_p(i,j,k),eps_p_p(i,j,k),temperature_p_p(i,j,k),y_e_p_p(i,j,k),& @@ -406,6 +432,7 @@ tau_p_p(i,j,k), rho_p_p(i,j,k), velx_p_p(i,j,k), vely_p_p(i,j,k), & velz_p_p(i,j,k), eps_p_p(i,j,k), press_p_p(i,j,k), w_lorentz_p_p(i,j,k),& temperature_p_p(i,j,k),y_e_p_p(i,j,k)) + y_e_con_p_p(i,j,k) = dens_p_p(i,j,k) * y_e_p_p(i,j,k) else call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho_p_p(i,j,k),eps_p_p(i,j,k),xtemp,xye,press_p_p(i,j,k),keyerr,anyerr)

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/ (Rev. 48)
by cott＠tapir.caltech.edu 10 Jun '11

10 Jun '11

User: cott Date: 2011/06/10 07:32 PM Modified: /trunk/src/ EOS_Omni_SingleVarCalls.F90 Log: * fix dpdrho issue File Changes: Directory: /trunk/src/ ====================== File [modified]: EOS_Omni_SingleVarCalls.F90 Delta lines: +3 -1 =================================================================== --- trunk/src/EOS_Omni_SingleVarCalls.F90 2011-04-21 04:46:37 UTC (rev 47) +++ trunk/src/EOS_Omni_SingleVarCalls.F90 2011-06-11 00:32:35 UTC (rev 48) @@ -203,8 +203,10 @@ press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) - dpdrhoe(i) = hybrid_dp_poly + hybrid_dp_th1 + hybrid_dp_th2 + dpdrhoe(i) = hybrid_dp_poly + max(0.0d0,hybrid_dp_th1 + hybrid_dp_th2) + + enddo case (4)

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/tov/ (Rev. 92)
by knarf＠cct.lsu.edu 10 Jun '11

10 Jun '11

User: knarf Date: 2011/06/10 05:52 PM Modified: /examples/tov/ tov.rpar Log: enable checkpointing File Changes: Directory: /examples/tov/ ========================= File [modified]: tov.rpar Delta lines: +12 -0 =================================================================== --- examples/tov/tov.rpar 2011-05-22 22:05:52 UTC (rev 91) +++ examples/tov/tov.rpar 2011-06-10 22:52:35 UTC (rev 92) @@ -235,6 +235,18 @@ ML_BSSN::ML_Ham ML_BSSN::ML_mom " + +IOHDF5::checkpoint = "yes" +IO::checkpoint_dir = \$parfile +IO::checkpoint_ID = "yes" +IO::checkpoint_every = 57600 +#IO::checkpoint_every_walltime_hours = 6.0 +IO::checkpoint_keep = 3 +IO::checkpoint_on_terminate = "yes" + +IO::recover = "autoprobe" +IO::recover_dir = \$parfile + PAR_EOF open(OUT, ">tov_${res}${type}.par");

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits