Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] WeylScal4/trunk/src/ (Rev. 81)
by roland.haas＠physics.gatech.edu 27 Sep '11

27 Sep '11

User: rhaas Date: 2011/09/27 07:26 PM Modified: /trunk/src/ Differencing.h, psis_calc_2nd.cc, psis_calc_4th.cc, psis_calc_Nth.cc Log: Regenerate code. File Changes: Directory: /trunk/src/ ====================== File [modified]: Differencing.h Delta lines: +0 -216 =================================================================== --- trunk/src/Differencing.h 2011-09-28 00:25:42 UTC (rev 80) +++ trunk/src/Differencing.h 2011-09-28 00:26:04 UTC (rev 81) @@ -862,219 +862,3 @@ } #endif -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus1(u) ((-(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]))*p1odx) -#else -# define PDplus1(u) (PDplus1_impl(u,p1odx,cdj,cdk)) -static CCTK_REAL PDplus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]))*p1odx; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus2(u) ((-(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]))*p1ody) -#else -# define PDplus2(u) (PDplus2_impl(u,p1ody,cdj,cdk)) -static CCTK_REAL PDplus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]))*p1ody; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus3(u) ((-(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]))*p1odz) -#else -# define PDplus3(u) (PDplus3_impl(u,p1odz,cdj,cdk)) -static CCTK_REAL PDplus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]))*p1odz; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDminus1(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(-1)+cdj*(0)+cdk*(0)]))*p1odx) -#else -# define PDminus1(u) (PDminus1_impl(u,p1odx,cdj,cdk)) -static CCTK_REAL PDminus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDminus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(-1)+cdj*(0)+cdk*(0)]))*p1odx; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDminus2(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(-1)+cdk*(0)]))*p1ody) -#else -# define PDminus2(u) (PDminus2_impl(u,p1ody,cdj,cdk)) -static CCTK_REAL PDminus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDminus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(-1)+cdk*(0)]))*p1ody; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDminus3(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(-1)]))*p1odz) -#else -# define PDminus3(u) (PDminus3_impl(u,p1odz,cdj,cdk)) -static CCTK_REAL PDminus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDminus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(-1)]))*p1odz; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus11(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - 2*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(2)+cdj*(0)+cdk*(0)]))*p1odx2) -#else -# define PDplus11(u) (PDplus11_impl(u,p1odx2,cdj,cdk)) -static CCTK_REAL PDplus11_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus11_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx2, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - 2*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(2)+cdj*(0)+cdk*(0)]))*p1odx2; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus12(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(1)+cdk*(0)]))*p1odxdy) -#else -# define PDplus12(u) (PDplus12_impl(u,p1odxdy,cdj,cdk)) -static CCTK_REAL PDplus12_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus12_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(1)+cdk*(0)]))*p1odxdy; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus13(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(1)]))*p1odxdz) -#else -# define PDplus13(u) (PDplus13_impl(u,p1odxdz,cdj,cdk)) -static CCTK_REAL PDplus13_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus13_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(1)]))*p1odxdz; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus21(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(1)+cdk*(0)]))*p1odxdy) -#else -# define PDplus21(u) (PDplus21_impl(u,p1odxdy,cdj,cdk)) -static CCTK_REAL PDplus21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(1)+cdk*(0)]))*p1odxdy; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus22(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - 2*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(2)+cdk*(0)]))*p1ody2) -#else -# define PDplus22(u) (PDplus22_impl(u,p1ody2,cdj,cdk)) -static CCTK_REAL PDplus22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody2, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - 2*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(2)+cdk*(0)]))*p1ody2; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus23(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(1)]))*p1odydz) -#else -# define PDplus23(u) (PDplus23_impl(u,p1odydz,cdj,cdk)) -static CCTK_REAL PDplus23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odydz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(1)]))*p1odydz; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus31(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(1)]))*p1odxdz) -#else -# define PDplus31(u) (PDplus31_impl(u,p1odxdz,cdj,cdk)) -static CCTK_REAL PDplus31_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus31_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(1)+cdj*(0)+cdk*(1)]))*p1odxdz; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus32(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(1)]))*p1odydz) -#else -# define PDplus32(u) (PDplus32_impl(u,p1odydz,cdj,cdk)) -static CCTK_REAL PDplus32_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus32_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odydz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) - (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(1)+cdk*(1)]))*p1odydz; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDplus33(u) (((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - 2*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(2)]))*p1odz2) -#else -# define PDplus33(u) (PDplus33_impl(u,p1odz2,cdj,cdk)) -static CCTK_REAL PDplus33_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDplus33_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz2, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return ((*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) - 2*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(1)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(2)]))*p1odz2; -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDonesided2nd1(u) ((3*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(2*dir(1))+cdj*(0)+cdk*(0)]) - 4*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(dir(1))+cdj*(0)+cdk*(0)]))*pm1o2dx*dir(1)) -#else -# define PDonesided2nd1(u) (PDonesided2nd1_impl(u,pm1o2dx,cdj,cdk)) -static CCTK_REAL PDonesided2nd1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDonesided2nd1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (3*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(2*dir(1))+cdj*(0)+cdk*(0)]) - 4*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(dir(1))+cdj*(0)+cdk*(0)]))*pm1o2dx*dir(1); -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDonesided2nd2(u) ((3*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(2*dir(2))+cdk*(0)]) - 4*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(dir(2))+cdk*(0)]))*pm1o2dy*dir(2)) -#else -# define PDonesided2nd2(u) (PDonesided2nd2_impl(u,pm1o2dy,cdj,cdk)) -static CCTK_REAL PDonesided2nd2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDonesided2nd2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (3*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(2*dir(2))+cdk*(0)]) - 4*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(dir(2))+cdk*(0)]))*pm1o2dy*dir(2); -} -#endif - -#ifndef KRANC_DIFF_FUNCTIONS -# define PDonesided2nd3(u) ((3*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(2*dir(3))]) - 4*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(dir(3))]))*pm1o2dz*dir(3)) -#else -# define PDonesided2nd3(u) (PDonesided2nd3_impl(u,pm1o2dz,cdj,cdk)) -static CCTK_REAL PDonesided2nd3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; -static CCTK_REAL PDonesided2nd3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk) -{ - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (3*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(0)]) + (*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(2*dir(3))]) - 4*(*(CCTK_REAL const*)&((char const*)(u))[cdi*(0)+cdj*(0)+cdk*(dir(3))]))*pm1o2dz*dir(3); -} -#endif - File [modified]: psis_calc_2nd.cc Delta lines: +0 -9 =================================================================== --- trunk/src/psis_calc_2nd.cc 2011-09-28 00:25:42 UTC (rev 80) +++ trunk/src/psis_calc_2nd.cc 2011-09-28 00:26:04 UTC (rev 81) @@ -132,21 +132,12 @@ CCTK_REAL const p1o840dx = 0.00119047619047619047619047619048*INV(dx); CCTK_REAL const p1o840dy = 0.00119047619047619047619047619048*INV(dy); CCTK_REAL const p1o840dz = 0.00119047619047619047619047619048*INV(dz); - CCTK_REAL const p1odx = INV(dx); CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdy = INV(dx)*INV(dy); - CCTK_REAL const p1odxdz = INV(dx)*INV(dz); - CCTK_REAL const p1ody = INV(dy); CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odydz = INV(dy)*INV(dz); - CCTK_REAL const p1odz = INV(dz); CCTK_REAL const p1odz2 = INV(SQR(dz)); CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); /* Loop over the grid points */ #pragma omp parallel File [modified]: psis_calc_4th.cc Delta lines: +0 -9 =================================================================== --- trunk/src/psis_calc_4th.cc 2011-09-28 00:25:42 UTC (rev 80) +++ trunk/src/psis_calc_4th.cc 2011-09-28 00:26:04 UTC (rev 81) @@ -132,21 +132,12 @@ CCTK_REAL const p1o840dx = 0.00119047619047619047619047619048*INV(dx); CCTK_REAL const p1o840dy = 0.00119047619047619047619047619048*INV(dy); CCTK_REAL const p1o840dz = 0.00119047619047619047619047619048*INV(dz); - CCTK_REAL const p1odx = INV(dx); CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdy = INV(dx)*INV(dy); - CCTK_REAL const p1odxdz = INV(dx)*INV(dz); - CCTK_REAL const p1ody = INV(dy); CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odydz = INV(dy)*INV(dz); - CCTK_REAL const p1odz = INV(dz); CCTK_REAL const p1odz2 = INV(SQR(dz)); CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); /* Loop over the grid points */ #pragma omp parallel File [modified]: psis_calc_Nth.cc Delta lines: +0 -9 =================================================================== --- trunk/src/psis_calc_Nth.cc 2011-09-28 00:25:42 UTC (rev 80) +++ trunk/src/psis_calc_Nth.cc 2011-09-28 00:26:04 UTC (rev 81) @@ -132,21 +132,12 @@ CCTK_REAL const p1o840dx = 0.00119047619047619047619047619048*INV(dx); CCTK_REAL const p1o840dy = 0.00119047619047619047619047619048*INV(dy); CCTK_REAL const p1o840dz = 0.00119047619047619047619047619048*INV(dz); - CCTK_REAL const p1odx = INV(dx); CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdy = INV(dx)*INV(dy); - CCTK_REAL const p1odxdz = INV(dx)*INV(dz); - CCTK_REAL const p1ody = INV(dy); CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odydz = INV(dy)*INV(dz); - CCTK_REAL const p1odz = INV(dz); CCTK_REAL const p1odz2 = INV(SQR(dz)); CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); /* Loop over the grid points */ #pragma omp parallel

