void
F77_NAME(dtrmm)(const char *side, const char *uplo, const char *transa,
const char *diag, const int *m, const int *n,
const double *alpha, const double *a, const int *lda,
double *b, const int *ldb)
{
DTRMM(side, uplo, transa, diag, m, n, alpha, a, lda, b, ldb);
}
void ProtoMol::Lapack::dtrmm(char *sideA, char *ulA, char *transA, char *diagA,
int *m, int *n, double *alpha, double *A,
int *lda, double *B, int *ldb) {
FAHCheckIn();
#if defined(HAVE_LAPACK)
dtrmm_(sideA, ulA, transA, diagA, m, n, alpha, A, lda, B, ldb);
#elif defined(HAVE_MKL_LAPACK)
DTRMM(sideA, ulA, transA, diagA, m, n, alpha, A, lda, B, ldb);
#else
THROW(std::string(__func__) + " not supported");
#endif
}
void TRMM<double>(const char* side , const char* uplo,
const char* transA, const char* diag,
const int m , const int n , const double& a,
const double* A , const int ldA ,
double* B , const int ldB)
{
ASSERT((side[0] == 'L') || (side[0] == 'R'));
ASSERT((uplo[0] == 'U') || (uplo[0] == 'L'));
ASSERT((transA[0] == 'N') || (transA[0] == 'T'));
ASSERT((diag[0] == 'N') || (diag[0] == 'U'));
DTRMM(F77_CHARACTER(side[0]) , F77_CHARACTER(uplo[0]),
F77_CHARACTER(transA[0]), F77_CHARACTER(diag[0]),
&m, &n, &a, A, &ldA, B, &ldB);
}
