Change the default vme == .false. to .true.

EPW/src/epw_readin.f90

@@ -421,7 +421,7 @@
               ! INTEGER :: modenum = 0. In QE 5, modenum variable does not exist
               ! anymore. Change the default EPW value to match the previous QE
               ! one. 
-  vme = .FALSE.
+  vme = .TRUE. ! Was false by default until EPW 5.1
   ephwrite = .FALSE.
   band_plot = .FALSE.
   nqsmear = 10

test-suite/epw_base/epw1.in

@@ -11,6 +11,7 @@
   epbread     = .false.
 
   etf_mem     = 0
+  vme         = .false.
 
   epwwrite    = .true.
   epwread     = .false.

test-suite/epw_base/epw10.in

@@ -11,6 +11,7 @@
   epbread     = .true.
 
   etf_mem     = 1
+  vme         = .false.
   use_ws      = .true.
 
   epwwrite    = .true.

test-suite/epw_base/epw2.in

@@ -12,6 +12,7 @@
   epbread     = .false.
 
   etf_mem     = 0
+  vme         = .false.
 
   epwwrite    = .false.
   epwread     = .true.

test-suite/epw_base/epw3.in

@@ -12,6 +12,7 @@
   epbread     = .false.
 
   etf_mem     = 0
+  vme         = .false.
 
   epwwrite    = .false.
   epwread     = .true.

test-suite/epw_base/epw4.in

@@ -14,6 +14,7 @@
   epbread     = .false.
 
   etf_mem     = 0
+  vme         = .false.
 
   epwwrite    = .false.
   epwread     = .true.

test-suite/epw_base/epw5.in

@@ -12,6 +12,7 @@
   epbread     = .false.
 
   etf_mem     = 0
+  vme         = .false.
 
   epwwrite    = .true.
   epwread     = .false.
@@ -30,7 +31,6 @@
   phonselfen  = .false.
   a2f         = .false.
 
-
   fsthick     = 15 ! eV 
   eptemp      = 300 ! K (same as PRB 76, 165108)
   degaussw    = 1.0 ! eV

test-suite/epw_base/epw6.in

@@ -12,6 +12,7 @@
   epbread     = .false.
 
   etf_mem     = 0
+  vme         = .false. 
 
   epwwrite    = .false.
   epwread     = .true.
@@ -36,7 +37,6 @@
   wmax_specfun = 20
   nw_specfun   = 20
 
-
   fsthick     = 15 ! eV 
   eptemp      = 300 ! K 
   degaussw    = 1.0 ! eV

test-suite/epw_base/epw7.in

@@ -12,6 +12,7 @@
   epbread     = .false.
 
   etf_mem     = 2
+  vme         = .false. 
 
   restart      = .true.
   restart_freq = 50
@@ -34,7 +35,6 @@
   phonselfen  = .false.
   a2f         = .false.
 
-
   fsthick     = 1.5 ! eV 
   eptemp      = 300 ! K (same as PRB 76, 165108)
   degaussw    = 0.3 ! eV

test-suite/epw_base/epw8.in

@@ -12,6 +12,7 @@
   epbread     = .false.
 
   etf_mem     = 2
+  vme         = .false.
 
   restart      = .true.
   restart_freq = 50      ! We test the restart feature
@@ -34,7 +35,6 @@
   phonselfen  = .false.
   a2f         = .false.
 
-
   fsthick     = 1.5 ! eV 
   eptemp      = 300 ! K (same as PRB 76, 165108)
   degaussw    = 0.3 ! eV

test-suite/epw_base/epw9.in

@@ -13,6 +13,7 @@
   epbread     = .false.
 
   etf_mem     = 1
+  vme         = .false.
 
   epwwrite    = .false.
   epwread     = .true.
@@ -38,7 +39,6 @@
   cumulant     = .true.
   bnd_cum      = 4
 
-
   fsthick     = 15 ! eV 
   eptemp      = 300 ! K (same as PRB 76, 165108)
   degaussw    = 1.0 ! eV

test-suite/epw_metal/epw1.in

@@ -12,6 +12,8 @@
   epwwrite    = .true.
   epwread     = .false.
 
+  vme         = .false.
+
   nbndsub     =  4
   nbndskip    =  5
 
