From fcf2ad9cbc43cb95338a5c0a0d96eb49bcfb1866 Mon Sep 17 00:00:00 2001 From: D Schneider Date: Sat, 13 Jun 2020 17:18:03 +0200 Subject: [PATCH] Remove make_vectorized and fix get_symmetric_gradient() --- include/deal.II/matrix_free/fe_evaluation.h | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h index 8379764156..cbde241e1f 100644 --- a/include/deal.II/matrix_free/fe_evaluation.h +++ b/include/deal.II/matrix_free/fe_evaluation.h @@ -3839,8 +3839,7 @@ FEEvaluationBase:: const auto cells = this->get_cell_ids(); // 2) actually gather values - VectorizedArrayType out = - make_vectorized_array(Number(1.)); + VectorizedArrayType out = Number(1.); for (unsigned int i = 0; i < VectorizedArrayType::size(); ++i) if (cells[i] != numbers::invalid_unsigned_int) out[i] = array[cells[i] / VectorizedArrayType::size()] @@ -6057,7 +6056,7 @@ inline DEAL_II_ALWAYS_INLINE SymmetricTensor<2, dim, VectorizedArrayType> // copy from generic function into dim-specialization function const Tensor<2, dim, VectorizedArrayType> grad = get_gradient(q_point); VectorizedArrayType symmetrized[(dim * dim + dim) / 2]; - VectorizedArrayType half = make_vectorized_array(0.5); + VectorizedArrayType half = Number(0.5); for (unsigned int d = 0; d < dim; ++d) symmetrized[d] = grad[d][d]; switch (dim) -- 2.39.5