// Inverse of matrix AA void lapack_zinvs(int nn, dcmplx *AA) { int *ipiv = NULL; ipiv = (int *) calloc(nn, sizeof(int)); check_mem(ipiv, "ipiv"); lapack_zgetrf(nn, nn, AA, ipiv); lapack_zgetri(nn, AA, ipiv); freeup(ipiv); return; error: if(ipiv) freeup(ipiv); abort(); }
void fzgetrf_(int* m, int* n, void* a, int* lda, int* ipiv, int* info) { lapack_zgetrf(*m, *n, a, *lda, ipiv, info); }