<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On 28 September 2017 at 19:54, Jacopo Simoni <span dir="ltr"><<a href="mailto:simonij@tcd.ie" target="_blank">simonij@tcd.ie</a>></span> wrote:<br><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 dir="ltr"><div><div><div><div><div></div><br><span class="gmail-m_858173876573312787gmail-im">r_ijk(:) = (i-1)/dffts%nr1x * a_1(:) + (j-1)/dffts%nr2x * a_2(:) + (k-1)/dffts%nr3x * a_3(:) ?<br></span></div><div><span class="gmail-m_858173876573312787gmail-im">with a_1, a_2, a_3 lattice parameters (I have a simple cubic cell), this is what I have found looking around in the documentation. In particular this is what it is written in the developer manual. The issue is that in the wfck2r program the real space grid does not appear explicitly so I am wondering if the previous expression for the wave function index is still usable.<br></span></div><span class="gmail-m_858173876573312787gmail-im"></span></div></div></div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></blockquote><div><br><br></div><div>Almost, you have this kind of code in add_efield, for the hard grid but applies jus tthe same to the smooth grid. Be careful that nr1x, nr2x and nr3x are the dimensions of the FFT which can be (but almost never are) bigger than the physical dimension of the cell, which are actually nr1, nr2 and nr3. So you should use the latter, and eventually throw away the points which fall outside the physical range.<br><br></div><div>Also, be careful when doit it in parallel, each CPU has a share of the planes of charge along Z, hence he box dimension along that direction is nr3, but the actual number of points on a CPU is dfftp%my_nr3p and the index of the first plane is dfftp%my_i0r3p.<br><br></div><div>hth<br></div><div><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 dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div></div></div></div></div></div></div></div></div></div></div></div></div></div></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><br clear="all"><br>-- <br><div class="gmail_signature">Lorenzo Paulatto - Paris<br></div>
</div></div>