109   if (diagonal_.rows() != n) {  in ComputeStep()
110     diagonal_.resize(n, 1);  in ComputeStep()
121   jacobian->SquaredColumnNorm(diagonal_.data());  in ComputeStep()
123     diagonal_[i] = min(max(diagonal_[i], min_diagonal_), max_diagonal_);  in ComputeStep()
125   diagonal_ = diagonal_.array().sqrt();  in ComputeStep()
176   gradient_.array() /= diagonal_.array();  in ComputeGradient()
188       (gradient_.array() / diagonal_.array()).matrix();  in ComputeCauchyPoint()
208     dogleg_step.array() /= diagonal_.array();  in ComputeTraditionalDoglegStep()
220     dogleg_step.array() /= diagonal_.array();  in ComputeTraditionalDoglegStep()
251   dogleg_step.array() /= diagonal_.array();  in ComputeTraditionalDoglegStep()
282     dogleg_step.array() /= diagonal_.array();  in ComputeSubspaceDoglegStep()
312     dogleg_step.array() /= diagonal_.array();  in ComputeSubspaceDoglegStep()
362   dogleg_step.array() /= diagonal_.array();  in ComputeSubspaceDoglegStep()
558     lm_diagonal_ = diagonal_ * std::sqrt(mu_);  in ComputeGaussNewtonStep()
607     gauss_newton_step_.array() *= -diagonal_.array();  in ComputeGaussNewtonStep()
706   tmp = (subspace_basis_.col(0).array() / diagonal_.array()).matrix();  in ComputeSubspaceModel()
708   tmp = (subspace_basis_.col(1).array() / diagonal_.array()).matrix();  in ComputeSubspaceModel()