Create Hessian with one sided differencing on physical domain boundaries
Create new Hessian operator with one sided differencing (towards domain center) on physical domain boundaries. Central differencing in the inner cells should remain as currently defined in meta_hess()
.