famfun <- function(fam, deriv, theta) {
    setfam(list(fam))
    result <- .C(C_aster_family,
        index = as.integer(1),
        deriv = as.integer(deriv),
        theta = as.double(theta),
        value = double(1))$value
    clearfam()
    return(result)
}