1 0

[svn:einsteintoolkit] WeylScal4/trunk/m/ (Rev. 80)
by roland.haas＠physics.gatech.edu 27 Sep '11

27 Sep '11

User: rhaas Date: 2011/09/27 07:25 PM Modified: /trunk/m/ WeylScal4.m Log: remove unused differential operator which causes compilation to fail File Changes: Directory: /trunk/m/ ==================== File [modified]: WeylScal4.m Delta lines: +1 -12 =================================================================== --- trunk/m/WeylScal4.m 2011-09-25 19:20:15 UTC (rev 79) +++ trunk/m/WeylScal4.m 2011-09-28 00:25:42 UTC (rev 80) @@ -28,18 +28,7 @@ PDstandard[i_] -> StandardCenteredDifferenceOperator[1,fdOrder/2,i], PDstandard[i_, i_] -> StandardCenteredDifferenceOperator[2,fdOrder/2,i], PDstandard[i_, j_] -> StandardCenteredDifferenceOperator[1,fdOrder/2,i] * - StandardCenteredDifferenceOperator[1,fdOrder/2,j], - - PDplus[i_] -> DPlus[i], - PDminus[i_] -> DMinus[i], - PDplus[i_,j_] -> DPlus[i] DPlus[j], - - PDonesided2nd[1] -> dir[1] (-shift[1]^(2 dir[1]) + 4 shift[1]^dir[1] - 3 )/ - (2 spacing[1]), - PDonesided2nd[2] -> dir[2] (-shift[2]^(2 dir[2]) + 4 shift[2]^dir[2] - 3 )/ - (2 spacing[2]), - PDonesided2nd[3] -> dir[3] (-shift[3]^(2 dir[3]) + 4 shift[3]^dir[3] - 3 )/ - (2 spacing[3]) + StandardCenteredDifferenceOperator[1,fdOrder/2,j] }; (*PD = PDstandard2nd;*)

