Use LU factorization for determinant calculation
Currently the linpack routines are used for calculating determinants.
I would prefer to use the LU factorization mechanism in LAPACK to more consistently get the determinent.
Once the LU is formed, one can calculate the determinant as:
call zgetrf(A,...)
det = product(diagonal(A))
the sign of det
should be corrected according to the number of pivots. If there is an odd number of pivots, then -det
otherwise +det
This would allow us to remove the linpack code all-together.
Edited by Nick R. Papior