Added: Unit test that provokes the mixed/periodicity failure.

Changed: Using EXPECT instead of ASSERT to catch all errors.

Author: kmokstad

src/SIM/Test/TestSIM.C

@@ -12,10 +12,11 @@
 
 #include "SIM2D.h"
 #include "SIM3D.h"
+#include "ASMs3D.h"
+#include "ASMmxBase.h"
 #include "IntegrandBase.h"
 
 #include "gtest/gtest.h"
-#include "tinyxml.h"
 
 
 template<class Dim> class TestProjectSIM : public Dim
@@ -44,7 +45,7 @@ private:
                          const std::vector<int>&) const
     {
       s.resize(1);
-      s(1) = X[0] + X[1] + X[2];
+      s(1) = X.sum();
       return true;
     }
 
@@ -56,9 +57,9 @@ private:
 class DummyIntegrand : public IntegrandBase {};
 
 
-TEST(TestSIM, UniqueBoundaryNodes)
+TEST(TestSIM2D, UniqueBoundaryNodes)
 {
-  SIM2D sim(new DummyIntegrand(),1);
+  SIM2D sim(new DummyIntegrand());
   ASSERT_TRUE(sim.read("src/SIM/Test/refdata/boundary_nodes.xinp"));
   ASSERT_TRUE(sim.preprocess());
 
@@ -74,91 +75,120 @@ TEST(TestSIM2D, ProjectSolution)
 {
   TestProjectSIM<SIM2D> sim({1});
 
-  Matrix ssol;
-  ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
+  Vector ssol, psol(sim.getNoDOFs());
+  ASSERT_TRUE(sim.project(ssol,psol));
 
   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);
+      EXPECT_FLOAT_EQ(ssol(n++), i + j);
 }
 
 
 TEST(TestSIM2D, ProjectSolutionMixed)
 {
   TestProjectSIM<SIM2D> sim({1,1});
 
-  Matrix ssol;
-  ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
+  Vector ssol, psol(sim.getNoDOFs());
+  ASSERT_TRUE(sim.project(ssol,psol));
 
   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);
+      EXPECT_FLOAT_EQ(ssol(n++), i/2.0 + j/2.0);
 }
 
 
 TEST(TestSIM3D, ProjectSolution)
 {
   TestProjectSIM<SIM3D> sim({1});
 
-  Matrix ssol;
-  ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
+  Vector ssol, psol(sim.getNoDOFs());
+  ASSERT_TRUE(sim.project(ssol,psol));
 
   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);
+        EXPECT_FLOAT_EQ(ssol(n++), i + j + k);
 }
 
 
 TEST(TestSIM3D, ProjectSolutionMixed)
 {
   TestProjectSIM<SIM3D> sim({1,1});
 
-  Matrix ssol;
-  ASSERT_TRUE(sim.project(ssol, Vector(sim.getNoDOFs())));
+  Vector ssol, psol(sim.getNoDOFs());
+  ASSERT_TRUE(sim.project(ssol,psol));
 
   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);
+        EXPECT_FLOAT_EQ(ssol(n++), i/2.0 + j/2.0 + k/2.0);
 }
 
 
-TEST(TestSIM, InjectPatchSolution)
+TEST(TestSIM2D, InjectPatchSolution)
 {
   TestProjectSIM<SIM2D> sim({1,1});
+  size_t n1 = sim.getNoNodes(true,1);
+  size_t n2 = sim.getNoNodes(true,2);
 
-  Vector sol(2*sim.getNoNodes(true, 1) + sim.getNoNodes(true, 2));
-  Vector lsol(2*sim.getNoNodes(true, 1));
+  Vector sol(2*n1 + n2);
+  Vector lsol(2*n1);
   size_t i, ofs;
-  for (i = 0; i < sim.getNoNodes(true,1); i++)
+  for (i = 0; i < n1; i++)
     lsol[2*i] = lsol[2*i+1] = i+1;
 
   ASSERT_TRUE(sim.addMixedMADOF(1, 2));
   sim.injectPatchSolution(sol, lsol, 0, 2, 1);
-  for (i = ofs = 0; i < sim.getNoNodes(true,1); i++, ofs += 2) {
-    ASSERT_FLOAT_EQ(sol[ofs], i+1);
-    ASSERT_FLOAT_EQ(sol[ofs+1], i+1);
+  for (i = ofs = 0; i < n1; i++) {
+    EXPECT_FLOAT_EQ(sol[ofs++], i+1);
+    EXPECT_FLOAT_EQ(sol[ofs++], i+1);
   }
-  for (i = 0; i < sim.getNoNodes(true,2); i++, ofs++)
-    ASSERT_FLOAT_EQ(sol[ofs], 0);
+  for (i = 0; i < n2; i++)
+    EXPECT_FLOAT_EQ(sol[ofs++], 0);
 
-  ASSERT_TRUE(sim.addMixedMADOF(2, 2));
-  Vector sol2(sim.getNoNodes(true, 1) + 2*sim.getNoNodes(true, 2));
-  Vector lsol2(2*sim.getNoNodes(true, 2));
-  for (i = 0; i < sim.getNoNodes(true,2); i++)
+  Vector sol2(n1 + 2*n2);
+  Vector lsol2(2*n2);
+  for (i = 0; i < n2; i++)
     lsol2[2*i] = lsol2[2*i+1] = i+1;
 
+  ASSERT_TRUE(sim.addMixedMADOF(2, 2));
   sim.injectPatchSolution(sol2, lsol2, 0, 2, 2);
-  for (i = ofs = 0; i < sim.getNoNodes(true,1); i++, ofs++)
-    ASSERT_FLOAT_EQ(sol2[ofs], 0);
-
-  for (i = 0; i < sim.getNoNodes(true,2); i++, ofs += 2) {
-    ASSERT_FLOAT_EQ(sol2[ofs], i+1);
-    ASSERT_FLOAT_EQ(sol2[ofs+1], i+1);
+  for (i = ofs = 0; i < n1; i++)
+    EXPECT_FLOAT_EQ(sol2[ofs++], 0);
+  for (i = 0; i < n2; i++) {
+    EXPECT_FLOAT_EQ(sol2[ofs++], i+1);
+    EXPECT_FLOAT_EQ(sol2[ofs++], i+1);
   }
 }
+
+
+class TestSIM3Dmx : public SIM3D
+{
+public:
+  TestSIM3Dmx() : SIM3D({3,1}) { myProblem = new DummyIntegrand(); }
+};
+
+
+TEST(TestSIM3D, Periodic)
+{
+  ASMmxBase::Type = ASMmxBase::REDUCED_CONT_RAISE_BASIS1;
+
+  TestSIM3Dmx sim;
+  ASSERT_TRUE(sim.createDefaultModel());
+
+  ASMs3D* pch = dynamic_cast<ASMs3D*>(sim.getPatch(1));
+  ASSERT_TRUE(pch != nullptr);
+  ASSERT_TRUE(pch->uniformRefine(0,1));
+  ASSERT_TRUE(pch->uniformRefine(1,1));
+  ASSERT_TRUE(pch->uniformRefine(2,1));
+  ASSERT_TRUE(sim.createFEMmodel());
+
+  pch->closeFaces(1);
+  pch->closeFaces(3);
+
+  ASSERT_TRUE(sim.preprocess());
+}