reverted grad_surf and slope_error to the original version with the case distinction between 'bp' and 'op'. This caused problems (funny stripes, the rms epsilon did not decrease).

modified original optimize_surface to be able to compare it to Paul's version

corrected ns_3d

minor changes in optimize_surface_exact and scale_fac