Provide a default implementation for VectorizedArray::dot_product().

diff --git a/include/deal.II/base/vectorization.h b/include/deal.II/base/vectorization.h
index 58ad2e09e76804e3a7dfc0ecc1191b72fccbe807..9511a853a1d4003e46f7e6a63daa0803b2999292 100644 (file)
--- a/include/deal.II/base/vectorization.h
+++ b/include/deal.II/base/vectorization.h
@@ -334,6 +334,25 @@ public:
     return VectorizedArrayIterator<const VectorizedArrayType>(
       static_cast<const VectorizedArrayType &>(*this), width);
   }
+
+  /**
+   * A default implementation for computing the dot product between two
+   * vectorized arrays. It first forms the lane-by-lane product between
+   * the current object and the argument, and then the across-lanes
+   * sum of the product. The function returns the kind of object that
+   * the VectorizedArray::sum() function returns.
+   *
+   * This function is inherited by all derived classes and provides
+   * the dot product to them unless they override it with their own
+   * implementation (presumably using a more efficient approach).
+   */
+  auto
+  dot_product(const VectorizedArrayType &v) const
+  {
+    VectorizedArrayType p = static_cast<const VectorizedArrayType &>(*this);
+    p *= v;
+    return p.sum();
+  }
 };