Added the gradient for what is now called the DeltaLeaf solver

The next change will parameterize the LinearMixedModel by the solver type.
The solver types will be renamed as PLSDiag, PLSOne, and PLSGeneral.