gaussapprox: Compute a Gaussian approximation to a density
In jmhewitt/dsdive: Discrete Space Models for Dive Trajectories of Animals

Computes a Gaussian approximation to the target density logf, and additionally return methods to aid in sampling from the approximation. The approximation is computed numerically via the stats::optim routine.

gaussapprox(
  logf,
  init,
  method = "BFGS",
  optim.output = FALSE,
  lower = -Inf,
  upper = Inf,
  control = list(fnscale = -1)
)

`logf`	log of the density for which an approximation should be built.
`init`	initial guess for mode of `logf`
`method`	method to be used in `stats::optim`
`optim.output`	`TRUE` to return output from `stats::optim`
`lower`	lower bound value to pass to `stats::optim`
`upper`	upper bound value to pass to `stats::optim`
`control`	control list to pass to `stats::optim`