If ecutwfc < 4*ecutrho (something that should not be done but sometimes is),

add warning message, check preventing ecutfock to become larger than ecutrho
make.depend in EPW updated

EPW/src/make.depend

@@ -186,16 +186,17 @@ ephwann_shuffle.o : grid.o
 ephwann_shuffle.o : indabs.o
 ephwann_shuffle.o : io_eliashberg.o
 ephwann_shuffle.o : io_epw.o
+ephwann_shuffle.o : io_selfen.o
 ephwann_shuffle.o : io_transport.o
 ephwann_shuffle.o : io_var.o
 ephwann_shuffle.o : low_lvl.o
-ephwann_shuffle.o : plot.o
 ephwann_shuffle.o : printing.o
 ephwann_shuffle.o : rigid_epw.o
 ephwann_shuffle.o : selfen.o
 ephwann_shuffle.o : spectral_func.o
 ephwann_shuffle.o : transport.o
 ephwann_shuffle.o : transport_iter.o
+ephwann_shuffle.o : utilities.o
 ephwann_shuffle.o : wan2bloch.o
 ephwann_shuffle.o : wigner.o
 ephwann_shuffle_mem.o : ../../Modules/cell_base.o
@@ -221,16 +222,17 @@ ephwann_shuffle_mem.o : grid.o
 ephwann_shuffle_mem.o : indabs.o
 ephwann_shuffle_mem.o : io_eliashberg.o
 ephwann_shuffle_mem.o : io_epw.o
+ephwann_shuffle_mem.o : io_selfen.o
 ephwann_shuffle_mem.o : io_transport.o
 ephwann_shuffle_mem.o : io_var.o
 ephwann_shuffle_mem.o : low_lvl.o
-ephwann_shuffle_mem.o : plot.o
 ephwann_shuffle_mem.o : printing.o
 ephwann_shuffle_mem.o : rigid_epw.o
 ephwann_shuffle_mem.o : selfen.o
 ephwann_shuffle_mem.o : spectral_func.o
 ephwann_shuffle_mem.o : transport.o
 ephwann_shuffle_mem.o : transport_iter.o
+ephwann_shuffle_mem.o : utilities.o
 ephwann_shuffle_mem.o : wan2bloch.o
 ephwann_shuffle_mem.o : wigner.o
 epw.o : ../../Modules/control_flags.o
@@ -394,6 +396,15 @@ io_epw.o : elph2.o
 io_epw.o : epwcom.o
 io_epw.o : io_var.o
 io_epw.o : low_lvl.o
+io_selfen.o : ../../Modules/io_files.o
+io_selfen.o : ../../Modules/io_global.o
+io_selfen.o : ../../Modules/kind.o
+io_selfen.o : ../../Modules/mp_world.o
+io_selfen.o : ../../UtilXlib/mp.o
+io_selfen.o : constants_epw.o
+io_selfen.o : elph2.o
+io_selfen.o : epwcom.o
+io_selfen.o : io_var.o
 io_transport.o : ../../Modules/cell_base.o
 io_transport.o : ../../Modules/io_files.o
 io_transport.o : ../../Modules/io_global.o
@@ -453,20 +464,6 @@ low_lvl.o : eliashbergcom.o
 low_lvl.o : elph2.o
 low_lvl.o : epwcom.o
 low_lvl.o : io_var.o
-plot.o : ../../Modules/cell_base.o
-plot.o : ../../Modules/io_files.o
-plot.o : ../../Modules/io_global.o
-plot.o : ../../Modules/kind.o
-plot.o : ../../Modules/mp_global.o
-plot.o : ../../Modules/mp_world.o
-plot.o : ../../PHonon/PH/phcom.o
-plot.o : ../../PW/src/pwcom.o
-plot.o : ../../UtilXlib/mp.o
-plot.o : constants_epw.o
-plot.o : elph2.o
-plot.o : epwcom.o
-plot.o : io_var.o
-plot.o : poolgathering.o
 poolgathering.o : ../../Modules/kind.o
 poolgathering.o : ../../Modules/mp_global.o
 poolgathering.o : ../../UtilXlib/mp.o
