<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Aug 20, 2014 at 3:52 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class=""><br>
On Aug 20, 2014, at 3:34 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
<br>
> The blurb should also answer the question, Why do we need another FFT library?<br>
<br>
</div>  Because it actually does multidimensional FFTs in parallel?<br></blockquote><div><br></div><div>For something this old and established, it might be important to say that it does blah</div><div>blah blah that Spiral, FFTW, etc.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  In theory, with this one could write very efficient parallel 3d Poisson solvers in PETSc for boxes, which is an important special case that PETSc does not currently support.</blockquote><div><br></div><div>Wouldn't you just use MG? I have other uses for FFT, including DFT codes where it makes a lot of sense.</div>
<div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
><br>
>    Matt<br>
><br>
><br>
> On Wed, Aug 20, 2014 at 1:53 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
><br>
> From: Dmitry Pekurovsky <a href="mailto:dmitry@sdsc.edu">dmitry@sdsc.edu</a><br>
> Date: August 12, 2014<br>
> Subject: Library for spectral transforms in 3D for parallel machines<br>
><br>
> P3DFFT is an open-source numerical library providing highly<br>
> scalable implementation of 3D spectral transforms namely Fast<br>
> Fourier Transform (FFT), with an option to combine it with<br>
> cosine/sine/Chebyshev/empty transform in the third dimension. (The<br>
> empty transform allows the user to substitute their own custom<br>
> transform in the third dimension. This can be useful in<br>
> applications such as inhomogeneous wall bounded turbulence.) P3DFFT<br>
> implements 2D domain decomposition which allows it to overcome a<br>
> scalability restriction inherent in 1D decomposition. This approach<br>
> has shown good scalability up to 131,072 cores.<br>
><br>
> A new version of P3DFFT 2.7.1 is now available. The project Home<br>
> Page is <a href="http://code/google.com/p/p3dfft" target="_blank">http://code/google.com/p/p3dfft</a> where instructions for<br>
> obtaining the source code are provided. Installation instructions<br>
> and a User Guide are also available.<br>
><br>
> P3DFFT features include real-to-complex and complex-to-real<br>
> transforms, in-place transforms, pruned transforms (with less than<br>
> full input or output), and multi-variable transforms. The package<br>
> includes example programs in Fortran and C. This is a project in<br>
> active development, with a user mailing list, a wiki page and a<br>
> version control system. P3DFFT is considered community software and<br>
> is being installed in public space at many supercomputer centers.<br>
> Contributions and feedback from users are welcome.<br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
> -- Norbert Wiener<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>