Computing the determinant of a C array

Question

I have written a program that generates a few N x N matrices for which I would like to compute their determinant. Related to this I have two questions

• What library is the best for doing so? I would like the fastest possible library since I have millions of such matrices.

• Are there any specifics I should take care of when casting the result to an integer? All matrices that I will generate have integer determinants and I would like to make sure that no rounding error skews the correct value of the determinant.

Edit. If possible provide an example of computing the determinant for the recommended library.

Answer 1

As for Matrix Libraries, it looks like that question is answered here:

Recommendations for a small c-based vector and matrix library

As for casting to an integer: if the determinant is not an integer, then you shouldn't be casting it to an integer, you should be using round, floor, or ceil to convert it in an acceptable way. These will probably give you integral values, but you will still need to cast them; however, you will now be able to do so without fear of losing any information.

Answer 2

You have GSL, but the choice really depends on your matrices. Are the matrices dense or sparse? Is N big or small? For small Ns you may find that coding the determinant yourself using Cramer's rule or Gauss elimination is faster, since most hight performance libraries focus on big matrices and their optimisations may introduce overhead on simple problems.

Answer 3

You can work wonders with matrices by blas and lapack. They are actually written in fortran and using them from "c" is kind of a tweak. but all in all they can crunch numbers in horrendous speed.

http://www.netlib.org/lapack/lug/node11.html