Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] www/about/members/ (Rev. 458)
by knarf＠cct.lsu.edu 30 Jan '11

30 Jan '11

User: knarf Date: 2011/01/30 09:37 PM Modified: /about/members/ members.txt Log: add Roberto De Pietri File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +3 -0 =================================================================== --- about/members/members.txt 2011-01-28 19:43:16 UTC (rev 457) +++ about/members/members.txt 2011-01-31 03:37:58 UTC (rev 458) @@ -58,6 +58,9 @@ Osaka University Luca Baiotti +Parma University + Roberto De Pietri + http://www.perimeterinstitute.ca/ Perimeter Institute David Rideout

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 22)
by tanja.bode＠physics.gatech.edu 30 Jan '11

30 Jan '11

User: tbode Date: 2011/01/30 04:54 PM Modified: / ET.tex Log: Analysis section draft. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +32 -4 =================================================================== --- ET.tex 2011-01-24 19:24:58 UTC (rev 21) +++ ET.tex 2011-01-30 22:54:11 UTC (rev 22) @@ -508,14 +508,42 @@ \subsection{Analysis} -\todo{1/2 page Tanja} +%\todo{1/2 page Tanja} -\paragraph{Masses} +Any simulated system may have many applicable quantities for analysis which +are most efficiently calculated during evolution. The Einstein Toolkit +currently includes modules to facilitate the calculation of the most common +analysis quantities: masses, horizons, and gravitational waves. In addition +to these, detailed below, the \codename{ADMConstraints} module calculates the +constraints required, but not imposed by, the formulation for assessing the +simulation's accuracy and stability. Finally, \codename{ADMAnalysis} calculates +several quantities often needed in analyzing the spacetime evolution: the trace +of the extrinsic curvature $K$, the determinant of the 3-metric, the spatial +metric and extrinsic curvature in spherical coordinates, and the spatial Ricci +tensor and trace. -\paragraph{Event and Apparent Horizons} +\subsubsection*{Horizons} For spacetimes which contain a \bh{,} the Toolkit +contains modules to find both the event and apparent horizons and their +properties. The \codename{EHFinder} module traces null surfaces backwards in +time through a pre-evolved spacetime. The Toolkit also includes two separate +apparent horizon finders: Thornburg's \codename{AHFinderDirect} and +\codename{AHFinder}. -\paragraph{Gravitational Waves} +\subsubsection*{Masses} In analyzing spacetimes, it is important to calculate +the masses and momenta of the system. Several Toolkit modules include the +capability of calculating several types of masses. The apparent horizon modules +include in their output the areal masses of the found horizons. The +\codename{Extract} module calculates ADM mass and momenta when extracting +gravitational waves. Additionally, \codename{HydroAnalysis} calculates finds +the center of mass and the point of maximum matter density of a hydronamic +field. +\subsubsection*{Gravitational Waves} Within the Toolkit are modules for extracting +gravitational waves using two formulations: the Moncrief formalism assuming +the spacetime approximates spherical symmetry \codename{Extract}) and the Weyl +scalar $\Psi_4$ constructed on a fiducial tetrad (\codename{WeylScal4}). The +Weyl scalars, as well as any other quantity, can be projected onto a basis +of spin-weighted spherical harmonics using the \codename{Multipole} module. \subsection{Relativity Tools} \todo{1/2 page Peter}

1 0

[svn:einsteintoolkit] www/toolkit/releases/ (Rev. 457)
by knarf＠cct.lsu.edu 28 Jan '11

28 Jan '11

User: knarf Date: 2011/01/28 01:43 PM Modified: /toolkit/releases/ index.php Log: fix typo File Changes: Directory: /toolkit/releases/ ============================= File [modified]: index.php Delta lines: +1 -1 =================================================================== --- toolkit/releases/index.php 2011-01-24 16:20:33 UTC (rev 456) +++ toolkit/releases/index.php 2011-01-28 19:43:16 UTC (rev 457) @@ -2,7 +2,7 @@ include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> <p>The most recent release of the Einstein Toolkit - is <tt>ET_2010_1</tt>, code named + is <tt>ET_2010_11</tt>, code named "<a href="http://en.wikipedia.org/wiki/Subrahmanyan_Chandrasekhar">Chandrasekhar</a>". We expect this version to remain unchanged, save for correcting serious problems that may be detected. This version should be your

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 35)
by schnetter＠cct.lsu.edu 26 Jan '11

