{"header":{"message":"AHFinder: Correct area calculation\n","description":null,"timestamp":"2020-04-05T19:25:18Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC","authors":[],"deps":{"hashes":[["K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC","AHFinder: Correct area calculation\n"],["AA",null],["UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC","AHFinder: Stronger tests\n"],["FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC","AHFinder: Works\n"],["AQLAEDIO5N7225NC2CE6YTGM2CW5JNHDNGMXCMQWJWRFE7DRS3LAC","AHFinder: Can now evaluate expansion\n"]],"deps":["UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC"],"known":["AQLAEDIO5N7225NC2CE6YTGM2CW5JNHDNGMXCMQWJWRFE7DRS3LAC","FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC"]},"hunks":[{"Replacement":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":69,"ins":{"up_context":[{"change":2,"pos":599}],"down_context":[{"change":2,"pos":711}],"flag":"BLOCK","start":0,"end":124,"inode":{"change":3,"pos":2243}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":599},"to":{"change":2,"start":599,"end":711},"introduced_by":2}],"inode":{"change":3,"pos":2243}},"ins_contents":[{"Add":{"text":"metric_t brill_lindquist_metric(const cGH *const cctkGH,","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:59"}},{"Add":{"text":" const coords_t &coords) {","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:123"}}],"del_contents":[{"Del":{"text":"metric_t calculate_metric(const cGH *const cctkGH,","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:652"}},{"Del":{"text":" const coords_t &coords) {","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:710"}}]}},{"Edit":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":74,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":776}],"down_context":[{"change":2,"pos":776}],"flag":"BLOCK","start":125,"end":205,"inode":{"change":3,"pos":2243}}},"contents":[{"Add":{"text":" const T x0 = 0.0;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:144"}},{"Add":{"text":" const T y0 = 0.0;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:164"}},{"Add":{"text":" const T z0 = 0.0;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:184"}},{"Add":{"text":" const T M = 1.0;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:203"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:204"}}]}},{"Edit":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":81,"atom":{"EdgeMap":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":836},"to":{"change":2,"start":836,"end":923},"introduced_by":2}],"inode":{"change":3,"pos":2243}}},"contents":[{"Del":{"text":" const T x0 = 0.0;","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:857"}},{"Del":{"text":" const T y0 = 0.0;","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:879"}},{"Del":{"text":" const T z0 = 0.0;","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:901"}},{"Del":{"text":" const T M = 1.0;","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:922"}}]}},{"Replacement":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":82,"ins":{"up_context":[{"change":2,"pos":949}],"down_context":[{"change":2,"pos":976}],"flag":"BLOCK","start":206,"end":266,"inode":{"change":3,"pos":2243}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":2,"pos":949},"to":{"change":2,"start":949,"end":976},"introduced_by":2}],"inode":{"change":3,"pos":2243}},"ins_contents":[{"Add":{"text":" using DT = dual; // spatial derivatives (x, y, z)","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:265"}}],"del_contents":[{"Del":{"text":" using DT = dual;","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:975"}}]}},{"Edit":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":397,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":11931}],"down_context":[{"change":3,"pos":17312}],"flag":"BLOCK","start":267,"end":285,"inode":{"change":3,"pos":2243}}},"contents":[{"Add":{"text":"#ifdef CCTK_DEBUG","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:284"}}]}},{"Edit":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":399,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":17373}],"down_context":[{"change":3,"pos":17373}],"flag":"BLOCK","start":286,"end":333,"inode":{"change":3,"pos":2243}}},"contents":[{"Add":{"text":" assert(fabs(s2 - 1) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:324"}},{"Add":{"text":"#endif","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:331"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:332"}}]}},{"Replacement":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":403,"ins":{"up_context":[{"change":3,"pos":17391}],"down_context":[{"change":3,"pos":17430}],"flag":"BLOCK","start":334,"end":3907,"inode":{"change":3,"pos":2243}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":17391},"to":{"change":3,"start":17391,"end":17430},"introduced_by":3}],"inode":{"change":3,"pos":2243}},"ins_contents":[{"Add":{"text":" // m1_a","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:347"}},{"Add":{"text":" // TODO: rename to m1^a","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:377"}},{"Add":{"text":" dual m1ur = h_t;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:403"}},{"Add":{"text":" dual m1ut = 1;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:427"}},{"Add":{"text":" dual z_m1up = 0;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:453"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:454"}},{"Add":{"text":" dual m1r = qrr * m1ur + qrt * m1ut + s_qrp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:516"}},{"Add":{"text":" dual m1t = qrt * m1ur + qtt * m1ut + s_qtp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:578"}},{"Add":{"text":" dual s_m1p = s_qrp * m1ur + s_qtp * m1ut + ss_qpp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:647"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:648"}},{"Add":{"text":" dual m12 = m1r * m1ur + m1t * m1ut + s_m1p * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:710"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:711"}},{"Add":{"text":" m1ur /= sqrt(m12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:736"}},{"Add":{"text":" m1ut /= sqrt(m12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:761"}},{"Add":{"text":" z_m1up /= sqrt(m12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:788"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:789"}},{"Add":{"text":" m1r = qrr * m1ur + qrt * m1ut + s_qrp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:843"}},{"Add":{"text":" m1t = qrt * m1ur + qtt * m1ut + s_qtp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:897"}},{"Add":{"text":" s_m1p = s_qrp * m1ur + s_qtp * m1ut + ss_qpp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:958"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:959"}},{"Add":{"text":" m12 = m1r * m1ur + m1t * m1ut + s_m1p * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1013"}},{"Add":{"text":" assert(fabs(m12 - 1) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1053"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1054"}},{"Add":{"text":" const dual sm1 = sr * m1ur + st * m1ut + s_sp * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1119"}},{"Add":{"text":" assert(fabs(sm1) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1155"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1156"}},{"Add":{"text":" // m2_a","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1170"}},{"Add":{"text":" dual s_m2ur = s_h_p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1200"}},{"Add":{"text":" dual s_m2ut = 0;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1226"}},{"Add":{"text":" dual m2up = 1;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1250"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1251"}},{"Add":{"text":" const dual s_m1m2 = s_m2ur * m1r + s_m2ut * m1t + m2up * s_m1p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1324"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1325"}},{"Add":{"text":" s_m2ur -= s_m1m2 * m1ur;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1356"}},{"Add":{"text":" s_m2ut -= s_m1m2 * m1ut;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1387"}},{"Add":{"text":" m2up -= s_m1m2 * z_m1up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1418"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1419"}},{"Add":{"text":" const dual s_m2r = qrr * s_m2ur + qrt * s_m2ut + s_qrp * m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1491"}},{"Add":{"text":" const dual s_m2t = qrt * s_m2ur + qtt * s_m2ut + s_qtp * m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1563"}},{"Add":{"text":" const dual ss_m2p = s_qrp * s_m2ur + s_qtp * s_m2ut + ss_qpp * m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1641"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1642"}},{"Add":{"text":" const dual ss_m22 = s_m2r * s_m2ur + s_m2t * s_m2ut + ss_m2p * m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1720"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1721"}},{"Add":{"text":" const dual m2ur = s_m2ur / sqrt(ss_m22);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1771"}},{"Add":{"text":" const dual m2ut = s_m2ut / sqrt(ss_m22);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1821"}},{"Add":{"text":" const dual z_m2up = m2up / sqrt(ss_m22);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1871"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1872"}},{"Add":{"text":" const dual m2r = qrr * m2ur + qrt * m2ut + s_qrp * z_m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:1940"}},{"Add":{"text":" const dual m2t = qrt * m2ur + qtt * m2ut + s_qtp * z_m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2008"}},{"Add":{"text":" const dual s_m2p = s_qrp * m2ur + s_qtp * m2ut + ss_qpp * z_m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2083"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2084"}},{"Add":{"text":" const dual m1m2 = m2ur * m1r + m2ut * m1t + z_m2up * s_m1p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2153"}},{"Add":{"text":" assert(fabs(m1m2) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2190"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2191"}},{"Add":{"text":" const dual m22 = m2r * m2ur + m2t * m2ut + s_m2p * z_m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2259"}},{"Add":{"text":" assert(fabs(m22 - 1) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2299"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2300"}},{"Add":{"text":" const dual sm2 = sr * m2ur + st * m2ut + s_sp * z_m2up;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2365"}},{"Add":{"text":" assert(fabs(sm2) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2401"}},{"Add":{"text":"#endif","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2408"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2409"}},{"Add":{"text":"#if 0","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2415"}},{"Add":{"text":" {","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2423"}},{"Add":{"text":" // m1_a","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2439"}},{"Add":{"text":" const dual m1ux = x_t + h_t * sin(theta) * cos(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2503"}},{"Add":{"text":" const dual m1uy = y_t + h_t * sin(theta) * sin(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2567"}},{"Add":{"text":" const dual m1uz = z_t + h_t * cos(theta);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2620"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2621"}},{"Add":{"text":" const dual m1x = gxx * m1ux + gxy * m1uy + gxz * m1uz;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2687"}},{"Add":{"text":" const dual m1y = gxy * m1ux + gyy * m1uy + gyz * m1uz;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2753"}},{"Add":{"text":" const dual m1z = gxz * m1ux + gyz * m1uy + gzz * m1uz;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2819"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2820"}},{"Add":{"text":" const dual m1r = m1x * x_r + m1y * y_r + m1z * z_r;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2883"}},{"Add":{"text":" const dual m1t = m1x * x_t + m1y * y_t + m1z * z_t;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:2946"}},{"Add":{"text":" const dual s_m1p = m1x * s_x_p + m1y * s_y_p + m1z * s_z_p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3017"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3018"}},{"Add":{"text":" const dual sm1 = sur * m1r + sut * m1t + z_sup * s_m1p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3085"}},{"Add":{"text":" assert(fabs(sm1) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3123"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3124"}},{"Add":{"text":" // m2_a","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3140"}},{"Add":{"text":" const dual s_m2ux = s_x_p + s_h_p * sin(theta) * cos(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3210"}},{"Add":{"text":" const dual s_m2uy = s_y_p + s_h_p * sin(theta) * sin(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3280"}},{"Add":{"text":" const dual s_m2uz = s_z_p + s_h_p * cos(theta);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3339"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3340"}},{"Add":{"text":" const dual s_m2x = gxx * s_m2ux + gxy * s_m2uy + gxz * s_m2uz;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3414"}},{"Add":{"text":" const dual s_m2y = gxy * s_m2ux + gyy * s_m2uy + gyz * s_m2uz;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3488"}},{"Add":{"text":" const dual s_m2z = gxz * s_m2ux + gyz * s_m2uy + gzz * s_m2uz;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3562"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3563"}},{"Add":{"text":" const dual s_m2r = s_m2x * x_r + s_m2y * y_r + s_m2z * z_r;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3634"}},{"Add":{"text":" const dual s_m2t = s_m2x * x_t + s_m2y * y_t + s_m2z * z_t;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3705"}},{"Add":{"text":" const dual ss_m2p = s_m2x * s_x_p + s_m2y * s_y_p + s_m2z * s_z_p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3783"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3784"}},{"Add":{"text":" const dual s_sm2 = sur * s_m2r + sut * s_m2t + z_sup * ss_m2p;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3858"}},{"Add":{"text":" assert(fabs(s_sm2) <= 1.0e-12);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3898"}},{"Add":{"text":" }","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3906"}}],"del_contents":[{"Del":{"text":" assert(fabs(s2 - 1) <= 1.0e-12);","pos":"FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC:17429"}}]}},{"Replacement":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":509,"ins":{"up_context":[{"change":3,"pos":13079}],"down_context":[{"change":3,"pos":13415}],"flag":"BLOCK","start":3908,"end":4397,"inode":{"change":3,"pos":2243}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":13079},"to":{"change":3,"start":17768,"end":17871},"introduced_by":3},{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":17871},"to":{"change":2,"start":2558,"end":2728},"introduced_by":2}],"inode":{"change":3,"pos":2243}},"ins_contents":[{"Add":{"text":" // Construct a 2-metric on the surface","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3952"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3953"}},{"Add":{"text":" const dual x_T =","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:3979"}},{"Add":{"text":" h_t * sin(theta) * cos(phi) + r * cos(theta) * cos(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4046"}},{"Add":{"text":" const dual y_T =","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4072"}},{"Add":{"text":" h_t * sin(theta) * sin(phi) + r * cos(theta) * sin(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4139"}},{"Add":{"text":" const dual z_T = h_t * cos(theta) - r * sin(theta);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4200"}},{"Add":{"text":" const dual s_x_P = s_h_p * sin(theta) * cos(phi) - r * sin(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4274"}},{"Add":{"text":" const dual s_y_P = s_h_p * sin(theta) * sin(phi) + r * cos(phi);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4348"}},{"Add":{"text":" const dual s_z_P = s_h_p * cos(theta);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4396"}}],"del_contents":[{"Del":{"text":" const T ss_detQ = (qtt * ss_qpp - pow(s_qtp, 2)).val;","pos":"FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC:17827"}},{"Del":{"text":" const T s_sqrt_detQ = sqrt(ss_detQ);","pos":"FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC:17870"}},{"Del":{"text":" // const T darea = sin(theta.val) * s_sqrt_detQ * geom.coord_dtheta(i, j)","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:2637"}},{"Del":{"text":" // *","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:2648"}},{"Del":{"text":" // geom.coord_dphi(i, j);","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:2696"}},{"Del":{"text":" const T A = s_sqrt_detQ;","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:2727"}}]}},{"Edit":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":520,"atom":{"NewVertex":{"up_context":[{"change":3,"pos":13416}],"down_context":[{"change":3,"pos":17996}],"flag":"BLOCK","start":4398,"end":6062,"inode":{"change":3,"pos":2243}}},"contents":[{"Add":{"text":" const dual qTT = 0 //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4444"}},{"Add":{"text":" + gxx * x_T * x_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4491"}},{"Add":{"text":" + gxy * x_T * y_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4538"}},{"Add":{"text":" + gxz * x_T * z_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4585"}},{"Add":{"text":" + gxy * y_T * x_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4632"}},{"Add":{"text":" + gyy * y_T * y_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4679"}},{"Add":{"text":" + gyz * y_T * z_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4726"}},{"Add":{"text":" + gxz * z_T * x_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4773"}},{"Add":{"text":" + gyz * z_T * y_T //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4820"}},{"Add":{"text":" + gzz * z_T * z_T;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4865"}},{"Add":{"text":" const dual s_qTP = 0 //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4916"}},{"Add":{"text":" + gxx * x_T * s_x_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:4967"}},{"Add":{"text":" + gxy * x_T * s_y_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5018"}},{"Add":{"text":" + gxz * x_T * s_z_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5069"}},{"Add":{"text":" + gxy * y_T * s_x_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5120"}},{"Add":{"text":" + gyy * y_T * s_y_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5171"}},{"Add":{"text":" + gyz * y_T * s_z_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5222"}},{"Add":{"text":" + gxz * z_T * s_x_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5273"}},{"Add":{"text":" + gyz * z_T * s_y_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5324"}},{"Add":{"text":" + gzz * z_T * s_z_P;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5373"}},{"Add":{"text":" const dual ss_qPP = 0 //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5427"}},{"Add":{"text":" + gxx * s_x_P * s_x_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5481"}},{"Add":{"text":" + gxy * s_x_P * s_y_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5535"}},{"Add":{"text":" + gxz * s_x_P * s_z_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5589"}},{"Add":{"text":" + gxy * s_y_P * s_x_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5643"}},{"Add":{"text":" + gyy * s_y_P * s_y_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5697"}},{"Add":{"text":" + gyz * s_y_P * s_z_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5751"}},{"Add":{"text":" + gxz * s_z_P * s_x_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5805"}},{"Add":{"text":" + gyz * s_z_P * s_y_P //","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5859"}},{"Add":{"text":" + gzz * s_z_P * s_z_P;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5911"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5912"}},{"Add":{"text":" const dual ss_detq2 = qTT * ss_qPP - pow(s_qTP, 2);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:5973"}},{"Add":{"text":" const dual s_sqrt_detq2 = sqrt(ss_detq2);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6024"}},{"Add":{"text":" const T A = s_sqrt_detq2.val;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6060"}},{"Add":{"text":"","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6061"}}]}},{"Replacement":{"inode":{"change":3,"pos":2243},"path":"AHFinder/src/ahfinder.cxx","line":856,"ins":{"up_context":[{"change":3,"pos":28213}],"down_context":[{"change":3,"pos":28271}],"flag":"BLOCK","start":6063,"end":6252,"inode":{"change":3,"pos":2243}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":28213},"to":{"change":2,"start":3425,"end":3558},"introduced_by":2}],"inode":{"change":3,"pos":2243}},"ins_contents":[{"Add":{"text":" // CCTK_VWARN(CCTK_WARN_ALERT, \"Using Brill-Lindquist metric\");","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6128"}},{"Add":{"text":" // const auto metric = brill_lindquist_metric(cctkGH, coords);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6193"}},{"Add":{"text":" const auto metric = interpolate_metric(cctkGH, coords);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6251"}}],"del_contents":[{"Del":{"text":"#warning \"TODO\"","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:3440"}},{"Del":{"text":" // const auto metric = interpolate_metric(cctkGH, coords);","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:3501"}},{"Del":{"text":" const auto metric = calculate_metric(cctkGH, coords);","pos":"UYXLP77T2LHOV62FHRAX637S2DNE275ZI2KNCWNOVJEX5KUNKOFAC:3557"}}]}},{"Replacement":{"inode":{"change":3,"pos":30913},"path":"AHFinder/src/discretization.hxx","line":36,"ins":{"up_context":[{"change":3,"pos":31599}],"down_context":[{"change":3,"pos":31647}],"flag":"BLOCK","start":6253,"end":6374,"inode":{"change":3,"pos":30913}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":31599},"to":{"change":3,"start":31599,"end":31647},"introduced_by":3}],"inode":{"change":3,"pos":30913}},"ins_contents":[{"Add":{"text":" const double theta = ssht_sampling_mw_t2theta(i, nmodes);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6314"}},{"Add":{"text":" assert(theta >= 0 && theta <= M_PI);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6355"}},{"Add":{"text":" return theta;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6373"}}],"del_contents":[{"Del":{"text":" return ssht_sampling_mw_t2theta(i, nmodes);","pos":"FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC:31646"}}]}},{"Replacement":{"inode":{"change":3,"pos":30913},"path":"AHFinder/src/discretization.hxx","line":44,"ins":{"up_context":[{"change":3,"pos":31771}],"down_context":[{"change":3,"pos":31817}],"flag":"BLOCK","start":6375,"end":6489,"inode":{"change":3,"pos":30913}},"del":{"edges":[{"previous":"BLOCK","flag":"BLOCK | DELETED","from":{"change":3,"pos":31771},"to":{"change":3,"start":31771,"end":31817},"introduced_by":3}],"inode":{"change":3,"pos":30913}},"ins_contents":[{"Add":{"text":" const double phi = ssht_sampling_mw_p2phi(j, nmodes);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6432"}},{"Add":{"text":" assert(phi >= 0 && phi < 2 * M_PI);","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6472"}},{"Add":{"text":" return phi;","pos":"K6AJPVHABM2UUWBFGLQG6XCXN2EARPR6MQ4FPXSQ3FKMMJP43F3QC:6488"}}],"del_contents":[{"Del":{"text":" return ssht_sampling_mw_p2phi(j, nmodes);","pos":"FRJS2N2AVBNG3DWJLR57OYJDGGQTSOQ5B5DC7U6L7V3FGRZ5PQCAC:31816"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review