From 907d5641d5e82205c052f53874278e961ae3ec05 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Mon, 30 Jan 2017 00:18:56 +0100
Subject: [PATCH] Testing integration point MeshProperties.
---
Tests/MeshLib/MeshProperties.cpp | 54 ++++++++++++++++++++++++++++++++
1 file changed, 54 insertions(+)
diff --git a/Tests/MeshLib/MeshProperties.cpp b/Tests/MeshLib/MeshProperties.cpp
index 13f9804f898..14b24478a14 100644
--- a/Tests/MeshLib/MeshProperties.cpp
+++ b/Tests/MeshLib/MeshProperties.cpp
@@ -14,6 +14,7 @@
#include "MeshLib/MeshGenerators/MeshGenerator.h"
#include "MeshLib/Mesh.h"
+#include "MeshLib/Elements/Element.h"
#include "MeshLib/PropertyVector.h"
class MeshLibProperties : public ::testing::Test
@@ -49,6 +50,59 @@ TEST_F(MeshLibProperties, PropertyVectorTestMetaData)
ASSERT_EQ(0u, p->size());
}
+TEST_F(MeshLibProperties, PropertyVectorTestIntegrationPoint)
+{
+ ASSERT_TRUE(mesh != nullptr);
+
+ int const n_integration_points = 4;
+ int const vector_length = 6;
+ std::string const prop_name("ip_field");
+
+ auto* const p_ptr = mesh->getProperties().createNewPropertyVector<double>(
+ prop_name, MeshLib::MeshItemType::IntegrationPoint, vector_length);
+ ASSERT_TRUE(p_ptr != nullptr);
+
+ auto& p = *p_ptr;
+ ASSERT_EQ(p.getPropertyName(), prop_name);
+ ASSERT_EQ(MeshLib::MeshItemType::IntegrationPoint, p.getMeshItemType());
+ ASSERT_EQ(vector_length, p.getNumberOfComponents());
+ ASSERT_EQ(0u, p.size());
+
+ // Fill the property vector with double data in following pattern:
+ // for element's id xx and integration point number yy the value is xx.yy.
+ p.resize(mesh->getNumberOfElements() * n_integration_points);
+
+ std::vector<std::size_t> offsets(mesh->getNumberOfElements());
+ std::size_t offset = 0; // last position in the property vector
+ for (auto const& e : mesh->getElements())
+ {
+ offsets[e->getID()] = offset;
+ for (int ip = 0; ip < n_integration_points; ++ip)
+ {
+ p[offset + ip] = e->getID() + ip * 0.01;
+ }
+ offset += n_integration_points;
+ }
+
+ // Check the values at each offset.
+
+ // Last element is checked after the for-loop.
+ std::size_t i = 0;
+ for (; i < offsets.size() - 1; ++i)
+ {
+ std::size_t const size = offsets[i + 1] - offsets[i];
+ ASSERT_EQ(n_integration_points, size);
+ for (int ip = 0; ip < n_integration_points; ++ip)
+ ASSERT_EQ(i + ip * 0.01, p[offsets[i] + ip]);
+ }
+ { // Last element
+ std::size_t const size = p.size() - offsets[i];
+ ASSERT_EQ(n_integration_points, size);
+ for (int ip = 0; ip < n_integration_points; ++ip)
+ ASSERT_EQ(i + ip * 0.01, p[offsets[i] + ip]);
+ }
+}
+
TEST_F(MeshLibProperties, AddDoubleProperties)
{
ASSERT_TRUE(mesh != nullptr);
--
GitLab