<div dir="ltr">This happens in the LR code that I develop. I tried to generate an example with QE and failed and now I realized what happened. According to the recent change I replaced all h_psiq with h_psi and s_psi calls. Because I wrote a wrapper, I accidently called it with h_psi(lda, <b>lda</b>, m, psi, hpsi). Because the second dimension was to large, it looped over components where igk_k is 0, so it is only a bug in my code and not in QE.<div><br></div><div>Might I suggest to add a check to h_psi/s_psi that igk_k(n) is nonzero to avoid someone else running into the same problem?</div><div><br></div><div>Martin</div><div><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jul 1, 2016 at 5:37 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Please provide a test where this happens<br>
<br>
Paolo<br>
<div><div class="h5"><br>
On Fri, Jul 1, 2016 at 6:14 PM, Martin Schlipf <<a href="mailto:martin.schlipf@gmail.com">martin.schlipf@gmail.com</a>> wrote:<br>
> I recently ran the code with array boundary checking and it crashed in<br>
> vloc_psi.f90. In particular the array igk_k has components that are 0 and<br>
> then accessing nls of these elements is not within the bounds of nls (which<br>
> is 1:ngms).<br>
><br>
> Martin<br>
><br>
> Martin Schlipf<br>
> Department of Materials<br>
> University of Oxford<br>
> Parks Road<br>
> Oxford OX1 3PH, UK<br>
><br>
</div></div>> _______________________________________________<br>
> Q-e-developers mailing list<br>
> <a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
> <a href="http://qe-forge.org/mailman/listinfo/q-e-developers" rel="noreferrer" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
><br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
--<br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</a><br>
_______________________________________________<br>
Q-e-developers mailing list<br>
<a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
<a href="http://qe-forge.org/mailman/listinfo/q-e-developers" rel="noreferrer" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
</font></span></blockquote></div><div><br></div>
</div></div></div>