Error: TridiagEigen: eigen decomposition failed

[alexpghayes]

Running into the following during some simulations, possibly related to #1.

library(RSpectra)

M <- new(
  "dgCMatrix",
  i = c(
    31L, 33L, 2L, 5L, 23L, 21L, 14L, 34L,
    20L, 39L, 11L, 1L, 22L, 14L, 19L, 7L, 6L, 2L, 20L, 24L, 28L,
    13L, 27L, 8L, 39L, 9L, 4L, 12L, 16L, 35L, 24L, 32L, 2L, 15L,
    3L, 36L, 18L, 2L, 38L, 26L, 29L, 10L, 3L, 2L, 0L, 30L, 37L, 25L,
    12L, 17L
  ), p = c(
    0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
    11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L,
    24L, 25L, 26L, 27L, 30L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L,
    40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L
  ), Dim = c(
    40L,
    47L
  ), Dimnames = list(NULL, NULL), x = c(
    1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1
  ), factors = list()
)

dim(M)
#> [1] 40 47
sum(M)
#> [1] 50
image(as.matrix(M))

M
#> 40 x 47 sparse Matrix of class "dgCMatrix"
#>                                                                                
#>  [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>  [2,] . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . .
#>  [3,] . . 1 . . . . . . . . . . . . . . 1 . . . . . . . . . . . 1 . . . . 1 . .
#>  [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . .
#>  [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . .
#>  [6,] . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>  [7,] . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . .
#>  [8,] . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . .
#>  [9,] . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . .
#> [10,] . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . .
#> [11,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [12,] . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [13,] . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . .
#> [14,] . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . .
#> [15,] . . . . . . 1 . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . .
#> [16,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . .
#> [17,] . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . .
#> [18,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [19,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . .
#> [20,] . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . .
#> [21,] . . . . . . . . 1 . . . . . . . . . 1 . . . . . . . . . . . . . . . . . .
#> [22,] . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [23,] . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . .
#> [24,] . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [25,] . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . 1 . . . . . . . .
#> [26,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [27,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
#> [28,] . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . .
#> [29,] . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . .
#> [30,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [31,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [32,] 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [33,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . .
#> [34,] . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [35,] . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [36,] . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . .
#> [37,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . .
#> [38,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [39,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 .
#> [40,] . . . . . . . . . 1 . . . . . . . . . . . . . . 1 . . . . . . . . . . . .
#>                          
#>  [1,] . . . . 1 . . . . .
#>  [2,] . . . . . . . . . .
#>  [3,] . . . 1 . . . . . .
#>  [4,] . . 1 . . . . . . .
#>  [5,] . . . . . . . . . .
#>  [6,] . . . . . . . . . .
#>  [7,] . . . . . . . . . .
#>  [8,] . . . . . . . . . .
#>  [9,] . . . . . . . . . .
#> [10,] . . . . . . . . . .
#> [11,] . 1 . . . . . . . .
#> [12,] . . . . . . . . . .
#> [13,] . . . . . . . . 1 .
#> [14,] . . . . . . . . . .
#> [15,] . . . . . . . . . .
#> [16,] . . . . . . . . . .
#> [17,] . . . . . . . . . .
#> [18,] . . . . . . . . . 1
#> [19,] . . . . . . . . . .
#> [20,] . . . . . . . . . .
#> [21,] . . . . . . . . . .
#> [22,] . . . . . . . . . .
#> [23,] . . . . . . . . . .
#> [24,] . . . . . . . . . .
#> [25,] . . . . . . . . . .
#> [26,] . . . . . . . 1 . .
#> [27,] . . . . . . . . . .
#> [28,] . . . . . . . . . .
#> [29,] . . . . . . . . . .
#> [30,] 1 . . . . . . . . .
#> [31,] . . . . . 1 . . . .
#> [32,] . . . . . . . . . .
#> [33,] . . . . . . . . . .
#> [34,] . . . . . . . . . .
#> [35,] . . . . . . . . . .
#> [36,] . . . . . . . . . .
#> [37,] . . . . . . . . . .
#> [38,] . . . . . . 1 . . .
#> [39,] . . . . . . . . . .
#> [40,] . . . . . . . . . .
svds(M, 2)
#> Error in fun(A, k, nu, nv, opts, mattype = "dgCMatrix"): TridiagEigen: eigen decomposition failed

^{Created on 2021-10-11 by the reprex package (v2.0.1.9000)}