Change CUDAWrappers::MatrixFree::point_type to Point

diff --git a/examples/step-64/step-64.cu b/examples/step-64/step-64.cu
index 26b742ce012c54ef3df4067d09a9af15dcc48c61..abc67f1696ea8ca5b86d2c711f6e17c2991166ca 100644 (file)
--- a/examples/step-64/step-64.cu
+++ b/examples/step-64/step-64.cu
@@ -100,7 +100,7 @@ namespace Step64
   {
     const unsigned int pos = CUDAWrappers::local_q_point_id<dim, double>(
       cell, gpu_data, n_dofs_1d, n_q_points);
-    const Tensor<1, dim> q_point =
+    const Point<dim> q_point =
       CUDAWrappers::get_quadrature_point<dim, double>(cell,
                                                       gpu_data,
                                                       n_dofs_1d);

diff --git a/include/deal.II/matrix_free/cuda_matrix_free.h b/include/deal.II/matrix_free/cuda_matrix_free.h
index b1d4f349df02f066db21b62e8e2a1a7b4884ea8e..304278e6e42834f82ea230b31fc42a9bab38dc00 100644 (file)
--- a/include/deal.II/matrix_free/cuda_matrix_free.h
+++ b/include/deal.II/matrix_free/cuda_matrix_free.h
@@ -78,8 +78,7 @@ namespace CUDAWrappers
   {
   public:
     using jacobian_type = Tensor<2, dim, Tensor<1, dim, Number>>;
-    // TODO this should really be a CUDAWrappers::Point
-    using point_type = Tensor<1, dim, Number>;
+    using point_type    = Point<dim, Number>;
 
     /**
      * Parallelization scheme used: parallel_in_elem (parallelism at the level

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 81a4e45d2ab11a5b2249f397525db2ddcba45b00..19e2f5a3aae004bb8c02938ec17c3fc42b9f4779 100644 (file)
--- a/include/deal.II/matrix_free/cuda_matrix_free.templates.h
+++ b/include/deal.II/matrix_free/cuda_matrix_free.templates.h
@@ -362,7 +362,7 @@ namespace CUDAWrappers
             transpose_in_place(q_points_host, n_cells, padding_length);
 
           alloc_and_copy(&data->q_points[color],
-                         ArrayView<const Tensor<1, dim, Number>>(
+                         ArrayView<const Point<dim, Number>>(
                            q_points_host.data(), q_points_host.size()),
                          n_cells * padding_length);
         }