<div dir="ltr">On Thu, May 12, 2016 at 4:14 PM, dario rocca <span dir="ltr"><<a href="mailto:roccad@gmail.com" target="_blank">roccad@gmail.com</a>></span> wrote:<br><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>In order to menage the g vectors for each k-point, you can use the arrays ngk (number of G vectors for each k-point) and igk (index of G corresponding to a given index of k+G; basically an index that allows you to identify the G vectors corresponding to a given k and order them).<br></div></div></div></div></div></blockquote><div><br></div><div style="font-family:arial,helvetica,sans-serif;font-size:small">just a note: "igk" is going to be replaced by another variable, "igk_k", that contains indices of k+G for all k-points. Indices "igk" used to be read from file in each loop over k-vectors, to save memory and computer time. This may have been a good idea many years ago, but it isn't any longer.<br><br></div><div style="font-family:arial,helvetica,sans-serif;font-size:small">Paolo<br></div></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div></div>