{"header":{"message":"ODESolvers: Implement RKF78 and DP87 integrators\n","description":null,"timestamp":"2020-06-27T22:00:54Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC","authors":[],"deps":{"hashes":[["P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC","ODESolvers: Implement RKF78 and DP87 integrators\n"],["AA",null],["75DPHH3IRB65ARCIQXSFGMU25R43IPMTZ2N4X54BBA75EJ2CX66QC","ODESolvers: Implement classic RK4\n"],["V37ACC7MAD57VQOKAT5FBO2ITNX73AWD2DUEIHMF7UPKV63YRTQAC","ODESolvers: Read RHS variables from interface.ccl\n"],["VLYBJZ6OP2XOVYCQIQOE2VHMX3DTWXYIPGZNYNME6F2H2EXXI2BQC","ODESolvers: Introduce ODESolvers_PostStep schedule group\n"],["3NWSBBMPHW3XOWLAMWI2AXX5VQDRHYT7CKGHL32RD6AUAUH67SVAC","ODESolvers: Reduce default verbosity\n"],["P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC","ODESolvers: Implement Euler and RK2 ODE integrators\n"],["2BSBB53DVGURXHAIBJQI2YO6V42HYBLMF5IULSCYMZQPN4EH7IMQC","ODESolvers: Check whether state vector and RHS contain valid data\n"],["W5IGLN4UQIR5MF5MIQT25JQZ5F2DGRRHJFLN3WDL3T6PLSWRS5DAC","New thorn ODESolvers\n"]],"deps":["75DPHH3IRB65ARCIQXSFGMU25R43IPMTZ2N4X54BBA75EJ2CX66QC"],"known":["P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC","V37ACC7MAD57VQOKAT5FBO2ITNX73AWD2DUEIHMF7UPKV63YRTQAC","2BSBB53DVGURXHAIBJQI2YO6V42HYBLMF5IULSCYMZQPN4EH7IMQC","VLYBJZ6OP2XOVYCQIQOE2VHMX3DTWXYIPGZNYNME6F2H2EXXI2BQC","3NWSBBMPHW3XOWLAMWI2AXX5VQDRHYT7CKGHL32RD6AUAUH67SVAC"]},"hunks":[{"Edit":{"inode":{"change":9,"pos":492},"path":"ODESolvers/param.ccl","line":5,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":207}],"down_context":[{"change":2,"pos":0}],"flag":"BLOCK","start":0,"end":70,"inode":{"change":9,"pos":492}}},"contents":[{"Add":{"text":" # See ","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:69"}}]}},{"Edit":{"inode":{"change":9,"pos":492},"path":"ODESolvers/param.ccl","line":12,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":137}],"down_context":[{"change":3,"pos":256}],"flag":"BLOCK","start":71,"end":148,"inode":{"change":9,"pos":492}}},"contents":[{"Add":{"text":" \"RKF78\" :: \"Runge-Kutta-Fehlberg 7(8)\"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:111"}},{"Add":{"text":" \"DP87\" :: \"Dormand & Prince 8(7)\"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:147"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":12,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":715}],"down_context":[{"change":3,"pos":715}],"flag":"BLOCK","start":149,"end":166,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":"#include ","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:165"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":14,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":737}],"down_context":[{"change":3,"pos":737}],"flag":"BLOCK","start":167,"end":185,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":"#include ","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:184"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":17,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":157}],"down_context":[{"change":3,"pos":755}],"flag":"BLOCK","start":186,"end":203,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":"#include ","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:202"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":42,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":1231}],"down_context":[{"change":3,"pos":1231}],"flag":"BLOCK","start":204,"end":270,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":" static void assign(const statecomp_t &y, const statecomp_t &x);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:269"}}]}},{"Replacement":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":69,"ins":{"up_context":[{"change":3,"pos":1564}],"down_context":[{"change":3,"pos":1612}],"flag":"BLOCK","start":271,"end":320,"inode":{"change":3,"pos":541}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":1564},"to":{"change":3,"start":1564,"end":1612},"introduced_by":3}],"inode":{"change":3,"pos":541}},"ins_contents":[{"Add":{"text":"// Copy state vector into newly allocated memory","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:319"}}],"del_contents":[{"Del":{"text":"// Copy state vector into newly allocate memory","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:1611"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":100,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":2149}],"down_context":[{"change":3,"pos":2149}],"flag":"BLOCK","start":321,"end":1144,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":"// y = x","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:329"}},{"Add":{"text":"void statecomp_t::assign(const statecomp_t &y, const statecomp_t &x) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:400"}},{"Add":{"text":" x.check_valid();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:419"}},{"Add":{"text":" y.check_valid();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:438"}},{"Add":{"text":" const size_t size = y.mfabs.size();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:476"}},{"Add":{"text":" assert(x.mfabs.size() == size);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:510"}},{"Add":{"text":" for (size_t n = 0; n < size; ++n) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:548"}},{"Add":{"text":" assert(x.mfabs.at(n)->nGrowVect() == y.mfabs.at(n)->nGrowVect());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:618"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:636"}},{"Add":{"text":" if (x.mfabs.at(n)->contains_nan())","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:675"}},{"Add":{"text":" CCTK_VERROR(\"statecomp_t::assign.x: Group %s contains nans\",","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:742"}},{"Add":{"text":" x.groupnames.at(n).c_str());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:789"}},{"Add":{"text":"#endif","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:796"}},{"Add":{"text":" amrex::MultiFab::Copy(*y.mfabs.at(n), *x.mfabs.at(n), 0, 0,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:860"}},{"Add":{"text":" y.mfabs.at(n)->nComp(), y.mfabs.at(n)->nGrowVect());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:939"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:957"}},{"Add":{"text":" if (y.mfabs.at(n)->contains_nan())","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:996"}},{"Add":{"text":" CCTK_VERROR(\"statecomp_t::assign.y: Group %s contains nans\",","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1063"}},{"Add":{"text":" y.groupnames.at(n).c_str());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1110"}},{"Add":{"text":"#endif","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1117"}},{"Add":{"text":" }","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1121"}},{"Add":{"text":" y.check_valid();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1140"}},{"Add":{"text":"}","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1142"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1143"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":441,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":1701}],"down_context":[{"change":3,"pos":4968}],"flag":"BLOCK","start":1145,"end":8220,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1145"}},{"Add":{"text":" } else if (CCTK_EQUALS(method, \"RKF78\")) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1190"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1191"}},{"Add":{"text":" typedef CCTK_REAL T;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1216"}},{"Add":{"text":" const auto R = [](T x, T y) { return x / y; };","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1267"}},{"Add":{"text":" const tuple > >, vector > tableau{","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1335"}},{"Add":{"text":" {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1345"}},{"Add":{"text":" {/* 1 */ 0, {}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1419"}},{"Add":{"text":" {/* 2 */ R(2, 27), {R(2, 27)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1493"}},{"Add":{"text":" {/* 3 */ R(1, 9), {R(1, 36), R(3, 36)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1567"}},{"Add":{"text":" {/* 4 */ R(1, 6), {R(1, 24), 0, R(3, 24)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1641"}},{"Add":{"text":" {/* 5 */ R(5, 12), {R(20, 48), 0, R(-75, 48), R(75, 48)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1715"}},{"Add":{"text":" {/* 6 */ R(1, 2), {R(1, 20), 0, 0, R(5, 20), R(4, 20)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1789"}},{"Add":{"text":" {/* 7 */ R(5, 6),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1819"}},{"Add":{"text":" {R(-25, 108), 0, 0, R(125, 108), R(-260, 108), R(250, 108)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1897"}},{"Add":{"text":" {/* 8 */ R(1, 6),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:1927"}},{"Add":{"text":" {R(31, 300), 0, 0, 0, R(61, 225), R(-2, 9), R(13, 900)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2001"}},{"Add":{"text":" {/* 9 */ R(2, 3),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2031"}},{"Add":{"text":" {2, 0, 0, R(-53, 6), R(704, 45), R(-107, 9), R(67, 90), 3}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2108"}},{"Add":{"text":" {/* 10 */ R(1, 3),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2139"}},{"Add":{"text":" {R(-91, 108), 0, 0, R(23, 108), R(-976, 135), R(311, 54),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2210"}},{"Add":{"text":" R(-19, 60), R(17, 6), R(-1, 12)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2262"}},{"Add":{"text":" {/* 11 */ 1,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2287"}},{"Add":{"text":" {R(2383, 4100), 0, 0, R(-341, 164), R(4496, 1025), R(-301, 82),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2364"}},{"Add":{"text":" R(2133, 4100), R(45, 82), R(45, 164), R(18, 41)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2432"}},{"Add":{"text":" // {/* 12 */ 0,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2513"}},{"Add":{"text":" // {R(3, 205), 0, 0, 0, 0, R(-6, 41), R(-3, 205), R(-3, 41), R(3,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2592"}},{"Add":{"text":" // 41),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2613"}},{"Add":{"text":" // R(6, 41)}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2645"}},{"Add":{"text":" // {/* 13 */ 1,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2673"}},{"Add":{"text":" // {R(-1777, 4100), 0, 0, R(-341, 164), R(4496, 1025), R(-289, 82),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2754"}},{"Add":{"text":" // R(2193, 4100), R(51, 82), R(33, 164), R(12, 41), 0, 1}}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2831"}},{"Add":{"text":" },","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2842"}},{"Add":{"text":" {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2852"}},{"Add":{"text":" R(41, 840), 0, 0, 0, 0, R(34, 105), R(9, 35), R(9, 35), R(9, 280),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2931"}},{"Add":{"text":" R(9, 280), R(41, 840),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2966"}},{"Add":{"text":" // 0,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:2984"}},{"Add":{"text":" // 0,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3002"}},{"Add":{"text":" }};","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3014"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3015"}},{"Add":{"text":" // Check Butcher tableau","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3044"}},{"Add":{"text":" const size_t nsteps = get<0>(tableau).size();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3094"}},{"Add":{"text":" {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3100"}},{"Add":{"text":" for (size_t step = 0; step < nsteps; ++step) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3153"}},{"Add":{"text":" // TODO: Could allow <=","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3185"}},{"Add":{"text":" assert(get<1>(get<0>(tableau).at(step)).size() == step);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3250"}},{"Add":{"text":" const auto &c = get<0>(get<0>(tableau).at(step));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3308"}},{"Add":{"text":" const auto &as = get<1>(get<0>(tableau).at(step));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3367"}},{"Add":{"text":" T x = 0;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3384"}},{"Add":{"text":" for (const auto &a : as)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3417"}},{"Add":{"text":" x += a;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3435"}},{"Add":{"text":" assert(fabs(x - c) <= 10 * numeric_limits::epsilon());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3501"}},{"Add":{"text":" }","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3509"}},{"Add":{"text":" // TODO: Could allow <=","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3539"}},{"Add":{"text":" assert(get<1>(tableau).size() == nsteps);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3587"}},{"Add":{"text":" const auto &bs = get<1>(tableau);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3627"}},{"Add":{"text":" T x = 0;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3642"}},{"Add":{"text":" for (const auto &b : bs)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3673"}},{"Add":{"text":" x += b;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3689"}},{"Add":{"text":" assert(fabs(x - 1) <= 10 * numeric_limits::epsilon());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3753"}},{"Add":{"text":" }","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3759"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3760"}},{"Add":{"text":" const auto old = var.copy();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3793"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3794"}},{"Add":{"text":" vector ks;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3822"}},{"Add":{"text":" ks.reserve(nsteps);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3846"}},{"Add":{"text":" // Special-case step 1 (why?)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3880"}},{"Add":{"text":" ks.push_back(rhs.copy());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3910"}},{"Add":{"text":" for (size_t step = 1; step < nsteps; ++step) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:3961"}},{"Add":{"text":" const auto &c = get<0>(get<0>(tableau).at(step));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4017"}},{"Add":{"text":" const auto &as = get<1>(get<0>(tableau).at(step));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4074"}},{"Add":{"text":" // Set current time","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4100"}},{"Add":{"text":" *const_cast(&cctkGH->cctk_time) = old_time + c * dt;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4172"}},{"Add":{"text":" // Add scaled RHS to state vector","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4212"}},{"Add":{"text":" statecomp_t::assign(var, old);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4249"}},{"Add":{"text":" for (size_t i = 0; i < as.size(); ++i)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4294"}},{"Add":{"text":" if (as.at(i) != 0)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4321"}},{"Add":{"text":" statecomp_t::axpy(var, as.at(i) * dt, ks.at(i));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4380"}},{"Add":{"text":" // TODO: Deallocate ks that are not needed any more","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4438"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_PostStep\");","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4494"}},{"Add":{"text":" if (verbose)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4513"}},{"Add":{"text":" CCTK_VINFO(\"Calculating RHS #2 at t=%g\", double(cctkGH->cctk_time));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4590"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_RHS\");","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4641"}},{"Add":{"text":" ks.push_back(rhs.copy());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4673"}},{"Add":{"text":" }","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4679"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4680"}},{"Add":{"text":" // Calculate new state vector","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4714"}},{"Add":{"text":" const auto &bs = get<1>(tableau);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4752"}},{"Add":{"text":" statecomp_t::assign(var, old);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4787"}},{"Add":{"text":" for (size_t i = 0; i < bs.size(); ++i)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4830"}},{"Add":{"text":" if (bs.at(i) != 0)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4855"}},{"Add":{"text":" statecomp_t::axpy(var, bs.at(i) * dt, ks.at(i));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4912"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4913"}},{"Add":{"text":" } else if (CCTK_EQUALS(method, \"DP87\")) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4957"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4958"}},{"Add":{"text":" typedef CCTK_REAL T;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:4983"}},{"Add":{"text":" const auto R = [](T x, T y) { return x / y; };","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5034"}},{"Add":{"text":" // These coefficients are taken from the Einstein Toolkit, thorn","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5103"}},{"Add":{"text":" // CactusNumerical/MoL, file RK87.c, written by Peter Diener,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5169"}},{"Add":{"text":" // following P. J. Prince and J. R. Dormand, Journal of","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5229"}},{"Add":{"text":" // Computational and Applied Mathematics, volume 7, no 1, 1981","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5296"}},{"Add":{"text":" const tuple >, vector > tableau{","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5353"}},{"Add":{"text":" {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5363"}},{"Add":{"text":" {/*1*/}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5422"}},{"Add":{"text":" {/*2*/ R(1, 18)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5481"}},{"Add":{"text":" {/*3*/ R(1, 48), R(1, 16)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5540"}},{"Add":{"text":" {/*4*/ R(1, 32), 0, R(3, 32)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5599"}},{"Add":{"text":" {/*5*/ R(5, 16), 0, -R(75, 64), R(75, 64)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5658"}},{"Add":{"text":" {/*6*/ R(3, 80), 0, 0, R(3, 16), R(3, 20)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5717"}},{"Add":{"text":" {/*7*/ R(29443841, 614563906), 0, 0, R(77736538, 692538347),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5790"}},{"Add":{"text":" -R(28693883, 1125000000), R(23124283, 1800000000)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5858"}},{"Add":{"text":" {/*8*/ R(16016141, 946692911), 0, 0, R(61564180, 158732637),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5931"}},{"Add":{"text":" R(22789713, 633445777), R(545815736, 2771057229),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:5994"}},{"Add":{"text":" -R(180193667, 1043307555)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6038"}},{"Add":{"text":" {/*9*/ R(39632708, 573591083), 0, 0, -R(433636366, 683701615),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6113"}},{"Add":{"text":" -R(421739975, 2616292301), R(100302831, 723423059),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6178"}},{"Add":{"text":" R(790204164, 839813087), R(800635310, 3783071287)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6246"}},{"Add":{"text":" {/*10*/ R(246121993, 1340847787), 0, 0,","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6298"}},{"Add":{"text":" -R(37695042795, 15268766246), -R(309121744, 1061227803),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6368"}},{"Add":{"text":" -R(12992083, 490766935), R(6005943493, 2108947869),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6433"}},{"Add":{"text":" R(393006217, 1396673457), R(123872331, 1001029789)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6502"}},{"Add":{"text":" {/*11*/ -R(1028468189, 846180014), 0, 0, R(8478235783, 508512852),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6581"}},{"Add":{"text":" R(1311729495, 1432422823), -R(10304129995, 1701304382),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6650"}},{"Add":{"text":" -R(48777925059, 3047939560), R(15336726248, 1032824649),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6720"}},{"Add":{"text":" -R(45442868181, 3398467696), R(3065993473, 597172653)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6792"}},{"Add":{"text":" {/*12*/ R(185892177, 718116043), 0, 0, -R(3185094517, 667107341),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6870"}},{"Add":{"text":" -R(477755414, 1098053517), -R(703635378, 230739211),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:6936"}},{"Add":{"text":" R(5731566787, 1027545527), R(5232866602, 850066563),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7002"}},{"Add":{"text":" -R(4093664535, 808688257), R(3962137247, 1805957418),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7069"}},{"Add":{"text":" R(65686358, 487910083)}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7110"}},{"Add":{"text":" {/*13*/ R(403863854, 491063109), 0, 0, -R(5068492393, 434740067),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7188"}},{"Add":{"text":" -R(411421997, 543043805), R(652783627, 914296604),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7252"}},{"Add":{"text":" R(11173962825, 925320556), -R(13158990841, 6184727034),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7321"}},{"Add":{"text":" R(3936647629, 1978049680), -R(160528059, 685178525),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7387"}},{"Add":{"text":" R(248638103, 1413531060), 0}, //","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7433"}},{"Add":{"text":" },","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7444"}},{"Add":{"text":" {R(14005451, 335480064), 0, 0, 0, 0, -R(59238493, 1068277825),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7515"}},{"Add":{"text":" R(181606767, 758867731), R(561292985, 797845732),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7574"}},{"Add":{"text":" -R(1041891430, 1371343529), R(760417239, 1151165299),","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7637"}},{"Add":{"text":" R(118820643, 751138087), -R(528747749, 2220607170), R(1, 4)}};","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7709"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7710"}},{"Add":{"text":" // Check Butcher tableau","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7739"}},{"Add":{"text":" const size_t nsteps = get<0>(tableau).size();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7789"}},{"Add":{"text":" {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7795"}},{"Add":{"text":" for (size_t step = 0; step < nsteps; ++step)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7846"}},{"Add":{"text":" // TODO: Could allow <=","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7878"}},{"Add":{"text":" assert(get<0>(tableau).at(step).size() == step);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7935"}},{"Add":{"text":" // TODO: Could allow <=","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:7965"}},{"Add":{"text":" assert(get<1>(tableau).size() == nsteps);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8013"}},{"Add":{"text":" const auto &bs = get<1>(tableau);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8053"}},{"Add":{"text":" T x = 0;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8068"}},{"Add":{"text":" for (const auto &b : bs)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8099"}},{"Add":{"text":" x += b;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8115"}},{"Add":{"text":" assert(fabs(x - 1) <= 10 * numeric_limits::epsilon());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8179"}},{"Add":{"text":" }","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8185"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8186"}},{"Add":{"text":" const auto old = var.copy();","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8219"}}]}},{"Edit":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":599,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":4969}],"down_context":[{"change":3,"pos":4969}],"flag":"BLOCK","start":8221,"end":9338,"inode":{"change":3,"pos":541}}},"contents":[{"Add":{"text":" vector ks;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8248"}},{"Add":{"text":" ks.reserve(nsteps);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8272"}},{"Add":{"text":" // Special-case step 1 (why?)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8306"}},{"Add":{"text":" ks.push_back(rhs.copy());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8336"}},{"Add":{"text":" for (size_t step = 1; step < nsteps; ++step) {","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8387"}},{"Add":{"text":" const auto &as = get<0>(tableau).at(step);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8436"}},{"Add":{"text":" T c = 0;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8451"}},{"Add":{"text":" for (const auto &a : as)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8482"}},{"Add":{"text":" c += a;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8498"}},{"Add":{"text":" // Set current time","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8524"}},{"Add":{"text":" *const_cast(&cctkGH->cctk_time) = old_time + c * dt;","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8596"}},{"Add":{"text":" // Add scaled RHS to state vector","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8636"}},{"Add":{"text":" statecomp_t::assign(var, old);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8673"}},{"Add":{"text":" for (size_t i = 0; i < as.size(); ++i)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8718"}},{"Add":{"text":" if (as.at(i) != 0)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8745"}},{"Add":{"text":" statecomp_t::axpy(var, as.at(i) * dt, ks.at(i));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8804"}},{"Add":{"text":" // TODO: Deallocate ks that are not needed any more","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8862"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_PostStep\");","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8918"}},{"Add":{"text":" if (verbose)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:8937"}},{"Add":{"text":" CCTK_VINFO(\"Calculating RHS #2 at t=%g\", double(cctkGH->cctk_time));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9014"}},{"Add":{"text":" CallScheduleGroup(cctkGH, \"ODESolvers_RHS\");","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9065"}},{"Add":{"text":" ks.push_back(rhs.copy());","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9097"}},{"Add":{"text":" }","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9103"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9104"}},{"Add":{"text":" // Calculate new state vector","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9138"}},{"Add":{"text":" const auto &bs = get<1>(tableau);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9176"}},{"Add":{"text":" statecomp_t::assign(var, old);","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9211"}},{"Add":{"text":" for (size_t i = 0; i < bs.size(); ++i)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9254"}},{"Add":{"text":" if (bs.at(i) != 0)","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9279"}},{"Add":{"text":" statecomp_t::axpy(var, bs.at(i) * dt, ks.at(i));","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9336"}},{"Add":{"text":"","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9337"}}]}},{"Replacement":{"inode":{"change":3,"pos":541},"path":"ODESolvers/src/solve.cxx","line":642,"ins":{"up_context":[{"change":3,"pos":5163}],"down_context":[{"change":3,"pos":5165}],"flag":"BLOCK","start":9339,"end":9364,"inode":{"change":3,"pos":541}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":5163},"to":{"change":3,"start":5163,"end":5165},"introduced_by":3}],"inode":{"change":3,"pos":541}},"ins_contents":[{"Add":{"text":"} // namespace ODESolver","pos":"P2MXPV6QGEX37BSXN4IKHNSS6TIA5T3TKV4CCP7TBCKLHC6L6SXAC:9363"}}],"del_contents":[{"Del":{"text":"}","pos":"P2M2LDDITHKC7HC6RMWQDTS5CT6CCH4LSIN37HQUA3NTC6SFNZOQC:5164"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review