Comment 1 for bug 1112763

If I define the dimensions of the periodic cell first, and then the radius of the big sphere as a multiple of the width of the periodic cell (which makes more sense, as I want the big body to be bigger than the periodic cell), I get something that I don't understand.

Using the script in the previous link with following lines to define the geometry:

# Periodic cell
cell0 = 1
cell1 = 10
cell2 = 10
# Spheres
DbigSphere = cell0*1.
radius_bigSphere= DbigSphere/2.
radius_sph = 1./10.

and varying " DbigSphere = cell0*NUMBER"
if NUMBER = 1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 3, 3.1 to 3.7 the interaction is not detected
if NUMBER = 1.8, 1.9, 2, 2.1 to 2.9, 3.8, 3.9, 4., 4.1 the interaction is detected
I am totally lost now. Plus, geometries different from this one give different values for which the interaction is or is not detected. So the value that works for one simulation, does not necessarily work for a different one.

This is a problem when I save and then reload my simulation with the body bigger than the period, where the big body is already in contact with some of the small bodies.
In that case some interactions are not detected (probably because their bounding boxes were exactly touching each other at the step I reloaded the test), and I cannot predict it. There are too many variables (such as width of the period, dimension of the big body, relative position of the big body with the small ones) that change and I do not have a clear understanding of how they are related to this bug.