{"header":{"message":"Poisson: New thorn\n","description":null,"timestamp":"2020-06-08T20:14:56Z","authors":[{"email":"schnetter@gmail.com","name":"Erik Schnetter"}]},"hash":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC","authors":[],"deps":{"hashes":[["3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC","Poisson: New thorn\n"],["AA",null],["ZDS34ZNQ5NNS2FVYGL5WIT7KMCDGKEQC2LSYA6AGP5YSRRZTRSPAC","Pipelines: Update thorn list\n"],["HNQLJR6ZWME6VBJ2Y7PSENLJPXC7INSS7NC2CKIWQAC776CQ74TQC","Pipelines: Build Cactus\n"],["W5IGLN4UQIR5MF5MIQT25JQZ5F2DGRRHJFLN3WDL3T6PLSWRS5DAC","New thorn ODESolvers\n"]],"deps":["ZDS34ZNQ5NNS2FVYGL5WIT7KMCDGKEQC2LSYA6AGP5YSRRZTRSPAC"],"known":[]},"hunks":[{"FileAdd":{"meta":{"basename":"Poisson","metadata":1005},"path":"Poisson","context":[{"change":1,"pos":0}],"contents":[]}},{"FileAdd":{"meta":{"basename":"README","metadata":420},"path":"Poisson/README","context":[{"change":0,"pos":10}],"contents":[{"Add":{"text":"Cactus Code Thorn Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:46"}},{"Add":{"text":"Author(s) : Erik Schnetter ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:98"}},{"Add":{"text":"Maintainer(s): Erik Schnetter ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:150"}},{"Add":{"text":"Licence : LGPL","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:170"}},{"Add":{"text":"--------------------------------------------------------------------------","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:245"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:246"}},{"Add":{"text":"1. Purpose","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:257"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:258"}},{"Add":{"text":"Solve the Poisson equation","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:285"}}]}},{"FileAdd":{"meta":{"basename":"configuration.ccl","metadata":420},"path":"Poisson/configuration.ccl","context":[{"change":0,"pos":10}],"contents":[{"Add":{"text":"# Configuration definitions for thorn Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:353"}}]}},{"FileAdd":{"meta":{"basename":"interface.ccl","metadata":420},"path":"Poisson/interface.ccl","context":[{"change":0,"pos":10}],"contents":[{"Add":{"text":"# Interface definition for thorn Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:412"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:413"}},{"Add":{"text":"IMPLEMENTS: Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:433"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:434"}},{"Add":{"text":"USES INCLUDE HEADER: loop.hxx","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:464"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:465"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:466"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:467"}},{"Add":{"text":"void FUNCTION SolvePoisson()","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:496"}},{"Add":{"text":"REQUIRES FUNCTION SolvePoisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:527"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:528"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:529"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:530"}},{"Add":{"text":"CCTK_REAL rhs TYPE=gf TAGS='index={0 0 0}' \"rhs\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:579"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:580"}},{"Add":{"text":"CCTK_REAL phi TYPE=gf TAGS='index={0 0 0}' \"phi\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:629"}},{"Add":{"text":"CCTK_REAL res TYPE=gf TAGS='index={0 0 0}' \"res\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:678"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:679"}},{"Add":{"text":"CCTK_REAL ires TYPE=gf TAGS='index={0 0 0}' \"ires\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:730"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:731"}},{"Add":{"text":"CCTK_REAL asol TYPE=gf TAGS='index={0 0 0}' \"asol\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:782"}},{"Add":{"text":"CCTK_REAL aerr TYPE=gf TAGS='index={0 0 0}' \"aerr\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:833"}}]}},{"FileAdd":{"meta":{"basename":"par","metadata":1005},"path":"Poisson/par","context":[{"change":0,"pos":10}],"contents":[]}},{"FileAdd":{"meta":{"basename":"poisson.par","metadata":420},"path":"Poisson/par/poisson.par","context":[{"change":0,"pos":841}],"contents":[{"Add":{"text":"ActiveThorns = \"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:873"}},{"Add":{"text":" CarpetX","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:885"}},{"Add":{"text":" Formaline","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:899"}},{"Add":{"text":" IOUtil","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:910"}},{"Add":{"text":" Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:922"}},{"Add":{"text":"\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:924"}},{"Add":{"text":" ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:926"}},{"Add":{"text":"$nlevels = 1","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:939"}},{"Add":{"text":"$ncells = 8","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:951"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:952"}},{"Add":{"text":"Cactus::cctk_show_schedule = no","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:984"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:985"}},{"Add":{"text":"Cactus::presync_mode = \"mixed-error\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1022"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1023"}},{"Add":{"text":"Cactus::cctk_itlast = 0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1047"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1048"}},{"Add":{"text":"CarpetX::verbose = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1071"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1072"}},{"Add":{"text":"CarpetX::xmin = -1.0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1093"}},{"Add":{"text":"CarpetX::ymin = -1.0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1114"}},{"Add":{"text":"CarpetX::zmin = -1.0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1135"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1136"}},{"Add":{"text":"CarpetX::xmax = +1.0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1157"}},{"Add":{"text":"CarpetX::ymax = +1.0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1178"}},{"Add":{"text":"CarpetX::zmax = +1.0","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1199"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1200"}},{"Add":{"text":"CarpetX::ncells_x = $ncells","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1228"}},{"Add":{"text":"CarpetX::ncells_y = $ncells","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1256"}},{"Add":{"text":"CarpetX::ncells_z = $ncells","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1284"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1285"}},{"Add":{"text":"CarpetX::periodic_x = no","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1310"}},{"Add":{"text":"CarpetX::periodic_y = no","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1335"}},{"Add":{"text":"CarpetX::periodic_z = no","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1360"}},{"Add":{"text":"CarpetX::periodic = no","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1383"}},{"Add":{"text":"# CarpetX::reflection_x = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1413"}},{"Add":{"text":"# CarpetX::reflection_y = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1443"}},{"Add":{"text":"# CarpetX::reflection_z = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1473"}},{"Add":{"text":"# CarpetX::reflection_upper_x = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1509"}},{"Add":{"text":"# CarpetX::reflection_upper_y = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1545"}},{"Add":{"text":"# CarpetX::reflection_upper_z = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1581"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1582"}},{"Add":{"text":"CarpetX::ghost_size = 1","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1606"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1607"}},{"Add":{"text":"CarpetX::max_num_levels = $nlevels","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1642"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1643"}},{"Add":{"text":"CarpetX::prolongation_type = \"ddf\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1678"}},{"Add":{"text":"CarpetX::prolongation_order = 1","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1710"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1711"}},{"Add":{"text":"IO::out_dir = $parfile","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1734"}},{"Add":{"text":"IO::out_every = 1","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1752"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1753"}},{"Add":{"text":"CarpetX::out_silo_vars = \"all\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1784"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1785"}},{"Add":{"text":"CarpetX::out_tsv = yes","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1808"}}]}},{"FileAdd":{"meta":{"basename":"param.ccl","metadata":420},"path":"Poisson/param.ccl","context":[{"change":0,"pos":10}],"contents":[{"Add":{"text":"# Parameter definitions for thorn Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1864"}}]}},{"FileAdd":{"meta":{"basename":"schedule.ccl","metadata":420},"path":"Poisson/schedule.ccl","context":[{"change":0,"pos":10}],"contents":[{"Add":{"text":"# Schedule definitions for thorn Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1922"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1923"}},{"Add":{"text":"SCHEDULE Poisson_setup AT initial","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1957"}},{"Add":{"text":"{","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1959"}},{"Add":{"text":" LANG: C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1969"}},{"Add":{"text":" WRITES: phi(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:1995"}},{"Add":{"text":" WRITES: rhs(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2021"}},{"Add":{"text":"} \"Set up Poisson equation\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2049"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2050"}},{"Add":{"text":"SCHEDULE Poisson_solve AT initial AFTER Poisson_setup","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2104"}},{"Add":{"text":"{","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2106"}},{"Add":{"text":" LANG: C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2116"}},{"Add":{"text":" OPTIONS: level","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2133"}},{"Add":{"text":" READS: phi(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2158"}},{"Add":{"text":" READS: rhs(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2183"}},{"Add":{"text":" WRITES: phi(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2209"}},{"Add":{"text":" WRITES: res(interior)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2233"}},{"Add":{"text":"} \"Solve Poisson equation (this overwrites the boundary conditions)\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2302"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2303"}},{"Add":{"text":"# SCHEDULE Poisson_fixup AT initial AFTER Poisson_solve","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2360"}},{"Add":{"text":"# {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2364"}},{"Add":{"text":"# LANG: C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2376"}},{"Add":{"text":"# WRITES: phi(boundary)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2402"}},{"Add":{"text":"# } \"Apply boundary conditions again\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2440"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2441"}},{"Add":{"text":"SCHEDULE Poisson_residual AT initial AFTER (Poisson_solve, Poisson_fixup)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2515"}},{"Add":{"text":"{","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2517"}},{"Add":{"text":" LANG: C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2527"}},{"Add":{"text":" READS: phi(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2552"}},{"Add":{"text":" WRITES: ires(interior)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2577"}},{"Add":{"text":" SYNC: ires","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2590"}},{"Add":{"text":"} \"Evaluate residual\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2612"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2613"}},{"Add":{"text":"SCHEDULE Poisson_residual_boundary AT initial AFTER Poisson_residual","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2682"}},{"Add":{"text":"{","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2684"}},{"Add":{"text":" LANG: C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2694"}},{"Add":{"text":" READS: phi(boundary)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2717"}},{"Add":{"text":" WRITES: ires(boundary)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2742"}},{"Add":{"text":"} \"Evaluate residual on boundary\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2776"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2777"}},{"Add":{"text":"SCHEDULE Poisson_error AT initial AFTER (Poisson_solve, Poisson_fixup)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2848"}},{"Add":{"text":"{","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2850"}},{"Add":{"text":" LANG: C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2860"}},{"Add":{"text":" READS: phi(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2885"}},{"Add":{"text":" WRITES: asol(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2912"}},{"Add":{"text":" WRITES: aerr(everywhere)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2939"}},{"Add":{"text":"} \"Analyse Poisson solution\"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:2968"}}]}},{"FileAdd":{"meta":{"basename":"src","metadata":1005},"path":"Poisson/src","context":[{"change":0,"pos":10}],"contents":[]}},{"FileAdd":{"meta":{"basename":"make.code.defn","metadata":420},"path":"Poisson/src/make.code.defn","context":[{"change":0,"pos":2976}],"contents":[{"Add":{"text":"# Main make.code.defn file for thorn Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3039"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3040"}},{"Add":{"text":"# Source files in this directory","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3073"}},{"Add":{"text":"SRCS = poisson.cxx","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3092"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3093"}},{"Add":{"text":"# Subdirectories containing source files","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3134"}},{"Add":{"text":"SUBDIRS =","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3144"}}]}},{"FileAdd":{"meta":{"basename":"poisson.cxx","metadata":420},"path":"Poisson/src/poisson.cxx","context":[{"change":0,"pos":2976}],"contents":[{"Add":{"text":"#include ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3180"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3181"}},{"Add":{"text":"#include ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3199"}},{"Add":{"text":"#include ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3235"}},{"Add":{"text":"#include ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3264"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3265"}},{"Add":{"text":"#include ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3284"}},{"Add":{"text":"#include ","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3301"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3302"}},{"Add":{"text":"namespace Poisson {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3322"}},{"Add":{"text":"using namespace Loop;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3344"}},{"Add":{"text":"using namespace std;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3365"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3366"}},{"Add":{"text":"////////////////////////////////////////////////////////////////////////////////","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3447"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3448"}},{"Add":{"text":"constexpr CCTK_REAL rmin = 0.1;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3480"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3481"}},{"Add":{"text":"constexpr CCTK_REAL A = 1.0;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3510"}},{"Add":{"text":"constexpr CCTK_REAL W = 0.1;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3539"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3540"}},{"Add":{"text":"// Gaussian RHS","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3556"}},{"Add":{"text":"template constexpr T frhs(T x, T y, T z) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3612"}},{"Add":{"text":" T r = sqrt(pow(x, 2) + pow(y, 2) + pow(z, 2));","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3661"}},{"Add":{"text":" return 4 * M_PI * A * exp(-pow(r / (sqrt(2) * W), 2));","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3718"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3720"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3721"}},{"Add":{"text":"// Solution for the above RHS, assuming homogeneous Dirichlet boundary","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3792"}},{"Add":{"text":"// conditions at infinity","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3818"}},{"Add":{"text":"template constexpr T fsol(T x, T y, T z) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3874"}},{"Add":{"text":" T r = sqrt(pow(x, 2) + pow(y, 2) + pow(z, 2));","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3923"}},{"Add":{"text":" T C = -A * pow(sqrt(2 * M_PI) * W, 3);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3964"}},{"Add":{"text":" if (r < 1.0e-12)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:3983"}},{"Add":{"text":" return C * 2 / (sqrt(2 * M_PI) * W);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4024"}},{"Add":{"text":" return C * erf(r / (sqrt(2) * W)) / r;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4065"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4067"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4068"}},{"Add":{"text":"template constexpr T fbnd(T x, T y, T z) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4124"}},{"Add":{"text":" T r = sqrt(pow(x, 2) + pow(y, 2) + pow(z, 2));","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4173"}},{"Add":{"text":" T alpha = -A * pow(sqrt(2 * M_PI) * W, 3);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4218"}},{"Add":{"text":" return alpha / fmax(rmin, r);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4250"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4252"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4253"}},{"Add":{"text":"// Robin boundary conditions:","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4283"}},{"Add":{"text":"// f(r) = \\alpha / r + C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4310"}},{"Add":{"text":"// f'(r) = - \\alpha / r^2","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4338"}},{"Add":{"text":"//","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4341"}},{"Add":{"text":"// r f(r) - r C = \\alpha","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4368"}},{"Add":{"text":"// -r^2 f'(r) = \\alpha","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4393"}},{"Add":{"text":"//","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4396"}},{"Add":{"text":"// r f(r) - r C = -r^2 f'(r)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4427"}},{"Add":{"text":"//","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4430"}},{"Add":{"text":"// f(r) + r f'(r) = C","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4454"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4455"}},{"Add":{"text":"////////////////////////////////////////////////////////////////////////////////","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4536"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4537"}},{"Add":{"text":"extern \"C\" void Poisson_setup(CCTK_ARGUMENTS) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4585"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_Poisson_setup;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4625"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4652"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4653"}},{"Add":{"text":" const GF3D phi_(cctkGH, phi);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4705"}},{"Add":{"text":" const GF3D rhs_(cctkGH, rhs);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4757"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4758"}},{"Add":{"text":" // Note: The boundary conditions are applied on the outermost layer","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4828"}},{"Add":{"text":" // of the interior points, not on the boundary points.","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4885"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4886"}},{"Add":{"text":" // Initialize all points (including ghosts) and set boundary conditions","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:4960"}},{"Add":{"text":" loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { phi_(p.I) = 0.0; });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5035"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5036"}},{"Add":{"text":" // Set boundary conditions","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5065"}},{"Add":{"text":" loop_int<0, 0, 0>(cctkGH, [&](const PointDesc &p) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5119"}},{"Add":{"text":" if (cctk_lbnd[0] + p.i == 1 || cctk_lbnd[0] + p.i == cctk_gsh[0] - 1 ||","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5195"}},{"Add":{"text":" cctk_lbnd[1] + p.j == 1 || cctk_lbnd[1] + p.j == cctk_gsh[1] - 1 ||","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5271"}},{"Add":{"text":" cctk_lbnd[2] + p.k == 1 || cctk_lbnd[2] + p.k == cctk_gsh[2] - 1)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5345"}},{"Add":{"text":" phi_(p.I) = fbnd(p.x, p.y, p.z);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5384"}},{"Add":{"text":" });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5390"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5391"}},{"Add":{"text":" // Set RHS","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5404"}},{"Add":{"text":" loop_all<0, 0, 0>(","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5425"}},{"Add":{"text":" cctkGH, [&](const PointDesc &p) { rhs_(p.I) = frhs(p.x, p.y, p.z); });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5502"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5504"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5505"}},{"Add":{"text":"extern \"C\" void Poisson_solve(CCTK_ARGUMENTS) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5553"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_Poisson_solve;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5593"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5620"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5621"}},{"Add":{"text":" SolvePoisson();","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5639"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5641"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5642"}},{"Add":{"text":"// extern \"C\" void Poisson_fixup(CCTK_ARGUMENTS) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5693"}},{"Add":{"text":"// DECLARE_CCTK_ARGUMENTS_Poisson_fixup;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5736"}},{"Add":{"text":"// DECLARE_CCTK_PARAMETERS;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5766"}},{"Add":{"text":"//","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5769"}},{"Add":{"text":"// const GF3D phi_(cctkGH, phi);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5824"}},{"Add":{"text":"//","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5827"}},{"Add":{"text":"// // Set boundary conditions","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5859"}},{"Add":{"text":"// loop_bnd<0, 0, 0>(","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5883"}},{"Add":{"text":"// cctkGH, [&](const PointDesc &p) { phi_(p.I) = fbnd(p.x, p.y, p.z); });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5963"}},{"Add":{"text":"// }","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5968"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:5969"}},{"Add":{"text":"extern \"C\" void Poisson_residual(CCTK_ARGUMENTS) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6020"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_Poisson_residual;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6063"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6090"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6091"}},{"Add":{"text":" const GF3D phi_(cctkGH, phi);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6149"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6150"}},{"Add":{"text":" const GF3D ires_(cctkGH, ires);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6204"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6205"}},{"Add":{"text":" loop_int<0, 0, 0>(cctkGH, [&](const PointDesc &p) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6259"}},{"Add":{"text":" if (cctk_lbnd[0] + p.i == 1 || cctk_lbnd[0] + p.i == cctk_gsh[0] - 1 ||","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6335"}},{"Add":{"text":" cctk_lbnd[1] + p.j == 1 || cctk_lbnd[1] + p.j == cctk_gsh[1] - 1 ||","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6411"}},{"Add":{"text":" cctk_lbnd[2] + p.k == 1 || cctk_lbnd[2] + p.k == cctk_gsh[2] - 1)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6485"}},{"Add":{"text":" ires_(p.I) = phi_(p.I) - fbnd(p.x, p.y, p.z);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6537"}},{"Add":{"text":" else","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6546"}},{"Add":{"text":" ires_(p.I) = (phi_(p.I - p.DI(0)) - 2 * phi_(p.I) + phi_(p.I + p.DI(0))) /","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6627"}},{"Add":{"text":" pow(p.DX[0], 2) +","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6668"}},{"Add":{"text":" (phi_(p.I - p.DI(1)) - 2 * phi_(p.I) + phi_(p.I + p.DI(1))) /","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6749"}},{"Add":{"text":" pow(p.DX[1], 2) +","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6790"}},{"Add":{"text":" (phi_(p.I - p.DI(2)) - 2 * phi_(p.I) + phi_(p.I + p.DI(2))) /","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6871"}},{"Add":{"text":" pow(p.DX[2], 2) -","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6912"}},{"Add":{"text":" frhs(p.x, p.y, p.z);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6952"}},{"Add":{"text":" });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6958"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6960"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:6961"}},{"Add":{"text":"extern \"C\" void Poisson_residual_boundary(CCTK_ARGUMENTS) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7021"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_Poisson_residual_boundary;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7073"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7100"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7101"}},{"Add":{"text":" const GF3D phi_(cctkGH, phi);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7159"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7160"}},{"Add":{"text":" const GF3D ires_(cctkGH, ires);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7214"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7215"}},{"Add":{"text":" loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7269"}},{"Add":{"text":" // ires_(p.I) = phi_(p.I) - fbnd(p.x, p.y, p.z);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7322"}},{"Add":{"text":" ires_(p.I) = 0.0;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7344"}},{"Add":{"text":" });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7350"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7352"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7353"}},{"Add":{"text":"extern \"C\" void Poisson_error(CCTK_ARGUMENTS) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7401"}},{"Add":{"text":" DECLARE_CCTK_ARGUMENTS_Poisson_error;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7441"}},{"Add":{"text":" DECLARE_CCTK_PARAMETERS;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7468"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7469"}},{"Add":{"text":" const GF3D phi_(cctkGH, phi);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7527"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7528"}},{"Add":{"text":" const GF3D asol_(cctkGH, asol);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7582"}},{"Add":{"text":" const GF3D aerr_(cctkGH, aerr);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7636"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7637"}},{"Add":{"text":" loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) {","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7691"}},{"Add":{"text":" if (cctk_lbnd[0] + p.i == 1 || cctk_lbnd[0] + p.i == cctk_gsh[0] - 1 ||","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7767"}},{"Add":{"text":" cctk_lbnd[1] + p.j == 1 || cctk_lbnd[1] + p.j == cctk_gsh[1] - 1 ||","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7843"}},{"Add":{"text":" cctk_lbnd[2] + p.k == 1 || cctk_lbnd[2] + p.k == cctk_gsh[2] - 1)","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7917"}},{"Add":{"text":" asol_(p.I) = 0.0;","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7941"}},{"Add":{"text":" else","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7950"}},{"Add":{"text":" asol_(p.I) = fsol(p.x, p.y, p.z);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:7990"}},{"Add":{"text":" aerr_(p.I) = phi_(p.I) - asol_(p.I);","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:8031"}},{"Add":{"text":" });","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:8037"}},{"Add":{"text":"}","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:8039"}},{"Add":{"text":"","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:8040"}},{"Add":{"text":"} // namespace Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:8063"}}]}},{"Edit":{"inode":{"change":4,"pos":11},"path":"azure-pipelines/carpetx.th","line":148,"atom":{"NewVertex":{"up_context":[{"change":5,"pos":786}],"down_context":[{"change":2,"pos":519}],"flag":"BLOCK","start":8065,"end":8081,"inode":{"change":4,"pos":11}}},"contents":[{"Add":{"text":"CarpetX/Poisson","pos":"3CP3B7DKC4YAZNE5PMOPZXWBCQMDXYWHZMZUUHDXE5QVKRW2JYCQC:8080"}}]}}],"owner":"eschnett","repo":"cactusamrex","login":null}

…

Dependencies

Change contents

Review