I have a determinant which I know is a square of an integer, but because its value is bigger than .Machine$integer.max
I used the mpfr package.
But I still have problems.
Here is the algorithm:
> a<- mpfr(sqrt(det(M)), precBits=512);a
1 'mpfr' number of precision 512 bits
[1] 430080000000001.1875
Could you please help me?
determinant()
function packaged with theRmpfr
package. See inside-r.org/packages/cran/Rmpfr/docs/determinant.mpfrMatrix. – bgoldst Apr 24 '15 at 9:58