compute determinant of complex matrix fortran90

Question

I need to compute the determinant of complex matrix which is symmetric. Size of matrix ranges from 500*500 to 2000*2000. Is there any subroutine for me to call? By the way, I use ifort to compile.

Answer 1

The easiest way would be to do an LU-decomposition as described here. I would suggest using LAPACK for this task...

This article has some code in C doing that for a real-valued symmetric matrix, so you need to exchange dspsv by zspsv to handle double-precision complex matrices.