OverflowError

[hendriksuess]

Hi,

when trying to calculate persistent homology up to dimension 5 for a point cloud of size 2000 via

ripser_parallel(X, coeff=2, maxdim=5, n_threads=30,thresh=0.255)

I am getting the following error message

---------------------------------------------------------------------------
OverflowError                             Traceback (most recent call last)
Cell In[11], line 1
----> 1 ripser_parallel(X, coeff=2, maxdim=5, n_threads=30,thresh=0.255)
File ~/.local/lib/python3.10/site-packages/gph/python/ripser_interface.py:603, in ripser_parallel(X, maxdim, thresh, coeff, metric, metric_params, nearest_neighbors_params, weights, weight_params, collapse_edges, n_threads, return_generators)
    600         use_sparse_computer = False
    602 if use_sparse_computer:
--> 603     res = _compute_ph_vr_sparse(
    604         np.asarray(row, dtype=np.int64, order="C"),
    605         np.asarray(col, dtype=np.int64, order="C"),
    606         np.asarray(data, dtype=np.float32, order="C"),
    607         n_points,
    608         maxdim,
    609         thresh,
    610         coeff,
    611         n_threads,
    612         return_generators
    613         )
    614 else:
    615     # Only consider upper diagonal
    616     diagonal = np.diagonal(dm).astype(np.float32)
File ~/.local/lib/python3.10/site-packages/gph/python/ripser_interface.py:51, in _compute_ph_vr_sparse(I, J, V, N, maxHomDim, thresh, coeff, n_threads, return_generators)
     48 def _compute_ph_vr_sparse(I, J, V, N, maxHomDim, thresh=-1, coeff=2,
     49                           n_threads=1, return_generators=False):
     50     if coeff == 2:
---> 51         ret = gph_ripser.rips_dm_sparse(I, J, V, I.size, N, coeff,
     52                                         maxHomDim, thresh, n_threads,
     53                                         return_generators)
     54     else:
     55         ret = gph_ripser_coeff.rips_dm_sparse(I, J, V, I.size, N, coeff,
     56                                               maxHomDim, thresh, n_threads,
     57                                               return_generators)
OverflowError: simplex index 9239466805310005056 in filtration is larger than maximum index 36028797018963967

Is it possible to increase the maximum for the number of simplices?