Currently, this package computes the partial derivative matrix A = drdy (residual, y, x, p) as part of the forward pass for use in the reverse pass. This approach appears to yield incorrect results ...