{"header":{"message":"CarpetX: Add linear solver\n","description":null,"timestamp":"2019-12-16T15:10:54Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC","authors":[],"deps":{"hashes":[["JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC","CarpetX: Add linear solver\n"],["AA",null],["HVV73EHEMEFVCPS75XZAFUGWOY7LBYOJOIDNHJUK4MB5QTS3CTJQC","CarpetX: Export function CallScheduleGroup\n"],["722HZ7UFINNE3YKSYKP2NHZ5XEG5QQLQHSKC7PREJZR3EX6RDYUAC","Import initial AMReX driver for Cactus\n"],["MSBBCXVGD3GRLE5KAI6BKAFRV7SQUWI2SNN43AJAUD3ISRCEXY6QC","Snapshot before transitioning from AmrMesh to AmrCore\n"],["GKKJ75HX2ERLVBZVE2CUB6T3J2SUT7R3UKEKTEYNOG43ZKX6X5MQC","The big rename: AMReX -> CarpetX\n"],["DGEGTDYVGRHEGNFDHCCVTD6AVLRJ2XHDXQPWR5MMR62V5NU2XGIQC","Add generic prolongation operators\n"]],"deps":["HVV73EHEMEFVCPS75XZAFUGWOY7LBYOJOIDNHJUK4MB5QTS3CTJQC"],"known":[]},"hunks":[{"Edit":{"inode":{"change":4,"pos":17976},"path":"CarpetX/interface.ccl","line":24,"atom":{"NewVertex":{"up_context":[{"change":5,"pos":2}],"down_context":[{"change":5,"pos":2}],"flag":"BLOCK","start":0,"end":97,"inode":{"change":4,"pos":17976}}},"contents":[{"Add":{"text":"void FUNCTION SolvePoisson()","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:28"}},{"Add":{"text":"PROVIDES FUNCTION SolvePoisson WITH CarpetX_SolvePoisson LANGUAGE C","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:96"}}]}},{"Edit":{"inode":{"change":4,"pos":17976},"path":"CarpetX/interface.ccl","line":28,"atom":{"NewVertex":{"up_context":[{"change":2,"pos":190}],"down_context":[{"change":5,"pos":3}],"flag":"BLOCK","start":98,"end":99,"inode":{"change":4,"pos":17976}}},"contents":[{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:98"}}]}},{"FileAdd":{"meta":{"basename":"linsolve.cxx","metadata":420},"path":"CarpetX/src/linsolve.cxx","context":[{"change":6,"pos":2677}],"contents":[{"Add":{"text":"#include \"driver.hxx\"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:137"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:138"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:156"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:184"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:213"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:214"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:233"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:257"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:292"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:293"}},{"Add":{"text":"#include ","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:311"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:312"}},{"Add":{"text":"namespace CarpetX {","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:332"}},{"Add":{"text":"using namespace amrex;","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:355"}},{"Add":{"text":"using namespace std;","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:376"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:377"}},{"Add":{"text":"extern \"C\" void CarpetX_SolvePoisson() {","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:418"}},{"Add":{"text":" // Create operator","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:439"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:440"}},{"Add":{"text":" Vector geoms(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:491"}},{"Add":{"text":" Vector grids(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:542"}},{"Add":{"text":" Vector dmaps(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:604"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level) {","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:675"}},{"Add":{"text":" geoms.at(level) = ghext->amrcore->Geom(level);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:726"}},{"Add":{"text":" grids.at(level) = ghext->amrcore->boxArray(level);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:781"}},{"Add":{"text":" dmaps.at(level) = ghext->amrcore->DistributionMap(level);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:843"}},{"Add":{"text":" }","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:847"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:848"}},{"Add":{"text":" MLNodeLaplacian mlnodelaplacian(geoms, grids, dmaps);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:904"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:905"}},{"Add":{"text":" mlnodelaplacian.setDomainBC(","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:936"}},{"Add":{"text":" {LinOpBCType::Periodic, LinOpBCType::Periodic, LinOpBCType::Periodic},","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1013"}},{"Add":{"text":" {LinOpBCType::Periodic, LinOpBCType::Periodic, LinOpBCType::Periodic});","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1091"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1092"}},{"Add":{"text":" vector sigmas(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1144"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level) {","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1215"}},{"Add":{"text":" auto &sigma = sigmas.at(level);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1251"}},{"Add":{"text":" sigma.define(ghext->amrcore->boxArray(level),","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1301"}},{"Add":{"text":" ghext->amrcore->DistributionMap(level), 1, 0);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1365"}},{"Add":{"text":" sigma.setVal(1.0);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1388"}},{"Add":{"text":" mlnodelaplacian.setSigma(level, sigma);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1432"}},{"Add":{"text":" }","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1436"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1437"}},{"Add":{"text":" mlnodelaplacian.setVerbose(10);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1471"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1472"}},{"Add":{"text":" // Create solver","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1491"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1492"}},{"Add":{"text":" MLMG mlmg(mlnodelaplacian);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1522"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1523"}},{"Add":{"text":" const int gi_rhs = CCTK_GroupIndex(\"Maxwell::dive\");","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1578"}},{"Add":{"text":" assert(gi_rhs >= 0);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1601"}},{"Add":{"text":" const int gi_sol = CCTK_GroupIndex(\"Maxwell::phi1\");","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1656"}},{"Add":{"text":" assert(gi_sol >= 0);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1679"}},{"Add":{"text":" const int gi_res = CCTK_GroupIndex(\"Maxwell::dive1\");","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1735"}},{"Add":{"text":" assert(gi_res >= 0);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1758"}},{"Add":{"text":" const int tl = 0;","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1778"}},{"Add":{"text":" const int vi = 0;","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1798"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1799"}},{"Add":{"text":" Vector ress(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1851"}},{"Add":{"text":" Vector sols(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1903"}},{"Add":{"text":" Vector rhss(ghext->leveldata.size());","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:1961"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level) {","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2032"}},{"Add":{"text":" const auto &restrict leveldata = ghext->leveldata.at(level);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2097"}},{"Add":{"text":" const auto &restrict groupdata_rhs = leveldata.groupdata.at(gi_rhs);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2170"}},{"Add":{"text":" rhss.at(level) = groupdata_rhs.mfab.at(tl).get();","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2224"}},{"Add":{"text":" const auto &restrict groupdata_sol = leveldata.groupdata.at(gi_sol);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2297"}},{"Add":{"text":" sols.at(level) = groupdata_sol.mfab.at(tl).get();","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2351"}},{"Add":{"text":" const auto &restrict groupdata_res = leveldata.groupdata.at(gi_res);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2424"}},{"Add":{"text":" ress.at(level) = groupdata_res.mfab.at(tl).get();","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2478"}},{"Add":{"text":" }","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2482"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2483"}},{"Add":{"text":" mlmg.setVerbose(10);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2506"}},{"Add":{"text":" mlmg.setBottomVerbose(10);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2535"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2536"}},{"Add":{"text":" // Initial guess","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2555"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2556"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level) {","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2627"}},{"Add":{"text":" sols.at(level)->setVal(0.0, vi);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2664"}},{"Add":{"text":" }","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2668"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2669"}},{"Add":{"text":" // Solve","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2680"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2681"}},{"Add":{"text":" mlmg.compResidual(ress, sols, rhss);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2720"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level)","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2789"}},{"Add":{"text":" CCTK_VINFO(\"norm_inf rhs[%d]: %g\", level,","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2835"}},{"Add":{"text":" double(rhss.at(level)->norminf(0, 0, false, true)));","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2903"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level)","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:2972"}},{"Add":{"text":" CCTK_VINFO(\"norm_inf sol[%d]: %g\", level,","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3018"}},{"Add":{"text":" double(sols.at(level)->norminf(0, 0, false, true)));","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3086"}},{"Add":{"text":" for (int level = 0; level < int(ghext->leveldata.size()); ++level)","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3155"}},{"Add":{"text":" CCTK_VINFO(\"norm_inf res[%d]: %g\", level,","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3201"}},{"Add":{"text":" double(ress.at(level)->norminf(0, 0, false, true)));","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3269"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3270"}},{"Add":{"text":" const CCTK_REAL rtol = 0.0;","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3300"}},{"Add":{"text":" const CCTK_REAL atol = 1.0e-12;","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3334"}},{"Add":{"text":" mlmg.solve(sols, rhss, rtol, atol);","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3372"}},{"Add":{"text":"}","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3374"}},{"Add":{"text":"","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3375"}},{"Add":{"text":"} // namespace CarpetX","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3398"}}]}},{"Edit":{"inode":{"change":4,"pos":15114},"path":"CarpetX/src/make.code.defn","line":8,"atom":{"NewVertex":{"up_context":[{"change":7,"pos":1026}],"down_context":[{"change":7,"pos":1055}],"flag":"BLOCK","start":3400,"end":3450,"inode":{"change":4,"pos":15114}}},"contents":[{"Add":{"text":" linsolve.cxx \\","pos":"JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC:3449"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review