Implements the second-order central difference approximation to the second-order space derivative.
For the boundary points fourth-order biased approximation is used.
Release Notes: