Finish tests on post-processings

Author: thelfer

bindings/python/tests/PostProcessingTest.py

@@ -31,5 +31,24 @@ def test_behaviour(self):
         self.assertTrue("PrincipalStrain" in postprocessings, "invalid post-processing names")
         pass
 
+    def test_behaviour_data(self):
+        eps = 1.e-14
+        e = numpy.asarray([1.3e-2, 1.2e-2, 1.4e-2, 0., 0., 0.], dtype=numpy.float64)
+        e2 = numpy.asarray([1.2e-2, 1.3e-2, 1.4e-2, 0., 0., 0.], dtype=numpy.float64)
+        h = mgis_bv.Hypothesis.Tridimensional
+        b = self.__get_behaviour(h)
+        m = mgis_bv.MaterialDataManager(b, 2)
+        mgis_bv.setMaterialProperty(m.s1, "YoungModulus", 150e9);
+        mgis_bv.setMaterialProperty(m.s1, "PoissonRatio", 0.3);
+        mgis_bv.setExternalStateVariable(m.s1, "Temperature", 293.15);
+        mgis_bv.update(m);
+        m.s1.gradients[0:] = e;
+        m.s1.gradients[1:] = e;
+        outputs =numpy.empty(shape = (2, 3), dtype=numpy.float64)
+        mgis_bv.executePostProcessing(outputs.reshape(6), m, "PrincipalStrain");
+        for i in range(0, 3):
+             self.assertTrue(abs(outputs[0,i] - e2[i]) < eps, "invalid output value")
+             self.assertTrue(abs(outputs[1,i] - e2[i]) < eps, "invalid output value")
+
 if __name__ == '__main__':
     unittest.main()

tests/PostProcessingTest.cxx

@@ -13,6 +13,7 @@
  */
 
 #include <cmath>
+#include <array>
 #include <cstdlib>
 #include <iostream>
 #include <string_view>
@@ -53,6 +54,63 @@ void check_behaviour(const mgis::behaviour::Behaviour& b,
   check(getVariableSize(o, h) == 3, "invalid post-processing output size");
 }  // end of check_behaviour
 
+void call_postprocessing(const mgis::behaviour::Behaviour& b){
+  using namespace mgis::behaviour;
+  constexpr auto e =
+      std::array<mgis::real, 6u>{1.3e-2, 1.2e-2, 1.4e-2, 0., 0., 0.};
+  constexpr auto e2 =
+      std::array<mgis::real, 6u>{1.2e-2, 1.3e-2, 1.4e-2, 0., 0., 0.};
+  constexpr auto eps = 10 * std::numeric_limits<mgis::real>::epsilon();
+  auto d = BehaviourData{b};
+  // initialize the states
+  setExternalStateVariable(d.s0, "Temperature", 293.15);
+  setExternalStateVariable(d.s1, "Temperature", 293.15);
+  //
+  for (mgis::size_type i = 0; i != 6; ++i) {
+    d.s1.gradients[i] = e[i];
+  }
+  //
+  auto outputs = allocatePostProcessingVariables(b, "PrincipalStrain");
+  if (!check(outputs.size() == 3u, "invalid outputs initialisation")) {
+    return;
+  }
+  auto v = make_view(d);
+  executePostProcessing(outputs, v, b, "PrincipalStrain");
+  for (mgis::size_type i = 0; i != 3; ++i) {
+    check(std::abs(outputs[i] - e2[i]) < eps, "invalid output value");
+  }
+} // end of call_postprocessing
+
+void call_postprocessing2(const mgis::behaviour::Behaviour& b){
+  using namespace mgis::behaviour;
+  constexpr auto e =
+      std::array<mgis::real, 6u>{1.3e-2, 1.2e-2, 1.4e-2, 0., 0., 0.};
+  constexpr auto e2 =
+      std::array<mgis::real, 6u>{1.2e-2, 1.3e-2, 1.4e-2, 0., 0., 0.};
+  constexpr auto eps = 10 * std::numeric_limits<mgis::real>::epsilon();
+  auto m = MaterialDataManager{b, 2u};
+  // initialize the states
+  setMaterialProperty(m.s1, "YoungModulus", 150e9);
+  setMaterialProperty(m.s1, "PoissonRatio", 0.3);
+  setExternalStateVariable(m.s1, "Temperature", 293.15);
+  update(m);
+  //
+  for (mgis::size_type i = 0; i != 6; ++i) {
+    m.s1.gradients[i] = e[i];
+    m.s1.gradients[6 + i] = e[i];
+  }
+  //
+  auto outputs = allocatePostProcessingVariables(m, "PrincipalStrain");
+  if (!check(outputs.size() == 6u, "invalid outputs initialisation")) {
+    return;
+  }
+  executePostProcessing(outputs, m, "PrincipalStrain");
+  for (mgis::size_type i = 0; i != 3; ++i) {
+    check(std::abs(outputs[i] - e2[i]) < eps, "invalid output value");
+    check(std::abs(outputs[3 + i] - e2[i]) < eps, "invalid output value");
+  }
+} // end of call_postprocessing2
+
 int main(const int argc, const char* const* argv) {
   using namespace mgis::behaviour;
   constexpr const auto h = Hypothesis::TRIDIMENSIONAL;
@@ -62,6 +120,8 @@ int main(const int argc, const char* const* argv) {
   try {
     const auto b = load(argv[1], argv[2], h);
     check_behaviour(b, h);
+    call_postprocessing(b);
+    call_postprocessing2(b);
   } catch (std::exception& e) {
     std::cerr << e.what() << std::endl;
     return EXIT_FAILURE;

tests/PostProcessingTest.mfront

@@ -31,5 +31,5 @@ nu.setGlossaryName("PoissonRatio");
 εᵖ.setEntryName("PrincipalStrain");
 //! compute the principal strain
 @PostProcessing PrincipalStrain {
-  εᵖ = eto.computeEigenValues();
+  εᵖ = eto.computeEigenValues(Stensor::ASCENDING);
 }