LCPQ/qmcchem
10
1
Fork 1
mirror of https://gitlab.com/scemama/qmcchem.git synced 2024-12-22 12:23:30 +01:00
Code Releases Activity

Unused variable

Browse Source
This commit is contained in:
Anthony Scemama 2017-09-05 17:39:11 +02:00
parent fed153e513
commit 0e2ee815f6
1 changed files with 0 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/TOOLS/invert.irp.f
View File

@ -43,7 +43,6 @@ subroutine invert_general(a,LDA,na,det_l)
integer :: ipiv(LDA) integer :: ipiv(LDA)
!DIR$ ATTRIBUTES ALIGN: $IRP_ALIGN :: ipiv !DIR$ ATTRIBUTES ALIGN: $IRP_ALIGN :: ipiv
integer :: lwork integer :: lwork
double precision :: f
integer :: i, j integer :: i, j
call dgetrf(na, na, a, LDA, ipiv, inf ) call dgetrf(na, na, a, LDA, ipiv, inf )
det_l = 1.d0 det_l = 1.d0
@ -74,7 +73,6 @@ subroutine sinvert(a,LDA,na,det_l)
integer :: ipiv(LDA) integer :: ipiv(LDA)
!DIR$ ATTRIBUTES ALIGN: $IRP_ALIGN :: ipiv !DIR$ ATTRIBUTES ALIGN: $IRP_ALIGN :: ipiv
integer :: lwork integer :: lwork
real :: f
integer :: i, j integer :: i, j
call sgetrf(na, na, a, LDA, ipiv, inf ) call sgetrf(na, na, a, LDA, ipiv, inf )
det_l = 1.d0 det_l = 1.d0
@ -113,7 +111,6 @@ subroutine invert2(a,LDA,na,det_l)
double precision :: det_l double precision :: det_l
double precision :: b(2,2) double precision :: b(2,2)
double precision :: f
b(1,1) = a(1,1) b(1,1) = a(1,1)
b(2,1) = a(2,1) b(2,1) = a(2,1)
b(1,2) = a(1,2) b(1,2) = a(1,2)
@ -134,7 +131,6 @@ subroutine invert3(a,LDA,na,det_l)
double precision :: b(4,3) double precision :: b(4,3)
!DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b
integer :: i integer :: i
double precision :: f
det_l = a(1,1)*(a(2,2)*a(3,3)-a(2,3)*a(3,2)) & 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)) & -a(1,2)*(a(2,1)*a(3,3)-a(2,3)*a(3,1)) &
+a(1,3)*(a(2,1)*a(3,2)-a(2,2)*a(3,1)) +a(1,3)*(a(2,1)*a(3,2)-a(2,2)*a(3,1))
@ -166,7 +162,6 @@ subroutine invert4(a,LDA,na,det_l)
double precision :: b(4,4) double precision :: b(4,4)
!DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b
integer :: i,j integer :: i,j
double precision :: f
det_l = a(1,1)*(a(2,2)*(a(3,3)*a(4,4)-a(3,4)*a(4,3)) & 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)) & -a(2,3)*(a(3,2)*a(4,4)-a(3,4)*a(4,2)) &
+a(2,4)*(a(3,2)*a(4,3)-a(3,3)*a(4,2))) & +a(2,4)*(a(3,2)*a(4,3)-a(3,3)*a(4,2))) &
@ -217,7 +212,6 @@ subroutine invert5(a,LDA,na,det_l)
double precision :: b(5,5) double precision :: b(5,5)
!DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: b
integer :: i,j integer :: i,j
double precision :: f
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)*( & 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)))- & 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)))- &
a(2,3)*(a(3,2)*(a(4,4)*a(5,5)-a(4,5)*a(5,4))-a(3,4)*(a(4,2)*a(5,5)- & a(2,3)*(a(3,2)*(a(4,4)*a(5,5)-a(4,5)*a(5,4))-a(3,4)*(a(4,2)*a(5,5)- &
@ -410,7 +404,6 @@ subroutine invert_update(a,LDA,na,det_l,b)
integer :: ipiv(LDA) integer :: ipiv(LDA)
!DIR$ ATTRIBUTES ALIGN: $IRP_ALIGN :: ipiv !DIR$ ATTRIBUTES ALIGN: $IRP_ALIGN :: ipiv
integer :: lwork integer :: lwork
double precision :: f
integer :: i, j integer :: i, j
double precision :: bold(LDA,na) double precision :: bold(LDA,na)
double precision :: ba(LDA,na) double precision :: ba(LDA,na)