{"header":{"message":"CarpetX: Correct tests for conservative interpolation\n","description":null,"timestamp":"2020-01-04T21:25:22Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC","authors":[],"deps":{"hashes":[["C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC","CarpetX: Correct tests for conservative interpolation\n"],["AA",null],["BRC6KOCKJIARRSZI3TJOKIK6XBSNYAS2IP6LJOPS2FPMLN57QKVQC","CarpetX: Use CCTK_isnan instead of isnan in some places\n"],["DGEGTDYVGRHEGNFDHCCVTD6AVLRJ2XHDXQPWR5MMR62V5NU2XGIQC","Add generic prolongation operators\n"],["DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC","Begin to correct conservative interpolation operators\n"]],"deps":["BRC6KOCKJIARRSZI3TJOKIK6XBSNYAS2IP6LJOPS2FPMLN57QKVQC"],"known":["DGEGTDYVGRHEGNFDHCCVTD6AVLRJ2XHDXQPWR5MMR62V5NU2XGIQC","DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC"]},"hunks":[{"Replacement":{"inode":{"change":3,"pos":1186},"path":"CarpetX/src/prolongate_3d_rf2.cxx","line":312,"ins":{"up_context":[{"change":3,"pos":2495}],"down_context":[{"change":3,"pos":2688}],"flag":"BLOCK","start":0,"end":229,"inode":{"change":3,"pos":1186}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":2495},"to":{"change":3,"start":2495,"end":2688},"introduced_by":3}],"inode":{"change":3,"pos":1186}},"ins_contents":[{"Add":{"text":" for (int order = 0; order <= ORDER; ++order) {","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:50"}},{"Add":{"text":" // const auto f{[&](T x) {","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:83"}},{"Add":{"text":" // return (order + 1) * (order == 0 ? T(1) : pow(x, order));","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:152"}},{"Add":{"text":" // }};","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:165"}},{"Add":{"text":" const auto fint{[&](T x) { return pow(x, order + 1); }};","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:228"}}],"del_contents":[{"Del":{"text":" for (int order = 0; order < ORDER; ++order) {","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:2544"}},{"Del":{"text":" auto f = [&](T x) { return order == 0 ? T(1) : pow(x, order); };","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:2615"}},{"Del":{"text":" auto fint = [&](T x) { return pow(x, order + 1) / (order + 1); };","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:2687"}}]}},{"Replacement":{"inode":{"change":3,"pos":1186},"path":"CarpetX/src/prolongate_3d_rf2.cxx","line":323,"ins":{"up_context":[{"change":3,"pos":2900}],"down_context":[{"change":3,"pos":2967}],"flag":"BLOCK","start":230,"end":449,"inode":{"change":3,"pos":1186}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":2900},"to":{"change":3,"start":2900,"end":2967},"introduced_by":3}],"inode":{"change":3,"pos":1186}},"ins_contents":[{"Add":{"text":" const T x = (i - i0) + CENTERING / T(2);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:280"}},{"Add":{"text":" // T y = f(x);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:305"}},{"Add":{"text":" const T dx = 1;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:331"}},{"Add":{"text":" const T xlo = x - dx / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:367"}},{"Add":{"text":" const T xhi = x + dx / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:403"}},{"Add":{"text":" const T y = fint(xhi) - fint(xlo);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:448"}}],"del_contents":[{"Del":{"text":" T x = (i - i0) + CENTERING / T(2);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:2944"}},{"Del":{"text":" T y = f(x);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:2966"}}]}},{"Replacement":{"inode":{"change":3,"pos":1186},"path":"CarpetX/src/prolongate_3d_rf2.cxx","line":339,"ins":{"up_context":[{"change":3,"pos":3320}],"down_context":[{"change":3,"pos":3454}],"flag":"BLOCK","start":450,"end":608,"inode":{"change":3,"pos":1186}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":3320},"to":{"change":3,"start":3320,"end":3454},"introduced_by":3}],"inode":{"change":3,"pos":1186}},"ins_contents":[{"Add":{"text":" const T dx = x1[1] - x1[0];","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:485"}},{"Add":{"text":" const T xlo = x1[0] - dx / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:523"}},{"Add":{"text":" const T xhi = x1[1] + dx / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:561"}},{"Add":{"text":" const T yint = fint(xhi) - fint(xlo);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:607"}}],"del_contents":[{"Del":{"text":" T dx = x1[1] - x1[0];","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3349"}},{"Del":{"text":" T xlo = x1[0] - dx / 2;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3381"}},{"Del":{"text":" T xhi = x1[1] + dx / 2;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3413"}},{"Del":{"text":" T yint = fint(xhi) - fint(xlo);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3453"}}]}},{"Replacement":{"inode":{"change":3,"pos":1186},"path":"CarpetX/src/prolongate_3d_rf2.cxx","line":345,"ins":{"up_context":[{"change":3,"pos":3518}],"down_context":[{"change":3,"pos":5008}],"flag":"BLOCK","start":609,"end":2600,"inode":{"change":3,"pos":1186}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":3518},"to":{"change":3,"start":3518,"end":4210},"introduced_by":3},{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":4210},"to":{"change":2,"start":133,"end":177},"introduced_by":2},{"previous":"","flag":"DELETED","from":{"change":2,"pos":177},"to":{"change":3,"start":4249,"end":5008},"introduced_by":2},{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":4249},"to":{"change":3,"start":4249,"end":5008},"introduced_by":3}],"inode":{"change":3,"pos":1186}},"ins_contents":[{"Add":{"text":" // Don't test this, the case (VC,CONS) should not be used","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:674"}},{"Add":{"text":" if (false) {","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:695"}},{"Add":{"text":" array xs, ys;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:729"}},{"Add":{"text":" xs[0] = xs[n + 2] = 0 / T(0);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:769"}},{"Add":{"text":" ys[0] = ys[n + 2] = 0 / T(0);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:809"}},{"Add":{"text":" const int i0 = n / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:841"}},{"Add":{"text":" for (int i = -1; i < n; ++i) {","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:882"}},{"Add":{"text":" const T x = (i - i0) + CENTERING / T(2);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:935"}},{"Add":{"text":" // T y = f(x);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:962"}},{"Add":{"text":" const T dx = 1;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:990"}},{"Add":{"text":" const T xlo = x - dx / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1028"}},{"Add":{"text":" const T xhi = x + dx / 2;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1066"}},{"Add":{"text":" const T y = fint(xhi) - fint(xlo);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1113"}},{"Add":{"text":" xs[i + 2] = x;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1140"}},{"Add":{"text":" ys[i + 2] = y;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1167"}},{"Add":{"text":" }","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1179"}},{"Add":{"text":" array x1;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1205"}},{"Add":{"text":" array y1;","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1231"}},{"Add":{"text":" for (int off = -1; off < 2; ++off) {","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1278"}},{"Add":{"text":" x1[off + 1] = CENTERING / T(4) + off / T(2);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1335"}},{"Add":{"text":" if (off < 0)","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1360"}},{"Add":{"text":" y1[off + 1] =","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1388"}},{"Add":{"text":" interp1d()(&ys[i0 + 1], 1, off + 2);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1467"}},{"Add":{"text":" else","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1484"}},{"Add":{"text":" y1[off + 1] =","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1512"}},{"Add":{"text":" interp1d()(&ys[i0 + 2], 1, off);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1587"}},{"Add":{"text":" assert(!CCTK_isnan(y1[off + 1]));","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1633"}},{"Add":{"text":" }","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1645"}},{"Add":{"text":" const T dx = x1[1] - x1[0];","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1683"}},{"Add":{"text":" const T xlo = x1[0];","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1714"}},{"Add":{"text":" const T xhi = x1[2];","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1745"}},{"Add":{"text":" const T yint = fint(xhi) - fint(xlo);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1793"}},{"Add":{"text":" if (!(y1[0] / 4 + y1[1] / 2 + y1[2] / 4 == ys[i0 + 2]) ||","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1861"}},{"Add":{"text":" !(y1[0] * dx / 2 + y1[1] * dx + y1[2] * dx / 2 == yint)) {","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:1934"}},{"Add":{"text":" cerr << \"settings: CENTERING=\" << CENTERING << \" ORDER=\" << ORDER","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2012"}},{"Add":{"text":" << \" order=\" << order << \"\\n\";","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2060"}},{"Add":{"text":" cerr << \"input:\\n\";","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2092"}},{"Add":{"text":" for (int i = -1; i < n; ++i)","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2133"}},{"Add":{"text":" cerr << \" xs=\" << xs[i + 2] << \" ys=\" << ys[i + 2] << \"\\n\";","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2208"}},{"Add":{"text":" cerr << \"output:\\n\";","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2241"}},{"Add":{"text":" for (int off = -1; off < 2; ++off)","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2288"}},{"Add":{"text":" cerr << \" x1=\" << x1[off + 1] << \" y1=\" << y1[off + 1] << \"\\n\";","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2367"}},{"Add":{"text":" cerr << \"xlo=\" << xlo << \" xhi=\" << xhi << \" yint=\" << yint << \"\\n\";","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2448"}},{"Add":{"text":" }","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2460"}},{"Add":{"text":" assert(y1[0] / 4 + y1[1] / 2 + y1[2] / 4 == ys[i0 + 2]);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2527"}},{"Add":{"text":" assert(y1[0] * dx / 2 + y1[1] * dx + y1[2] * dx / 2 == yint);","pos":"C7Y67X3DJBDM4QW3NAVSVJF6GSSLD466S4FFJSIKUEJANLPC6FIAC:2599"}}],"del_contents":[{"Del":{"text":" array xs, ys;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3549"}},{"Del":{"text":" xs[0] = xs[n + 2] = 0 / T(0);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3587"}},{"Del":{"text":" ys[0] = ys[n + 2] = 0 / T(0);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3625"}},{"Del":{"text":" const int i0 = n / 2;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3655"}},{"Del":{"text":" for (int i = -1; i < n; ++i) {","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3694"}},{"Del":{"text":" T x = (i - i0) + CENTERING / T(2);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3739"}},{"Del":{"text":" T y = f(x);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3761"}},{"Del":{"text":" xs[i + 2] = x;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3786"}},{"Del":{"text":" ys[i + 2] = y;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3811"}},{"Del":{"text":" }","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3821"}},{"Del":{"text":" array x1;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3845"}},{"Del":{"text":" array y1;","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3869"}},{"Del":{"text":" for (int off = -1; off < 2; ++off) {","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3914"}},{"Del":{"text":" x1[off + 1] = CENTERING / T(4) + off / T(2);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3969"}},{"Del":{"text":" if (off < 0)","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:3992"}},{"Del":{"text":" y1[off + 1] =","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4018"}},{"Del":{"text":" interp1d()(&ys[i0 + 1], 1, off + 2);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4095"}},{"Del":{"text":" else","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4110"}},{"Del":{"text":" y1[off + 1] =","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4136"}},{"Del":{"text":" interp1d()(&ys[i0 + 2], 1, off);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4209"}},{"Del":{"text":" assert(!CCTK_isnan(y1[off + 1]));","pos":"BRC6KOCKJIARRSZI3TJOKIK6XBSNYAS2IP6LJOPS2FPMLN57QKVQC:176"}},{"Del":{"text":" }","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4258"}},{"Del":{"text":" T dx = x1[1] - x1[0];","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4288"}},{"Del":{"text":" T xlo = x1[0];","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4311"}},{"Del":{"text":" T xhi = x1[2];","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4334"}},{"Del":{"text":" T yint = fint(xhi) - fint(xlo);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4374"}},{"Del":{"text":" if (!(y1[0] / 4 + y1[1] / 2 + y1[2] / 4 == ys[i0 + 2]) ||","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4440"}},{"Del":{"text":" !(y1[0] * dx / 2 + y1[1] * dx + y1[2] * dx / 2 == yint)) {","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4511"}},{"Del":{"text":" cerr << \"settings: CENTERING=\" << CENTERING << \" ORDER=\" << ORDER","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4587"}},{"Del":{"text":" << \" order=\" << order << \"\\n\";","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4633"}},{"Del":{"text":" cerr << \"input:\\n\";","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4663"}},{"Del":{"text":" for (int i = -1; i < n; ++i)","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4702"}},{"Del":{"text":" cerr << \" xs=\" << xs[i + 2] << \" ys=\" << ys[i + 2] << \"\\n\";","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4775"}},{"Del":{"text":" cerr << \"output:\\n\";","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4806"}},{"Del":{"text":" for (int off = -1; off < 2; ++off)","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4851"}},{"Del":{"text":" cerr << \" x1=\" << x1[off + 1] << \" y1=\" << y1[off + 1] << \"\\n\";","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:4928"}},{"Del":{"text":" cerr << \"xlo=\" << xlo << \" xhi=\" << xhi << \" yint=\" << yint << \"\\n\";","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:5007"}}]}},{"Edit":{"inode":{"change":3,"pos":1186},"path":"CarpetX/src/prolongate_3d_rf2.cxx","line":392,"atom":{"EdgeMap":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":5018},"to":{"change":3,"start":5018,"end":5153},"introduced_by":3}],"inode":{"change":3,"pos":1186}}},"contents":[{"Del":{"text":" assert(y1[0] / 4 + y1[1] / 2 + y1[2] / 4 == ys[i0 + 2]);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:5082"}},{"Del":{"text":" assert(y1[0] * dx / 2 + y1[1] * dx + y1[2] * dx / 2 == yint);","pos":"DQZRKU4B6C3WWWXFMBBSVKYD5AXDGUMHHPACVWKERGKB5WLHZGYAC:5152"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review