@@ -40,7 +42,6 @@
   phonselfen  = .true.
   a2f         = .true.
 
-
   fsthick     = 6 ! eV
   eptemp      = 300 ! K
   degaussw    = 0.1  ! eV

test-suite/epw_metal/epw2.in

@@ -12,6 +12,8 @@
   epwwrite    = .true.
   epwread     = .false.
 
+  vme         = .false. 
+
   nbndsub     =  4
   nbndskip    =  5
 
@@ -43,7 +45,6 @@
   lifc        = .true.
   asr_typ     = 'crystal'
 
-
   fsthick     = 6 ! eV
   eptemp      = 300 ! K
   degaussw    = 0.1  ! eV

test-suite/epw_metal/epw3.in

@@ -0,0 +1,75 @@
+--
+&inputepw
+  prefix      = 'pb'
+  amass(1)    = 207.2
+  outdir      = './'
+
+  elph        = .true.
+  kmaps       = .true.
+  epbwrite    = .false.
+  epbread     = .false.
+
+  epwwrite    = .false.
+  epwread     = .true.
+
+  vme         = .false.
+
+  nbndsub     =  4
+  nbndskip    =  5
+
+  wannierize  = .false.
+  num_iter    = 300
+  dis_win_max = 21
+  dis_win_min = -3
+  dis_froz_min= -3
+  dis_froz_max= 13.5
+  proj(1)     = 'Pb:sp3'
+
+  wdata(1) = 'bands_plot = .true.'
+  wdata(2) = 'begin kpoint_path'
+  wdata(3) = 'G 0.00 0.00 0.00 X 0.00 0.50 0.50'
+  wdata(4) = 'X 0.00 0.50 0.50 W 0.25 0.50 0.75'
+  wdata(5) = 'W 0.25 0.50 0.75 L 0.50 0.50 0.50'
+  wdata(6) = 'L 0.50 0.50 0.50 K 0.375 0.375 0.75'
+  wdata(7) = 'K 0.375 0.375 0.75 G 0.00 0.00 0.00'
+  wdata(8) = 'G 0.00 0.00 0.00 L 0.50 0.50 0.50'
+  wdata(9) = 'end kpoint_path'
+  wdata(10) = 'bands_plot_format = gnuplot'
+
+  iverbosity  = 0
+
+  elecselfen  = .false.
+  phonselfen  = .true.
+  a2f         = .false.
+
+  lifc        = .true.
+  asr_typ     = 'crystal'
+
+  fsthick     = 6 ! eV
+  eptemp      = 300 ! K
+  degaussw    = 0.1  ! eV
+  degaussq    = 0.05 ! meV
+
+  dvscf_dir   = './save'
+
+  nkf1        = 6
+  nkf2        = 6
+  nkf3        = 6
+
+  nqf1        = 6
+  nqf2        = 6
+  nqf3        = 6
+
+  nk1         = 3
+  nk2         = 3
+  nk3         = 3
+
+  nq1         = 3
+  nq2         = 3
+  nq3         = 3
+ /
+      4 cartesian
+  0.000000000000   0.000000000000   0.000000000000 
+ -0.333333333333   0.333333333333  -0.333333333333 
+  0.000000000000   0.666666666667   0.000000000000 
+  0.666666666667   0.000000000000   0.666666666667 

test-suite/epw_mob/epw1.in

@@ -14,6 +14,7 @@
   epwread     = .false.
 
   etf_mem     = 1
+  vme         = .false. 
 
   lifc        = .true.
   asr_typ     = 'crystal'

test-suite/epw_mob/epw2.in

@@ -32,6 +32,7 @@
   asr_typ     = 'crystal'
 
   etf_mem     = 1
+  vme         = .false.
 
   nbndsub     =  16
   nbndskip    =  0

test-suite/epw_mob/epw3.in

@@ -32,6 +32,7 @@
   asr_typ     = 'crystal'
 
   etf_mem     = 1
+  vme         = .false. 
 
   nbndsub     =  16
   nbndskip    =  0

test-suite/epw_mob/epw4.in

@@ -27,6 +27,7 @@
   asr_typ     = 'crystal'
 
   etf_mem     = 1
+  vme         = .false. 
 
   nbndsub     =  16
   nbndskip    =  0

test-suite/epw_mob_polar/epw1.in

