PRD implementation:
- Calculate particle-constant part of PRD in separate function
- Main calculation in calc_accel
- Print out particle radius and add it to the configuration_values struct
- Add CMake option for PRD