R/GetFullJacobian.r

In zarquon42b/RvtkStatismo: Integrating statismo and R using the vtkStandardMeshRepresenter

#' get Jacobian for all model points or a subset
#'
#' get Jacobian for all model points or a subset
#'
#' @param model statistical shape model of class 'pPCA'           
#' @param pt integer vector specifying which points to use
#' @return
#' returns a matrix with rows n: (n+2) containing the Jacobian for the n-th point
#' @seealso \code{\link{GetJacobian}}
#' @export
#' @rdname GetFullJacobian 
setGeneric("GetFullJacobian", function(model,pt) {
    standardGeneric("GetFullJacobian")
})

#' @rdname GetFullJacobian
setMethod("GetFullJacobian", signature(model="pPCA"), function(model,pt) {
    if (missing(pt))
        pt <- GetDomainPoints(model)
    else
        pt <- GetDomainPoints(model)[pt,]
    # lse
      
    out <- .Call("GetFullJacobian",model,pt)
    return(out)
})