<div dir="ltr">Dear Paolo and Lorenzo,<div><br></div><div style>Thank you for such a detailed explanation!</div><div style><br></div><div style>Best regards,</div><div style>Andrei</div></div><div class="gmail_extra"><br><br>
<div class="gmail_quote">On Fri, May 17, 2013 at 3:59 AM, 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">
<div class="im">On Fri, 2013-05-17 at 01:18 +0800, Fang Liu wrote:<br>
<br>
> That is, I use a 25*25*25 3D FFT mesh. The original 3d array g3d =<br>
> (1.0_DP, -1.0_DP).<br>
><br>
> First I do invfft and then fwfft. Finally I don't get the original 3d<br>
> array, but a 3d array g3dnew in which<br>
> g3dnew(1),g3dnew(2),...,g3dnew(3025)= (1.0_DP, -1.0_DP)<br>
> while g3dnew(3026),....,g3dnew(15625)=0.0.<br>
<br>
> But if I do fwfft first and then invfft, I can finally get the<br>
> original 3d array.<br>
<br>
</div>The 3D FFTs in QE are tailored for the "dual-space" algorithm<br>
as used in plane-wave codes. If you specify "Wave", the code<br>
will assume that there are no G components beyond a sphere<br>
of radius G^2 < Ecut. For instance, when performing the G->R<br>
FT along z, many vectors will be ignored, because they are<br>
assumed to contain zeros. See Lorenzo's post for a more detailed<br>
explanation.<br>
<div class="HOEnZb"><div class="h5"><br>
P.<br>
<br>
--<br>
 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">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</a><br>
<br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div>