{"header":{"message":"HydroToyCarpetX: Experiment with Kokkos\n","description":null,"timestamp":"2019-12-16T15:17:01Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC","authors":[],"deps":{"hashes":[["4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC","HydroToyCarpetX: Experiment with Kokkos\n"],["AA",null],["CCU5R34POG34UUEP3KF5GAWXNSUX55DXQKNNZCFWI5EKULTTZY6AC","HydroToyCarpetX: Don't require NSIMD (yet?)\n"],["GKKJ75HX2ERLVBZVE2CUB6T3J2SUT7R3UKEKTEYNOG43ZKX6X5MQC","The big rename: AMReX -> CarpetX\n"],["RZOODNAYNAGAXHNLGMDDPYVIBNFJBDYMWO6I3C3J5LEZCMAQUZRAC","HydroToyCarpetX: Vectorize functions\n"]],"deps":["CCU5R34POG34UUEP3KF5GAWXNSUX55DXQKNNZCFWI5EKULTTZY6AC"],"known":[]},"hunks":[{"Edit":{"inode":{"change":4,"pos":4614},"path":"HydroToyCarpetX/configuration.ccl","line":3,"atom":{"NewVertex":{"up_context":[{"change":5,"pos":1}],"down_context":[{"change":2,"pos":0}],"flag":"BLOCK","start":0,"end":18,"inode":{"change":4,"pos":4614}}},"contents":[{"Add":{"text":"# REQUIRES Kokkos","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:17"}}]}},{"FileAdd":{"meta":{"basename":"evolve-kokkos.cxx","metadata":420},"path":"HydroToyCarpetX/src/evolve-kokkos.cxx","context":[{"change":4,"pos":9280}],"contents":[{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:59"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:60"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:87"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:88"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:109"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:110"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:128"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:164"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:193"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:194"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:215"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:232"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:251"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:270"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:271"}},{"Add":{"text":"namespace HydroToyCarpetX {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:299"}},{"Add":{"text":"using namespace std;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:320"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:321"}},{"Add":{"text":"constexpr int dim = 3;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:344"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:345"}},{"Add":{"text":"extern \"C\" void HydroToyCarpetX_Evolve(CCTK_ARGUMENTS) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:402"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_HydroToyCarpetX_Evolve;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:451"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:478"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:479"}},{"Add":{"text":" const CCTK_REAL dx = CCTK_DELTA_SPACE(0);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:523"}},{"Add":{"text":" const CCTK_REAL dy = CCTK_DELTA_SPACE(1);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:567"}},{"Add":{"text":" const CCTK_REAL dz = CCTK_DELTA_SPACE(2);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:611"}},{"Add":{"text":" const CCTK_REAL dV1 = 1 / (dx * dy * dz);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:655"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:656"}},{"Add":{"text":" const Loop::vect di{{1, 0, 0}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:700"}},{"Add":{"text":" const Loop::vect dj{{0, 1, 0}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:744"}},{"Add":{"text":" const Loop::vect dk{{0, 0, 1}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:788"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:789"}},{"Add":{"text":" const Loop::vect lsh{{cctk_lsh[0], cctk_lsh[1], cctk_lsh[2]}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:864"}},{"Add":{"text":" const Loop::vect nghostzones{","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:906"}},{"Add":{"text":" {cctk_nghostzones[0], cctk_nghostzones[1], cctk_nghostzones[2]}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:978"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:979"}},{"Add":{"text":" array tmin_, tmax_;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1016"}},{"Add":{"text":" GetTileExtent(cctkGH, tmin_.data(), tmax_.data());","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1069"}},{"Add":{"text":" const Loop::vect tmin(tmin_);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1111"}},{"Add":{"text":" const Loop::vect tmax(tmax_);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1153"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1154"}},{"Add":{"text":" const Loop::vect imin = max(tmin, nghostzones);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1214"}},{"Add":{"text":" const Loop::vect imax = min(tmax, lsh - nghostzones);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1280"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1281"}},{"Add":{"text":" // regular, cell centred variables","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1318"}},{"Add":{"text":" const Loop::vect ash{{cctk_ash[0], cctk_ash[1], cctk_ash[2]}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1393"}},{"Add":{"text":" // variables without ghosts","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1423"}},{"Add":{"text":" const Loop::vect ash_ng = ash - 2 * nghostzones;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1484"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1485"}},{"Add":{"text":" // fluxes :face-centred without ghosts","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1526"}},{"Add":{"text":" const Loop::vect ash_fx = ash_ng + di;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1577"}},{"Add":{"text":" const Loop::vect ash_fy = ash_ng + dj;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1628"}},{"Add":{"text":" const Loop::vect ash_fz = ash_ng + dk;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1679"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1680"}},{"Add":{"text":" const auto mkstr{[](const Loop::vect &ash) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1737"}},{"Add":{"text":" Loop::vect str;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1773"}},{"Add":{"text":" ptrdiff_t s = 1;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1794"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1830"}},{"Add":{"text":" str[d] = s;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1848"}},{"Add":{"text":" s *= ash[d];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1867"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1873"}},{"Add":{"text":" return str;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1889"}},{"Add":{"text":" }};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1895"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1896"}},{"Add":{"text":" const Loop::vect str = mkstr(ash);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:1949"}},{"Add":{"text":" const Loop::vect str_fx = mkstr(ash_fx);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2008"}},{"Add":{"text":" const Loop::vect str_fy = mkstr(ash_fy);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2067"}},{"Add":{"text":" const Loop::vect str_fz = mkstr(ash_fz);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2126"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2127"}},{"Add":{"text":" constexpr ptrdiff_t off = 0;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2158"}},{"Add":{"text":" const ptrdiff_t off_fx = -str_fx[0] - str_fx[1] - str_fx[2];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2221"}},{"Add":{"text":" const ptrdiff_t off_fy = -str_fy[0] - str_fy[1] - str_fy[2];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2284"}},{"Add":{"text":" const ptrdiff_t off_fz = -str_fz[0] - str_fz[1] - str_fz[2];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2347"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2348"}},{"Add":{"text":" typedef vectype CCTK_VEC_REAL;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2392"}},{"Add":{"text":" constexpr int VS = vecprops::size();","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2442"}},{"Add":{"text":" // assert((imax[0] - imin[0]) % VS == 0);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2486"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2487"}},{"Add":{"text":" const auto vec_dV1 = CCTK_VEC_REAL::set1(dV1);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2536"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2537"}},{"Add":{"text":" // Transport","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2552"}},{"Add":{"text":" // dt rho + d_i (rho vel^i) = 0","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2586"}},{"Add":{"text":" // dt mom_j + d_i (mom_j vel^i) = 0","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2624"}},{"Add":{"text":" // dt etot + d_i (etot vel^i) = 0","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2660"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2661"}},{"Add":{"text":" const auto calcupdate{[&](auto &u_, const auto &u_p_, const auto &fx_,","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2734"}},{"Add":{"text":" const auto &fy_, const auto &fz_) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2798"}},{"Add":{"text":" for (int k = imin[2]; k < imax[2]; ++k) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2844"}},{"Add":{"text":" for (int j = imin[1]; j < imax[1]; ++j) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2892"}},{"Add":{"text":" ptrdiff_t idx_i0 = off + str[1] * j + str[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:2950"}},{"Add":{"text":" ptrdiff_t idx_fx_i0 = off_fx + str_fx[1] * j + str_fx[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3020"}},{"Add":{"text":" ptrdiff_t idx_fy_i0 = off_fy + str_fy[1] * j + str_fy[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3090"}},{"Add":{"text":" ptrdiff_t idx_fz_i0 = off_fz + str_fz[1] * j + str_fz[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3160"}},{"Add":{"text":" for (int i = imin[0]; i < imax[0]; i += VS) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3214"}},{"Add":{"text":" ptrdiff_t idx = idx_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3252"}},{"Add":{"text":" ptrdiff_t idx_fx = idx_fx_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3296"}},{"Add":{"text":" ptrdiff_t idx_fy = idx_fy_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3340"}},{"Add":{"text":" ptrdiff_t idx_fz = idx_fz_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3384"}},{"Add":{"text":" auto u_p = CCTK_VEC_REAL::loadu(u_p_[idx]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3438"}},{"Add":{"text":" auto fx1 = CCTK_VEC_REAL::loadu(fx_[idx_fx + str_fx[0]]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3506"}},{"Add":{"text":" auto fx0 = CCTK_VEC_REAL::loadu(fx_[idx_fx]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3562"}},{"Add":{"text":" auto fy1 = CCTK_VEC_REAL::loadu(fy_[idx_fy + str_fy[1]]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3630"}},{"Add":{"text":" auto fy0 = CCTK_VEC_REAL::loadu(fy_[idx_fy]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3686"}},{"Add":{"text":" auto fz1 = CCTK_VEC_REAL::loadu(fz_[idx_fz + str_fz[2]]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3754"}},{"Add":{"text":" auto fz0 = CCTK_VEC_REAL::loadu(fz_[idx_fz]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3810"}},{"Add":{"text":" auto u = u_p - vec_dV1 * ((fx1 - fx0) + (fy1 - fy0) + (fz1 - fz0));","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3888"}},{"Add":{"text":" u.storeu_partial(u_[idx], i, i, imax[0]);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3940"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3950"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3958"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3964"}},{"Add":{"text":" }};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3970"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:3971"}},{"Add":{"text":" calcupdate(rho, rho_p, fxrho, fyrho, fzrho);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4018"}},{"Add":{"text":" calcupdate(momx, momx_p, fxmomx, fymomx, fzmomx);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4070"}},{"Add":{"text":" calcupdate(momy, momy_p, fxmomy, fymomy, fzmomy);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4122"}},{"Add":{"text":" calcupdate(momz, momz_p, fxmomz, fymomz, fzmomz);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4174"}},{"Add":{"text":" calcupdate(etot, etot_p, fxetot, fyetot, fzetot);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4226"}},{"Add":{"text":"}","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4228"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4229"}},{"Add":{"text":"} // namespace HydroToyCarpetX","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4260"}}]}},{"FileAdd":{"meta":{"basename":"evolve-openmp.cxx","metadata":420},"path":"HydroToyCarpetX/src/evolve-openmp.cxx","context":[{"change":4,"pos":9280}],"contents":[{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4302"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4303"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4321"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4357"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4386"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4387"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4408"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4425"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4444"}},{"Add":{"text":"#include ","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4463"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4464"}},{"Add":{"text":"namespace HydroToyCarpetX {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4492"}},{"Add":{"text":"using namespace std;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4513"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4514"}},{"Add":{"text":"constexpr int dim = 3;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4537"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4538"}},{"Add":{"text":"extern \"C\" void HydroToyCarpetX_Evolve(CCTK_ARGUMENTS) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4595"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_HydroToyCarpetX_Evolve;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4644"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4671"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4672"}},{"Add":{"text":" const CCTK_REAL dx = CCTK_DELTA_SPACE(0);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4716"}},{"Add":{"text":" const CCTK_REAL dy = CCTK_DELTA_SPACE(1);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4760"}},{"Add":{"text":" const CCTK_REAL dz = CCTK_DELTA_SPACE(2);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4804"}},{"Add":{"text":" const CCTK_REAL dV1 = 1 / (dx * dy * dz);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4848"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4849"}},{"Add":{"text":" const Loop::vect di{{1, 0, 0}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4893"}},{"Add":{"text":" const Loop::vect dj{{0, 1, 0}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4937"}},{"Add":{"text":" const Loop::vect dk{{0, 0, 1}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4981"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:4982"}},{"Add":{"text":" const Loop::vect lsh{{cctk_lsh[0], cctk_lsh[1], cctk_lsh[2]}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5057"}},{"Add":{"text":" const Loop::vect nghostzones{","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5099"}},{"Add":{"text":" {cctk_nghostzones[0], cctk_nghostzones[1], cctk_nghostzones[2]}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5171"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5172"}},{"Add":{"text":" array tmin_, tmax_;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5209"}},{"Add":{"text":" GetTileExtent(cctkGH, tmin_.data(), tmax_.data());","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5262"}},{"Add":{"text":" const Loop::vect tmin(tmin_);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5304"}},{"Add":{"text":" const Loop::vect tmax(tmax_);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5346"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5347"}},{"Add":{"text":" const Loop::vect imin = max(tmin, nghostzones);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5407"}},{"Add":{"text":" const Loop::vect imax = min(tmax, lsh - nghostzones);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5473"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5474"}},{"Add":{"text":" // regular, cell centred variables","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5511"}},{"Add":{"text":" const Loop::vect ash{{cctk_ash[0], cctk_ash[1], cctk_ash[2]}};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5586"}},{"Add":{"text":" // variables without ghosts","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5616"}},{"Add":{"text":" const Loop::vect ash_ng = ash - 2 * nghostzones;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5677"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5678"}},{"Add":{"text":" // fluxes :face-centred without ghosts","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5719"}},{"Add":{"text":" const Loop::vect ash_fx = ash_ng + di;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5770"}},{"Add":{"text":" const Loop::vect ash_fy = ash_ng + dj;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5821"}},{"Add":{"text":" const Loop::vect ash_fz = ash_ng + dk;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5872"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5873"}},{"Add":{"text":" const auto mkstr{[](const Loop::vect &ash) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5930"}},{"Add":{"text":" Loop::vect str;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5966"}},{"Add":{"text":" ptrdiff_t s = 1;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:5987"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6023"}},{"Add":{"text":" str[d] = s;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6041"}},{"Add":{"text":" s *= ash[d];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6060"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6066"}},{"Add":{"text":" return str;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6082"}},{"Add":{"text":" }};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6088"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6089"}},{"Add":{"text":" const Loop::vect str = mkstr(ash);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6142"}},{"Add":{"text":" const Loop::vect str_fx = mkstr(ash_fx);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6201"}},{"Add":{"text":" const Loop::vect str_fy = mkstr(ash_fy);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6260"}},{"Add":{"text":" const Loop::vect str_fz = mkstr(ash_fz);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6319"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6320"}},{"Add":{"text":" constexpr ptrdiff_t off = 0;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6351"}},{"Add":{"text":" const ptrdiff_t off_fx = -str_fx[0] - str_fx[1] - str_fx[2];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6414"}},{"Add":{"text":" const ptrdiff_t off_fy = -str_fy[0] - str_fy[1] - str_fy[2];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6477"}},{"Add":{"text":" const ptrdiff_t off_fz = -str_fz[0] - str_fz[1] - str_fz[2];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6540"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6541"}},{"Add":{"text":" constexpr int VS = 4; // TODO: choose automatically","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6595"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6596"}},{"Add":{"text":" // Transport","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6611"}},{"Add":{"text":" // dt rho + d_i (rho vel^i) = 0","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6645"}},{"Add":{"text":" // dt mom_j + d_i (mom_j vel^i) = 0","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6683"}},{"Add":{"text":" // dt etot + d_i (etot vel^i) = 0","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6719"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6720"}},{"Add":{"text":" const auto calcupdate{[&](auto &u_, const auto &u_p_, const auto &fx_,","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6793"}},{"Add":{"text":" const auto &fy_, const auto &fz_) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6857"}},{"Add":{"text":" for (int k = imin[2]; k < imax[2]; ++k) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6903"}},{"Add":{"text":" for (int j = imin[1]; j < imax[1]; ++j) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:6951"}},{"Add":{"text":" ptrdiff_t idx_i0 = off + str[1] * j + str[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7009"}},{"Add":{"text":" ptrdiff_t idx_fx_i0 = off_fx + str_fx[1] * j + str_fx[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7079"}},{"Add":{"text":" ptrdiff_t idx_fy_i0 = off_fy + str_fy[1] * j + str_fy[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7149"}},{"Add":{"text":" ptrdiff_t idx_fz_i0 = off_fz + str_fz[1] * j + str_fz[2] * k;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7219"}},{"Add":{"text":" int ni = (imax[0] - imin[0] + VS - 1) / VS * VS;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7276"}},{"Add":{"text":" assert(ni % VS == 0);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7306"}},{"Add":{"text":"#pragma omp simd simdlen(VS)","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7335"}},{"Add":{"text":" // for (int i = imin[0]; i < imax[0]; ++i) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7388"}},{"Add":{"text":" for (int i0 = 0; i0 < ni; ++i0) {","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7430"}},{"Add":{"text":" int i = imin[0] + i0;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7462"}},{"Add":{"text":" ptrdiff_t idx = idx_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7500"}},{"Add":{"text":" ptrdiff_t idx_fx = idx_fx_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7544"}},{"Add":{"text":" ptrdiff_t idx_fy = idx_fy_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7588"}},{"Add":{"text":" ptrdiff_t idx_fz = idx_fz_i0 + i;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7632"}},{"Add":{"text":" auto u_p = u_p_[idx];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7664"}},{"Add":{"text":" auto fx1 = fx_[idx_fx + str_fx[0]];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7710"}},{"Add":{"text":" auto fx0 = fx_[idx_fx];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7744"}},{"Add":{"text":" auto fy1 = fy_[idx_fy + str_fy[1]];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7790"}},{"Add":{"text":" auto fy0 = fy_[idx_fy];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7824"}},{"Add":{"text":" auto fz1 = fz_[idx_fz + str_fz[2]];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7870"}},{"Add":{"text":" auto fz0 = fz_[idx_fz];","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7904"}},{"Add":{"text":" auto u = u_p - dV1 * ((fx1 - fx0) + (fy1 - fy0) + (fz1 - fz0));","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:7978"}},{"Add":{"text":" if (i < imax[0])","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8005"}},{"Add":{"text":" u_[idx] = u;","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8030"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8040"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8048"}},{"Add":{"text":" }","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8054"}},{"Add":{"text":" }};","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8060"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8061"}},{"Add":{"text":" calcupdate(rho, rho_p, fxrho, fyrho, fzrho);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8108"}},{"Add":{"text":" calcupdate(momx, momx_p, fxmomx, fymomx, fzmomx);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8160"}},{"Add":{"text":" calcupdate(momy, momy_p, fxmomy, fymomy, fzmomy);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8212"}},{"Add":{"text":" calcupdate(momz, momz_p, fxmomz, fymomz, fzmomz);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8264"}},{"Add":{"text":" calcupdate(etot, etot_p, fxetot, fyetot, fzetot);","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8316"}},{"Add":{"text":"}","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8318"}},{"Add":{"text":"","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8319"}},{"Add":{"text":"} // namespace HydroToyCarpetX","pos":"4GFYOLM5UKLGVSDJD6CDHEA23QBUIE6X67REX7C7PW5RUJLKZHVAC:8350"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review