R/doClpConstr.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"doClpConstr" <-
  function (X, clp_ind, clpCon, clpequ, num_clpequ, usecompnames0, 
            usecompnamesequ)  
  {   
    
    #add equality constraints
    if(num_clpequ > 0) {
      for(i in 1:num_clpequ) {
        X[, clpCon$clpMod[clp_ind, i] ] <- 
          X[, clpCon$clpMod[clp_ind, i] ] + 
          X[, clpCon$clpRem[clp_ind,i] ] * clpequ[i]
      }
    }
    ## column indices that are removed because of 0 constraints
    if(usecompnames0)
      rem0 <- colnames(X)[ which(clpCon$clp0mat[clp_ind, ] > 0)]
    else 
      rem0 <- which(clpCon$clp0mat[clp_ind, ] > 0)
    ## column indices that are removed because of equ constraints
    if(usecompnamesequ)
      remE <- colnames(X)[clpCon$clpRem[clp_ind, ]]
    else 
      remE <- clpCon$clpRem[clp_ind, ]
    rem <- append(rem0, remE) 
    
    if(sum(rem) > 0 )
      X <- as.matrix(X[, - rem])
    
    
    X
  }