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.