Matrix rows are already sorted in AMGCL, do not sort them twice !

BrunoLevy · Dec 16, 2024 · 9ea057a · 9ea057a
1 parent a42e3d3
commit 9ea057a
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 2 deletions.
diff --git a/src/lib/geogram/NL/nl_amgcl.cpp b/src/lib/geogram/NL/nl_amgcl.cpp
@@ -858,8 +858,6 @@ template <class Backend> NLboolean nlSolveAMGCL_generic() {
         size_t(n), (rowptr_t*)M->rowptr, (colind_t *)M->colind, M->val
     );
 
-    amgcl::backend::sort_rows(*M_amgcl);
-
     GEO::Stopwatch* Wbuild = new GEO::Stopwatch("AMGCL build", ctxt->verbose);
 
     if(nlExtensionIsInitialized_CUDA()) {

diff --git a/src/lib/geogram/NL/nl_cuda.c b/src/lib/geogram/NL/nl_cuda.c
@@ -1718,6 +1718,7 @@ static NLCUDASparseMatrix* CreateCUDASlicesFromCRSMatrixSlices(
 /**
  * \brief Gets the amount of memory required to store a matrix in CUDA
  * \param[in] M a pointer to a CRS matrix
+ * \param[in] with_buffer true if one needs to allocate a buffer for X and Y
  * \return the required GPU RAM to store \p M on the GPU
  */
 static size_t nlCUDAMatrixNeededMem(NLCRSMatrix* M, NLboolean with_buffer) {