The scalar λk, i.e. damping factor, controls both the magnitude and
direction of pk at each major iteration. If the reduction of F is rapid, a
smaller λk can be used, but if an iteration gives insufficient reduction
of F, λk can be increased and pk tends to move towards a more steep
descent direction. In this work the initial damping factor λ0 is given as
0.001.