R/gpGradient.R

In alkalait/gptk: Gaussian Processes Tool-Kit

gpGradient <- function(params, model) {

  model = gpExpandParam(model, params)
  ## only necessary parts from gpLogLikeGradients implemented!

#   if (model$optimiseBeta && !model$fixInducing)
#     browser()

  g = - gpLogLikeGradients(model)

  return (g)
}