WH_2d_perf function

2D Whittaker-Henderson Smoothing (Maximum Likelihood, Generalized Fellner-Schall update)