Add more fHS-WCA tests

This adds some more cpp tests for the fHS-WCA potential.
The implementation wrapped is compared to an independent basic cpp implementation.
A function to plot the fHS-WCA pair-potential is added to the interaction struct.
Also the energy is made strictly non-negative to avoid numerical problems.