Test periodic bc

This tests the periodic boundary conditions which are generated in space_getsid() and used in runner_dopair1_density.

An 8x8x8 mesh is generated that is periodic and 26 pair interactions are computed for the following cells:

• All 8 corner cells
• The cell at the centre of each face
• The cell at the centre of each edge

These checks are performed for both unperturbed and perturbed cells.

Added 24 commits:

• e23dc2dd...bd229c7a - 18 commits from branch master
• 6842e754 - Also use the 'nearest()' function call in the brute-force version of the force pairs.
• 2545e4ab - Merge branch 'master' into testPeriodicBC
• a80b5fc5 - Make the difffloat script use a new line in the files to set a limit below which…
• c1a698ca - Also update the other tests' tolerance files to include the new lower limit lines.
• e15d6dbf - Also use the DOPAIR branching function for gizmo's extra loop
• 9c30c862 - Code formatting

Title changed from [WIP] Test periodic bc to Test periodic bc

Added 1 commit:

• 2f7f2aea - Corrected mistake in the periodic boundary test that prevented compilation when …

@jwillis: are you sure the test is comprehensive ? My recollection is that the python script only checks the first n^3 entries so 216 here and not what comes next. So only one of the directions is effectively being checked here.

Ah, I thought that difffloat.py would check the whole file?

So I would just need to change the number of lines to check in the script? From n^3 to n^3 * 26, as there are 26 cells to test?

But that would break all the other tests though, I'll have a think

Added 1 commit:

• 473ba476 - Check all particles in output files when testing the periodic boundary conditions.

Will that work?

Possibly. I'll try it out. Might change the code a bit to analyse each direction separately.

Added 1 commit:

• b2692841 - Output the position of the cell in the output of the periodic boundary condition…

Here is a version that works based on your changes. I am still unhappy with the low level of accuracy but will look into this at a later stage I think.

Added 1 commit:

• 5fbf74b3 - Code formatting.

Added 1 commit:

• 7dfb9465 - The test should also work in the non-optimized case.

Status changed to merged

mentioned in commit c9d95e88