Class Field2DInterpToTest

tests/unit/mesh/test_interpolation.cxx:276–337

}

class Field2DInterpToTest : public ::testing::Test {
protected:
  Field2DInterpToTest() : input(mesh) {
    input = 0.;
    input(2, 2) = 2.;
    input(1, 2) = 1.8;
    input(0, 2) = 1.7;
    input(3, 2) = 1.3;
    input(4, 2) = 1.5;
    input(2, 1) = 3.8;
    input(2, 0) = 3.7;
    input(2, 3) = 3.3;
    input(2, 4) = 3.5;
  }

  static void SetUpTestCase() {
    WithQuietOutput quiet_info{output_info};
    WithQuietOutput quiet_warn{output_warn};
    delete mesh;
    mesh = new FakeMesh(nx, ny, nz);
    mesh->StaggerGrids = true;
    mesh->xstart = 2;
    mesh->ystart = 2;
    mesh->xend = nx - 3;
    mesh->yend = ny - 3;

    mesh->createDefaultRegions();

    // We need Coordinates so a parallel transform is available as
    // FieldFactory::create3D wants to un-field-align the result
    for (const auto& location :
         std::list<CELL_LOC>{CELL_CENTRE, CELL_XLOW, CELL_YLOW, CELL_ZLOW}) {

      static_cast<FakeMesh*>(mesh)->setCoordinates(nullptr, location);
      static_cast<FakeMesh*>(mesh)->setCoordinates(
          std::make_shared<Coordinates>(
              mesh, Field2D{1.0, mesh}, Field2D{1.0, mesh}, BoutReal{1.0},
              Field2D{1.0, mesh}, Field2D{0.0, mesh}, Field2D{1.0, mesh},
              Field2D{1.0, mesh}, Field2D{1.0, mesh}, Field2D{0.0, mesh},
              Field2D{0.0, mesh}, Field2D{0.0, mesh}, Field2D{1.0, mesh},
              Field2D{1.0, mesh}, Field2D{1.0, mesh}, Field2D{0.0, mesh},
              Field2D{0.0, mesh}, Field2D{0.0, mesh}, Field2D{0.0, mesh},
              Field2D{0.0, mesh}),
          location);
      // No call to Coordinates::geometry() needed here
      mesh->getCoordinates(location)->setParallelTransform(
          bout::utils::make_unique<ParallelTransformIdentity>(*mesh));
    }
  }

  static void TearDownTestCase() {
    delete mesh;
    mesh = nullptr;
  }

  Field2D input;

public: