A structured diagonal Hessian approximation method with evaluation complexity analysis for nonlinear least squares