@@ -15,7 +15,7 @@
   epwread     = .false.
 
   etf_mem     = 1
-  vme         = .false.
+  vme         = .true.
 
   lifc        = .false.
   asr_typ     = 'simple'

test-suite/epw_mob_polar/epw2.in

@@ -34,7 +34,7 @@
   selecqread   = .false.
 
   etf_mem     = 1
-  vme         = .false.
+  vme         = .true.
 
   lifc        = .false.
   asr_typ     = 'simple'

test-suite/epw_mob_polar/epw3.in

@@ -34,7 +34,7 @@
   selecqread   = .false.
 
   etf_mem     = 1
-  vme         = .false.
+  vme         = .true.
 
   lifc        = .false.
   asr_typ     = 'simple'

test-suite/epw_pl/epw1.in

@@ -12,6 +12,8 @@
   epwwrite    = .true.
   epwread     = .false.
 
+  vme         = .false. 
+
   nbndsub     =  8
   nbndskip    =  0
 
@@ -43,7 +45,6 @@
    
   band_plot   = .true.
 
-
   fsthick     = 3.0
   eptemp      = 20.0
   degaussw    = 0.01

test-suite/epw_polar/epw1.in

@@ -14,6 +14,7 @@
   epwread     = .false.
 
   etf_mem     = 1
+  vme         = .false. 
   lpolar      = .true.
   lphase      = .true. 
 
@@ -40,7 +41,6 @@
   phonselfen  = .false.
   a2f         = .false.
 
-
   fsthick     = 2.0 ! eV
   eptemp      = 300 ! K
   degaussw    = 0.1 ! eV

test-suite/epw_polar/epw2.in

@@ -14,6 +14,7 @@
   epwread     = .false.
 
   etf_mem     = 2
+  vme         = .false. 
   lpolar      = .true.
   lphase      = .true.
 
@@ -40,7 +41,6 @@
   phonselfen  = .false.
   a2f         = .false.
 
-
   fsthick     = 2.0 ! eV
   eptemp      = 300 ! K
   degaussw    = 0.1 ! eV

test-suite/epw_polar/epw3.in

@@ -14,6 +14,7 @@
   epwread     = .true.
 
   etf_mem     = 1
+  vme         = .false.
   lpolar      = .true.
   lphase      = .true.
 
@@ -47,7 +48,6 @@
   phonselfen  = .false.
   a2f         = .false.
 
-
   fsthick     = 2.0 ! eV
   eptemp      = 300 ! K
   degaussw    = 0.1 ! eV

test-suite/epw_soc/epw1.in

@@ -12,6 +12,8 @@
   epwwrite    = .true.
   epwread     = .false.
 
+  vme         = .false.
+
   nbndsub     =  8
   nbndskip    =  10
 

test-suite/epw_soc/epw2.in

@@ -12,6 +12,8 @@
   epwwrite    = .true.
   epwread     = .false.
 
+  vme         = .false.
+
   nbndsub     =  8
   nbndskip    =  10
 
@@ -43,7 +45,6 @@
   phonselfen  = .true.
   a2f         = .true.
 
-
   fsthick     = 20 ! eV
   eptemp      = 500 ! K
   degaussw    = 0.2  ! eV

test-suite/epw_super/epw1.in

