Computation of the gradient
First we note that:
Since, (S-1)T = S which follows from S being a unitary transform.
Variation of H(X) due to a variation dX = B1/2dx of the model state
where X0 = Xb + Gx0 . H’(X0) is the Jacobian matrix of H evaluated at X0.