From 60356cdbbb04531e0076b27a5d9f3b1519262e5f Mon Sep 17 00:00:00 2001
From: Anthony Scemama <scemama@irsamc.ups-tlse.fr>
Date: Sun, 12 Dec 2021 11:20:11 +0100
Subject: [PATCH] Removed DIR$

---
 org/qmckl_blas.org | 15 +++------------
 1 file changed, 3 insertions(+), 12 deletions(-)

diff --git a/org/qmckl_blas.org b/org/qmckl_blas.org
index 279fd80..e560095 100644
--- a/org/qmckl_blas.org
+++ b/org/qmckl_blas.org
@@ -434,24 +434,18 @@ integer function qmckl_adjoint_f(context, ma, na, LDA, A, det_l) &
 
  select case (na)
   case default
-!DIR$ forceinline
     print *," TODO: Implement general adjoint"
     stop 0
   case (5)
-!DIR$ forceinline
     call adjoint5(a,LDA,na,det_l)
   case (4)
-!DIR$ forceinline
     call adjoint4(a,LDA,na,det_l)
     
   case (3)
-!DIR$ forceinline
     call adjoint3(a,LDA,na,det_l)
   case (2)
-!DIR$ forceinline
     call adjoint2(a,LDA,na,det_l)
   case (1)
-!DIR$ forceinline
     call adjoint1(a,LDA,na,det_l)
   case (0)
     det_l=1.d0
@@ -668,7 +662,6 @@ subroutine cofactor3(a,LDA,na,det_l)
  integer*8, intent(in)             :: na
  double precision, intent(inout) :: det_l
  double precision :: b(4,3)
- !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b
  integer :: i
  det_l = a(1,1)*(a(2,2)*a(3,3)-a(2,3)*a(3,2)) &
         -a(1,2)*(a(2,1)*a(3,3)-a(2,3)*a(3,1)) &
@@ -677,7 +670,7 @@ subroutine cofactor3(a,LDA,na,det_l)
   b(i,1) = a(i,1)
   b(i,2) = a(i,2)
   b(i,3) = a(i,3)
- enddo
+ end do
  a(1,1) =  b(2,2)*b(3,3) - b(2,3)*b(3,2)
  a(2,1) =  b(2,3)*b(3,1) - b(2,1)*b(3,3)
  a(3,1) =  b(2,1)*b(3,2) - b(2,2)*b(3,1)
@@ -699,7 +692,6 @@ subroutine cofactor4(a,LDA,na,det_l)
  integer*8, intent(in)             :: na
  double precision, intent(inout) :: det_l
  double precision :: b(4,4)
- !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b
  integer :: i,j
  det_l =  a(1,1)*(a(2,2)*(a(3,3)*a(4,4)-a(3,4)*a(4,3))  &
                  -a(2,3)*(a(3,2)*a(4,4)-a(3,4)*a(4,2))  &
@@ -718,7 +710,7 @@ subroutine cofactor4(a,LDA,na,det_l)
   b(2,i) = a(2,i)
   b(3,i) = a(3,i)
   b(4,i) = a(4,i)
- enddo
+ end do
  
  a(1,1) =  b(2,2)*(b(3,3)*b(4,4)-b(3,4)*b(4,3))-b(2,3)*(b(3,2)*b(4,4)-b(3,4)*b(4,2))+b(2,4)*(b(3,2)*b(4,3)-b(3,3)*b(4,2))
  a(2,1) = -b(2,1)*(b(3,3)*b(4,4)-b(3,4)*b(4,3))+b(2,3)*(b(3,1)*b(4,4)-b(3,4)*b(4,1))-b(2,4)*(b(3,1)*b(4,3)-b(3,3)*b(4,1))
@@ -749,7 +741,6 @@ subroutine cofactor5(a,LDA,na,det_l)
  integer*8, intent(in)             :: na
  double precision, intent(inout) :: det_l
  double precision :: b(5,5)
- !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b
  integer :: i,j
  det_l = a(1,1)*(a(2,2)*(a(3,3)*(a(4,4)*a(5,5)-a(4,5)*a(5,4))-a(3,4)*( &
  a(4,3)*a(5,5)-a(4,5)*a(5,3))+a(3,5)*(a(4,3)*a(5,4)-a(4,4)*a(5,3)))- &
@@ -795,7 +786,7 @@ subroutine cofactor5(a,LDA,na,det_l)
   b(3,i) = a(3,i)
   b(4,i) = a(4,i)
   b(5,i) = a(5,i)
- enddo
+ end do
  
  a(1,1) = &
  (b(2,2)*(b(3,3)*(b(4,4)*b(5,5)-b(4,5)*b(5,4))-b(3,4)*(b(4,3)*b(5,5)-b(4,5)*b(5,3))+b(3,5)*(b(4,3)*b(5,4)-b(4,4)*b(5,3)))-b(2,3)* &