Fixed a small bug
Switched to 2nd order derivatives
Tweaked boundary condition a bit
Commented out debug message
Added cuda_runtime to avoid flycheck warnings
Reduced precision to allow for 2nd order derivative