Commit e8aab95b authored by Dmitry Yu. Naumov's avatar Dmitry Yu. Naumov

Extrapolation; Ignore empty int. point values.

Used in LIE for the stress/strain output, where the fracture
(lower order) elements should not contribute to the extrapolation.
parent 7819b3b2
......@@ -150,6 +150,12 @@ void LocalLinearLeastSquaresExtrapolator::extrapolateElement(
extrapolatables.getIntegrationPointValues(
element_index, t, x, dof_table, _integration_point_values_cache);
// Empty vector means to ignore the values and not to change the counts.
if (integration_point_values.empty())
{
return;
}
auto const& N_0 = extrapolatables.getShapeMatrix(element_index, 0);
auto const num_nodes = static_cast<unsigned>(N_0.cols());
auto const num_values =
......
......@@ -222,15 +222,14 @@ private:
std::vector<double> const& getIntPtSigma(
std::vector<double>& cache, std::size_t const /*component*/) const
{
cache.resize(_ip_data.size());
cache.resize(0);
return cache;
}
std::vector<double> const& getIntPtEpsilon(
std::vector<double>& cache, std::size_t const /*component*/) const
{
cache.resize(_ip_data.size());
cache.resize(0);
return cache;
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment