R/test_singularity.R

In imabc: Incremental Mixture Approximate Bayesian Computation (IMABC)

test_singularity <- function(x, tol) {
  # Initialize return value to FALSE
  singular <- FALSE

  # Get eigen values
  e_value <- eigen(x, symmetric = TRUE, only.values = TRUE)$values
  if (any(e_value < tol)) { singular <- TRUE }

  return(singular)
}