grHs: Computes the grad and hessian of 'recl'.
In kohleth/spcreml: Spatial Modelling using Composite REML

Computes the grad and hessian of recl.

1 2	grHs(theta, dS, coordsMat, Y, Xmat, pointsj, cov.model = "exp", gr = TRUE, hs = TRUE, ...)

`theta`	named vector of covariance parameter. Must (at least) contain named elements sigma2 and phi. Can also contain nugget and kappa.
`dS`	a list of functions to compute the first derivative of Sigma (dSigma/dtheta).
`coordsMat`	a N by D coordinate matrix. D being the spatial dimension, typically 2 or 3.
`Y`	the response vector
`Xmat`	the design matrix (not the data.frame)
`cov.model`	character string that specifies the spatial covariance model. See `?geoR::cov.spatial`
`gr`	TRUE/FALSE. Whether the grad vector should be returned.
`hs`	TRUE/FALSE. Whether the Hessian matrix should be returned.
`poinstj`	A list of row numbers that indicates which row belongs to block j (block j is itself a union of block k and block l). See Details.