TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2024-06-30 00:44:52 +02:00
Code Issues Projects Releases Wiki Activity

Removed DIR$

Browse Source
This commit is contained in:
Anthony Scemama 2021-12-12 11:20:11 +01:00
parent 3f9daae3d3
commit 347b2153b5
1 changed files with 3 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
org/qmckl_blas.org
View File

@ -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)* &