@@ -0,0 +1,87 @@
+--
+&inputepw
+  prefix      = 'MgB2',
+  amass(1)    = 24.305,
+  amass(2)    = 10.811
+  outdir      = './'
+
+  ep_coupling = .true.
+  elph        = .true.
+  kmaps       = .true.
+  epbwrite    = .false.
+  epbread     = .false.
+  epwwrite    = .false.
+  epwread     = .true.
+
+  etf_mem     = 0
+  vme         = .false. 
+
+  nbndsub     =  5,
+  nbndskip    =  0
+
+  wannierize  = .false.
+  num_iter    = 500
+  dis_froz_max= 8.8
+  proj(1)     = 'B:pz'
+  proj(2)     = 'f=0.5,1.0,0.5:s'
+  proj(3)     = 'f=0.0,0.5,0.5:s'
+  proj(4)     = 'f=0.5,0.5,0.5:s'
+
+  iverbosity  = 0
+
+  eps_acustic = 5.0    ! Lowest boundary for the 
+  ephwrite    = .true. ! Writes .ephmat files used when wliasberg = .true.
+
+  fsthick     = 20.0  ! eV
+  eptemp      = 300  ! K
+  degaussw    = 0.10 ! eV
+  degaussq    = 0.5 ! meV
+  nqstep      = 500
+
+  eliashberg  = .true.
+
+  liso = .true.
+  limag = .true.
+  lpade = .true.
+  lacon = .true.
+
+  conv_thr_iaxis = 1.0d-2
+  conv_thr_racon = 1.0d-2
+
+  wscut = 0.5   ! eV   Upper limit over frequency integration/summation in the Elisashberg eq
+
+  nstemp   = 1
+  tempsmin = 15.00
+  tempsmax = 20.00
+
+  nsiter   = 500
+
+  muc     = 0.16
+
+  dvscf_dir   = './save'
+ 
+  nk1         = 3
+  nk2         = 3
+  nk3         = 3
+
+  nq1         = 3
+  nq2         = 3
+  nq3         = 3
+ 
+  mp_mesh_k = .true.
+  nkf1        = 6
+  nkf2        = 6
+  nkf3        = 6
+
+  nqf1        = 6
+  nqf2        = 6
+  nqf3        = 6
+ /
+  6 cartesian
+          0.000000000   0.000000000   0.000000000 
+          0.000000000   0.000000000   0.291867841 
+          0.000000000   0.384900179   0.000000000 
+          0.000000000   0.384900179   0.291867841 
+          0.333333333   0.577350269   0.000000000  
+          0.333333333   0.577350269   0.291867841 
+

test-suite/epw_super/epw2.in

@@ -10,11 +10,11 @@
   kmaps       = .false.
   epbwrite    = .true.
   epbread     = .false.
-
-  epwwrite = .true.
-  epwread  = .false.
+  epwwrite    = .true.
+  epwread     = .false.
 
   etf_mem     = 0
+  vme         = .false.
 
   nbndsub     =  5,
   nbndskip    =  0
@@ -29,7 +29,6 @@
 
   iverbosity  = 0
 
-
   eps_acustic = 5.0    ! Lowest boundary for the 
   ephwrite    = .true. ! Writes .ephmat files used when wliasberg = .true.
 

test-suite/epw_trev/epw1.in

@@ -9,10 +9,11 @@
   kmaps       = .false.
   epbwrite    = .true.
   epbread     = .false.
-
   epwwrite    = .true.
   epwread     = .false.
 
+  vme         = .false. 
+
   nbndsub     =  4
   nbndskip    =  0
 
@@ -29,7 +30,6 @@
   phonselfen  = .true.
   a2f         = .true.
 
-
   fsthick     = 15 ! eV
   eptemp      = 2000 ! K
   degaussw    = 1.0 ! eV

test-suite/jobconfig

@@ -64,7 +64,7 @@ inputs_args = ('scf.in', '1'), ('ph.in', '2'), ('scf_epw.in', '1'), ('nscf_epw.i
 
 [epw_trev/]
 program = EPW
-inputs_args = ('scf.in', '1'), ('ph.in', '2'), ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw.in', '3')
+inputs_args = ('scf.in', '1'), ('ph.in', '2'), ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw1.in', '3')
 
 [epw_trev_uspp/]
 program = EPW
@@ -72,15 +72,15 @@ inputs_args = ('scf.in', '1'), ('ph.in', '2'), ('scf.in', '1'), ('nscf.in', '1')
 
 [epw_metal/]
 program = EPW
-inputs_args = ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw.in', '3'), ('epw2.in', '3')
+inputs_args = ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw1.in', '3'), ('epw2.in', '3')
 
 [epw_soc/]
 program = EPW
-inputs_args = ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw.in', '3'), ('epw2.in', '3')
+inputs_args = ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw1.in', '3'), ('epw2.in', '3')
 
 [epw_super/]
 program = EPW
-inputs_args = ('scf.in', '1'), ('ph.in', '2'), ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw_iso.in', '3'), ('epw_aniso.in', '3')
+inputs_args = ('scf.in', '1'), ('ph.in', '2'), ('scf_epw.in', '1'), ('nscf_epw.in', '1'), ('epw1.in', '3'), ('epw2.in', '3')
 
 [epw_polar/]
 program = EPW