Fails with smoothing radius <0.4

[NH89]

Noticed by Shafeeq: The program fails when reducing smoothing radius below 0.4 in my computer.

[NH89]

I will investigate this, once I have the C++ branch working. My first suspicion is that it may be 'correct' behaviour - i.e. if particles are created wider apart than the outer limit of their smoothing radius, then they are separate (like grains of sand) and should move like a granular fluid.

[eshafeeqe]

Dear Nick,

I am not telling the mathematical or theoretical problem here.

The system is crashing while smoothing radius <0.4. and crashing is happening in kernel insert_particles_in_bin_array

The code in the kernel given below,

kernel void insert_particles_in_bin_array (PSO_ARGS, global uint * backup_prefix_sum) {
    USE_FIELD_FIRST_VALUE(n, uint) USE_FIELD(gridcell, uint) USE_FIELD(cellparticles, uint)

    size_t i = get_global_id(0);

    if (i < n)
    {
        
        cellparticles[atomic_inc(backup_prefix_sum + gridcell[i])] = i;
    }

}

I am suspecting that there could be some memory access violation happening here if the smoothing radius is less than 0.4.