The internal helper validate_arr_cub() ensures that a numeric array has dimensions that are acceptable to the cubical Ripser source code. However, these constraints in the source code ($2000 \times 1000$ for dim2, $512^3$ for dim3, and $64^4$ for dim4) seem arbitrary. It may be appropriate to warn the user before proceeding with the computation on a large array, but i believe these hard constraints should be lifted. This can be done easily in the dim2 case and with caution in the dim3 and dim4 cases.
@rrrlw does this seem reasonable to you? I'd also like to contact the authors to ask their reasons.
The internal helper$2000 \times 1000$ for $512^3$ for $64^4$ for
validate_arr_cub()ensures that a numeric array has dimensions that are acceptable to the cubical Ripser source code. However, these constraints in the source code (dim2,dim3, anddim4) seem arbitrary. It may be appropriate to warn the user before proceeding with the computation on a large array, but i believe these hard constraints should be lifted. This can be done easily in thedim2case and with caution in thedim3anddim4cases.@rrrlw does this seem reasonable to you? I'd also like to contact the authors to ask their reasons.