R/getAnisotropy.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"getAnisotropy" <- function(model) {
  ani <- model@anispec
  if (length(ani$angle) > 0) {
    # non-associative, or associative w/same func. per comp
    if (length(ani$anifunc) == 1) {
      ani$anifunc <- rep(
        ani$anifunc,
        max(model@ncolc)
      )
    }
    model@ncolc <- model@ncolc + length(ani$super)
    ani$calcani <- TRUE
    if (length(ani$useparperp) == 0) {
      ani$useparperp <- FALSE
    }
  } else {
    ani$calcani <- FALSE
    ani$useparperp <- FALSE
  }
  if (length(ani$super) > 0) {
    model@cohcol <- model@cohcol + length(ani$super)
  }
  model@anispec <- ani
  model
}