<div dir="ltr">Hi Marton,<div><br></div><div>Yes, I am aware of that. I think there is a issue with nspin=2, nspin=1 is totally fine.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 6, 2017 at 11:15 PM, Marton <span dir="ltr"><<a href="mailto:vormar@gmail.com" target="_blank">vormar@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="auto"><div dir="ltr">Hi Jia,<div><br></div><div>I don't think there is a bug there. If you go a few lines above in the same file, you find, for example:</div><div><span style="font-variant-ligatures:no-common-ligatures;font-family:menlo;font-size:11px"><br></span></div><div><span style="font-variant-ligatures:no-common-ligatures;font-family:menlo;font-size:11px"> </span><span style="font-variant-ligatures:no-common-ligatures;font-family:menlo;font-size:11px">IF ( nspin == 1) full_occ = 2.0d0</span></div><div><br></div><div>which means that <font face="monospace, monospace">full_occ</font> takes care of the occupations.<br></div><div><br></div>HTH<div><br><div>Marton Voros</div><div><br></div><div>--</div><div>Aneesur Rahman Fellow</div><div>Materials Science Division</div><div>Argonne National Laboratory</div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Feb 6, 2017 at 3:27 PM, Jia Chen <span dir="ltr"><<a href="mailto:jiachenchem@gmail.com" target="_blank">jiachenchem@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Dear All,<div><br></div><div>I am suspicious that a bud exit in epsilon.f90 qe-6.0. The check of nbnd:</div><div>==============================<wbr>=========================</div><div>IF ( REAL(nbnd,DP)*full_occ <= nelec ) CALL errore('epsilon', 'bad band number', 1)<br></div><div>==============================<wbr>=========================</div><div>I think it should be:</div><div>REAL(nbnd,DP)*full_occ <= nelec/2.0<br></div><div><br></div><div>Let me know if I missed something...</div><div><br></div><div>Cheers</div><span class="m_7557928156303730596m_-2977406803739256980m_1995206756938415840m_-227488833213356m_6358309315408987227HOEnZb"><font color="#888888"><div>Jia</div></font></span></div>
<br></div></div>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listi<wbr>nfo/pw_forum</a><br></blockquote></div><br></div>
</div>
<br>______________________________<wbr>_________________<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/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div>