1 0

[svn:einsteintoolkit] SphericalHarmonicRecon/trunk/test/ (Rev. 8)
by yosef＠astro.rit.edu 27 Sep '11

27 Sep '11

User: zlochower Date: 2011/09/27 03:31 PM Modified: /trunk/test/ regression_test.par Log: incorrect data file location fixed File Changes: Directory: /trunk/test/ ======================= File [modified]: regression_test.par Delta lines: +2 -2 =================================================================== --- trunk/test/regression_test.par 2011-09-06 18:39:29 UTC (rev 7) +++ trunk/test/regression_test.par 2011-09-27 20:31:09 UTC (rev 8) @@ -33,8 +33,8 @@ #Override extraction parameters--------------- SphericalHarmonicRecon::time_derivative_in_file = "yes" -#SphericalHarmonicRecon::metric_data_filename = "../../../arrangements/PITTNullCode/SphericalHarmonicRecon/test/metric_test_case.h5" -SphericalHarmonicRecon::metric_data_filename = "metric_test_case.h5" +SphericalHarmonicRecon::metric_data_filename = "../../../arrangements/PITTNullCode/SphericalHarmonicRecon/test/metric_test_case.h5" +#SphericalHarmonicRecon::metric_data_filename = "metric_test_case.h5" #Extraction data------------------------

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 619)
by roland.haas＠physics.gatech.edu 27 Sep '11

27 Sep '11