26 Jan '11

User: eschnett Date: 2011/01/26 10:15 AM Added: /trunk/src/nuc_eos/ linterp.F Removed: /trunk/src/nuc_eos/ temp.f Log: Rename file, stage 2. Change suffix from .f to .F so that Cactus expands macros. File Changes: Directory: /trunk/src/nuc_eos/ ============================== File [added]: linterp.F Delta lines: +134 -0 =================================================================== --- trunk/src/nuc_eos/linterp.F (rev 0) +++ trunk/src/nuc_eos/linterp.F 2011-01-26 16:15:44 UTC (rev 35) @@ -0,0 +1,134 @@ +#include "cctk.h" + + SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, + . d1, d2, d3 ) +c + implicit none +c +c--------------------------------------------------------------------- +c +c purpose: interpolation of a function of three variables in an +c equidistant(!!!) table. +c +c method: 8-point Lagrange linear interpolation formula +c +c x input vector of first variable +c y input vector of second variable +c z input vector of third variable +c +c f output vector of interpolated function values +c +c kt vector length of input and output vectors +c +c ft 3d array of tabulated function values +c nx x-dimension of table +c ny y-dimension of table +c nz z-dimension of table +c xt vector of x-coordinates of table +c yt vector of y-coordinates of table +c zt vector of z-coordinates of table +c +c d1 centered derivative of ft with respect to x +c d2 centered derivative of ft with respect to y +c d3 centered derivative of ft with respect to z +c Note that d? only make sense when intp3d is called with kt=1 +c--------------------------------------------------------------------- +c +c + +c + integer kt,nx,ny,nz,ktx + double precision x(kt),y(kt),z(kt),f(kt) + double precision xt(nx),yt(ny),zt(nz) + double precision ft(nx,ny,nz) + double precision d1,d2,d3 +c +c + PARAMETER (ktx = 400) + double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), + & a1(ktx), a2(ktx), a3(ktx), a4(ktx), + & a5(ktx), a6(ktx), a7(ktx), a8(ktx) + + double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi + integer n,ix,iy,iz + + IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') +c +c +c------ determine spacing parameters of (equidistant!!!) table +c + dx = (xt(nx) - xt(1)) / FLOAT(nx-1) + dy = (yt(ny) - yt(1)) / FLOAT(ny-1) + dz = (zt(nz) - zt(1)) / FLOAT(nz-1) +c + dxi = 1. / dx + dyi = 1. / dy + dzi = 1. / dz +c + dxyi = dxi * dyi + dxzi = dxi * dzi + dyzi = dyi * dzi +c + dxyzi = dxi * dyi * dzi +c +c +c------- loop over all points to be interpolated +c + DO n = 1, kt +c +c------- determine location in (equidistant!!!) table +c + ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) + iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) + iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) +c + ix = MAX( 2, MIN( ix, nx ) ) + iy = MAX( 2, MIN( iy, ny ) ) + iz = MAX( 2, MIN( iz, nz ) ) +c +c write(*,*) iy-1,iy,iy+1 +c +c------- set-up auxiliary arrays for Lagrange interpolation +c + delx(n) = xt(ix) - x(n) + dely(n) = yt(iy) - y(n) + delz(n) = zt(iz) - z(n) +c + fh(n,1) = ft(ix , iy , iz ) + fh(n,2) = ft(ix-1, iy , iz ) + fh(n,3) = ft(ix , iy-1, iz ) + fh(n,4) = ft(ix , iy , iz-1) + fh(n,5) = ft(ix-1, iy-1, iz ) + fh(n,6) = ft(ix-1, iy , iz-1) + fh(n,7) = ft(ix , iy-1, iz-1) + fh(n,8) = ft(ix-1, iy-1, iz-1) +c +c------ set up coefficients of the interpolation polynomial and +c evaluate function values +c + a1(n) = fh(n,1) + a2(n) = dxi * ( fh(n,2) - fh(n,1) ) + a3(n) = dyi * ( fh(n,3) - fh(n,1) ) + a4(n) = dzi * ( fh(n,4) - fh(n,1) ) + a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) + a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) + a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) + a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + + & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) +c + d1 = -a2(n) + d2 = -a3(n) + d3 = -a4(n) + f(n) = a1(n) + a2(n) * delx(n) + & + a3(n) * dely(n) + & + a4(n) * delz(n) + & + a5(n) * delx(n) * dely(n) + & + a6(n) * delx(n) * delz(n) + & + a7(n) * dely(n) * delz(n) + & + a8(n) * delx(n) * dely(n) * delz(n) +c + ENDDO +c + RETURN + END + File [removed]: temp.f Delta lines: +0 -134 =================================================================== --- trunk/src/nuc_eos/temp.f 2011-01-26 16:15:13 UTC (rev 34) +++ trunk/src/nuc_eos/temp.f 2011-01-26 16:15:44 UTC (rev 35) @@ -1,134 +0,0 @@ -#include "cctk.h" - - SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, - . d1, d2, d3 ) -c - implicit none -c -c--------------------------------------------------------------------- -c -c purpose: interpolation of a function of three variables in an -c equidistant(!!!) table. -c -c method: 8-point Lagrange linear interpolation formula -c -c x input vector of first variable -c y input vector of second variable -c z input vector of third variable -c -c f output vector of interpolated function values -c -c kt vector length of input and output vectors -c -c ft 3d array of tabulated function values -c nx x-dimension of table -c ny y-dimension of table -c nz z-dimension of table -c xt vector of x-coordinates of table -c yt vector of y-coordinates of table -c zt vector of z-coordinates of table -c -c d1 centered derivative of ft with respect to x -c d2 centered derivative of ft with respect to y -c d3 centered derivative of ft with respect to z -c Note that d? only make sense when intp3d is called with kt=1 -c--------------------------------------------------------------------- -c -c - -c - integer kt,nx,ny,nz,ktx - double precision x(kt),y(kt),z(kt),f(kt) - double precision xt(nx),yt(ny),zt(nz) - double precision ft(nx,ny,nz) - double precision d1,d2,d3 -c -c - PARAMETER (ktx = 400) - double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), - & a1(ktx), a2(ktx), a3(ktx), a4(ktx), - & a5(ktx), a6(ktx), a7(ktx), a8(ktx) - - double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi - integer n,ix,iy,iz - - IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') -c -c -c------ determine spacing parameters of (equidistant!!!) table -c - dx = (xt(nx) - xt(1)) / FLOAT(nx-1) - dy = (yt(ny) - yt(1)) / FLOAT(ny-1) - dz = (zt(nz) - zt(1)) / FLOAT(nz-1) -c - dxi = 1. / dx - dyi = 1. / dy - dzi = 1. / dz -c - dxyi = dxi * dyi - dxzi = dxi * dzi - dyzi = dyi * dzi -c - dxyzi = dxi * dyi * dzi -c -c -c------- loop over all points to be interpolated -c - DO n = 1, kt -c -c------- determine location in (equidistant!!!) table -c - ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) - iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) - iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) -c - ix = MAX( 2, MIN( ix, nx ) ) - iy = MAX( 2, MIN( iy, ny ) ) - iz = MAX( 2, MIN( iz, nz ) ) -c -c write(*,*) iy-1,iy,iy+1 -c -c------- set-up auxiliary arrays for Lagrange interpolation -c - delx(n) = xt(ix) - x(n) - dely(n) = yt(iy) - y(n) - delz(n) = zt(iz) - z(n) -c - fh(n,1) = ft(ix , iy , iz ) - fh(n,2) = ft(ix-1, iy , iz ) - fh(n,3) = ft(ix , iy-1, iz ) - fh(n,4) = ft(ix , iy , iz-1) - fh(n,5) = ft(ix-1, iy-1, iz ) - fh(n,6) = ft(ix-1, iy , iz-1) - fh(n,7) = ft(ix , iy-1, iz-1) - fh(n,8) = ft(ix-1, iy-1, iz-1) -c -c------ set up coefficients of the interpolation polynomial and -c evaluate function values -c - a1(n) = fh(n,1) - a2(n) = dxi * ( fh(n,2) - fh(n,1) ) - a3(n) = dyi * ( fh(n,3) - fh(n,1) ) - a4(n) = dzi * ( fh(n,4) - fh(n,1) ) - a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) - a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) - a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) - a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + - & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) -c - d1 = -a2(n) - d2 = -a3(n) - d3 = -a4(n) - f(n) = a1(n) + a2(n) * delx(n) - & + a3(n) * dely(n) - & + a4(n) * delz(n) - & + a5(n) * delx(n) * dely(n) - & + a6(n) * delx(n) * delz(n) - & + a7(n) * dely(n) * delz(n) - & + a8(n) * delx(n) * dely(n) * delz(n) -c - ENDDO -c - RETURN - END -

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 34)
by schnetter＠cct.lsu.edu 26 Jan '11

