OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

add some tests for SIM[2|3]D::projectSolution

Browse Source
This commit is contained in:
Arne Morten Kvarving 2016-09-02 10:14:16 +02:00 committed by Knut Morten Okstad
parent ff49abc34f
commit ea76fc6cf9
1 changed files with 102 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

102
src/SIM/Test/TestSIM.C
View File

@ -10,9 +10,45 @@
//!
//==============================================================================
#include "IntegrandBase.h"
#include "IFEM.h"
#include "SIM2D.h"
#include "SIM3D.h"
#include "gtest/gtest.h"
#include "tinyxml.h"
template<class Dim>
class TestProjectSIM : public Dim {
public:
TestProjectSIM(const SIMbase::CharVec& nf) : Dim(nf)
{
Dim::myProblem = new TestProjectIntegrand(Dim::dimension);
}
private:
class TestProjectIntegrand : public IntegrandBase {
public:
TestProjectIntegrand(int dim) : IntegrandBase(dim) {}
//! \brief Evaluates the secondary solution at a result point.
//! \param[out] s The solution field values at current point
//! \param[in] fe Finite element data at current point
//! \param[in] X Cartesian coordinates of current point
//! \param[in] MNPC Nodal point correspondance for the basis function values
bool evalSol(Vector& s, const FiniteElement& fe, const Vec3& X,
const std::vector<int>& MNPC) const override
{
s.resize(1);
s(1) = X[0] + X[1] + X[2];
return true;
}
size_t getNoFields(int = 2) const override { return 1; }
};
};
TEST(TestSIM, UniqueBoundaryNodes)
{
@ -27,3 +63,69 @@ TEST(TestSIM, UniqueBoundaryNodes)
std::sort(vec.begin(), vec.end());
ASSERT_TRUE(std::unique(vec.begin(), vec.end()) == vec.end());
}
TEST(TestSIM2D, ProjectSolution)
{
TestProjectSIM<SIM2D> sim({1});
sim.createDefaultModel();
ASSERT_TRUE(sim.preprocess());
Matrix ssol;
ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
size_t n = 1;
for (size_t j = 0; j < 2; ++j)
for (size_t i = 0; i < 2; ++i)
ASSERT_FLOAT_EQ(ssol(1, n++), i + j);
}
TEST(TestSIM2D, ProjectSolutionMixed)
{
TestProjectSIM<SIM2D> sim({1,1});
sim.createDefaultModel();
ASSERT_TRUE(sim.preprocess());
Matrix ssol;
ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
size_t n = 1;
for (size_t j = 0; j < 3; ++j)
for (size_t i = 0; i < 3; ++i)
ASSERT_FLOAT_EQ(ssol(1, n++), i/2.0 + j/2.0);
}
TEST(TestSIM3D, ProjectSolution)
{
TestProjectSIM<SIM3D> sim({1});
sim.createDefaultModel();
ASSERT_TRUE(sim.preprocess());
Matrix ssol;
ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
size_t n = 1;
for (size_t k = 0; k < 2; ++k)
for (size_t j = 0; j < 2; ++j)
for (size_t i = 0; i < 2; ++i)
ASSERT_FLOAT_EQ(ssol(1, n++), i + j + k);
}
TEST(TestSIM3D, ProjectSolutionMixed)
{
TestProjectSIM<SIM3D> sim({1,1});
sim.createDefaultModel();
ASSERT_TRUE(sim.preprocess());
Matrix ssol;
ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
size_t n = 1;
for (size_t k = 0; k < 3; ++k)
for (size_t j = 0; j < 3; ++j)
for (size_t i = 0; i < 3; ++i)
ASSERT_FLOAT_EQ(ssol(1, n++), i/2.0 + j/2.0 + k/2.0);
}