diff --git a/org/qmckl_blas.org b/org/qmckl_blas.org
index 81db1b7..1a50529 100644
--- a/org/qmckl_blas.org
+++ b/org/qmckl_blas.org
@@ -161,19 +161,19 @@ integer function qmckl_dgemm_f(context, TransA, TransB, m, n, k, alpha, A, LDA,
   endif
 
    if (TransA) then
-     if (alpha == 1.d0 && beta == 0.d0) then
+     if (alpha .eq. 1.0d0 .and. beta .eq. 0.0d0) 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 .eq. 1.0d0 .and. beta .eq. 0.0d0) 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 .eq. 1.0d0 .and. beta .eq. 0.0d0) then
        C = matmul(A,B)
      else
        C = beta*C + alpha*matmul(A,B)