Infrastructure for Kolmogorov-Smirnov tests, (tested!) methods for sampling with replacement, and changes to include guards

RandBLAS/base.hh

@@ -27,8 +27,7 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_base_hh
-#define randblas_base_hh
+#pragma once
 
 /// @file
 
@@ -54,9 +53,18 @@
 /// code common across the project
 namespace RandBLAS {
 
+/**
+ * Stores stride information for a matrix represented as a buffer.
+ * The intended semantics for a buffer "A" and the conceptualized
+ * matrix "mat(A)" are 
+ * 
+ *  mat(A)[i, j] == A[i * inter_row_stride + j * inter_col_stride].
+ * 
+ * for all (i, j) within the bounds of mat(A).
+ */
 struct stride_64t {
-    int64_t inter_row_stride;
-    int64_t inter_col_stride;
+    int64_t inter_row_stride; // step down a column
+    int64_t inter_col_stride; // step along a row
 };
 
 inline stride_64t layout_to_strides(blas::Layout layout, int64_t ldim) {
@@ -212,5 +220,3 @@ std::ostream &operator<<(
 }
 
 } // end namespace RandBLAS::base
-
-#endif

RandBLAS/config.h.in

@@ -1,5 +1,4 @@
-#ifndef RandBLAS_config_h
-#define RandBLAS_config_h
+#pragma once
 
 #define RandBLAS_FULL_VERSION "@RandBLAS_FULL_VERSION@"
 #define RandBLAS_VERSION_MAJOR @RandBLAS_VERSION_MAJOR@
@@ -53,5 +52,3 @@
 //
 //   if you are linking to OpenMP.
 //
-
-#endif

RandBLAS/dense_skops.hh

@@ -26,9 +26,7 @@
 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 // POSSIBILITY OF SUCH DAMAGE.
 //
-
-#ifndef randblas_dense_hh
-#define randblas_dense_hh
+#pragma once
 
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
@@ -640,5 +638,3 @@ RNGState<RNG> fill_dense(
 }
 
 }  // end namespace RandBLAS
-
-#endif

RandBLAS/exceptions.hh

@@ -29,9 +29,6 @@
 
 #pragma once
 
-#ifndef RandBLAS_EXCEPTIONS_HH
-#define RandBLAS_EXCEPTIONS_HH
-
 #include <exception>
 #include <cstdarg>
 #include <string>
@@ -164,5 +161,3 @@ inline void abort_if( bool cond, const char* func,  const char* format, ... ) {
 #endif
 
 } // namespace RandBLAS::exceptions
-
-#endif

RandBLAS/random_gen.hh

@@ -27,8 +27,7 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef random_gen_hh
-#define random_gen_hh
+#pragma once
 
 /// @file
 
@@ -178,4 +177,3 @@ struct uneg11
 
 } // end of namespace r123ext
 
-#endif

RandBLAS/skge.hh

@@ -26,9 +26,7 @@
 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 // POSSIBILITY OF SUCH DAMAGE.
 //
-
-#ifndef randblas_skge_hh
-#define randblas_skge_hh
+#pragma once
 
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
@@ -1231,6 +1229,3 @@ inline void sketch_general(
 };
 
 }  // end namespace RandBLAS
-
-
-#endif

RandBLAS/sksy.hh

@@ -27,8 +27,7 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sksy_hh
-#define randblas_sksy_hh
+#pragma once
 
 #include "RandBLAS/util.hh"
 #include "RandBLAS/base.hh"
@@ -538,4 +537,3 @@ inline void sketch_symmetric(
 }
 
 } // end namespace RandBLAS
-#endif

RandBLAS/skve.hh

@@ -26,9 +26,7 @@
 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 // POSSIBILITY OF SUCH DAMAGE.
 //
-
-#ifndef randblas_skve_hh
-#define randblas_skve_hh
+#pragma once
 
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
@@ -260,4 +258,3 @@ inline void sketch_vector(
 }
 
 }  // end namespace RandBLAS
-#endif

RandBLAS/sparse_data/base.hh

@@ -26,9 +26,7 @@
 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 // POSSIBILITY OF SUCH DAMAGE.
 //
-
-#ifndef randblas_sparse_data_hh
-#define randblas_sparse_data_hh
+#pragma once
 
 #include "RandBLAS/config.h"
 #include "RandBLAS/base.hh"
