Does Collider's performance depend on constitutive law?

Hi all,

as the title reveals we are struggling with the question of how the apllied contact law is affecting the computational time of the insertionSortCollider engine. We ran different simulations and made the observation that when using the Cundall model, the engine took about 4.09E10 (44.58%) ms to compute. Hertz-Mindlin only needed 3.32E09 (4.92%) ms. What is the reason for that? As far as I've understood, the engine creates bounding-boxes to detect mutual interactions of all body types. Shouldn't that work out regardless of the contact model? Documentation says (Page 18):

'[...] we want to compute bounds for all bodies in the simulation; suppose we want bound of type axis-aligned bounding box. Since the exact algorithm is different depending on particular shape, we need to provide functors for handling all specific cases'.

We haven't changed any of the bodies within different simualtions, so we would expect the algorithm to remain the same...
In our setup, we are currently working with ~ 45k spherical bodies and a specimen based on ~ 40k facets. If usefull, I can provide you with the code.

Thank you everyone in advance :)