Get Omega^-1 and derivatives
rxSymInvChol(
invObjOrMatrix,
theta = NULL,
type = "cholOmegaInv",
thetaNumber = 0L
)Object for inverse-type calculations. If this is a matrix,
setup the object for inversion by rxSymInvCholCreate with the default arguments and return
a reactive s3 object. Otherwise, use the inversion object to calculate the requested derivative/inverse.
Thetas to be used for calculation. If missing (NULL), a
special s3 class is created and returned to access Omega^1
objects as needed and cache them based on the theta that is
used.
The type of object. Currently the following types are supported:
cholOmegaInv gives the
Cholesky decomposition of the Omega Inverse matrix.
omegaInv gives the Omega Inverse matrix.
d(omegaInv) gives the d(Omega^-1) withe respect to the
theta parameter specified in thetaNumber.
d(D) gives the d(diagonal(Omega^-1)) with respect to
the theta parameter specified in the thetaNumber
parameter
For types d(omegaInv) and d(D),
the theta number that the derivative is taken against. This
must be positive from 1 to the number of thetas defining the
Omega matrix.
Matrix based on parameters or environment with all the matrixes calculated in variables omega, omegaInv, dOmega, dOmegaInv.