Dear Prof. Hütter,<div><br /></div><div>Thanks for the suggestion I'll switch to FFTW3 then, but I do still have one question about using FFTW3. Why do I get the error
<span style="color: rgb(80, 0, 80);">"the FFT in the x direction is not allowed n01 dimension 175" when switching from FFTSG/EXTENDED to FFTW3/EXTENDED, while I get an "</span><span style="color: rgb(80, 0, 80);">Index to radix array not found.</span><span style="color: rgb(80, 0, 80);">" error when switching of EXTENDED_FFT_LENGHTS?</span></div><div><font color="#500050"><br /></font></div><div><font color="#500050">Kind regards,</font></div><div><font color="#500050">Léon<br /></font><div><br /></div><div><br /><br /></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 4 October 2023 at 09:35:36 UTC+2 Jürg Hutter wrote:<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>there is an upper limit for the FFT length for the FFTSG library (1024).
<br>I would switch to FFTW3 where no such limit exists.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Sent: Tuesday, October 3, 2023 2:06 PM
<br>To: cp2k
<br>Subject: [CP2K:19284] FFT problems with large unit cell
<br>
<br>Dear all,
<br>
<br>I am running some tests with long chain adsorbates on a TiO2 slab and I'm running into some problems with the FFT libraries. I initially had the FFTSG set as preferred library, because the manual stated that "The use of PREFERRED_FFT_LIBRARY FFTSG is required" when using the WAVELET solver. However, I just read in this post<<a href="https://groups.google.com/g/cp2k/c/QEveWI7UqV0" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/g/cp2k/c/QEveWI7UqV0&source=gmail&ust=1696492184602000&usg=AOvVaw0dPVWSCGau9ZdKyGpIhUsm">https://groups.google.com/g/cp2k/c/QEveWI7UqV0</a>> from Prof. Hütter that you can use the FFTW3 library in combination with the WAVELET solver and I already read several times that using this library is generally faster compared to the internal FFTSG one. I also already had to enable the EXTENDED_FFT_LENGTHS option for calculations with smaller (10x10x30) cells as I ran into some "Index to radix array not found."-errors.
<br>
<br>Now, with these long adsorbates I have a cell of around 10 x 15 x 60 Å and the combination of FFTSG/EXTENDED_FFT_LENGTHS again gives "Index to radix array not found."-errors, while switching to FFTW3/EXTENDED_FFT_LENGHTS says that "the FFT in the x direction is not allowed n01 dimension 175".
<br>
<br>So I have two questions:
<br>1) Is it possible to calculate such a large cell without getting these FFT errors or is it just too large?
<br>2) Should I switch to FFTW3/EXTENDED_FFT_LENGTHS to speed up my calculations or is it fine to continue with FFTSG/EXTENDED_FFT_LENGTHS?
<br>
<br>Hopefully someone can help me out with this problem.
<br>
<br>All the best,
<br>Léon
<br>
<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">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/38650355-0f22-4606-82bb-16a670bc8619n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/38650355-0f22-4606-82bb-16a670bc8619n%2540googlegroups.com&source=gmail&ust=1696492184602000&usg=AOvVaw1AOeItcWi0tR95DDuFTFpz">https://groups.google.com/d/msgid/cp2k/38650355-0f22-4606-82bb-16a670bc8619n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/38650355-0f22-4606-82bb-16a670bc8619n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/38650355-0f22-4606-82bb-16a670bc8619n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1696492184602000&usg=AOvVaw1hUVjGICax7K8rTCL4-0SJ">https://groups.google.com/d/msgid/cp2k/38650355-0f22-4606-82bb-16a670bc8619n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></blockquote></div>
<p></p>
-- <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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/acbeac79-87b6-4a21-8dfe-e4b978507655n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/acbeac79-87b6-4a21-8dfe-e4b978507655n%40googlegroups.com</a>.<br />