R/calcB.R

"calcB" <-
  function(k) {
    lenk <- length(k)
    m <- rep(0, lenk * lenk)
    m <- as.matrix(.C("r_calcB",
      m = as.double(m), as.double(k),
      as.integer(lenk), PACKAGE = "TIMP"
    )$m)
    dim(m) <- c(lenk, lenk)
    m
  }
glotaran/TIMP documentation built on Feb. 11, 2023, 4:11 a.m.