R/fullKF.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

"fullKF" <-
  function(theta, kinscal, kmat, jvec, fixedkmat = FALSE,
           kinscalspecial = list(), kinscalspecialspec = list(),
           nocolsums = FALSE) {
    # fullKF

    dimk <- nrow(kmat)
    A <- matrix(nrow = dimk, ncol = dimk)
    K <- fillK(
      theta, kinscal, kmat, fixedkmat, kinscalspecial,
      kinscalspecialspec, nocolsums
    )
    eigenlijk <- eigen(K, only.values = FALSE)
    V <- eigenlijk$vectors
    gamma <- solve(V) %*% jvec
    for (j in 1:dimk) {
      A[j, ] <- V[, j] * gamma[j]
    }
    list(A = A, values = eigenlijk$values)
  }