@@ -193,7 +191,3 @@ namespace RandBLAS {
     using RandBLAS::sparse_data::IndexBase;
     using RandBLAS::sparse_data::SparseMatrix;
 }
-
-
-
-#endif

RandBLAS/sparse_data/conversions.hh

@@ -26,9 +26,8 @@
 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 // POSSIBILITY OF SUCH DAMAGE.
 //
+#pragma once
 
-#ifndef randblas_sparse_data_conversions
-#define randblas_sparse_data_conversions
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -209,5 +208,3 @@ void reindex_inplace(COOMatrix<T> &A, IndexBase desired) {
 }
 
 } // end namespace RandBLAS::sparse_data::conversions
-
-#endif

RandBLAS/sparse_data/coo_matrix.hh

@@ -27,8 +27,8 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_coo
-#define randblas_sparse_data_coo
+#pragma once
+
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -409,4 +409,3 @@ void coo_to_dense(const COOMatrix<T> &spmat, Layout layout, T *mat) {
 
 } // end namespace RandBLAS::sparse_data::coo
 
-#endif

RandBLAS/sparse_data/coo_spmm_impl.hh

@@ -27,8 +27,8 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_coo_multiply
-#define randblas_sparse_data_coo_multiply
+#pragma once
+
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -163,5 +163,3 @@ static void apply_coo_left_jki_p11(
 
 
 } // end namespace
-
-#endif

RandBLAS/sparse_data/csc_matrix.hh

@@ -27,8 +27,8 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_csc
-#define randblas_sparse_data_csc
+#pragma once
+
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -246,7 +246,4 @@ void dense_to_csc(Layout layout, T* mat, T abs_tol, CSCMatrix<T> &spmat) {
     return;
 }
 
-
-}
-
-#endif
+} // end namespace RandBLAS::sparse_data::csc

RandBLAS/sparse_data/csc_spmm_impl.hh

@@ -27,8 +27,7 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_csc_multiply
-#define randblas_sparse_data_csc_multiply
+#pragma once
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -213,5 +212,4 @@ static void apply_csc_left_kib_rowmajor_1p1(
     return;
 }
 
-}
-#endif
+} // end namespace RandBLAS::sparse_data::csc

RandBLAS/sparse_data/csr_matrix.hh

@@ -27,8 +27,8 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_csr
-#define randblas_sparse_data_csr
+#pragma once
+
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -262,4 +262,3 @@ void dense_to_csr(Layout layout, T* mat, T abs_tol, CSRMatrix<T> &spmat) {
 
 
 } // end namespace RandBLAS::sparse_data::csr
-#endif

RandBLAS/sparse_data/csr_spmm_impl.hh

@@ -27,8 +27,8 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_csr_multiply
-#define randblas_sparse_data_csr_multiply
+#pragma once
+
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -156,5 +156,3 @@ static void apply_csr_left_ikb_rowmajor(
 }
 
 } // end namespace RandBLAS::sparse_data::csr
-
-#endif

RandBLAS/sparse_data/sksp.hh

@@ -27,8 +27,7 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sksp_hh
-#define randblas_sksp_hh
+#pragma once
 
 #include "RandBLAS/base.hh"
 #include "RandBLAS/dense_skops.hh"
@@ -625,6 +624,3 @@ inline void sketch_sparse(
 }
 
 }  // end namespace RandBLAS
-
-
-#endif

RandBLAS/sparse_data/spmm_dispatch.hh

@@ -27,8 +27,8 @@
 // POSSIBILITY OF SUCH DAMAGE.
 //
 
-#ifndef randblas_sparse_data_spmm_dispatch
-#define randblas_sparse_data_spmm_dispatch
+#pragma once
+
 #include "RandBLAS/base.hh"
 #include "RandBLAS/exceptions.hh"
 #include "RandBLAS/sparse_data/base.hh"
@@ -384,5 +384,3 @@ inline void spmm(blas::Layout layout, blas::Op opA, blas::Op opB, int64_t m, int
 }
 
 }
-
-#endif

RandBLAS/sparse_skops.hh

@@ -26,9 +26,7 @@
 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 // POSSIBILITY OF SUCH DAMAGE.
 //
-
-#ifndef randblas_sparse_skops_hh
-#define randblas_sparse_skops_hh
+#pragma once
 
 #include "RandBLAS/config.h"
 #include "RandBLAS/base.hh"
@@ -483,5 +481,3 @@ static auto transpose(SKOP const &S) {
 }
 
 } // end namespace RandBLAS::sparse
-
-#endif