{"header":{"message":"CarpetX: Support higher interpolation orders\n","description":null,"timestamp":"2020-02-27T21:03:02Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC","authors":[],"deps":{"hashes":[["IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC","CarpetX: Support higher interpolation orders\n"],["AA",null],["TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC","CarpetX: Add interpolator\n"],["T56YRLRZ34J5KRCRLR7MQMHA6Z3LIEG67V3YESZU5E5S3MCMFOUQC","CarpetX: Add TODO item\n"],["722HZ7UFINNE3YKSYKP2NHZ5XEG5QQLQHSKC7PREJZR3EX6RDYUAC","Import initial AMReX driver for Cactus\n"],["N3HG6SWTZQ4EVDYZYNKLDPVCSGYVGN5PK37ULA2ZV7CIRYHUUSGQC","Add conservative prolongation operators\n"],["GQVQJCNQNO2KD7ZMC7RESCUAMUAP7OED6CTA6SYLZKQGXKXZ6T3QC","Implement refluxing\n"]],"deps":["TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC"],"known":["T56YRLRZ34J5KRCRLR7MQMHA6Z3LIEG67V3YESZU5E5S3MCMFOUQC"]},"hunks":[{"Edit":{"inode":{"change":5,"pos":17918},"path":"CarpetX/param.ccl","line":201,"atom":{"NewVertex":{"up_context":[{"change":6,"pos":1}],"down_context":[{"change":7,"pos":0}],"flag":"BLOCK","start":0,"end":86,"inode":{"change":5,"pos":17918}}},"contents":[{"Add":{"text":"CCTK_INT interpolation_order \"Interpolation order\" STEERABLE=always","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:67"}},{"Add":{"text":"{","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:69"}},{"Add":{"text":" 0:* :: \"\"","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:81"}},{"Add":{"text":"} 1","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:85"}}]}},{"Edit":{"inode":{"change":5,"pos":17918},"path":"CarpetX/param.ccl","line":207,"atom":{"NewVertex":{"up_context":[{"change":6,"pos":72}],"down_context":[{"change":7,"pos":1}],"flag":"BLOCK","start":87,"end":88,"inode":{"change":5,"pos":17918}}},"contents":[{"Add":{"text":"","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:87"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":25,"atom":{"EdgeMap":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":788},"to":{"change":2,"start":788,"end":914},"introduced_by":2}],"inode":{"change":2,"pos":373}}},"contents":[{"Del":{"text":"// Calculate 10^n","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:805"}},{"Del":{"text":"constexpr inline int iexp10(int n) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:842"}},{"Del":{"text":" int r = 1;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:855"}},{"Del":{"text":" while (n > 0) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:873"}},{"Del":{"text":" r *= 10;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:886"}},{"Del":{"text":" --n;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:895"}},{"Del":{"text":" }","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:899"}},{"Del":{"text":" return r;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:911"}},{"Del":{"text":"}","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:913"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":48,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":1662}],"down_context":[{"change":2,"pos":1662}],"flag":"BLOCK","start":89,"end":376,"inode":{"change":2,"pos":373}}},"contents":[{"Add":{"text":" case 0: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:102"}},{"Add":{"text":" const T x = di[dir];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:129"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:147"}},{"Add":{"text":" assert(fabs(x) <= T(0.5));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:180"}},{"Add":{"text":"#endif","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:187"}},{"Add":{"text":" const T y0 = interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:235"}},{"Add":{"text":" switch (derivs[dir]) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:264"}},{"Add":{"text":" case 0:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:278"}},{"Add":{"text":" return y0;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:297"}},{"Add":{"text":" case 1:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:311"}},{"Add":{"text":" return 0;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:329"}},{"Add":{"text":" case 2:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:343"}},{"Add":{"text":" return 0;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:361"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:369"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:375"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":64,"ins":{"up_context":[{"change":2,"pos":1676}],"down_context":[{"change":2,"pos":1765}],"flag":"BLOCK","start":377,"end":568,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":1676},"to":{"change":2,"start":1676,"end":1765},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" const T x = di[dir];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:403"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:421"}},{"Add":{"text":" assert(x >= T(0) && x <= T(1));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:459"}},{"Add":{"text":"#endif","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:466"}},{"Add":{"text":" const T y0 = interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:514"}},{"Add":{"text":" const T y1 = interpolate(i + DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:567"}}],"del_contents":[{"Del":{"text":" T y0 = interpolate(i, di);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:1717"}},{"Del":{"text":" T y1 = interpolate(i + DI, di);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:1764"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":72,"ins":{"up_context":[{"change":2,"pos":1808}],"down_context":[{"change":2,"pos":1858}],"flag":"BLOCK","start":569,"end":607,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":1808},"to":{"change":2,"start":1808,"end":1858},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" return (1 - x) * y0 + x * y1;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:606"}}],"del_contents":[{"Del":{"text":" return (1 - di[dir]) * y0 + di[dir] * y1;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:1857"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":75,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":1909}],"down_context":[{"change":2,"pos":1909}],"flag":"BLOCK","start":608,"end":640,"inode":{"change":2,"pos":373}}},"contents":[{"Add":{"text":" case 2:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:621"}},{"Add":{"text":" return 0;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:639"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":80,"ins":{"up_context":[{"change":2,"pos":1937}],"down_context":[{"change":2,"pos":2077}],"flag":"BLOCK","start":641,"end":899,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":1937},"to":{"change":2,"start":1937,"end":2077},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" const T x = di[dir] - order / T(2);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:682"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:700"}},{"Add":{"text":" assert(fabs(x) <= T(0.5));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:733"}},{"Add":{"text":"#endif","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:740"}},{"Add":{"text":" const T y0 = interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:788"}},{"Add":{"text":" const T y1 = interpolate(i + DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:841"}},{"Add":{"text":" const T y2 = interpolate(i + 2 * DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:898"}}],"del_contents":[{"Del":{"text":" T y0 = interpolate(i, di);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:1978"}},{"Del":{"text":" T y1 = interpolate(i + DI, di);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2025"}},{"Del":{"text":" T y2 = interpolate(i + 2 * DI, di);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2076"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":89,"ins":{"up_context":[{"change":2,"pos":2120}],"down_context":[{"change":2,"pos":2276}],"flag":"BLOCK","start":900,"end":1058,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":2120},"to":{"change":2,"start":2120,"end":2276},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" return (-1 / T(2) * x + 1 / T(2) * pow(x, 2)) * y0 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:960"}},{"Add":{"text":" (1 - pow(x, 2)) * y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:998"}},{"Add":{"text":" (1 / T(2) * x + 1 / T(2) * pow(x, 2)) * y2;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1057"}}],"del_contents":[{"Del":{"text":" return (di[dir] - 1) * (di[dir] - 2) / 2 * y0 +","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2175"}},{"Del":{"text":" di[dir] * (di[dir] - 2) / -1 * y1 +","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2226"}},{"Del":{"text":" di[dir] * (di[dir] - 1) / 2 * y2;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2275"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":93,"ins":{"up_context":[{"change":2,"pos":2290}],"down_context":[{"change":2,"pos":2411}],"flag":"BLOCK","start":1059,"end":1134,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":2290},"to":{"change":2,"start":2290,"end":2411},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" return ((-1 / T(2) + x) * y0 - 2 * x * y1 + (1 / T(2) + x) * y2) /","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1133"}}],"del_contents":[{"Del":{"text":" return ((2 * di[dir] - 3) / 2 * y0 + (2 * di[dir] - 2) / -1 * y1 +","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2364"}},{"Del":{"text":" (2 * di[dir] - 1) / 2 * y2) /","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2410"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":97,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":2502}],"down_context":[{"change":2,"pos":2502}],"flag":"BLOCK","start":1135,"end":2552,"inode":{"change":2,"pos":373}}},"contents":[{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1142"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1148"}},{"Add":{"text":" case 3: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1162"}},{"Add":{"text":" const T x = di[dir] - order / T(2);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1204"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1222"}},{"Add":{"text":" assert(fabs(x) <= T(0.5));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1255"}},{"Add":{"text":"#endif","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1262"}},{"Add":{"text":" const T y0 = interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1310"}},{"Add":{"text":" const T y1 = interpolate(i + DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1363"}},{"Add":{"text":" const T y2 = interpolate(i + 2 * DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1420"}},{"Add":{"text":" const T y3 = interpolate(i + 3 * DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1477"}},{"Add":{"text":" switch (derivs[dir]) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1506"}},{"Add":{"text":" case 0:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1520"}},{"Add":{"text":" return (-1 / T(16) + 1 / T(24) * x + 1 / T(4) * pow(x, 2) -","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1588"}},{"Add":{"text":" 1 / T(6) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1628"}},{"Add":{"text":" y0 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1652"}},{"Add":{"text":" (9 / T(16) - 9 / T(8) * x - 1 / T(4) * pow(x, 2) +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1718"}},{"Add":{"text":" 1 / T(2) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1758"}},{"Add":{"text":" y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1782"}},{"Add":{"text":" (9 / T(16) + 9 / T(8) * x - 1 / T(4) * pow(x, 2) -","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1848"}},{"Add":{"text":" 1 / T(2) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1888"}},{"Add":{"text":" y2 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1912"}},{"Add":{"text":" (-1 / T(16) - 1 / T(24) * x + 1 / T(4) * pow(x, 2) +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:1980"}},{"Add":{"text":" 1 / T(6) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2020"}},{"Add":{"text":" y3;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2043"}},{"Add":{"text":" case 1:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2057"}},{"Add":{"text":" return ((1 / T(24) + 1 / T(2) * x - 1 / T(2) * pow(x, 2)) * y0 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2130"}},{"Add":{"text":" (-9 / T(8) - 1 / T(2) * x + 3 / T(2) * pow(x, 2)) * y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2203"}},{"Add":{"text":" (9 / T(8) - 1 / T(2) * x - 3 / T(2) * pow(x, 2)) * y2 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2275"}},{"Add":{"text":" (-1 / T(24) + 1 / T(2) * x + 1 / T(2) * pow(x, 2)) * y3) /","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2350"}},{"Add":{"text":" dx[dir];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2374"}},{"Add":{"text":" case 2:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2388"}},{"Add":{"text":" return ((1 / T(2) - x) * y0 + (-1 / T(2) + 3 * x) * y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2453"}},{"Add":{"text":" (-1 / T(2) - 3 * x) * y2 + (1 / T(2) + x) * y3) /","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2519"}},{"Add":{"text":" pow(dx[dir], 2);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2551"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":134,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":2516}],"down_context":[{"change":2,"pos":2516}],"flag":"BLOCK","start":2553,"end":4613,"inode":{"change":2,"pos":373}}},"contents":[{"Add":{"text":" case 4: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2566"}},{"Add":{"text":" const T x = di[dir] - order / T(2);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2608"}},{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2626"}},{"Add":{"text":" assert(fabs(x) <= T(0.5));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2659"}},{"Add":{"text":"#endif","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2666"}},{"Add":{"text":" const T y0 = interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2714"}},{"Add":{"text":" const T y1 = interpolate(i + DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2767"}},{"Add":{"text":" const T y2 = interpolate(i + 2 * DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2824"}},{"Add":{"text":" const T y3 = interpolate(i + 3 * DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2881"}},{"Add":{"text":" const T y4 = interpolate(i + 4 * DI, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2938"}},{"Add":{"text":" switch (derivs[dir]) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2967"}},{"Add":{"text":" case 0:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:2981"}},{"Add":{"text":" return (1 / T(12) * x - 1 / T(24) * pow(x, 2) - 1 / T(12) * pow(x, 3) +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3061"}},{"Add":{"text":" 1 / T(24) * pow(x, 4)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3102"}},{"Add":{"text":" y0 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3126"}},{"Add":{"text":" (-2 / T(3) * x + 2 / T(3) * pow(x, 2) + 1 / T(6) * pow(x, 3) -","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3204"}},{"Add":{"text":" 1 / T(6) * pow(x, 4)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3244"}},{"Add":{"text":" y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3268"}},{"Add":{"text":" (1 - 5 / T(4) * pow(x, 2) + 1 / T(4) * pow(x, 4)) * y2 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3340"}},{"Add":{"text":" (2 / T(3) * x + 2 / T(3) * pow(x, 2) - 1 / T(6) * pow(x, 3) -","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3417"}},{"Add":{"text":" 1 / T(6) * pow(x, 4)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3457"}},{"Add":{"text":" y3 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3481"}},{"Add":{"text":" (-1 / T(12) * x - 1 / T(24) * pow(x, 2) + 1 / T(12) * pow(x, 3) +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3562"}},{"Add":{"text":" 1 / T(24) * pow(x, 4)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3603"}},{"Add":{"text":" y4;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3626"}},{"Add":{"text":" case 1:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3640"}},{"Add":{"text":" return ((1 / T(12) - 1 / T(12) * x - 1 / T(4) * pow(x, 2) +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3708"}},{"Add":{"text":" 1 / T(6) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3749"}},{"Add":{"text":" y0 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3774"}},{"Add":{"text":" (-2 / T(3) + 4 / T(3) * x + 1 / T(2) * pow(x, 2) -","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3841"}},{"Add":{"text":" 2 / T(3) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3882"}},{"Add":{"text":" y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3907"}},{"Add":{"text":" (-5 / T(2) * x + pow(x, 3)) * y2 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:3958"}},{"Add":{"text":" (2 / T(3) + 4 / T(3) * x - 1 / T(2) * pow(x, 2) -","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4024"}},{"Add":{"text":" 2 / T(3) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4065"}},{"Add":{"text":" y3 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4090"}},{"Add":{"text":" (-1 / T(12) - 1 / T(12) * x + 1 / T(4) * pow(x, 2) +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4159"}},{"Add":{"text":" 1 / T(6) * pow(x, 3)) *","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4200"}},{"Add":{"text":" y4) /","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4226"}},{"Add":{"text":" dx[dir];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4250"}},{"Add":{"text":" case 2:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4264"}},{"Add":{"text":" return ((-1 / T(12) - 1 / T(2) * x + 1 / T(2) * pow(x, 2)) * y0 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4338"}},{"Add":{"text":" (4 / T(3) + x - 2 * pow(x, 2)) * y1 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4392"}},{"Add":{"text":" (-5 / T(2) + 3 * pow(x, 2)) * y2 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4443"}},{"Add":{"text":" (4 / T(3) - x - 2 * pow(x, 2)) * y3 +","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4497"}},{"Add":{"text":" (-1 / T(12) + 1 / T(2) * x + 1 / T(2) * pow(x, 2)) * y4) /","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4572"}},{"Add":{"text":" pow(dx[dir], 2);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4604"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4612"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":183,"ins":{"up_context":[{"change":2,"pos":2522}],"down_context":[{"change":2,"pos":2537}],"flag":"BLOCK","start":4614,"end":4650,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":2522},"to":{"change":2,"start":2522,"end":2537},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" default:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4626"}},{"Add":{"text":" assert(0);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4643"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4649"}}],"del_contents":[{"Del":{"text":" assert(0);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:2536"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":200,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":3208}],"down_context":[{"change":2,"pos":3208}],"flag":"BLOCK","start":4651,"end":4678,"inode":{"change":2,"pos":373}}},"contents":[{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4677"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":204,"atom":{"EdgeMap":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":3315},"to":{"change":2,"start":3315,"end":3343},"introduced_by":2}],"inode":{"change":2,"pos":373}}},"contents":[{"Del":{"text":" constexpr int order = 2;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:3341"}},{"Del":{"text":"","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:3342"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":273,"atom":{"EdgeMap":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":5730},"to":{"change":2,"start":5730,"end":5768},"introduced_by":2}],"inode":{"change":2,"pos":373}}},"contents":[{"Del":{"text":" const int op = operations[v];","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:5767"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":274,"ins":{"up_context":[{"change":2,"pos":5799}],"down_context":[{"change":2,"pos":5882}],"flag":"BLOCK","start":4679,"end":4913,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":5799},"to":{"change":2,"start":5799,"end":5882},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" int op = operations[v];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4710"}},{"Add":{"text":" while (op > 0) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4735"}},{"Add":{"text":" const int dir = op % 10 - 1;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4774"}},{"Add":{"text":" if (dir >= 0) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4800"}},{"Add":{"text":" assert(dir >= 0 && dir < dim);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4843"}},{"Add":{"text":" ++derivs[dir];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4870"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4882"}},{"Add":{"text":" op /= 10;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4902"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4912"}}],"del_contents":[{"Del":{"text":" for (int d = 0; d < dim; ++d)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:5836"}},{"Del":{"text":" derivs[d] = (op / iexp10(d)) % 10;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:5881"}}]}},{"Edit":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":286,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":5957}],"down_context":[{"change":2,"pos":5957}],"flag":"BLOCK","start":4914,"end":5615,"inode":{"change":2,"pos":373}}},"contents":[{"Add":{"text":" switch (interpolation_order) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4952"}},{"Add":{"text":" case 0: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:4970"}},{"Add":{"text":" constexpr int order = 0;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5005"}},{"Add":{"text":"#pragma omp simd","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5022"}},{"Add":{"text":" for (int n = 0; n < np; ++n) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5063"}},{"Add":{"text":" vect i;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5093"}},{"Add":{"text":" vect di;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5130"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5174"}},{"Add":{"text":" CCTK_REAL x = particles[n].pos(d);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5223"}},{"Add":{"text":" CCTK_REAL ri = (x - x0[d]) / dx[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5273"}},{"Add":{"text":" i[d] = lrint(ri - (order / CCTK_REAL(2)));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5330"}},{"Add":{"text":" di[d] = ri - i[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5363"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5377"}},{"Add":{"text":" const interpolator interp{vars, vi, derivs, dx};","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5456"}},{"Add":{"text":" varresult.at(n) = interp.interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5522"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5534"}},{"Add":{"text":" break;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5551"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5561"}},{"Add":{"text":" case 1: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5579"}},{"Add":{"text":" constexpr int order = 1;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5614"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":307,"ins":{"up_context":[{"change":2,"pos":5974}],"down_context":[{"change":2,"pos":6299}],"flag":"BLOCK","start":5616,"end":7334,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":5974},"to":{"change":2,"start":5974,"end":6299},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" for (int n = 0; n < np; ++n) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5656"}},{"Add":{"text":" vect i;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5686"}},{"Add":{"text":" vect di;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5723"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5767"}},{"Add":{"text":" CCTK_REAL x = particles[n].pos(d);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5816"}},{"Add":{"text":" CCTK_REAL ri = (x - x0[d]) / dx[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5866"}},{"Add":{"text":" i[d] = lrint(ri - (order / CCTK_REAL(2)));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5923"}},{"Add":{"text":" di[d] = ri - i[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5956"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:5970"}},{"Add":{"text":" const interpolator interp{vars, vi, derivs, dx};","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6049"}},{"Add":{"text":" varresult.at(n) = interp.interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6115"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6127"}},{"Add":{"text":" break;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6144"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6154"}},{"Add":{"text":" case 2: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6172"}},{"Add":{"text":" constexpr int order = 2;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6207"}},{"Add":{"text":"#pragma omp simd","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6224"}},{"Add":{"text":" for (int n = 0; n < np; ++n) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6265"}},{"Add":{"text":" vect i;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6295"}},{"Add":{"text":" vect di;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6332"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6376"}},{"Add":{"text":" CCTK_REAL x = particles[n].pos(d);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6425"}},{"Add":{"text":" CCTK_REAL ri = (x - x0[d]) / dx[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6475"}},{"Add":{"text":" i[d] = lrint(ri - (order / CCTK_REAL(2)));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6532"}},{"Add":{"text":" di[d] = ri - i[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6565"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6579"}},{"Add":{"text":" const interpolator interp{vars, vi, derivs, dx};","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6658"}},{"Add":{"text":" varresult.at(n) = interp.interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6724"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6736"}},{"Add":{"text":" break;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6753"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6763"}},{"Add":{"text":" case 3: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6781"}},{"Add":{"text":" constexpr int order = 3;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6816"}},{"Add":{"text":"#pragma omp simd","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6833"}},{"Add":{"text":" for (int n = 0; n < np; ++n) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6874"}},{"Add":{"text":" vect i;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6904"}},{"Add":{"text":" vect di;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6941"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:6985"}},{"Add":{"text":" CCTK_REAL x = particles[n].pos(d);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7034"}},{"Add":{"text":" CCTK_REAL ri = (x - x0[d]) / dx[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7084"}},{"Add":{"text":" i[d] = lrint(ri - (order / CCTK_REAL(2)));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7141"}},{"Add":{"text":" di[d] = ri - i[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7174"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7188"}},{"Add":{"text":" const interpolator interp{vars, vi, derivs, dx};","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7267"}},{"Add":{"text":" varresult.at(n) = interp.interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7333"}}],"del_contents":[{"Del":{"text":" for (int n = 0; n < np; ++n) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6012"}},{"Del":{"text":" vect i;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6040"}},{"Del":{"text":" vect di;","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6075"}},{"Del":{"text":" for (int d = 0; d < dim; ++d) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6117"}},{"Del":{"text":" CCTK_REAL x = particles[n].pos(d);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6164"}},{"Del":{"text":" CCTK_REAL ri = (x - x0[d]) / dx[d];","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6212"}},{"Del":{"text":" i[d] = lrint(ri - (order / CCTK_REAL(2)));","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6267"}},{"Del":{"text":" di[d] = ri - i[d];","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6298"}}]}},{"Replacement":{"inode":{"change":2,"pos":373},"path":"CarpetX/src/interpolate.cxx","line":353,"ins":{"up_context":[{"change":2,"pos":6311}],"down_context":[{"change":2,"pos":7751}],"flag":"BLOCK","start":7335,"end":8009,"inode":{"change":2,"pos":373}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":6311},"to":{"change":2,"start":6311,"end":7751},"introduced_by":2}],"inode":{"change":2,"pos":373}},"ins_contents":[{"Add":{"text":" break;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7351"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7361"}},{"Add":{"text":" case 4: {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7379"}},{"Add":{"text":" constexpr int order = 4;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7414"}},{"Add":{"text":"#pragma omp simd","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7431"}},{"Add":{"text":" for (int n = 0; n < np; ++n) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7472"}},{"Add":{"text":" vect i;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7502"}},{"Add":{"text":" vect di;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7539"}},{"Add":{"text":" for (int d = 0; d < dim; ++d) {","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7583"}},{"Add":{"text":" CCTK_REAL x = particles[n].pos(d);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7632"}},{"Add":{"text":" CCTK_REAL ri = (x - x0[d]) / dx[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7682"}},{"Add":{"text":" i[d] = lrint(ri - (order / CCTK_REAL(2)));","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7739"}},{"Add":{"text":" di[d] = ri - i[d];","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7772"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7786"}},{"Add":{"text":" const interpolator interp{vars, vi, derivs, dx};","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7865"}},{"Add":{"text":" varresult.at(n) = interp.interpolate(i, di);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7931"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7943"}},{"Add":{"text":" break;","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7960"}},{"Add":{"text":" }","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7970"}},{"Add":{"text":" default:","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:7987"}},{"Add":{"text":" assert(0);","pos":"IC2NW4EHRU42VTEF6RNOC26DVB3BKSA6JLGMO6CUXGIHHOZ3KTCQC:8008"}}],"del_contents":[{"Del":{"text":"","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6311"}},{"Del":{"text":" // const vect DI{1, 0, 0};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6358"}},{"Del":{"text":" // const vect DJ{0, 1, 0};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6405"}},{"Del":{"text":" // const vect DK{0, 0, 1};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6452"}},{"Del":{"text":" // const auto interp0{[&](const vect i) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6514"}},{"Del":{"text":" // return vars(i[0], i[1], i[2], vi);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6564"}},{"Del":{"text":" // }};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6581"}},{"Del":{"text":" // const auto interp1{[&](const vect i) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6643"}},{"Del":{"text":" // if (derivs[0] == 0)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6678"}},{"Del":{"text":" // return (1 - di[0]) * interp0(i) + di[0] * interp0(i + DI);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6754"}},{"Del":{"text":" // if (derivs[0] == 1)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6789"}},{"Del":{"text":" // return (-interp0(i) + interp0(i + DI)) / dx[0];","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6854"}},{"Del":{"text":" // return 0 * interp0(i);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6892"}},{"Del":{"text":" // }};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6909"}},{"Del":{"text":" // const auto interp2{[&](const vect i) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:6971"}},{"Del":{"text":" // if (derivs[1] == 0)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7006"}},{"Del":{"text":" // return (1 - di[1]) * interp1(i) + di[1] * interp1(i + DJ);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7082"}},{"Del":{"text":" // if (derivs[1] == 1)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7117"}},{"Del":{"text":" // return (-interp1(i) + interp1(i + DJ)) / dx[1];","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7182"}},{"Del":{"text":" // return 0 * interp1(i);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7220"}},{"Del":{"text":" // }};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7237"}},{"Del":{"text":" // const auto interp3{[&](const vect i) {","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7299"}},{"Del":{"text":" // if (derivs[2] == 0)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7334"}},{"Del":{"text":" // return (1 - di[2]) * interp2(i) + di[2] * interp2(i + DK);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7410"}},{"Del":{"text":" // if (derivs[2] == 1)","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7445"}},{"Del":{"text":" // return (-interp2(i) + interp2(i + DK)) / dx[2];","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7510"}},{"Del":{"text":" // return 0 * interp2(i);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7548"}},{"Del":{"text":" // }};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7565"}},{"Del":{"text":" // varresult.at(n) = interp3(i);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7608"}},{"Del":{"text":"","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7609"}},{"Del":{"text":" const interpolator interp{vars, vi, derivs, dx};","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7686"}},{"Del":{"text":" varresult.at(n) = interp.interpolate(i, di);","pos":"TUDUMVD5MTJJJGVYVOCALOOKOUYNB7LOFJRCYCZDVGN536JAO2OQC:7750"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review