 partial derivatives
of the likelihood with respect to all model variables) must also be done with
the same accuracy as the likelihood computation itself.Fournier developed a protocol for writing code to compute the required
derivatives based on the chain rule of differential calculus.This
protocol is very similar to the suite of methods that came to be known
as ``reverse mode automatic differentiation" 
.