Handled similar edge cases that might lead to nullpointer issues.

Author: pcarbo

DESCRIPTION

@@ -1,7 +1,7 @@
 Encoding: UTF-8
 Type: Package
 Package: fastTopics
-Version: 0.7-35
+Version: 0.7-36
 Date: 2025-12-21
 Title: Fast Algorithms for Fitting Topic Models and Non-Negative
     Matrix Factorizations to Count Data

R/fit_poisson_nmf.R

@@ -455,7 +455,7 @@ fit_poisson_nmf <- function (X, k, fit0, numiter = 100,
     cat(sprintf("Running at most %d %s updates, %s extrapolation ",
                 numiter,method.text,
                 ifelse(control$extrapolate,"with","without")))
-    cat("(fastTopics 0.7-35).\n")
+    cat("(fastTopics 0.7-36).\n")
   }
 
   # INITIALIZE ESTIMATES

src/ccd.cpp

@@ -29,7 +29,9 @@ inline void ccd_update_factor_sparse (const sp_mat& V, const mat& W,
 				      unsigned int j, double e) {
   vec          v = nonzeros(V.col(j));
   unsigned int n = v.n_elem;
-  uvec         i(n);
+  if (n == 0)
+    return;
+  uvec i(n);
   getcolnonzeros(V,i,j);
   H.col(j) = ccd_kl_update(W.rows(i),sumw,v,H.col(j),e);
 }

src/pnmfem.cpp

@@ -25,7 +25,9 @@ inline vec pnmfem_update_factor_sparse (const sp_mat& X, const mat& F,
   vec          x = nonzeros(X.col(j));
   unsigned int n = x.n_elem;
   vec          f = F.col(j);
-  uvec         i(n);
+  if (n == 0)
+    return;
+  uvec i(n);
   getcolnonzeros(X,i,j);
   poismixem(L1,u,x,i,f,numiter);
   return f;