@@ -485,6 +482,7 @@ printing.o : division.o
 printing.o : eliashbergcom.o
 printing.o : elph2.o
 printing.o : epwcom.o
+printing.o : io_var.o
 printing.o : kinds_epw.o
 printing.o : poolgathering.o
 pw2wan2epw.o : ../../FFTXlib/fft_interfaces.o
@@ -548,10 +546,11 @@ selfen.o : ../../UtilXlib/mp.o
 selfen.o : constants_epw.o
 selfen.o : elph2.o
 selfen.o : epwcom.o
-selfen.o : io_transport.o
+selfen.o : io_selfen.o
 selfen.o : io_var.o
 selfen.o : poolgathering.o
 spectral_func.o : ../../Modules/cell_base.o
+spectral_func.o : ../../Modules/io_files.o
 spectral_func.o : ../../Modules/io_global.o
 spectral_func.o : ../../Modules/kind.o
 spectral_func.o : ../../Modules/mp_global.o
@@ -562,7 +561,7 @@ spectral_func.o : ../../UtilXlib/mp.o
 spectral_func.o : constants_epw.o
 spectral_func.o : elph2.o
 spectral_func.o : epwcom.o
-spectral_func.o : io_transport.o
+spectral_func.o : io_selfen.o
 spectral_func.o : io_var.o
 spectral_func.o : poolgathering.o
 spectral_func.o : selfen.o
@@ -660,10 +659,17 @@ transport_iter.o : kinds_epw.o
 transport_iter.o : poolgathering.o
 transport_iter.o : printing.o
 transport_iter.o : wigner.o
+utilities.o : ../../Modules/cell_base.o
 utilities.o : ../../Modules/io_global.o
 utilities.o : ../../Modules/kind.o
+utilities.o : ../../Modules/mp_global.o
+utilities.o : ../../Modules/noncol.o
+utilities.o : ../../PHonon/PH/phcom.o
+utilities.o : ../../PW/src/pwcom.o
+utilities.o : ../../UtilXlib/mp.o
 utilities.o : constants_epw.o
 utilities.o : eliashbergcom.o
+utilities.o : elph2.o
 utilities.o : epwcom.o
 wan2bloch.o : ../../Modules/cell_base.o
 wan2bloch.o : ../../Modules/io_files.o

PW/src/input.f90

@@ -1510,7 +1510,7 @@ SUBROUTINE iosys()
   !
   IF(ecutfock <= 0.0_DP) THEN
      ! default case
-     ecutfock_ = 4.0_DP*ecutwfc
+     ecutfock_ = MIN ( ecutrho, 4.0_DP*ecutwfc)
   ELSE
      IF(ecutfock < ecutwfc .OR. ecutfock > ecutrho) CALL errore('iosys', &
           'ecutfock can not be < ecutwfc or > ecutrho!', 1) 
@@ -1603,6 +1603,7 @@ SUBROUTINE set_cutoff ( ecutwfc_in, ecutrho_in, ecutwfc_pp, ecutrho_pp )
   USE gvecs, ONLY : dual
   USE gvect, ONLY : ecutrho
   USE gvecw, ONLY : ecutwfc
+  USE constants, ONLY : eps8
   !
   IMPLICIT NONE
   REAL(dp), INTENT(INOUT) :: ecutwfc_in, ecutrho_in
@@ -1632,6 +1633,8 @@ SUBROUTINE set_cutoff ( ecutwfc_in, ecutrho_in, ecutwfc_pp, ecutrho_pp )
   ecutrho_in = ecutrho
   dual = ecutrho / ecutwfc
   IF ( dual <= 1.0_dp ) CALL errore( 'set_cutoff', 'ecutrho <= ecutwfc?!?', 1 )
+  IF ( dual < 4.0_dp - eps8 ) CALL infomsg( 'set_cutoff', &
+          'ecutrho < 4*ecutwfc, are you sure?' )
   !
 END SUBROUTINE set_cutoff
 !