Esempio n. 1
0
// 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();
}
Esempio n. 2
0
void fzgetrf_(int* m, int* n, void* a, int* lda, int* ipiv, int* info)
{
	lapack_zgetrf(*m, *n, a, *lda, ipiv, info);
}