[REVIEW] Spectral clustering functionality in RAFT.

CHANGELOG.md

@@ -1,6 +1,7 @@
 # RAFT 0.15.0 (Date TBD)
 
 ## New Features
+- PR #12: Spectral clustering.
 - PR #7: Migrating cuml comms -> raft comms_t
 - PR #15: add exception based error handling macros
 

cpp/CMakeLists.txt

@@ -204,6 +204,7 @@ set(RAFT_LINK_LIBRARIES
   ${CUDA_cusolver_LIBRARY}
   ${CUDA_CUDART_LIBRARY}
   ${CUDA_cusparse_LIBRARY}
+  ${CUDA_curand_LIBRARY}
   rmm)
 
 set(RAFT_LINK_DIRECTORIES "${RMM_LIBRARY}")
@@ -224,7 +225,10 @@ if(BUILD_RAFT_TESTS)
     test/handle.cpp
     test/mr/device/buffer.cpp
     test/mr/host/buffer.cpp
-    test/test.cpp)
+    test/test.cpp
+    test/spectral_matrix.cu
+    test/eigen_solvers.cu
+    test/cluster_solvers.cu)
 
   target_include_directories(test_raft
     PRIVATE

cpp/include/raft/linalg/cublas_wrappers.h

@@ -578,5 +578,51 @@ inline cublasStatus_t cublasdot(cublasHandle_t handle, int n, const double *x,
 }
 /** @} */
 
+/**
+ * @defgroup setpointermode cublas set pointer mode method
+ * @{
+ */
+// no T dependency...
+// template <typename T>
+// cublasStatus_t cublassetpointermode(  // NOLINT
+//                                         cublasHandle_t  handle,
+//                                         cublasPointerMode_t mode,
+//                                         cudaStream_t stream);
+
+// template<>
+inline cublasStatus_t cublassetpointermode(cublasHandle_t handle,
+                                           cublasPointerMode_t mode,
+                                           cudaStream_t stream) {
+  CUBLAS_CHECK(cublasSetStream(handle, stream));
+  return cublasSetPointerMode(handle, mode);
+}
+/** @} */
+
+/**
+ * @defgroup scal cublas dot calls
+ * @{
+ */
+template <typename T>
+cublasStatus_t cublasscal(cublasHandle_t handle, int n, const T *alpha, T *x,
+                          int incx, cudaStream_t stream);
+
+template <>
+inline cublasStatus_t cublasscal(cublasHandle_t handle, int n,
+                                 const float *alpha, float *x, int incx,
+                                 cudaStream_t stream) {
+  CUBLAS_CHECK(cublasSetStream(handle, stream));
+  return cublasSscal(handle, n, alpha, x, incx);
+}
+
+template <>
+inline cublasStatus_t cublasscal(cublasHandle_t handle, int n,
+                                 const double *alpha, double *x, int incx,
+                                 cudaStream_t stream) {
+  CUBLAS_CHECK(cublasSetStream(handle, stream));
+  return cublasDscal(handle, n, alpha, x, incx);
+}
+
+/** @} */
+
 }  // namespace linalg
 }  // namespace raft