Created a branching function for DOPAIR2 to only call the vectorised function when vectorisation is enabled and call the serial version when corner cells interact. Moved some debug statements inside the branching functions.

Pad particle caches with sensible values for the smoothing length to prevent floating point exceptions.

runner_dopair2_force_vec now runs correctly without looping over all particles. max_index have been populated using the max between hi and hj_max.

pairs_all_force now only updates active particles. Removed error message in shuffle_particles for arrays of length 1.

Fixed bug where max_index_i could be count_j instead of count_j - 1. Also fixed bug where an exit iteration of 0 prevents the inner loop from executing. For dopair2.

Added a test which interacts a pair of cells with only one particle each. Disabled error in shuffle_particles.

Fixed bug where max_index_i could be count_j instead of count_j - 1. Also fixed bug where an exit iteration of 0 prevents the inner loop from executing.

Take the maximum between hi and hj_max when working out which particles to skip. This ensures all interactions are found. Sedov blast now runs correctly with runner_dopair2_force_vec.

Pad cache with better fake particles so that they won't be in range of real particles. Sedov blast now runs to completion but timesteps don't match serial output.

Also send dx_max_part and dx_max_gpart alongside ti_end_min at the end of a time-step when running with MPI

Implements the solution to #367. Instead of just sending the time-step we send the distances moved by any part and gpart in that cell. 

I have also done some cleaning up around the packing and unpacking and improved the documentation.

See merge request !434

Correct vec neighbour find

Fixes neighbour find bug in the vectorised version of the density interactions.

See merge request !435

Added an extra test case to testActivePair.sh to include the RNG seed that creates a vectorised neighbour find bug.

Also send dx_max_part and dx_max_gpart alongside ti_end_min at the end of a time-step when running with MPI.

Fixes to the gravity code

This fixes #363. The long-range gravity task has to consider all the cells that were assigned to it at the rebuild time. 

Additional improvements:
 - Drift the multipoles while activating the tasks,
 - Avoid FPEs in the FFT task,
 - Adjust the number of tasks per top-level required to run.

See merge request !430

Small cosmetic improvements to the documentation of the FMM code.



See merge request !432

Added debugging checks to verify that the particles are in the correct frame after the shift in DOPAIR1 and DOPAIR2. The check takes into account directions and drift since last rebuild.

Corrected the same check in the vectorized versions of the pair interactions. This was not using the correct criterion for the maximal displacement any particle can have had since the last rebuild.

See merge request !431