swh:1:snp:81eadaa089e8253d8469bcef66aa332632c6c669
Tip revision: 3f1c6b45c918375fc2b65bb2a7d6e3590979cb61 authored by HwB on 12 December 2011, 00:00:00 UTC
version 0.9.1
version 0.9.1
Tip revision: 3f1c6b4
inv.R
##
## i n v . R matrix inverse
##
inv <- function(a) {
if (length(a) == 0) return(matrix(0, nrow=0, ncol=0))
if ((!is.numeric(a) && !is.complex(a)) || !is.matrix(a))
stop("Argument 'a' must be a numeric or complex matrix.")
if (nrow(a) != ncol(a))
stop("Matrix 'a' must be square.")
e <- try(b <- solve(a), silent=TRUE)
if (class(e) == "try-error") {
warning("Matrix appears to be singular.")
b <- rep(Inf, length(a))
dim(b) <- dim(a)
}
return(b)
}