[Mplapack-devel] FORTRAN interface for MPACK

[Fletcher, J. P <j.p...@as...>]

// mpack_qd_cpp.cpp

// Interface qd_complex sqrt function.

// John Fletcher   J.P...@as...  April 2011

#include <mblas_qd.h>
#include <mlapack_qd.h>

// The function bodies to be called from FORTRAN have to be extern "C".
// They also need a trailing _

extern "C" {


  qd_complex qdc_sqrt_(const qd_complex &z)
  {
    return sqrt(z);
  }

  void rgetrf_(const mpackint &m, const mpackint &n, qd_real * A, const mpackint &lda, mpackint *ipiv, mpackint *info )
  {
    Rgetrf(m,n,A,lda,ipiv,info);
  }

  void rgetrs_(const char *trans, const mpackint &n, const mpackint &nrhs, qd_real * A, const mpackint &lda, mpackint *ipiv, qd_real *B, const mpackint &ldb, mpackint *info )
  {
    Rgetrs(trans,n,nrhs,A,lda,ipiv,B,ldb,info);
  }

  void cgetrf_(const mpackint &m, const mpackint &n, qd_complex * A, const mpackint &lda, mpackint *ipiv, mpackint *info )
  {
    Cgetrf(m,n,A,lda,ipiv,info);
  }

  void cgetrs_(const char *trans, const mpackint &n, const mpackint &nrhs, qd_complex * A, const mpackint &lda, mpackint *ipiv, qd_complex *B, const mpackint &ldb, mpackint *info )
  {
    Cgetrs(trans,n,nrhs,A,lda,ipiv,B,ldb,info);
  }

} /* extern "C" */