<div>Hi!</div><div><br></div><div>Excellent. Thank you for the reply!</div><div><br></div><div>Follow-up question: Some of the slabs I will be considering are quite large, so I will need to set EXTENDED_FFT_LENGTHS to True. This leads to a shift in total energy, at least with FFTW3, as expected per the documentation. Could this lead to inconsistencies in results with the wavelet solver using FFTW3?</div><div><br></div><div>All the best,</div><div>Nicklas<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">torsdag 11 mars 2021 kl. 12:59:58 UTC+1 skrev jgh:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>I think the manual is not correct in this case. The wavelet solver
<br>uses its own FFT library. The wavelet FFT library is the same as the
<br>FFTSG library in CP2K and they share the allowed FFT roots. 
<br>That is the reason behind the statement in the manual.
<br>However, we restricted the allowed roots when using the FFTW3 to more
<br>or less the same as in FFTSG (not exactly). But I think
<br>they will only diverge for very large grids.
<br>So, you can use FFTW3 in CP2K with the wavelet solver. Should there
<br>be an inconsistency of grid sizes, this will be checked in the wavelet
<br>solver and the code will stop with an error message.
<br>
<br>regards
<br>
<br>Juerg Hutter
<br>--------------------------------------------------------------
<br>Juerg Hutter                         Phone : <a href="tel:+41%2044%20635%2044%2091" value="+41446354491" target="_blank" rel="nofollow">++41 44 635 4491</a>
<br>Institut für Chemie C                FAX   : <a href="tel:+41%2044%20635%2068%2038" value="+41446356838" target="_blank" rel="nofollow">++41 44 635 6838</a>
<br>Universität Zürich                   E-mail: <a href data-email-masked rel="nofollow">h...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>---------------------------------------------------------------
<br>
<br>-----<a href data-email-masked rel="nofollow">c...@googlegroups.com</a> wrote: -----
<br>To: "cp2k" <<a href data-email-masked rel="nofollow">c...@googlegroups.com</a>>
<br>From: "Nicklas Österbacka" 
<br>Sent by: <a href data-email-masked rel="nofollow">c...@googlegroups.com</a>
<br>Date: 03/10/2021 11:42AM
<br>Subject: [CP2K:14920] FFT library for wavelet Poisson solver
<br>
<br>Hi!
<br>
<br>I'm working towards doing some calculations on adsorption energies on slabs using CP2K 6.1 and would prefer to use non-periodic boundary conditions in the surface normal direction. As I have little experience with aperiodic DFT I ran some small tests to compare with periodic calculations.
<br>
<br>The wavelet Poisson solver requires the use of the FFTSG FFT library according to the manual, but using FFTW3 instead yields virtually the same result for the geometry optimization of an oxygen molecule in a 20x20x20 Å box with no periodicity in any direction. (I realize that the box is much larger than necessary.) The energy difference between the two solvers is about 2e-13 Ha, while the difference in O-O bond length is about 3e-9 Å.
<br>
<br>As the calculation using FFTW3 ran about 20 % faster it seems preferable to using FFTSG, but I want to make sure that CP2K isn't doing anything unexpected when combining FFTW3 with the wavelet Poisson solver. Is the manual incorrect on the FFTSG requirement, or is it perhaps only required for XZ periodicity?
<br>
<br>Any insight would be appreciated!
<br>Nicklas
<br>  
<br>  -- 
<br> You received this message because you are subscribed to the Google Groups "cp2k" group.
<br> To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>.
<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/5efc1d76-1d28-4d02-bbde-eecc987f4e20n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=sv&q=https://groups.google.com/d/msgid/cp2k/5efc1d76-1d28-4d02-bbde-eecc987f4e20n%2540googlegroups.com&source=gmail&ust=1615567669113000&usg=AFQjCNHBKv0LRRpi1Y2VJKfbqykg9EhMHg">https://groups.google.com/d/msgid/cp2k/5efc1d76-1d28-4d02-bbde-eecc987f4e20n%40googlegroups.com</a>.
<br> 
<br></blockquote></div>