Fix a bug in a CUDA vector constructor

The constructor did not set the values in the vector to zero

diff --git a/include/deal.II/lac/cuda_vector.h b/include/deal.II/lac/cuda_vector.h
index aa407fff13f420aabaadc5ea85c3181b560aed9f..488138ccfbdc5d4fba7127e4270817b368212b56 100644 (file)
--- a/include/deal.II/lac/cuda_vector.h
+++ b/include/deal.II/lac/cuda_vector.h
@@ -27,7 +27,6 @@
 DEAL_II_NAMESPACE_OPEN
 
 class CommunicationPatternBase;
-class IndexSet;
 template <typename Number> class ReadWriteVector;
 
 namespace LinearAlgebra

diff --git a/source/lac/cuda_vector.cu b/source/lac/cuda_vector.cu
index 41501d896889d5f4486d284f7b43b13c0cde469b..36b2be442b29024a8d1409e1ae6da7270dd9e4a6 100644 (file)
--- a/source/lac/cuda_vector.cu
+++ b/source/lac/cuda_vector.cu
@@ -503,11 +503,10 @@ namespace LinearAlgebra
     template <typename Number>
     Vector<Number>::Vector(const size_type n)
       :
-      n_elements(n)
+      val(nullptr),
+      n_elements(0)
     {
-      // Allocate the memory
-      cudaError_t error_code = cudaMalloc(&val, n_elements*sizeof(Number));
-      AssertCuda(error_code);
+      reinit(n, false);
     }
 
 
@@ -542,7 +541,7 @@ namespace LinearAlgebra
         }
       else
         {
-          if (n_elements != n)
+          if ((n_elements != n) && (val != nullptr))
             {
               cudaError_t error_code = cudaFree(val);
               AssertCuda(error_code);