diff --git a/examples/laplacexy/laplace_perp/square/BOUT.inp b/examples/laplacexy/laplace_perp/square/BOUT.inp index 95ecb7119f..38fb175d90 100644 --- a/examples/laplacexy/laplace_perp/square/BOUT.inp +++ b/examples/laplacexy/laplace_perp/square/BOUT.inp @@ -1,5 +1,5 @@ -input = sin(pi*x - y) +input_field = sin(pi*x - y) non_uniform = true # Include corrections to second derivatives diff --git a/examples/laplacexy/laplace_perp/test.cxx b/examples/laplacexy/laplace_perp/test.cxx index 1312b005da..75577bb8b7 100644 --- a/examples/laplacexy/laplace_perp/test.cxx +++ b/examples/laplacexy/laplace_perp/test.cxx @@ -1,4 +1,5 @@ #include +#include #include #include @@ -10,11 +11,15 @@ int main(int argc, char** argv) { BoutInitialise(argc, argv); /////////////////////////////////////// - bool calc_metric; - calc_metric = Options::root()["calc_metric"].withDefault(false); + const bool calc_metric = Options::root()["calc_metric"].withDefault(false); if (calc_metric) { // Read metric tensor - Field2D Rxy, Btxy, Bpxy, B0, hthe, I; + Field2D Rxy; + Field2D Btxy; + Field2D Bpxy; + Field2D B0; + Field2D hthe; + Field2D I; mesh->get(Rxy, "Rxy"); // m mesh->get(Btxy, "Btxy"); // T mesh->get(Bpxy, "Bpxy"); // T @@ -47,13 +52,14 @@ int main(int argc, char** argv) { /////////////////////////////////////// // Read an analytic input - Field2D input = FieldFactory::get()->create2D("input", Options::getRoot(), mesh); + const Field2D input = + FieldFactory::get()->create2D("input_field", Options::getRoot(), mesh); // Create a LaplaceXY solver - LaplaceXY* laplacexy = new LaplaceXY(mesh); + LaplaceXY laplacexy{mesh}; // Solve, using 0.0 as starting guess - Field2D solved = laplacexy->solve(input, 0.0); + Field2D solved = laplacexy.solve(input, 0.0); // Need to communicate guard cells mesh->communicate(solved); diff --git a/examples/laplacexy/laplace_perp/torus/BOUT.inp b/examples/laplacexy/laplace_perp/torus/BOUT.inp index 365294174f..4109b44bab 100644 --- a/examples/laplacexy/laplace_perp/torus/BOUT.inp +++ b/examples/laplacexy/laplace_perp/torus/BOUT.inp @@ -1,5 +1,5 @@ -input = sin(pi*x - y) +input_field = sin(pi*x - y) calc_metric = true # Read Rxy, Bpxy etc and calculate metric non_uniform = true # Include corrections to second derivatives diff --git a/tests/MMS/laplace/data/BOUT.inp b/tests/MMS/laplace/data/BOUT.inp index 1965b1e8f7..5b3cc37b09 100644 --- a/tests/MMS/laplace/data/BOUT.inp +++ b/tests/MMS/laplace/data/BOUT.inp @@ -5,7 +5,7 @@ MXG = 1 solution = sin(pi*x) -input = -pi^2*sin(pi*x) +input_field = -pi^2*sin(pi*x) [mesh] diff --git a/tests/MMS/laplace/laplace.cxx b/tests/MMS/laplace/laplace.cxx index 54dbaaba67..fbcdee355c 100644 --- a/tests/MMS/laplace/laplace.cxx +++ b/tests/MMS/laplace/laplace.cxx @@ -1,16 +1,21 @@ #include - #include +#include #include #include +#include + +#include using bout::globals::mesh; +using namespace std::string_literals; int main(int argc, char** argv) { int init_err = BoutInitialise(argc, argv); if (init_err < 0) { return 0; - } else if (init_err > 0) { + } + if (init_err > 0) { return init_err; } @@ -31,16 +36,14 @@ int main(int argc, char** argv) { FieldFactory fact(mesh); - std::shared_ptr gen = fact.parse("input"); - output << "GEN = " << gen->str() << endl; - - Field3D input = fact.create3D("input"); - - Field3D result = lap->solve(input); - - Field3D solution = fact.create3D("solution"); + const auto input_name = "input_field"s; + const auto gen = fact.parse(input_name); + output.write("GEN = {}\n", gen->str()); - Field3D error = result - solution; + const Field3D input = fact.create3D(input_name); + const Field3D result = lap->solve(input); + const Field3D solution = fact.create3D("solution"); + const Field3D error = result - solution; Options dump; dump["input"] = input; diff --git a/tests/MMS/spatial/d2dx2/data/BOUT.inp b/tests/MMS/spatial/d2dx2/data/BOUT.inp index 5c1deaf0af..cdc0d96139 100644 --- a/tests/MMS/spatial/d2dx2/data/BOUT.inp +++ b/tests/MMS/spatial/d2dx2/data/BOUT.inp @@ -5,7 +5,7 @@ MXG = 1 MYG = 1 # No guard cells in Y -input = sin(0.5*pi*x) +input_field = sin(0.5*pi*x) solution = -sin(0.5*pi*x) * 0.25*pi*pi [mesh] diff --git a/tests/MMS/spatial/d2dx2/test_d2dx2.cxx b/tests/MMS/spatial/d2dx2/test_d2dx2.cxx index 6b35b719e1..7f7127cea9 100644 --- a/tests/MMS/spatial/d2dx2/test_d2dx2.cxx +++ b/tests/MMS/spatial/d2dx2/test_d2dx2.cxx @@ -4,7 +4,9 @@ #include #include +#include #include +#include using bout::globals::mesh; @@ -12,8 +14,9 @@ int main(int argc, char** argv) { BoutInitialise(argc, argv); - Field3D input = FieldFactory::get()->create3D("input", Options::getRoot(), mesh); - Field3D solution = FieldFactory::get()->create3D("solution", Options::getRoot(), mesh); + Field3D input = FieldFactory::get()->create3D("input_field", Options::getRoot(), mesh); + const Field3D solution = + FieldFactory::get()->create3D("solution", Options::getRoot(), mesh); // At this point the boundary cells are set to the analytic solution input.setBoundary("bndry"); @@ -21,7 +24,7 @@ int main(int argc, char** argv) { // Boundaries of input now set using extrapolation around mid-point boundary - Field3D result = D2DX2(input); + const Field3D result = D2DX2(input); // At this point result is not set in the boundary cells Options dump; diff --git a/tests/MMS/spatial/d2dz2/data/BOUT.inp b/tests/MMS/spatial/d2dz2/data/BOUT.inp index 6618d6f574..640411c112 100644 --- a/tests/MMS/spatial/d2dz2/data/BOUT.inp +++ b/tests/MMS/spatial/d2dz2/data/BOUT.inp @@ -5,7 +5,7 @@ zperiod = 1 MXG = 0 # No guard cells in X MYG = 1 # No guard cells in Y -input = sin(z) +input_field = sin(z) solution = -sin(z) [mesh] diff --git a/tests/MMS/spatial/d2dz2/test_d2dz2.cxx b/tests/MMS/spatial/d2dz2/test_d2dz2.cxx index 0b7aa2cde4..68a12b062f 100644 --- a/tests/MMS/spatial/d2dz2/test_d2dz2.cxx +++ b/tests/MMS/spatial/d2dz2/test_d2dz2.cxx @@ -4,7 +4,9 @@ #include #include +#include #include +#include using bout::globals::mesh; @@ -12,8 +14,10 @@ int main(int argc, char** argv) { BoutInitialise(argc, argv); - Field3D input = FieldFactory::get()->create3D("input", Options::getRoot(), mesh); - Field3D solution = FieldFactory::get()->create3D("solution", Options::getRoot(), mesh); + const Field3D input = + FieldFactory::get()->create3D("input_field", Options::getRoot(), mesh); + const Field3D solution = + FieldFactory::get()->create3D("solution", Options::getRoot(), mesh); Field3D result = D2DZ2(input);