void
F77_NAME(dtrsm)(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)
{
DTRSM(side, uplo, transa, diag, m, n, alpha, a, lda, b, ldb);
}
void ProtoMol::Lapack::dtrsm(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)
dtrsm_(sideA, ulA, transA, diagA, m, n, alpha, A, lda, B, ldb);
#elif defined(HAVE_MKL_LAPACK)
DTRSM(sideA, ulA, transA, diagA, m, n, alpha, A, lda, B, ldb);
#else
THROW(std::string(__func__) + " not supported");
#endif
}
void IpBlasDtrsm(bool trans, Index ndim, Index nrhs, Number alpha,
const Number* A, Index ldA, Number* B, Index ldB)
{
ipfint M=ndim, N=nrhs, LDA=ldA, LDB=ldB;
char SIDE = 'L';
char UPLO = 'L';
char TRANSA;
if (trans) {
TRANSA = 'T';
}
else {
TRANSA = 'N';
}
char DIAG = 'N';
DTRSM(&SIDE, &UPLO, &TRANSA, &DIAG, &M, &N,
&alpha, A, &LDA, B, &LDB, 1, 1, 1, 1);
}
