<div dir="ltr"><div>There is actually a stupid bug in the computation of the maximum number of plane waves (only for hybrid functionals and ultrasoft PPs). Please replace routine exx_n_plane_waves in file PW/src/exx.f90 with the attached one, report if it doesn't work<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 1, 2015 at 9:52 PM, Xiaoming Wang <span dir="ltr"><<a href="mailto:xw213@rci.rutgers.edu" target="_blank">xw213@rci.rutgers.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div link="#0563C1" vlink="#954F72" lang="EN-US"><div><p class="MsoNormal">Hello qe users and developers,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I am trying the new feature as in QE-5.2.0, which implement USPP with hybrid functional calculations. However, as for my calculations, I found errors in output which depend on number of cores I used.  The error reads “ Error in routine gk_sort (1):  array gk out-of-bounds” when I use 16 cores as ‘mpirun –np 16 pw.x’, but it works for 4 cores. However, for NCPP, I found no such problems, at least for my system. Does anyone know what the problem is? And the following is my input:<u></u><u></u></p><p class="MsoNormal">----------------------------------------------------------------------------------------------------------<u></u><u></u></p><p class="MsoNormal">&control<u></u><u></u></p><p class="MsoNormal">    calculation = 'scf'<u></u><u></u></p><p class="MsoNormal">    restart_mode='from_scratch',<u></u><u></u></p><p class="MsoNormal">    prefix='bp',<u></u><u></u></p><p class="MsoNormal">    tstress = .true.<u></u><u></u></p><p class="MsoNormal">    tprnfor = .true.<u></u><u></u></p><p class="MsoNormal">    etot_conv_thr = 1.0D-6<u></u><u></u></p><p class="MsoNormal">    forc_conv_thr = 1.0D-4<u></u><u></u></p><p class="MsoNormal">    pseudo_dir = ' /pseudo',<u></u><u></u></p><p class="MsoNormal">    outdir='./tmp'<u></u><u></u></p><p class="MsoNormal">    wf_collect=.true.<u></u><u></u></p><p class="MsoNormal"> /<u></u><u></u></p><p class="MsoNormal"> &system<u></u><u></u></p><p class="MsoNormal">    ibrav=8<u></u><u></u></p><p class="MsoNormal">    a=3.33826<u></u><u></u></p><p class="MsoNormal">    b=4.59076<u></u><u></u></p><p class="MsoNormal">    c=13.0<u></u><u></u></p><p class="MsoNormal">    nat=4<u></u><u></u></p><p class="MsoNormal">    ntyp=1<u></u><u></u></p><p class="MsoNormal">    ecutwfc=30<u></u><u></u></p><p class="MsoNormal">    ecutrho=200<u></u><u></u></p><p class="MsoNormal">    occupations='fixed'<u></u><u></u></p><p class="MsoNormal">    nosym=.true.<u></u><u></u></p><p class="MsoNormal">    noinv = .true.<u></u><u></u></p><p class="MsoNormal">    input_dft = 'HSE'<u></u><u></u></p><p class="MsoNormal">    nqx1=5, nqx2=4, nqx3=1<u></u><u></u></p><p class="MsoNormal"> /<u></u><u></u></p><p class="MsoNormal"> &electrons<u></u><u></u></p><p class="MsoNormal">    diagonalization='david'<u></u><u></u></p><p class="MsoNormal">    mixing_mode = 'plain'<u></u><u></u></p><p class="MsoNormal">    mixing_beta = 0.7<u></u><u></u></p><p class="MsoNormal">    conv_thr =  1.0d-8<u></u><u></u></p><p class="MsoNormal"> /<u></u><u></u></p><p class="MsoNormal">ATOMIC_SPECIES<u></u><u></u></p><p class="MsoNormal"> P  30.9738  p_pbe_v1.uspp.F.UPF<u></u><u></u></p><p class="MsoNormal">K_POINTS automatic<u></u><u></u></p><p class="MsoNormal">10 8 1 0 0 0<u></u><u></u></p><p class="MsoNormal">ATOMIC_POSITIONS (crystal)<u></u><u></u></p><p class="MsoNormal">P       0.500000      0.914274      0.002447<u></u><u></u></p><p class="MsoNormal">P       0.000000      0.585727      0.002447<u></u><u></u></p><p class="MsoNormal">P       0.000000      0.414273      0.167618<u></u><u></u></p><p class="MsoNormal">P       0.500000      0.085726      0.167617<u></u><u></u></p><p class="MsoNormal">---------------------------------------------------------------------------------------------------------------------<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Best,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Xiaoming Wang<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">IAMDN, Rutgers<u></u><u></u></p></div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div>