From d54d64ded85b33ccdb6ec1b3211c8e91b866c2e4 Mon Sep 17 00:00:00 2001
From: Anthony Scemama <scemama@irsamc.ups-tlse.fr>
Date: Tue, 6 Jun 2017 00:26:47 +0200
Subject: [PATCH] Trying to fix travis

---
 plugins/Hartree_Fock/diagonalize_fock.irp.f | 55 +++++++++------------
 1 file changed, 24 insertions(+), 31 deletions(-)

diff --git a/plugins/Hartree_Fock/diagonalize_fock.irp.f b/plugins/Hartree_Fock/diagonalize_fock.irp.f
index 3d5000b5..ff67a4e4 100644
--- a/plugins/Hartree_Fock/diagonalize_fock.irp.f
+++ b/plugins/Hartree_Fock/diagonalize_fock.irp.f
@@ -55,49 +55,42 @@
        lwork = 1+6*n + 2*n*n
        liwork = 3 + 5*n
        
-       allocate(work(lwork))
+!       allocate(work(lwork))
 !       allocate(iwork(liwork) )
-       
-       lwork = -1
-       liwork = -1
-       
+!       
+!       lwork = -1
+!       liwork = -1
+!       
 !       call dsyevd( 'V', 'U', mo_tot_num, F,                         &
 !           size(F,1), diagonal_Fock_matrix_mo,                       &
 !           work, lwork, iwork, liwork, info)
-       
-       call dsyev( 'V', 'U', mo_tot_num, F,                         &
-           size(F,1), diagonal_Fock_matrix_mo,                       &
-           work, lwork, info)
-       
-       if (info /= 0) then
-         print *,  irp_here//' failed : ', info
-         stop 1
-       endif
-       lwork = int(work(1))
+!       
+!       if (info /= 0) then
+!         print *,  irp_here//' failed : ', info
+!         stop 1
+!       endif
+!       lwork = int(work(1))
 !       liwork = iwork(1)
 !       deallocate(iwork)
-       deallocate(work)
+!       deallocate(work)
 
        allocate(work(lwork))
-!       allocate(iwork(liwork) )
-!       call dsyevd( 'V', 'U', mo_tot_num, F,                         &
-!           size(F,1), diagonal_Fock_matrix_mo,                       &
-!           work, lwork, iwork, liwork, info)
-!       deallocate(iwork)
-
-       call dsyev( 'V', 'U', mo_tot_num, F,                         &
+       allocate(iwork(liwork) )
+       call dsyevd( 'V', 'U', mo_tot_num, F,                         &
            size(F,1), diagonal_Fock_matrix_mo,                       &
-           work, lwork, info)
+           work, lwork, iwork, liwork, info)
+       deallocate(iwork)
+
 
        if (info /= 0) then
-!         call dsyev( 'V', 'U', mo_tot_num, F,                        &
-!             size(F,1), diagonal_Fock_matrix_mo,                     &
-!             work, lwork, info)
-!         
-!         if (info /= 0) then
-           print *,  irp_here//' failed : ', info
+         call dsyev( 'V', 'U', mo_tot_num, F,                        &
+             size(F,1), diagonal_Fock_matrix_mo,                     &
+             work, lwork, info)
+
+         if (info /= 0) then
+           print *,  irp_here//' DSYEV failed : ', info
            stop 1
-!         endif
+         endif
        endif
 
        call dgemm('N','N',ao_num,mo_tot_num,mo_tot_num, 1.d0, &