User: rhaas Date: 2011/09/27 12:38 PM Modified: /release-info/ run_tests_py Log: specify number processors per node in run_tests_py File Changes: Directory: /release-info/ ========================= File [modified]: run_tests_py Delta lines: +2 -2 =================================================================== --- release-info/run_tests_py 2011-09-27 16:50:30 UTC (rev 618) +++ release-info/run_tests_py 2011-09-27 17:38:50 UTC (rev 619) @@ -19,8 +19,8 @@ export CCTK_TESTSUITE_RUN_PROCESSORS=2 export OMP_NUM_THREADS=2;; esac - simfactory/bin/sim create-submit --testsuite ${HOSTNAME}__${CCTK_TESTSUITE_RUN_PROCESSORS}_${OMP_NUM_THREADS} --procs ${CCTK_TESTSUITE_RUN_PROCESSORS} --num-threads ${OMP_NUM_THREADS} --walltime 3:00:00 + simfactory/bin/sim create-submit --testsuite ${HOSTNAME}__${CCTK_TESTSUITE_RUN_PROCESSORS}_${OMP_NUM_THREADS} --procs ${CCTK_TESTSUITE_RUN_PROCESSORS} --num-threads ${OMP_NUM_THREADS} ---ppn-used ${CCTK_TESTSUITE_RUN_PROCESSORS} -walltime 3:00:00 # additional runs hostname | grep ^numrel && - ./simfactory/bin/sim create-submit --configuration gcc --testsuite ${HOSTNAME}-gcc__${CCTK_TESTSUITE_RUN_PROCESSORS}_${OMP_NUM_THREADS} --procs ${CCTK_TESTSUITE_RUN_PROCESSORS} --num-threads ${OMP_NUM_THREADS} --walltime 3:00:00 + ./simfactory/bin/sim create-submit --configuration gcc --testsuite ${HOSTNAME}-gcc__${CCTK_TESTSUITE_RUN_PROCESSORS}_${OMP_NUM_THREADS} --procs ${CCTK_TESTSUITE_RUN_PROCESSORS} --num-threads ${OMP_NUM_THREADS} --ppn-used ${CCTK_TESTSUITE_RUN_PROCESSORS} --walltime 3:00:00 done

1 0

[svn:einsteintoolkit] www/ (Rev. 618)
by roland.haas＠physics.gatech.edu 27 Sep '11

27 Sep '11

User: rhaas Date: 2011/09/27 11:50 AM Added: /release-info/ collect_results_py Log: add script to collect test results out of simfactories ouptut directories Directory Changes: Directory: /svn:executable/ =========================== + * File Changes: Directory: /release-info/ ========================= File [added]: collect_results_py Delta lines: +9 -0 =================================================================== --- release-info/collect_results_py (rev 0) +++ release-info/collect_results_py 2011-09-27 16:50:30 UTC (rev 618) @@ -0,0 +1,9 @@ +#!/bin/bash + +set -e # Abort on errors + +cd ET_2011_10 + +for sim in $(./simfactory/bin/sim list-simulations --name-only) ; do + cp $(./simfactory/bin/sim get-output-dir $sim | awk '/Output directory:/{print $3}')/TEST/sim/summary.log ../$sim.log +done Property changes on: release-info/collect_results_py ___________________________________________________________________

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 617)
by roland.haas＠physics.gatech.edu 27 Sep '11

27 Sep '11

User: rhaas Date: 2011/09/27 11:46 AM Added: /release-info/ bethe-intel11.1__1_1.log, bethe-intel11.1__2_1.log, bethe-intel11.1__2_2.log, zwicky-intel11.1__1_1.log, zwicky-intel11.1__2_1.log, zwicky-intel11.1__2_2.log Removed: /release-info/ bethe.tapir.caltech.edu__1_1.log, bethe.tapir.caltech.edu__2_1.log, bethe.tapir.caltech.edu__2_2.log Modified: /release-info/ run_tests_py Log: add results for zwicky the local Caltech cluster File Changes: Directory: /release-info/ ========================= File [added]: bethe-intel11.1__1_1.log Delta lines: <Unavailable> File [added]: bethe-intel11.1__2_1.log Delta lines: <Unavailable> File [added]: bethe-intel11.1__2_2.log Delta lines: <Unavailable> File [removed]: bethe.tapir.caltech.edu__1_1.log Delta lines: <Unavailable> File [removed]: bethe.tapir.caltech.edu__2_1.log Delta lines: <Unavailable> File [removed]: bethe.tapir.caltech.edu__2_2.log Delta lines: <Unavailable> File [modified]: run_tests_py Delta lines: <Unavailable> File [added]: zwicky-intel11.1__1_1.log Delta lines: <Unavailable> File [added]: zwicky-intel11.1__2_1.log Delta lines: <Unavailable> File [added]: zwicky-intel11.1__2_2.log Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 113)
by knarf＠cct.lsu.edu 26 Sep '11

26 Sep '11