26 Jan '11

User: eschnett Date: 2011/01/26 10:15 AM Added: /trunk/src/nuc_eos/ temp.f Removed: /trunk/src/nuc_eos/ linterp.f Log: Rename file, stage 1 File Changes: Directory: /trunk/src/nuc_eos/ ============================== File [removed]: linterp.f Delta lines: +0 -134 =================================================================== --- trunk/src/nuc_eos/linterp.f 2011-01-26 16:14:37 UTC (rev 33) +++ trunk/src/nuc_eos/linterp.f 2011-01-26 16:15:13 UTC (rev 34) @@ -1,134 +0,0 @@ -#include "cctk.h" - - SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, - . d1, d2, d3 ) -c - implicit none -c -c--------------------------------------------------------------------- -c -c purpose: interpolation of a function of three variables in an -c equidistant(!!!) table. -c -c method: 8-point Lagrange linear interpolation formula -c -c x input vector of first variable -c y input vector of second variable -c z input vector of third variable -c -c f output vector of interpolated function values -c -c kt vector length of input and output vectors -c -c ft 3d array of tabulated function values -c nx x-dimension of table -c ny y-dimension of table -c nz z-dimension of table -c xt vector of x-coordinates of table -c yt vector of y-coordinates of table -c zt vector of z-coordinates of table -c -c d1 centered derivative of ft with respect to x -c d2 centered derivative of ft with respect to y -c d3 centered derivative of ft with respect to z -c Note that d? only make sense when intp3d is called with kt=1 -c--------------------------------------------------------------------- -c -c - -c - integer kt,nx,ny,nz,ktx - double precision x(kt),y(kt),z(kt),f(kt) - double precision xt(nx),yt(ny),zt(nz) - double precision ft(nx,ny,nz) - double precision d1,d2,d3 -c -c - PARAMETER (ktx = 400) - double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), - & a1(ktx), a2(ktx), a3(ktx), a4(ktx), - & a5(ktx), a6(ktx), a7(ktx), a8(ktx) - - double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi - integer n,ix,iy,iz - - IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') -c -c -c------ determine spacing parameters of (equidistant!!!) table -c - dx = (xt(nx) - xt(1)) / FLOAT(nx-1) - dy = (yt(ny) - yt(1)) / FLOAT(ny-1) - dz = (zt(nz) - zt(1)) / FLOAT(nz-1) -c - dxi = 1. / dx - dyi = 1. / dy - dzi = 1. / dz -c - dxyi = dxi * dyi - dxzi = dxi * dzi - dyzi = dyi * dzi -c - dxyzi = dxi * dyi * dzi -c -c -c------- loop over all points to be interpolated -c - DO n = 1, kt -c -c------- determine location in (equidistant!!!) table -c - ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) - iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) - iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) -c - ix = MAX( 2, MIN( ix, nx ) ) - iy = MAX( 2, MIN( iy, ny ) ) - iz = MAX( 2, MIN( iz, nz ) ) -c -c write(*,*) iy-1,iy,iy+1 -c -c------- set-up auxiliary arrays for Lagrange interpolation -c - delx(n) = xt(ix) - x(n) - dely(n) = yt(iy) - y(n) - delz(n) = zt(iz) - z(n) -c - fh(n,1) = ft(ix , iy , iz ) - fh(n,2) = ft(ix-1, iy , iz ) - fh(n,3) = ft(ix , iy-1, iz ) - fh(n,4) = ft(ix , iy , iz-1) - fh(n,5) = ft(ix-1, iy-1, iz ) - fh(n,6) = ft(ix-1, iy , iz-1) - fh(n,7) = ft(ix , iy-1, iz-1) - fh(n,8) = ft(ix-1, iy-1, iz-1) -c -c------ set up coefficients of the interpolation polynomial and -c evaluate function values -c - a1(n) = fh(n,1) - a2(n) = dxi * ( fh(n,2) - fh(n,1) ) - a3(n) = dyi * ( fh(n,3) - fh(n,1) ) - a4(n) = dzi * ( fh(n,4) - fh(n,1) ) - a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) - a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) - a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) - a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + - & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) -c - d1 = -a2(n) - d2 = -a3(n) - d3 = -a4(n) - f(n) = a1(n) + a2(n) * delx(n) - & + a3(n) * dely(n) - & + a4(n) * delz(n) - & + a5(n) * delx(n) * dely(n) - & + a6(n) * delx(n) * delz(n) - & + a7(n) * dely(n) * delz(n) - & + a8(n) * delx(n) * dely(n) * delz(n) -c - ENDDO -c - RETURN - END - File [added]: temp.f Delta lines: +134 -0 =================================================================== --- trunk/src/nuc_eos/temp.f (rev 0) +++ trunk/src/nuc_eos/temp.f 2011-01-26 16:15:13 UTC (rev 34) @@ -0,0 +1,134 @@ +#include "cctk.h" + + SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, + . d1, d2, d3 ) +c + implicit none +c +c--------------------------------------------------------------------- +c +c purpose: interpolation of a function of three variables in an +c equidistant(!!!) table. +c +c method: 8-point Lagrange linear interpolation formula +c +c x input vector of first variable +c y input vector of second variable +c z input vector of third variable +c +c f output vector of interpolated function values +c +c kt vector length of input and output vectors +c +c ft 3d array of tabulated function values +c nx x-dimension of table +c ny y-dimension of table +c nz z-dimension of table +c xt vector of x-coordinates of table +c yt vector of y-coordinates of table +c zt vector of z-coordinates of table +c +c d1 centered derivative of ft with respect to x +c d2 centered derivative of ft with respect to y +c d3 centered derivative of ft with respect to z +c Note that d? only make sense when intp3d is called with kt=1 +c--------------------------------------------------------------------- +c +c + +c + integer kt,nx,ny,nz,ktx + double precision x(kt),y(kt),z(kt),f(kt) + double precision xt(nx),yt(ny),zt(nz) + double precision ft(nx,ny,nz) + double precision d1,d2,d3 +c +c + PARAMETER (ktx = 400) + double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), + & a1(ktx), a2(ktx), a3(ktx), a4(ktx), + & a5(ktx), a6(ktx), a7(ktx), a8(ktx) + + double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi + integer n,ix,iy,iz + + IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') +c +c +c------ determine spacing parameters of (equidistant!!!) table +c + dx = (xt(nx) - xt(1)) / FLOAT(nx-1) + dy = (yt(ny) - yt(1)) / FLOAT(ny-1) + dz = (zt(nz) - zt(1)) / FLOAT(nz-1) +c + dxi = 1. / dx + dyi = 1. / dy + dzi = 1. / dz +c + dxyi = dxi * dyi + dxzi = dxi * dzi + dyzi = dyi * dzi +c + dxyzi = dxi * dyi * dzi +c +c +c------- loop over all points to be interpolated +c + DO n = 1, kt +c +c------- determine location in (equidistant!!!) table +c + ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) + iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) + iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) +c + ix = MAX( 2, MIN( ix, nx ) ) + iy = MAX( 2, MIN( iy, ny ) ) + iz = MAX( 2, MIN( iz, nz ) ) +c +c write(*,*) iy-1,iy,iy+1 +c +c------- set-up auxiliary arrays for Lagrange interpolation +c + delx(n) = xt(ix) - x(n) + dely(n) = yt(iy) - y(n) + delz(n) = zt(iz) - z(n) +c + fh(n,1) = ft(ix , iy , iz ) + fh(n,2) = ft(ix-1, iy , iz ) + fh(n,3) = ft(ix , iy-1, iz ) + fh(n,4) = ft(ix , iy , iz-1) + fh(n,5) = ft(ix-1, iy-1, iz ) + fh(n,6) = ft(ix-1, iy , iz-1) + fh(n,7) = ft(ix , iy-1, iz-1) + fh(n,8) = ft(ix-1, iy-1, iz-1) +c +c------ set up coefficients of the interpolation polynomial and +c evaluate function values +c + a1(n) = fh(n,1) + a2(n) = dxi * ( fh(n,2) - fh(n,1) ) + a3(n) = dyi * ( fh(n,3) - fh(n,1) ) + a4(n) = dzi * ( fh(n,4) - fh(n,1) ) + a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) + a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) + a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) + a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + + & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) +c + d1 = -a2(n) + d2 = -a3(n) + d3 = -a4(n) + f(n) = a1(n) + a2(n) * delx(n) + & + a3(n) * dely(n) + & + a4(n) * delz(n) + & + a5(n) * delx(n) * dely(n) + & + a6(n) * delx(n) * delz(n) + & + a7(n) * dely(n) * delz(n) + & + a8(n) * delx(n) * dely(n) * delz(n) +c + ENDDO +c + RETURN + END +

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 33)
by schnetter＠cct.lsu.edu 26 Jan '11

