tmp_index is shorter than to_index sometimes, especially in small and degenerate blockmodels about fastrg HOT 2 CLOSED

library(fastRG)
#> Loading required package: Matrix

set.seed(7)

n <- 10
k <- 10
pi <- rep(1, k) / k

B <- diag(rep(0.5, k))

sbm <- sbm(n = n, pi = pi, B = B)

sample_sparse(sbm)
#> Error in to[to_index] <- to_tmp[tmp_index]: replacement has length zero

> sample_sparse(sbm)
10 x 10 sparse Matrix of class "dsCMatrix"

[1,] . . . . . . . . . .
[2,] . . . . . . . . . .
[3,] . . . . . . . . . .
[4,] . . . . 3 1 . . . .
[5,] . . . 3 1 . . . . .
[6,] . . . 1 . . . . . .
[7,] . . . . . . . . . .
[8,] . . . . . . . . . .
[9,] . . . . . . . . . .
[10,] . . . . . . . . . .
Warning message:
  In to[to_index] <- to_tmp[tmp_index] :
  number of items to replace is not a multiple of replacement length