From ee531d17b1c2751a92dbb098f5cfe7064578d455 Mon Sep 17 00:00:00 2001
From: "Dmitry Yu. Naumov" <github@naumov.de>
Date: Mon, 1 Oct 2018 12:39:57 +0200
Subject: [PATCH] [NL] Assert scalar quantity and element type.
---
NumLib/Function/Interpolation.h | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/NumLib/Function/Interpolation.h b/NumLib/Function/Interpolation.h
index 876ac7d8ee0..cc9f5a4fa3e 100644
--- a/NumLib/Function/Interpolation.h
+++ b/NumLib/Function/Interpolation.h
@@ -88,8 +88,8 @@ void shapeFunctionInterpolate(
interpolated_values...);
}
-/// Interpolates \c node_values given in lower order element nodes (e.g. the
-/// base nodes) to higher order element's nodes (e.g. quadratic nodes) and
+/// Interpolates scalar \c node_values given in lower order element nodes (e.g.
+/// the base nodes) to higher order element's nodes (e.g. quadratic nodes) and
/// writes the result into the global property vector.
///
/// The base nodes' values are copied. For each higher order node the shape
@@ -102,6 +102,9 @@ void interpolateToHigherOrderNodes(
Eigen::MatrixBase<EigenMatrixType> const& node_values,
MeshLib::PropertyVector<double>& interpolated_values_global_vector)
{
+ assert(dynamic_cast<HigherOrderMeshElementType const*>(&element));
+ assert(node_values.cols() == 1); // Scalar quantity only.
+
using SF = LowerOrderShapeFunction;
using ShapeMatricesType = ShapeMatrixPolicyType<SF, GlobalDim>;
using ShapeMatrices = typename ShapeMatricesType::ShapeMatrices;
--
GitLab