26 Jan '11

User: eschnett Date: 2011/01/26 10:14 AM Modified: /trunk/src/nuc_eos/ linterp.f, linterp_many.F90, make.code.defn, nuc_eos.F90, readtable.F90 Log: Replace stop commands with calls to CCTK_WARN, since stop with a string argument is not supported by all Fortran compilers. File Changes: Directory: /trunk/src/nuc_eos/ ============================== File [modified]: linterp.f Delta lines: +3 -1 =================================================================== --- trunk/src/nuc_eos/linterp.f 2010-12-22 21:44:39 UTC (rev 32) +++ trunk/src/nuc_eos/linterp.f 2011-01-26 16:14:37 UTC (rev 33) @@ -1,3 +1,5 @@ +#include "cctk.h" + SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, . d1, d2, d3 ) c @@ -50,7 +52,7 @@ double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi integer n,ix,iy,iz - IF (kt .GT. ktx) STOP'***KTX**' + IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') c c c------ determine spacing parameters of (equidistant!!!) table File [modified]: linterp_many.F90 Delta lines: +4 -2 =================================================================== --- trunk/src/nuc_eos/linterp_many.F90 2010-12-22 21:44:39 UTC (rev 32) +++ trunk/src/nuc_eos/linterp_many.F90 2011-01-26 16:14:37 UTC (rev 33) @@ -1,3 +1,5 @@ +#include "cctk.h" + SUBROUTINE intp3d_many ( x, y, z, f, kt, ft, nx, ny, nz, nvars, xt, yt, zt) ! implicit none @@ -43,7 +45,7 @@ real*8 dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi integer n,ix,iy,iz - IF (kt .GT. ktx) STOP'***KTX**' + IF (kt .GT. ktx) call CCTK_WARN(0, '***KTX**') ! ! !------ determine spacing parameters of (equidistant!!!) table @@ -65,7 +67,7 @@ ! !------- loop over all points to be interpolated ! - dO n = 1, kt + do n = 1, kt ! !------- determine location in (equidistant!!!) table ! File [modified]: make.code.defn Delta lines: +1 -1 =================================================================== --- trunk/src/nuc_eos/make.code.defn 2010-12-22 21:44:39 UTC (rev 32) +++ trunk/src/nuc_eos/make.code.defn 2011-01-26 16:14:37 UTC (rev 33) @@ -1,6 +1,6 @@ SRCS = eosmodule.F90 nuc_eos.F90 bisection.F90 \ findtemp.F90 findrho.F90 linterp_many.F90 readtable.F90 \ - linterp.f + linterp.F SUBDIRS = File [modified]: nuc_eos.F90 Delta lines: +16 -13 =================================================================== --- trunk/src/nuc_eos/nuc_eos.F90 2010-12-22 21:44:39 UTC (rev 32) +++ trunk/src/nuc_eos/nuc_eos.F90 2011-01-26 16:14:37 UTC (rev 33) @@ -17,6 +17,9 @@ ! 1 -> coming in with temperature ! 2 -> coming in with entropy ! + +#include "cctk.h" + subroutine nuc_eos_full(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& xdpderho,xdpdrhoe,xxa,xxh,xxn,xxp,xabar,xzbar,xmu_e,xmu_n,xmu_p, & xmuhat,keytemp,keyerr,rfeps) @@ -45,11 +48,11 @@ keyerrr = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN(0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin) then - stop "nuc_eos: rho < rhomin" + call CCTK_WARN (0, "nuc_eos: rho < rhomin") endif if(xye.gt.eos_yemax) then @@ -64,11 +67,11 @@ if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then - stop "nuc_eos: temp < tempmin" + call CCTK_WARN (0, "nuc_eos: temp < tempmin") endif endif @@ -84,7 +87,7 @@ !need to find temperature based on xeps call findtemp(lr,lt,y,leps,keyerrt,rfeps) if(keyerrt.ne.0) then - stop "Did not find temperature" + call CCTK_WARN (0, "Did not find temperature") endif xtemp = 10.0d0**lt @@ -206,7 +209,7 @@ keyerrr = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN (0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin*1.2d0) then @@ -227,7 +230,7 @@ if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then @@ -329,7 +332,7 @@ keyerrt = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN (0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin*1.2d0) then @@ -349,7 +352,7 @@ if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then @@ -361,7 +364,7 @@ keyerr = 0 if(keytemp.gt.1) then - stop "eos_nuc_press does not support keytemp other than 0 and 1" + call CCTK_WARN (0, "eos_nuc_press does not support keytemp other than 0 and 1") endif lr = log10(xrho) @@ -408,7 +411,7 @@ keyerrt = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN (0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin*1.2d0) then @@ -428,7 +431,7 @@ if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then @@ -455,7 +458,7 @@ xtemp = 10.0d0**lt elseif(keytemp.gt.1) then - stop "eos_nuc_press does not support keytemp > 1" + call CCTK_WARN (0, "eos_nuc_press does not support keytemp > 1") endif ! have temperature, proceed: File [modified]: readtable.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/nuc_eos/readtable.F90 2010-12-22 21:44:39 UTC (rev 32) +++ trunk/src/nuc_eos/readtable.F90 2011-01-26 16:14:37 UTC (rev 33) @@ -42,7 +42,7 @@ call h5dclose_f(dset_id,error) if(error.ne.0) then - stop "Could not read EOS table file" + call CCTK_WARN (0, "Could not read EOS table file") endif dims1(1)=1 @@ -51,7 +51,7 @@ call h5dclose_f(dset_id,error) if(error.ne.0) then - stop "Could not read EOS table file" + call CCTK_WARN (0, "Could not read EOS table file") endif dims1(1)=1 @@ -60,7 +60,7 @@ call h5dclose_f(dset_id,error) if(error.ne.0) then - stop "Could not read EOS table file" + call CCTK_WARN (0, "Could not read EOS table file") endif ! write(message,"(a25,1P3i5)") "We have nrho ntemp nye: ", nrho,ntemp,nye @@ -213,7 +213,7 @@ accerr=accerr+error if(accerr.ne.0) then - stop "Problem reading EOS table file" + call CCTK_WARN (0, "Problem reading EOS table file") endif call h5fclose_f (file_id,error)