User: knarf Date: 2011/09/26 07:51 PM Modified: /trunk/ einsteintoolkit.th Log: sort alphabetically, remove duplicate thorns File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.th Delta lines: +14 -21 =================================================================== --- trunk/einsteintoolkit.th 2011-09-20 17:29:26 UTC (rev 112) +++ trunk/einsteintoolkit.th 2011-09-27 00:51:07 UTC (rev 113) @@ -213,8 +213,21 @@ EinsteinUtils/SetMask_SphericalSurface EinsteinUtils/TGRtensor +PITTNullCode/NullConstr +PITTNullCode/NullDecomp +PITTNullCode/NullEvolve +PITTNullCode/NullExact +PITTNullCode/NullGrid +PITTNullCode/NullInterp +PITTNullCode/NullNews +PITTNullCode/NullPsiInt +PITTNullCode/NullSHRExtract +PITTNullCode/NullVars +PITTNullCode/SphericalHarmonicDecomp +PITTNullCode/SphericalHarmonicRecon + # Additional Cactus thorns !TARGET = $ARR !TYPE = svn @@ -222,8 +235,8 @@ !URL = http://svn.cactuscode.org/projects/$1/$2/trunk !CHECKOUT = ExternalLibraries/BLAS -ExternalLibraries/GSL ExternalLibraries/FFTW3 +ExternalLibraries/GSL ExternalLibraries/HDF5 ExternalLibraries/LAPACK ExternalLibraries/libjpeg @@ -318,26 +331,6 @@ -# CCE -!TARGET = $ARR -!TYPE = svn -!AUTH_URL = https://svn.einsteintoolkit.org/cactus/$1/$2/trunk -!URL = http://svn.einsteintoolkit.org/cactus/$1/$2/trunk -!CHECKOUT = -PITTNullCode/NullConstr -PITTNullCode/NullDecomp -PITTNullCode/NullEvolve -PITTNullCode/NullExact -PITTNullCode/NullGrid -PITTNullCode/NullInterp -PITTNullCode/NullNews -PITTNullCode/NullPsiInt -PITTNullCode/NullSHRExtract -PITTNullCode/NullVars -PITTNullCode/SphericalHarmonicRecon - - - # Private thorns !TARGET = $ARR !TYPE = ignore

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 616)
by roland.haas＠physics.gatech.edu 26 Sep '11

26 Sep '11

User: rhaas Date: 2011/09/26 05:24 PM Added: /release-info/ bethe.tapir.caltech.edu__1_1.log, bethe.tapir.caltech.edu__2_1.log, bethe.tapir.caltech.edu__2_2.log, compile_tests_py, run_tests_py Log: add results for bethe add initial implementation of Frank's run_tests and compile_tests for simfactory/py (both end in _py). File Changes: Directory: /release-info/ ========================= File [added]: bethe.tapir.caltech.edu__1_1.log Delta lines: <Unavailable> File [added]: bethe.tapir.caltech.edu__2_1.log Delta lines: <Unavailable> File [added]: bethe.tapir.caltech.edu__2_2.log Delta lines: <Unavailable> File [added]: compile_tests_py Delta lines: <Unavailable> File [added]: run_tests_py Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] ADMBase/trunk/ (Rev. 63)
by bcmsma＠astro.rit.edu 26 Sep '11

26 Sep '11

User: bmundim Date: 2011/09/26 05:02 PM Modified: /trunk/ param.ccl Log: Change default value of initial_shift parameter to "zero" according to Trac Ticket #565. File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +1 -1 =================================================================== --- trunk/param.ccl 2011-08-17 22:07:20 UTC (rev 62) +++ trunk/param.ccl 2011-09-26 22:02:25 UTC (rev 63) @@ -17,7 +17,7 @@ { "none" :: "Shift is inactive" "zero" :: "Shift is zero" -} "none" +} "zero" KEYWORD initial_dtlapse "Initial dtlapse value" {

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 615)
by knarf＠cct.lsu.edu 26 Sep '11

26 Sep '11

User: knarf Date: 2011/09/26 02:12 PM Removed: /release-info/ bluedrop__1.log, bluedrop__1_1.log, bluedrop__2.log, bluedrop__2_1.log, bluedrop__2_2.log Log: remove bluedrop from current machine list File Changes: Directory: /release-info/ ========================= File [removed]: bluedrop__1.log Delta lines: <Unavailable> File [removed]: bluedrop__1_1.log Delta lines: <Unavailable> File [removed]: bluedrop__2.log Delta lines: <Unavailable> File [removed]: bluedrop__2_1.log Delta lines: <Unavailable> File [removed]: bluedrop__2_2.log Delta lines: <Unavailable>

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits