<div dir="ltr">Dear Lorenzo and others,<div><br></div><div>I think there is a bug in the above subroutine when either nqx1, nqx2 and nqx3, or all of them are greater than 1. The example does not work with the original version of the code because a lot of k-points with coordinates (0,0,0 are generated as shown below:</div><div><br></div><div>        k(    1) = (  -0.2500000   0.2500000   0.2500000), wk =   0.5000000<br>        k(    2) = (   0.2500000  -0.2500000   0.7500000), wk =   1.5000000<br>        k(    3) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000<br>        k(    4) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000<br>        k(    5) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000<br>        k(    6) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000<br>        k(    7) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000<br>        k(    8) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000<br>        k(    9) = (   0.5000000   0.5000000   0.5000000), wk =   0.0000000<br></div><div><br></div><div>In the above list, k-points #3 to #7 are all (0, 0,0). Because of this, when running the bands.x code, it complains and stops.<br></div><div><br></div><div> I was working with QE-6.8. The same problem I think persists also in 7.0.<br></div><div> </div><div>In line 89 of the subroutine, the value of the counter 'iq' seems to be wrongly initialized. With this the array elements are wrongly assigned in xk(:,ik). The correct one should be:</div><div><br></div><div>iq=nkstot_old</div><div><br></div><div>Further the array elements of 'wk', corresponding to the additional k-points also need to be set to zero. This can be done right after line 98 by adding the following line:</div><div>wk(iq)=0.0d0</div><div><br></div><div>The modified subroutine is attached.</div><div><br></div><div>With regards,</div><div>Prasenjit<br></div><div><br></div><div><div><br></div><div>With regards,</div><div>Prasenjit</div>-- <br><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr">PRASENJIT GHOSH,<br>IISER Pune,<br>Dr. Homi Bhabha Road, Pashan<br>Pune, Maharashtra 411008, India<br><br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790</div></div></div></div>