{"header":{"message":"ODESolvers: Implement Euler and RK2 ODE integrators\n","description":null,"timestamp":"2019-11-11T18:25:06Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC","authors":[],"deps":{"hashes":[["P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC","ODESolvers: Implement Euler and RK2 ODE integrators\n"],["AA",null],["W5IGLN4UQIR5MF5MIQT25JQZ5F2DGRRHJFLN3WDL3T6PLSWRS5DAC","New thorn ODESolvers\n"]],"deps":["W5IGLN4UQIR5MF5MIQT25JQZ5F2DGRRHJFLN3WDL3T6PLSWRS5DAC"],"known":[]},"hunks":[{"Edit":{"inode":{"change":2,"pos":343},"path":"ODEsolvers/configuration.ccl","line":2,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":393}],"down_context":[],"flag":"BLOCK","start":0,"end":24,"inode":{"change":2,"pos":343}}},"contents":[{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:0"}},{"Add":{"text":"REQUIRES AMReX CarpetX","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:23"}}]}},{"Edit":{"inode":{"change":2,"pos":410},"path":"ODEsolvers/interface.ccl","line":4,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":479}],"down_context":[],"flag":"BLOCK","start":25,"end":151,"inode":{"change":2,"pos":410}}},"contents":[{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:25"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:26"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:27"}},{"Add":{"text":"void FUNCTION CallScheduleGroup(","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:60"}},{"Add":{"text":" CCTK_POINTER IN cctkGH,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:86"}},{"Add":{"text":" CCTK_STRING IN groupname)","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:114"}},{"Add":{"text":"REQUIRES FUNCTION CallScheduleGroup","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:150"}}]}},{"Edit":{"inode":{"change":2,"pos":492},"path":"ODEsolvers/param.ccl","line":2,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":538}],"down_context":[],"flag":"BLOCK","start":152,"end":264,"inode":{"change":2,"pos":492}}},"contents":[{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:152"}},{"Add":{"text":"KEYWORD method \"ODE solver method\" STEERABLE=always","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:204"}},{"Add":{"text":"{","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:206"}},{"Add":{"text":" \"Constant\" :: \"\"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:225"}},{"Add":{"text":" \"Euler\" :: \"\"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:241"}},{"Add":{"text":" \"RK2\" :: \"\"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:255"}},{"Add":{"text":"} \"RK2\"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:263"}}]}},{"Edit":{"inode":{"change":2,"pos":554},"path":"ODEsolvers/schedule.ccl","line":2,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":599}],"down_context":[],"flag":"BLOCK","start":265,"end":510,"inode":{"change":2,"pos":554}}},"contents":[{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:265"}},{"Add":{"text":"SCHEDULE ODESolvers_solve AT evol","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:299"}},{"Add":{"text":"{","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:301"}},{"Add":{"text":" LANG: C","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:311"}},{"Add":{"text":" OPTIONS: level","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:328"}},{"Add":{"text":"} \"Solve ODEs loop\"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:348"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:349"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:350"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:351"}},{"Add":{"text":"# This group is not scheduled explicitly here; it will be scheduled dynamically","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:431"}},{"Add":{"text":"# from ODESolvers_solve","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:456"}},{"Add":{"text":"SCHEDULE GROUP ODESolvers_RHS","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:486"}},{"Add":{"text":"{","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:488"}},{"Add":{"text":"} \"Evaluate ODE RHS\"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:509"}}]}},{"Replacement":{"inode":{"change":2,"pos":624},"path":"ODEsolvers/src/make.code.defn","line":4,"ins":{"up_context":[{"change":2,"pos":707}],"down_context":[{"change":2,"pos":714}],"flag":"BLOCK","start":511,"end":528,"inode":{"change":2,"pos":624}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":707},"to":{"change":2,"start":707,"end":714},"introduced_by":2}],"inode":{"change":2,"pos":624}},"ins_contents":[{"Add":{"text":"SRCS = solve.cxx","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:527"}}],"del_contents":[{"Del":{"text":"SRCS =","pos":"W5IGLN4UQIR5MF5MIQT25JQZ5F2DGRRHJFLN3WDL3T6PLSWRS5DAC:713"}}]}},{"FileAdd":{"meta":{"basename":"solve.cxx","metadata":420},"path":"ODEsolvers/src/solve.cxx","context":[{"change":2,"pos":606}],"contents":[{"Add":{"text":"#include <../../CarpetX/src/driver.hxx>","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:581"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:582"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:600"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:629"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:665"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:666"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:694"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:695"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:714"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:736"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:754"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:773"}},{"Add":{"text":"#include ","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:791"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:792"}},{"Add":{"text":"namespace ODESolver {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:814"}},{"Add":{"text":"using namespace amrex;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:837"}},{"Add":{"text":"using namespace std;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:858"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:859"}},{"Add":{"text":"////////////////////////////////////////////////////////////////////////////////","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:940"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:941"}},{"Add":{"text":"// A state vector component, with mfabs for each level, group, and variable","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1017"}},{"Add":{"text":"struct statecomp_t {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1038"}},{"Add":{"text":" vector mfabs;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1066"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1067"}},{"Add":{"text":"private:","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1076"}},{"Add":{"text":" // These will be automatically freed when this object is deallocated","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1147"}},{"Add":{"text":" vector > owned_stuff;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1192"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1193"}},{"Add":{"text":"public:","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1201"}},{"Add":{"text":" statecomp_t copy() const;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1229"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1230"}},{"Add":{"text":" static void axpy(const statecomp_t &y, CCTK_REAL alpha, const statecomp_t &x);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1311"}},{"Add":{"text":" static void lincomb(const statecomp_t &z, CCTK_REAL alpha,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1372"}},{"Add":{"text":" const statecomp_t &x, CCTK_REAL beta,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1432"}},{"Add":{"text":" const statecomp_t &y);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1477"}},{"Add":{"text":"};","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1480"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1481"}},{"Add":{"text":"////////////////////////////////////////////////////////////////////////////////","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1562"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1563"}},{"Add":{"text":"// Copy state vector into newly allocate memory","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1611"}},{"Add":{"text":"statecomp_t statecomp_t::copy() const {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1651"}},{"Add":{"text":" const size_t size = mfabs.size();","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1687"}},{"Add":{"text":" statecomp_t result;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1709"}},{"Add":{"text":" result.mfabs.reserve(size);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1739"}},{"Add":{"text":" result.owned_stuff.reserve(size);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1775"}},{"Add":{"text":" for (size_t n = 0; n < size; ++n) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1813"}},{"Add":{"text":" const auto &x = mfabs.at(n);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1846"}},{"Add":{"text":" auto y = make_unique(x->boxArray(), x->DistributionMap(),","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1918"}},{"Add":{"text":" x->nComp(), x->nGrowVect());","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1982"}},{"Add":{"text":" MultiFab::Copy(*y, *x, 0, 0, y->nComp(), y->nGrowVect());","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2044"}},{"Add":{"text":" result.mfabs.push_back(y.get());","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2081"}},{"Add":{"text":" result.owned_stuff.push_back(move(y));","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2124"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2128"}},{"Add":{"text":" return result;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2145"}},{"Add":{"text":"}","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2147"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2148"}},{"Add":{"text":"// y += alpha * x","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2166"}},{"Add":{"text":"void statecomp_t::axpy(const statecomp_t &y, CCTK_REAL alpha,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2228"}},{"Add":{"text":" const statecomp_t &x) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2275"}},{"Add":{"text":" const size_t size = y.mfabs.size();","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2313"}},{"Add":{"text":" assert(x.mfabs.size() == size);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2347"}},{"Add":{"text":" for (size_t n = 0; n < size; ++n) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2385"}},{"Add":{"text":" // TODO: Update Saxpy call when using AMReX 19.12 or later","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2448"}},{"Add":{"text":" const auto nghosts = y.mfabs.at(n)->nGrowVect();","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2501"}},{"Add":{"text":" for (int d = 0; d < 3; ++d)","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2533"}},{"Add":{"text":" assert(nghosts[d] == nghosts[0]);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2573"}},{"Add":{"text":" MultiFab::Saxpy(*y.mfabs.at(n), alpha, *x.mfabs.at(n), 0, 0,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2638"}},{"Add":{"text":" y.mfabs.at(n)->nComp(), nghosts[0]);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2695"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2699"}},{"Add":{"text":"}","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2701"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2702"}},{"Add":{"text":"// z = alpha * x + beta * y","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2730"}},{"Add":{"text":"void statecomp_t::lincomb(const statecomp_t &z, CCTK_REAL alpha,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2795"}},{"Add":{"text":" const statecomp_t &x, CCTK_REAL beta,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2859"}},{"Add":{"text":" const statecomp_t &y) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2909"}},{"Add":{"text":" const size_t size = z.mfabs.size();","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2947"}},{"Add":{"text":" assert(x.mfabs.size() == size);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:2981"}},{"Add":{"text":" assert(y.mfabs.size() == size);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3015"}},{"Add":{"text":" for (size_t n = 0; n < size; ++n) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3053"}},{"Add":{"text":" // TODO: Update Lincomb call when using AMReX 19.12 or later","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3118"}},{"Add":{"text":" const auto nghosts = z.mfabs.at(n)->nGrowVect();","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3171"}},{"Add":{"text":" for (int d = 0; d < 3; ++d)","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3203"}},{"Add":{"text":" assert(nghosts[d] == nghosts[0]);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3243"}},{"Add":{"text":" MultiFab::LinComb(*z.mfabs.at(n), alpha, *x.mfabs.at(n), 0, beta,","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3313"}},{"Add":{"text":" *y.mfabs.at(n), 0, 0, z.mfabs.at(n)->nComp(), nghosts[0]);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3394"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3398"}},{"Add":{"text":"}","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3400"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3401"}},{"Add":{"text":"////////////////////////////////////////////////////////////////////////////////","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3482"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3483"}},{"Add":{"text":"extern \"C\" void ODESolvers_solve(CCTK_ARGUMENTS) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3534"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_ODESolvers_solve;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3577"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3604"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3605"}},{"Add":{"text":" const CCTK_REAL dt = CCTK_DELTA_TIME;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3645"}},{"Add":{"text":" const int tl = 0;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3665"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3666"}},{"Add":{"text":" statecomp_t var, rhs;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3690"}},{"Add":{"text":" for (const auto &leveldata : CarpetX::ghext->leveldata) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3750"}},{"Add":{"text":" for (const auto &groupdata : leveldata.groupdata) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3806"}},{"Add":{"text":" if (groupdata.groupindex == CCTK_GroupIndex(\"TestODESolvers::state\")) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3884"}},{"Add":{"text":" const int rhs_gi = CCTK_GroupIndex(\"TestODESolvers::rhs\");","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:3951"}},{"Add":{"text":" const auto &rhs_groupdata = leveldata.groupdata.at(rhs_gi);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4019"}},{"Add":{"text":" for (int vi = 0; vi < groupdata.numvars; ++vi) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4076"}},{"Add":{"text":" var.mfabs.push_back(groupdata.mfab.at(tl).get());","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4136"}},{"Add":{"text":" rhs.mfabs.push_back(rhs_groupdata.mfab.at(tl).get());","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4200"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4210"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4218"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4224"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4228"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4229"}},{"Add":{"text":" const CCTK_REAL saved_time = cctkGH->cctk_time;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4279"}},{"Add":{"text":" *const_cast(&cctkGH->cctk_time) -= dt;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4333"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4334"}},{"Add":{"text":" if (CCTK_EQUALS(method, \"constant\")) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4375"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4376"}},{"Add":{"text":" // do nothing","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4394"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4395"}},{"Add":{"text":" } else if (CCTK_EQUALS(method, \"Euler\")) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4440"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4441"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_RHS\");","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4490"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4491"}},{"Add":{"text":" // Add scaled RHS to state vector","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4529"}},{"Add":{"text":" statecomp_t::axpy(var, dt, rhs);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4566"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4567"}},{"Add":{"text":" } else if (CCTK_EQUALS(method, \"RK2\")) {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4610"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4611"}},{"Add":{"text":" const auto old = var.copy();","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4644"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4645"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_RHS\");","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4694"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4695"}},{"Add":{"text":" // Add scaled RHS to state vector","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4733"}},{"Add":{"text":" statecomp_t::axpy(var, dt / 2, rhs);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4774"}},{"Add":{"text":" *const_cast(&cctkGH->cctk_time) += dt / 2;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4834"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4835"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_RHS\");","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4884"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4885"}},{"Add":{"text":" // Calculate new state vector","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4919"}},{"Add":{"text":" statecomp_t::lincomb(var, 1, old, dt, rhs);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4967"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4968"}},{"Add":{"text":" } else {","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4979"}},{"Add":{"text":" assert(0);","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4994"}},{"Add":{"text":" }","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4998"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:4999"}},{"Add":{"text":" // Reset current time","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5023"}},{"Add":{"text":" *const_cast(&cctkGH->cctk_time) = saved_time;","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5084"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5085"}},{"Add":{"text":" // TODO: Update time here, and not during time level cycling in the driver","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5162"}},{"Add":{"text":"}","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5164"}},{"Add":{"text":"","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5165"}},{"Add":{"text":"} // namespace ODESolver","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5190"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review