diff --git a/org/qmckl_blas.org b/org/qmckl_blas.org index 33856fb..3d59d7b 100644 --- a/org/qmckl_blas.org +++ b/org/qmckl_blas.org @@ -146,35 +146,35 @@ integer function qmckl_dgemm_f(context, TransA, TransB, m, n, k, alpha, A, LDA, return endif - if (LDA_2 .ne. m) then + if (LDA_2 /= m) then info = QMCKL_INVALID_ARG_9 return endif - if (LDB_2 .ne. k) then + if (LDB_2 /= k) then info = QMCKL_INVALID_ARG_10 return endif - if (LDC .ne. m) then + if (LDC /= m) then info = QMCKL_INVALID_ARG_13 return endif if (TransA) then - if (alpha == 1.d0 && beta == 0.d0) then + if (alpha == 1.d0 .and. beta == 0.d0) then C = matmul(AT,B) else C = beta*C + alpha*matmul(AT,B) endif else if (TransB) then - if (alpha == 1.d0 && beta == 0.d0) then + if (alpha == 1.d0 .and. beta == 0.d0) then C = matmul(A,BT) else C = beta*C + alpha*matmul(A,BT) endif else - if (alpha == 1.d0 && beta == 0.d0) then + if (alpha == 1.d0 .and. beta == 0.d0) then C = matmul(A,B) else C = beta*C + alpha*matmul(A,B)