1 0

[svn:einsteintoolkit] AHFinderDirect/trunk/ (Rev. 1554)
by knarf＠cct.lsu.edu 25 Jan '11

25 Jan '11

User: knarf Date: 2011/01/25 03:48 PM Modified: /trunk/ param.ccl, schedule.ccl Log: Change the number of timelevels of the ahmask to match the number which is used for the metric. Setting this to three unconditionally leads to problems when running with Carpet and prolongation_order_time=1 because all GFs in Carpet are expected to have the same number of timelevels: prolongation_order_time+1. This is not the ideal solution. That would eliminate also variables ala metric_timelevels, because they directly depend on prolongation_order_time and should not have to be set (correctly) in a parameter file. This could be done automatically. However, this patch for now uses metric_timelevels, in order to get AHFinderDirect working quickly. File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +1 -0 =================================================================== --- trunk/param.ccl 2010-11-21 03:28:57 UTC (rev 1553) +++ trunk/param.ccl 2011-01-25 21:48:16 UTC (rev 1554) @@ -20,6 +20,7 @@ # note our access is read-only shares: ADMBase USES KEYWORD metric_type +USES CCTK_INT metric_timelevels # we need to look at SpaceMask::use_mask in order to find out about the # excision region File [modified]: schedule.ccl Delta lines: +13 -1 =================================================================== --- trunk/schedule.ccl 2010-11-21 03:28:57 UTC (rev 1553) +++ trunk/schedule.ccl 2011-01-25 21:48:16 UTC (rev 1554) @@ -2,8 +2,20 @@ # $Header$ storage: ah_radius ah_origin ah_centroid ah_flags -storage: ahmask[3] +if (metric_timelevels == 1) +{ + storage: ahmask[1] +} +else if (metric_timelevels == 2) +{ + storage: ahmask[2] +} +else if (metric_timelevels == 3) +{ + storage: ahmask[3] +} + # # setup #

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 21)
by bcmsma＠astro.rit.edu 24 Jan '11

