TimeGridsTest.cpp

Go to the documentation of this file.

#ifndef PRECICE_NO_MPI
 
#include <Eigen/Core>
#include <Eigen/Dense>
 
#include "acceleration/SharedPointer.hpp"
#include "acceleration/test/helper.hpp"
#include "cplscheme/CouplingData.hpp"
#include "cplscheme/SharedPointer.hpp"
#include "mesh/Mesh.hpp"
 
#include "testing/TestContext.hpp"
#include "testing/Testing.hpp"
#include "time/TimeGrids.hpp"
 
BOOST_AUTO_TEST_SUITE(TimeTests)
 
using namespace precice;
using namespace precice::time;
 
using precice::testing::makeCouplingData;
 
struct TimeGridsTestsFixture {
  using DataMap = std::map<int, cplscheme::PtrCouplingData>;
};
 
BOOST_FIXTURE_TEST_SUITE(TimeGridsTest, TimeGridsTestsFixture)
 
PRECICE_TEST_SETUP(1_rank)
BOOST_AUTO_TEST_CASE(TestMoveAndScaleTimeGrids)
{
  PRECICE_TEST();
  std::vector<int> dataIDs;
  dataIDs.push_back(0);
  mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
 
  mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
 
  // init displacements and waveform
 
  double timeWindowStart = 0.0;
  double dt              = 1.0;
  double timeWindowEnd   = timeWindowStart + dt;
 
  displacements->values().resize(4);
  displacements->values() << 1.0, 1.0, 1.0, 1.0;
  displacements->setSampleAtTime(timeWindowStart, displacements->sample());
  displacements->setSampleAtTime(timeWindowStart + 0.5 * dt, displacements->sample());
  displacements->setSampleAtTime(timeWindowEnd, displacements->sample());
 
  cplscheme::PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, true);
 
  DataMap data;
  data.insert(std::pair<int, cplscheme::PtrCouplingData>(0, dpcd));
 
  time::TimeGrids fullTimeGrid(data, dataIDs, false);
  time::TimeGrids reducedTimeGrid(data, dataIDs, true);
 
  Eigen::VectorXd storedFullTimeGrid    = fullTimeGrid.getTimeGridAfter(dataIDs[0], timeWindowStart);
  Eigen::VectorXd storedReducedTimeGrid = reducedTimeGrid.getTimeGridAfter(dataIDs[0], timeWindowStart);

  BOOST_TEST(storedFullTimeGrid.size() == 2);
  BOOST_TEST(testing::equals(storedFullTimeGrid(0), 0.5));
  BOOST_TEST(testing::equals(storedFullTimeGrid(1), 1.0));
 
  BOOST_TEST(storedReducedTimeGrid.size() == 1);
  BOOST_TEST(testing::equals(storedReducedTimeGrid(0), 1.0));
 
  // Create a second data pointer with the contentes of the second time window
  mesh::PtrData newdisplacements(new mesh::Data("dvalues", -1, 1));
 
  // move to next window
  timeWindowStart = timeWindowEnd;
  dt              = 2 * dt;
  timeWindowEnd   = timeWindowStart + dt;
 
  // init newdisplacements
  newdisplacements->values().resize(4);
  newdisplacements->values() << 1.0, 1.0, 1.0, 1.0;
  newdisplacements->setSampleAtTime(timeWindowStart, newdisplacements->sample());
  newdisplacements->setSampleAtTime(timeWindowEnd, newdisplacements->sample());
 
  cplscheme::PtrCouplingData newdpcd = makeCouplingData(newdisplacements, dummyMesh, true);
 
  DataMap newdata;
  newdata.insert(std::pair<int, cplscheme::PtrCouplingData>(0, newdpcd));
 
  // move the time grid to the new time window
  fullTimeGrid.moveTimeGridToNewWindow(newdata);
  reducedTimeGrid.moveTimeGridToNewWindow(newdata);
 
  storedFullTimeGrid    = fullTimeGrid.getTimeGridAfter(dataIDs[0], timeWindowStart);
  storedReducedTimeGrid = reducedTimeGrid.getTimeGridAfter(dataIDs[0], timeWindowStart);

  BOOST_TEST(storedFullTimeGrid.size() == 2);
  BOOST_TEST(testing::equals(storedFullTimeGrid(0), 2.0));
  BOOST_TEST(testing::equals(storedFullTimeGrid(1), 3.0));
 
  BOOST_TEST(storedReducedTimeGrid.size() == 1);
  BOOST_TEST(testing::equals(storedReducedTimeGrid(0), 3.0));
}
 
BOOST_AUTO_TEST_SUITE_END()
BOOST_AUTO_TEST_SUITE_END()
#endif