R/scaleConList.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"scaleConList" <-
  function(conList, specList) {
    maxDim <- max(unlist(lapply(conList, ncol)))
    listProd <- vector("list", maxDim)
    for (i in 1:length(conList)) {
      for (j in 1:maxDim) {
        if (j <= ncol(conList[[i]])) {
          dtemp <- abs(conList[[i]][, j] %*% t(specList[[i]][, j]))
          listProd[[j]] <- append(listProd[[j]], dtemp)
        }
      }
    }
    compMax <- unlist(lapply(listProd, max))
    maxall <- max(compMax)
    factScale <- sqrt(compMax / maxall)
    for (i in 1:length(conList)) {
      conListMax <- apply(conList[[i]], 2, max)
      conList[[i]] <- t(t(conList[[i]]) * (factScale / conListMax))
    }
    conList
  }