TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2025-01-03 10:06:09 +01:00
Code Issues Projects Releases Wiki Activity

Protected qmckldgemm stuff with ifdefs.

Browse Source
This commit is contained in:
v1j4y 2022-08-09 16:11:40 +02:00
parent 46a26c7028
commit e822694a56
1 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
org/qmckl_blas.org
View File

@ -889,7 +889,9 @@ integer function qmckl_dgemm_f(context, TransA, TransB, &
m, n, k, alpha, A, LDA, B, LDB, beta, C, LDC) & m, n, k, alpha, A, LDA, B, LDB, beta, C, LDC) &
result(info) result(info)
use qmckl use qmckl
#ifdef HAVE_LIBQMCKLDGEMM
use qmckl_dgemm_tiled_module use qmckl_dgemm_tiled_module
#endif
implicit none implicit none
integer(qmckl_context), intent(in) :: context integer(qmckl_context), intent(in) :: context
character , intent(in) :: TransA, TransB character , intent(in) :: TransA, TransB
@ -943,9 +945,7 @@ integer function qmckl_dgemm_f(context, TransA, TransB, &
return return
endif endif
!call dgemm(transA, transB, int(m,4), int(n,4), int(k,4), & #ifdef HAVE_LIBQMCKLDGEMM
! alpha, A, int(LDA,4), B, int(LDB,4), beta, C, int(LDC,4))
! Copy A to A1 ! Copy A to A1
allocate(A1(k,m)) allocate(A1(k,m))
do j=1,m do j=1,m
@ -984,6 +984,11 @@ integer function qmckl_dgemm_f(context, TransA, TransB, &
end do end do
deallocate(A1,B1,C1) deallocate(A1,B1,C1)
#else
call dgemm(transA, transB, int(m,4), int(n,4), int(k,4), &
alpha, A, int(LDA,4), B, int(LDB,4), beta, C, int(LDC,4))
#endif
end function qmckl_dgemm_f end function qmckl_dgemm_f
#+end_src #+end_src