From 17af4ffed9278e0aa704bb1d00149d70db1195b5 Mon Sep 17 00:00:00 2001 From: Daniel Arndt Date: Thu, 26 Sep 2019 00:12:54 -0400 Subject: [PATCH] Use alias --- source/numerics/solution_transfer.cc | 16 +++++++--------- 1 file changed, 7 insertions(+), 9 deletions(-) diff --git a/source/numerics/solution_transfer.cc b/source/numerics/solution_transfer.cc index ede47089dc..2c8609409d 100644 --- a/source/numerics/solution_transfer.cc +++ b/source/numerics/solution_transfer.cc @@ -518,23 +518,21 @@ SolutionTransfer::interpolate( { const unsigned int old_index = pointerstruct->second.active_fe_index; + const FullMatrix &interpolation_matrix = + interpolation_hp(active_fe_index, old_index); // The interpolation matrix might be empty when using // FE_Nothing. - if (interpolation_hp(active_fe_index, old_index).empty()) + if (interpolation_matrix.empty()) tmp.reinit(dofs_per_cell, false); else { tmp.reinit(dofs_per_cell, true); const unsigned int old_index = pointerstruct->second.active_fe_index; - AssertDimension( - (*valuesptr)[j].size(), - interpolation_hp(active_fe_index, old_index).n()); - AssertDimension( - tmp.size(), - interpolation_hp(active_fe_index, old_index).m()); - interpolation_hp(active_fe_index, old_index) - .vmult(tmp, (*valuesptr)[j]); + AssertDimension((*valuesptr)[j].size(), + interpolation_matrix.n()); + AssertDimension(tmp.size(), interpolation_matrix.m()); + interpolation_matrix.vmult(tmp, (*valuesptr)[j]); } data = &tmp; } -- 2.39.5