R/is.included.R

is.included <-
function(B, A, tol1=sqrt(.Machine$double.eps), tol2=sqrt(.Machine$double.eps)) {
  return(frob(B-A%*%ginv(A, tol=tol1)%*%B) / frob(B) < tol2)
}
marutter/lmreg documentation built on May 8, 2019, 1:51 p.m.