R/msnCP.dev.grad.out.prod.R

In MAINT.Data: Model and Analyse Interval Data

msnCP.dev.grad.out.prd <- function(param, y, grpind, Config, n=ifelse(is.matrix(y),nrow(y),1),
  p=ifelse(is.matrix(y),ncol(y),length(y)), k=1, inoptm=TRUE, 		     	
  limlnk2=log(1e8), trace=FALSE, c2tol=1e-3, ldRtol=-500, beta0tol=1e-6, ...)
{
  if (!all(is.finite(param))) {
    if (inoptm)  {
      return(rep(0.,length(param)))
    } else {
      return(NULL)
    }
  }
  if (!is.matrix(y)) { dim(y) <- c(1,p) }

  grad <- .Call( "msnCP_grad_out_prod", param, y, grpind, Config, n, p, k, 
    limlnk2, trace, c2tol, ldRtol, beta0tol, .Machine$double.eps, PACKAGE = "MAINT.Data" )

  if (!all(is.finite(grad))) { 
    if (inoptm)  {
      return(rep(0.,length(param)))
    } else {
      return(NULL)
    }
  }

  grad
}