From: Daniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Date: Sat, 18 May 2019 17:55:49 +0000 (+0200)
Subject: Change CUDAWrappers::MatrixFree::point_type to Point
X-Git-Tag: v9.1.0-rc3~1^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F8253%2Fhead;p=dealii.git

Change CUDAWrappers::MatrixFree::point_type to Point
---

diff --git a/examples/step-64/step-64.cu b/examples/step-64/step-64.cu
index 26b742ce01..abc67f1696 100644
--- 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 b1d4f349df..304278e6e4 100644
--- 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 81a4e45d2a..19e2f5a3aa 100644
--- 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);
         }