<div dir="ltr">On Mon, Jun 22, 2015 at 11:29 AM, Jiqiang Li <span dir="ltr"><<a href="mailto:jqli14@fudan.edu.cn" target="_blank">jqli14@fudan.edu.cn</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div> <br></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div>(1)Do the elements of array for total local potential (variable <b>vrs</b>) or charge density (variable <span style="background-color:rgb(255,255,255)"><b>rho</b></span>) in real space keep the same order in different  calculation when the volume and shape of supercell remain unchanged?</div></div></blockquote><div><br></div><div>Not sure what you mean, but the logic of real-space arrays is always the same:<br> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div>(2) If do, does the element in array <b>vrs</b>(nr1*nr2*(n3-1)+nr1*(n2-1)+n1) represent the position [x(n1), y(n2),z(n3)] in real space?</div></div></blockquote><div><br></div><div>in the serial case: index (nr1x*nr2x*(n3-1)+nr1x*(n2-1)+n1) for an array of dimension (nr1x*nr2x*nr3x) corresponds to index (n1,n2,n3) of the three-dimensional array of dimensions (nr1x,nr2x,nr3x).<br>Index (nr1x*nr2x*(n3-1)+nr1x*(n2-1)+n1) correspond to real-space point n1*at(:,1)+n2*at(:,1)+n3*at(:,3).<br></div><div>IN the parallel case, arrays are split across processors along the third dimension (nr3) but the logic is the same<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div>(3) Is it correct that \sum_i {<b>rho</b>(i)}*volume/N is not exactly equal to the total charge?  N is the size of array <b>rho</b>.</div></div></blockquote><div><br></div><div>No, it must be exactly equal to the total charge, and N=nr1*nr2*nr3, which may differ from the size of the array rho<br> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div>(4)What's the difference between variables <b>nr1</b> and <b>nr1x</b>?</div>
<div></div></div></blockquote><div><br></div><div>see above: one is the true dimension of the FFT transform, one is the physical dimension of the array. They may differ due to efficiency reasons.<br><br></div><div>Paolo<br><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div>Thank you for your consideration!</div><div><br></div><div>Best regards!</div><span class=""><font color="#888888"><div>Jiqiang Li</div><hr style="width:210px;min-height:1px" color="#b5c4df" align="left" size="1">
<div><span></span></div>
</font></span></div><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>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div></div>