from gudhi.sklearn.cubical_persistence import CubicalPersistence
cp = CubicalPersistence(homology_dimensions=[0])
cp.fit_transform([1])
segfaults. This can easily happen if you have a 1d array X and call fit_transform(X) instead of fit_transform([X]).