From d2f3fe075cf88b5c27d83361081f55c30298b645 Mon Sep 17 00:00:00 2001 From: Daniel Arndt Date: Wed, 24 May 2023 16:27:53 -0400 Subject: [PATCH] Remove inv_jacobians instead --- include/deal.II/matrix_free/cuda_matrix_free.templates.h | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/include/deal.II/matrix_free/cuda_matrix_free.templates.h b/include/deal.II/matrix_free/cuda_matrix_free.templates.h index fdefd46e6a..d6ef02b34f 100644 --- a/include/deal.II/matrix_free/cuda_matrix_free.templates.h +++ b/include/deal.II/matrix_free/cuda_matrix_free.templates.h @@ -260,13 +260,11 @@ namespace CUDAWrappers if (update_flags & update_gradients) { - const std::vector> &inv_jacobians = - fe_values.get_inverse_jacobians(); for (unsigned int i = 0; i < q_points_per_cell; ++i) for (unsigned int d = 0; d < dim; ++d) for (unsigned int e = 0; e < dim; ++e) inv_jacobian_host(cell_id, i, d, e) = - inv_jacobians[i][d][e]; + fe_values.inverse_jacobian(i)[d][e]; } } -- 2.39.5