Fix a bug in CUDAWrappers::MatrixFree::cell_loop when different partitioners are...

diff --git a/doc/news/changes/minor/20200610Turcksin b/doc/news/changes/minor/20200610Turcksin
new file mode 100644 (file)
index 0000000..9eabd04
--- /dev/null
+++ b/doc/news/changes/minor/20200610Turcksin
@@ -0,0 +1,5 @@
+Fixed: Fix a bug where CUDAWrappers::MatrixFree::cell_loop() would set the
+destination vector to zero if the partitioner of the MatrixFree object was
+different from the partitioner of the source or destination vector.
+<br>
+(Bruno Turcksin, 2020/06/10)

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 c7e81aa35f95701273fd50113f697ad287477d51..ba9f6283ec6c3d13ea9eebe5c912e723fdebc171 100644 (file)
--- a/include/deal.II/matrix_free/cuda_matrix_free.templates.h
+++ b/include/deal.II/matrix_free/cuda_matrix_free.templates.h
@@ -1230,6 +1230,7 @@ namespace CUDAWrappers
         LinearAlgebra::distributed::Vector<Number, MemorySpace::CUDA>
           ghosted_dst(ghosted_src);
         ghosted_src = src;
+        ghosted_dst = dst;
 
         // Execute the loop on the cells
         for (unsigned int i = 0; i < n_colors; ++i)