the true Hessian of the objective function, while the method based on the TD error term uses Gauss-Newton iterations and the PD error based optimization is a variant of the steepest descent method. The extension of SDM to surface fitting will also be discussed.