[T] Added IP data setter unit tests

Tests/ProcessLib/TestIPDataAccess.cpp

@@ -34,7 +34,7 @@ struct ProcessLib_IPDataAccess : ::testing::Test
 
         constexpr int off_diag_size = dim == 2 ? 1 : 3;
 
-        static constexpr std::size_t num_int_pts = 10;
+        constexpr std::size_t num_int_pts = 10;
 
         std::vector<IPData<dim>> ip_data(num_int_pts);
 
@@ -52,6 +52,24 @@ struct ProcessLib_IPDataAccess : ::testing::Test
         return ip_data;
     }
 
+    static std::vector<IPData<Dim::value>> getIPDataNaNs()
+    {
+        using KV = typename IPData<dim>::KV;
+
+        constexpr std::size_t num_int_pts = 10;
+        constexpr double nan = std::numeric_limits<double>::quiet_NaN();
+
+        std::vector<IPData<dim>> ip_data(num_int_pts);
+
+        for (std::size_t i = 0; i < num_int_pts; ++i)
+        {
+            ip_data[i].kelvin = KV::Constant(nan);
+            ip_data[i].scalar = nan;
+        }
+
+        return ip_data;
+    }
+
     static std::vector<double> getScalarData()
     {
         return {100, 101, 102, 103, 104, 105, 106, 107, 108, 109};
@@ -162,3 +180,49 @@ TYPED_TEST(ProcessLib_IPDataAccess, GetKelvinVectorDataTransposedOrder)
                 testing::Pointwise(testing::DoubleEq(),
                                    this->getKVDataTransposedOrder()));
 }
+
+TYPED_TEST(ProcessLib_IPDataAccess, SetScalarData)
+{
+    constexpr int dim = TypeParam::value;
+
+    auto ip_data = this->getIPDataNaNs();
+
+    auto const cache = this->getScalarData();
+
+    auto const num_read = ProcessLib::setIntegrationPointScalarData(
+        &cache[0], ip_data, &IPData<dim>::scalar);
+
+    ASSERT_EQ(ip_data.size(), num_read);
+
+    auto const ip_data_expected = this->getIPData();
+
+    for (std::size_t i = 0; i < ip_data_expected.size(); ++i)
+    {
+        EXPECT_DOUBLE_EQ(ip_data_expected[i].scalar, ip_data[i].scalar)
+            << "Values at integration point " << i << " differ.";
+    }
+}
+
+TYPED_TEST(ProcessLib_IPDataAccess, SetKelvinVectorData)
+{
+    constexpr int dim = TypeParam::value;
+
+    auto ip_data = this->getIPDataNaNs();
+
+    auto const cache = this->getKVDataTransposedOrder();
+
+    auto const num_read = ProcessLib::setIntegrationPointKelvinVectorData<dim>(
+        &cache[0], ip_data, &IPData<dim>::kelvin);
+
+    ASSERT_EQ(ip_data.size(), num_read);
+
+    auto const ip_data_expected = this->getIPData();
+
+    for (std::size_t i = 0; i < ip_data_expected.size(); ++i)
+    {
+        EXPECT_THAT(
+            ip_data[i].kelvin,
+            testing::Pointwise(testing::DoubleEq(), ip_data_expected[i].kelvin))
+            << "Values at integration point " << i << " differ.";
+    }
+}