Update fitting_surface_pdm.cpp

Thomas, you forgot to switch to the sparse representation of the LHS vector in a few (last 5) add* functions.