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