Implement CUDAWrappers::Vector::norm_sqr

diff --git a/include/deal.II/lac/cuda_vector.h b/include/deal.II/lac/cuda_vector.h
index 9c73d8d45d75cd43872dfd4a68b0db73cb70f54c..b6dfb083f024297c935856555d59e9272141d383 100644 (file)
--- a/include/deal.II/lac/cuda_vector.h
+++ b/include/deal.II/lac/cuda_vector.h
@@ -45,8 +45,9 @@ namespace LinearAlgebra
      *
      * @note Only float and double are supported.
      *
-     * @ingroup CUDAWrappers Vectors
-     * @author Karl Ljungkvist, Bruno Turcksin, 2016
+     * @see CUDAWrappers
+     * @ingroup Vectors
+     * @author Karl Ljungkvist, Bruno Turcksin, Daniel Arndt, 2016, 2018
      */
     template <typename Number>
     class Vector : public VectorSpaceVector<Number>
@@ -233,6 +234,12 @@ namespace LinearAlgebra
       virtual real_type
       l2_norm() const override;
 
+      /**
+       * Return the square of the $l_2$-norm.
+       */
+      real_type
+      norm_sqr() const;
+
       /**
        * Return the maximum norm of the vector (i.e., the maximum absolute
        * value among all entries and among all processors).

diff --git a/source/lac/cuda_vector.cu b/source/lac/cuda_vector.cu
index 9ced53b14ed3c41468b2b800ae8b9515cbf9ade5..0cb54b675b05559b700a3496320d8272eadb3432 100644 (file)
--- a/source/lac/cuda_vector.cu
+++ b/source/lac/cuda_vector.cu
@@ -597,7 +597,16 @@ namespace LinearAlgebra
     typename Vector<Number>::real_type
     Vector<Number>::l2_norm() const
     {
-      return std::sqrt((*this) * (*this));
+      return std::sqrt(norm_sqr());
+    }
+
+
+
+    template <typename Number>
+    typename Vector<Number>::real_type
+    Vector<Number>::norm_sqr() const
+    {
+      return (*this) * (*this);
     }