24 Jan '11

User: bmundim Date: 2011/01/24 01:24 PM Modified: / ET.tex Log: Add RIT affiliation address. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +4 -1 =================================================================== --- ET.tex 2011-01-24 19:15:09 UTC (rev 20) +++ ET.tex 2011-01-24 19:24:58 UTC (rev 21) @@ -78,6 +78,10 @@ \thanks{{}$^{(3)}$ Center for Computation \& Technology, Louisiana State University, Baton Rouge, LA 70803, USA} +\thanks{{}$^{(4)}$ Center for Computational Relativity \& Gravitation, + School of Mathematical Sciences, + Rochester Institute of Technology, Rochester, New York 14623, USA} + % TODO: Already add a few affiliations which we know are going % to be used @@ -413,7 +417,6 @@ \subsection{Initial Data} \todo{1/2 page Josh, Bruno} - The Einstein Toolkit contains many modules used to generate initial data for GR simulations, including both vacuum and hydrodynamical configurations. These include modules used primarily for testing out various components of

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 71)
by bcmsma＠astro.rit.edu 24 Jan '11

24 Jan '11

User: bmundim Date: 2011/01/24 01:15 PM Modified: /trunk/ einsteintoolkit.bib Log: Fix a few reference keys. File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: +11 -0 =================================================================== --- trunk/einsteintoolkit.bib 2011-01-24 19:01:49 UTC (rev 70) +++ trunk/einsteintoolkit.bib 2011-01-24 19:15:35 UTC (rev 71) @@ -230,6 +230,17 @@ SLACcitation = "%%CITATION = 0711.1097;%%" } +@Article{Brill:1959zz, + author = "Brill, Dieter R.", + title = "{On the positive definite mass of the Bondi-Weber-Wheeler + time-symmetric gravitational waves}", + journal = "Annals Phys.", + volume = "7", + year = "1959", + pages = "466-483", + SLACcitation = "%%CITATION = APNYA,7,466;%%" +} + @Article{Brown:2007pg, requested-for = "EinsteinInitialData/NoExcision", author = "David Brown and Olivier Sarbach and Erik Schnetter

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 20)
by bcmsma＠astro.rit.edu 24 Jan '11

