Added the 4th order boundary conditions
Working CPU version
Reenable GPU divergence calculation
Working boundary routine, but no divergence free field yet...