<div>Thank you Prasenjit, I tested it and it seems correct. I also fixeda couple of bugs in plotbands.</div><br><div>Do you want to do a pull request yourself for the file ? (In that case I'ld ask you to remove the debug lines) or should I take care to do it?</div><br><div>thank you again</div><br><div><signature id="initial">Dr. Lorenzo Paulatto<br>IdR @ IMPMC - CNRS UMR 7590 & Sorbonne Université<br>phone: +33 (0)1 442 79822 / skype: paulatz<br><a href="http://www.impmc.upmc.fr/~paulatto/">http://www.impmc.upmc.fr/~paulatto/</a> - <a href="https://anharmonic.github.io/">https://anharmonic.github.io/</a><br>23-24/423 B115, 4 place Jussieu 75252 Paris CX 05</signature></div><div class="gmail_quote_attribution">On Feb 15 2022, at 7:29 am, Prasenjit Ghosh <prasenjit.jnc@gmail.com> wrote:</div><blockquote><div><div>Dear Lorenzo and others,</div><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><div>        k(    1) = (  -0.2500000   0.2500000   0.2500000), wk =   0.5000000</div><div>        k(    2) = (   0.2500000  -0.2500000   0.7500000), wk =   1.5000000</div><div>        k(    3) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000</div><div>        k(    4) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000</div><div>        k(    5) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000</div><div>        k(    6) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000</div><div>        k(    7) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000</div><div>        k(    8) = (   0.0000000   0.0000000   0.0000000), wk =   0.0000000</div><div>        k(    9) = (   0.5000000   0.5000000   0.5000000), wk =   0.0000000</div></div><div><br></div><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.</div></div><div><br></div><div><div>I was working with QE-6.8. The same problem I think persists also in 7.0.</div></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><div>Prasenjit</div></div><div><br></div><div><div><br></div><div>With regards,</div><div>Prasenjit</div><div>--</div><div><div><div>PRASENJIT GHOSH,</div><div>IISER Pune,</div><div>Dr. Homi Bhabha Road, Pashan</div><div>Pune, Maharashtra 411008, India</div><br><div>Phone: +91 (20) 2590 8203</div><div>Fax: +91 (20) 2589 9790</div></div></div></div></div><div>_______________________________________________</div><div>Quantum ESPRESSO is supported by MaX (www.max-centre.eu)</div><div>users mailing list users@lists.quantum-espresso.org</div><div>https://lists.quantum-espresso.org/mailman/listinfo/users</div></blockquote><img class="mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border:0; width:0; height:0;" src="https://link.getmailspring.com/open/2F55062C-7DC8-47DD-9EFF-D793ED75DA16@getmailspring.com?me=ca1a67e4&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D">