Very interesting article. As a follow-on, starting with version 13.1 of SAS/IML there is a new built-in function LOGABSDET, which computes the log of the absolute value of the determinant of an arbitrary matrix. You can use this function when the magnitude of a determinant is too large to fit on your machine.