24 Jan '11

User: bmundim Date: 2011/01/24 01:15 PM Modified: / ET.tex Log: Fix a few reference keys. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +3 -3 =================================================================== --- ET.tex 2011-01-24 18:01:31 UTC (rev 19) +++ ET.tex 2011-01-24 19:15:09 UTC (rev 20) @@ -425,9 +425,9 @@ numerical software packages that are included in the toolkit as External libraries. The {\tt twopunctures} module \cite{Ansorg:2004ds}, commonly used in numerical relativity to generate binary black hole data, invokes -the GNU Scientific Library [GSL; \cite{Galassi:2009}]. Several modules +the GNU Scientific Library [GSL; \cite{GSL:web,Galassi:2009}]. Several modules have also been implemented to read in datafiles generated by the {\tt Lorene -code} \cite{Lorene:web,GGTMB}, including the BHBH, BHNS, and NSNS data +code} \cite{Lorene:web,Gourgoulhon:2000nn}, including the BHBH, BHNS, and NSNS data made publicly available through the Lorene website. Scheduling of initial data routines generally follows a standard format. @@ -469,7 +469,7 @@ \end{enumerate} \item Vacuum gravitational wave configurations: \begin{enumerate} -\item {\tt IDBrillData}: A Brill wave spacetime \cite{Brill:1959}. +\item {\tt IDBrillData}: A Brill wave spacetime \cite{Brill:1959zz}. \item {\tt IDLinearWaves}: A spacetime containing a linear gravitational wave. \end{enumerate} \item Black Hole configurations:

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits