36 #include <visp3/core/vpConfig.h>
37 #include <visp3/core/vpMatrix.h>
39 #ifndef DOXYGEN_SHOULD_SKIP_THIS
41 #if defined(VISP_HAVE_LAPACK) && !defined(VISP_HAVE_LAPACK_BUILT_IN)
44 extern "C" void dgemm_(
char *transa,
char *transb, integer *M, integer *N, integer *K,
double *alpha,
double *a,
45 integer *lda,
double *b, integer *ldb,
double *beta,
double *c, integer *ldc);
47 extern "C" void dgemv_(
char *trans, integer *M, integer *N,
double *alpha,
double *a, integer *lda,
double *x,
48 integer *incx,
double *beta,
double *y, integer *incy);
50 void vpMatrix::blas_dgemm(
char trans_a,
char trans_b,
const int M_,
const int N_,
const int K_,
double alpha,
51 double *a_data,
const int lda_,
double *b_data,
const int ldb_,
double beta,
double *c_data,
54 integer M = (integer)M_, K = (integer)K_, N = (integer)N_;
55 integer lda = (integer)lda_, ldb = (integer)ldb_, ldc = (integer)ldc_;
57 dgemm_(&trans_a, &trans_b, &M, &N, &K, &alpha, a_data, &lda, b_data, &ldb, &beta, c_data, &ldc);
60 void vpMatrix::blas_dgemv(
char trans,
const int M_,
const int N_,
double alpha,
double *a_data,
const int lda_,
61 double *x_data,
const int incx_,
double beta,
double *y_data,
const int incy_)
63 integer M = (integer)M_, N = (integer)N_;
64 integer lda = (integer)lda_, incx = (integer)incx_, incy = (integer)incy_;
66 dgemv_(&trans, &M, &N, &alpha, a_data, &lda, x_data, &incx, &beta, y_data, &incy);
72 void dummy_vpMatrix_blas(){};
75 #endif // #ifndef DOXYGEN_SHOULD_SKIP_THIS