From cf3550b6b7c2b1a7a6bd06bafc1e1ac898e84c1a Mon Sep 17 00:00:00 2001
From: v1j4y <vijay.gopal.c@gmail.com>
Date: Wed, 15 Sep 2021 11:55:45 +0200
Subject: [PATCH] Fixed some bugs, now compiles.

---
 org/qmckl_blas.org | 44 +++++++++++++++++++++++---------------------
 1 file changed, 23 insertions(+), 21 deletions(-)

diff --git a/org/qmckl_blas.org b/org/qmckl_blas.org
index 72e9cb5..fd7e34e 100644
--- a/org/qmckl_blas.org
+++ b/org/qmckl_blas.org
@@ -41,7 +41,7 @@ int main() {
    | double        | beta     | in  | Array containing the $n \times m$ matrix $B$ |
    | double        | C[][ldc] | out | Array containing the $n \times m$ matrix $B$ |
    | int64_t       | ldc      | in  | Leading dimension of array ~B~               |
-
+   
 *** Requirements
 
     - ~context~ is not ~QMCKL_NULL_CONTEXT~
@@ -86,14 +86,15 @@ integer function qmckl_dgemm_f(context, TransA, TransB, m, n, k, alpha, A, LDA,
   use qmckl
   implicit none
   integer(qmckl_context)  , intent(in)  :: context
-  logical*8  , intent(in)  :: TransA, TransN
+  logical*8  , intent(in)  :: TransA, TransB
   integer*8  , intent(in)  :: m, n, k
+  real*8     , intent(in)  :: alpha, beta
   integer*8  , intent(in)  :: lda
-  real*8     , intent(in)  :: A(lda,*)
+  real*8     , intent(in)  :: A(m,n)
   integer*8  , intent(in)  :: ldb
-  real*8     , intent(in)  :: B(ldb,*)
+  real*8     , intent(in)  :: B(n,k)
   integer*8  , intent(in)  :: ldc
-  real*8     , intent(out) :: C(ldc,*)
+  real*8     , intent(out) :: C(m,n)
 
   integer*8 :: i,j
 
@@ -158,13 +159,13 @@ end function qmckl_dgemm_f
       integer (c_int64_t) , intent(in)  , value :: n
       integer (c_int64_t) , intent(in)  , value :: k
       real    (c_double ) , intent(in)  , value :: alpha
-      real    (c_double ) , intent(in)          :: A(lda,*)
       integer (c_int64_t) , intent(in)  , value :: lda
-      real    (c_double ) , intent(in)          :: B(ldb,*)
+      real    (c_double ) , intent(in)          :: A(lda,*)
       integer (c_int64_t) , intent(in)  , value :: ldb
+      real    (c_double ) , intent(in)          :: B(ldb,*)
       real    (c_double ) , intent(in)  , value :: beta
-      real    (c_double ) , intent(out)         :: C(ldc,*)
       integer (c_int64_t) , intent(in)  , value :: ldc
+      real    (c_double ) , intent(out)         :: C(ldc,*)
 
       integer(c_int32_t), external :: qmckl_dgemm_f
       info = qmckl_dgemm_f &
@@ -180,26 +181,26 @@ end function qmckl_dgemm_f
     #+begin_src f90 :tangle (eval fh_func) :comments org :exports none
     interface
       integer(c_int32_t) function qmckl_dgemm &
-          (context, TransA, TransB, m, n, k, alpha, A, lda, B, ldb, beat, C, ldc) &
+          (context, TransA, TransB, m, n, k, alpha, A, lda, B, ldb, beta, C, ldc) &
           bind(C)
         use, intrinsic :: iso_c_binding
         import
         implicit none
 
         integer (c_int64_t) , intent(in)  , value :: context
-        logical*8 (c_bool) , intent(in)  , value :: TransA
-        logical*8 (c_bool) , intent(in)  , value :: TransB
+        logical*8          , intent(in)  , value :: TransA
+        logical*8          , intent(in)  , value :: TransB
         integer (c_int64_t) , intent(in)  , value :: m
         integer (c_int64_t) , intent(in)  , value :: n
         integer (c_int64_t) , intent(in)  , value :: k
         real    (c_double ) , intent(in)          :: alpha
-        real    (c_double ) , intent(in)          :: A(lda,*)
         integer (c_int64_t) , intent(in)  , value :: lda
+        real    (c_double ) , intent(in)          :: A(lda,*)
+        integer (c_int64_t) , intent(in)  , value :: ldb
         real    (c_double ) , intent(in)          :: B(ldb,*)
-        integer (c_int64_t) , intent(in)  , value :: ldb
         real    (c_double ) , intent(in)          :: beta
-        real    (c_double ) , intent(out)         :: C(ldb,*)
-        integer (c_int64_t) , intent(in)  , value :: ldb
+        integer (c_int64_t) , intent(in)  , value :: ldc
+        real    (c_double ) , intent(out)         :: C(ldc,*)
 
       end function qmckl_dgemm
     end interface
@@ -214,15 +215,15 @@ integer(qmckl_exit_code) function test_qmckl_dgemm(context) bind(C)
 
   double precision, allocatable :: A(:,:), B(:,:), C(:,:), D(:,:)
   integer*8                     :: m, n, k, LDA, LDB, LDC
-  integer*8                     :: i,j
+  integer*8                     :: i,j,l
   logical*8                     :: TransA, TransB
-  double precision              :: x
+  double precision              :: x, alpha, beta
 
   TransA = .False.
   TransB = .False.
-  m = 5
-  k = 4
-  n = 6
+  m = 5_8
+  k = 4_8
+  n = 6_8
   LDA = m
   LDB = k
   LDC = m
@@ -255,11 +256,12 @@ integer(qmckl_exit_code) function test_qmckl_dgemm(context) bind(C)
   do j=1,m
      do i=l,n
        do l=1,k
-          D(i,j) += D(i,j) + (A(i,k)*B(k,j))
+          D(i,j) = D(i,j) + A(i,k)*B(k,j)
        end do
        x = x + (D(i,j) - C(i,j))**2
      end do
   end do
+  print *,"DABS(X)= ",dabs(x)
 
   if (dabs(x) <= 1.d-15) then
      test_qmckl_dgemm = QMCKL_SUCCESS