multivator (version 1.1-4)

ipd: Positive definite matrices

Description

Is a matrix symmetric positive-definite?

Usage

ipd(mat)

Arguments

mat

A matrix

Value

Returns either TRUE or FALSE, printing a diagnostic message

Examples

data(mtoys)
stopifnot(ipd(crossprod(matrix(rnorm(30),10))))
stopifnot(ipd(M(toy_mhp)))