Wiki examples for BLAS2 functions are added (kokkos#2122)

Some small additional change the the function headers themselves
to add some missing header file inclusions.

Applying clang-format

Removing constexpr since it won't happen before some work in Core.

blas/src/KokkosBlas2_ger.hpp

@@ -17,6 +17,8 @@
 #ifndef KOKKOSBLAS2_GER_HPP_
 #define KOKKOSBLAS2_GER_HPP_
 
+#include "KokkosKernels_helpers.hpp"
+
 #include <KokkosBlas2_ger_spec.hpp>
 
 namespace KokkosBlas {
@@ -42,15 +44,6 @@ template <class ExecutionSpace, class XViewType, class YViewType,
 void ger(const ExecutionSpace& space, const char trans[],
          const typename AViewType::const_value_type& alpha, const XViewType& x,
          const YViewType& y, const AViewType& A) {
-  static_assert(
-      Kokkos::SpaceAccessibility<typename AViewType::memory_space,
-                                 typename XViewType::memory_space>::assignable,
-      "AViewType memory space must be assignable from XViewType");
-  static_assert(
-      Kokkos::SpaceAccessibility<typename AViewType::memory_space,
-                                 typename YViewType::memory_space>::assignable,
-      "AViewType memory space must be assignable from YViewType");
-
   static_assert(
       Kokkos::SpaceAccessibility<ExecutionSpace,
                                  typename AViewType::memory_space>::accessible,

blas/src/KokkosBlas2_syr.hpp

@@ -17,6 +17,8 @@
 #ifndef KOKKOSBLAS2_SYR_HPP_
 #define KOKKOSBLAS2_SYR_HPP_
 
+#include "KokkosKernels_helpers.hpp"
+
 #include <KokkosBlas2_syr_spec.hpp>
 
 namespace KokkosBlas {
@@ -64,11 +66,6 @@ template <class ExecutionSpace, class XViewType, class AViewType>
 void syr(const ExecutionSpace& space, const char trans[], const char uplo[],
          const typename AViewType::const_value_type& alpha, const XViewType& x,
          const AViewType& A) {
-  static_assert(
-      Kokkos::SpaceAccessibility<typename AViewType::memory_space,
-                                 typename XViewType::memory_space>::assignable,
-      "AViewType memory space must be assignable from XViewType");
-
   static_assert(
       Kokkos::SpaceAccessibility<ExecutionSpace,
                                  typename AViewType::memory_space>::accessible,

blas/src/KokkosBlas2_syr2.hpp

@@ -17,6 +17,8 @@
 #ifndef KOKKOSBLAS2_SYR2_HPP_
 #define KOKKOSBLAS2_SYR2_HPP_
 
+#include "KokkosKernels_helpers.hpp"
+
 #include <KokkosBlas2_syr2_spec.hpp>
 #include <sstream>
 

example/wiki/CMakeLists.txt

@@ -1,2 +1,3 @@
+ADD_SUBDIRECTORY(blas)
 ADD_SUBDIRECTORY(sparse)
 ADD_SUBDIRECTORY(graph)

example/wiki/blas/CMakeLists.txt

@@ -0,0 +1,19 @@
+KOKKOSKERNELS_INCLUDE_DIRECTORIES(${CMAKE_CURRENT_BINARY_DIR})
+KOKKOSKERNELS_INCLUDE_DIRECTORIES(${CMAKE_CURRENT_SOURCE_DIR})
+
+KOKKOSKERNELS_INCLUDE_DIRECTORIES(${CMAKE_CURRENT_SOURCE_DIR}/../../../../test_common)
+
+KOKKOSKERNELS_ADD_EXECUTABLE_AND_TEST(
+  wiki_blas2_ger
+  SOURCES KokkosBlas2_wiki_ger.cpp
+  )
+
+KOKKOSKERNELS_ADD_EXECUTABLE_AND_TEST(
+  wiki_blas2_syr
+  SOURCES KokkosBlas2_wiki_syr.cpp
+  )
+
+KOKKOSKERNELS_ADD_EXECUTABLE_AND_TEST(
+  wiki_blas2_syr2
+  SOURCES KokkosBlas2_wiki_syr2.cpp
+  )

example/wiki/blas/KokkosBlas2_wiki_ger.cpp

@@ -0,0 +1,23 @@
+#include <Kokkos_Core.hpp>
+#include <KokkosBlas2_ger.hpp>
+
+int main(int argc, char* argv[]) {
+  Kokkos::initialize(argc, argv);
+  {
+    constexpr int M = 5;
+    constexpr int N = 4;
+
+    Kokkos::View<double**> A("A", M, N);
+    Kokkos::View<double*> x("X", M);
+    Kokkos::View<double*> y("Y", N);
+
+    Kokkos::deep_copy(A, 1.0);
+    Kokkos::deep_copy(x, 3.0);
+    Kokkos::deep_copy(y, 1.3);
+
+    const double alpha = Kokkos::ArithTraits<double>::one();
+
+    KokkosBlas::ger("T", alpha, x, y, A);
+  }
+  Kokkos::finalize();
+}

example/wiki/blas/KokkosBlas2_wiki_syr.cpp

@@ -0,0 +1,20 @@
+#include <Kokkos_Core.hpp>
+#include <KokkosBlas2_syr.hpp>
+
+int main(int argc, char* argv[]) {
+  Kokkos::initialize(argc, argv);
+  {
+    constexpr int M = 5;
+
+    Kokkos::View<double**> A("A", M, M);
+    Kokkos::View<double*> x("X", M);
+
+    Kokkos::deep_copy(A, 1.0);
+    Kokkos::deep_copy(x, 3.0);
+
+    const double alpha = double(1.0);
+
+    KokkosBlas::syr("T", "U", alpha, x, A);
+  }
+  Kokkos::finalize();
+}

example/wiki/blas/KokkosBlas2_wiki_syr2.cpp

@@ -0,0 +1,22 @@
+#include <Kokkos_Core.hpp>
+#include <KokkosBlas2_syr2.hpp>
+
+int main(int argc, char* argv[]) {
+  Kokkos::initialize(argc, argv);
+  {
+    constexpr int M = 5;
+
+    Kokkos::View<double**> A("A", M, M);
+    Kokkos::View<double*> x("X", M);
+    Kokkos::View<double*> y("Y", M);
+
+    Kokkos::deep_copy(A, 1.0);
+    Kokkos::deep_copy(x, 3.0);
+    Kokkos::deep_copy(y, 1.3);
+
+    const double alpha = double(1.0);
+
+    KokkosBlas::syr2("T", "U", alpha, x, y, A);
+  }
+  Kokkos::finalize();
+}