<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><br class=""></div><div class=""> PETSc developers,</div><div class=""><br class=""></div> We are getting MatCreateFFT questions sometimes, does anyone have a pot of money to pay someone (or plenty of free time) to refactor that stuff to be usable? I dislike even suggesting people use the current CF as it stands.<div class=""><br class=""></div><div class=""> Barry</div><div class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">Begin forwarded message:</div><br class="Apple-interchange-newline"><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">From: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class="">Krzysztof Gawarecki <<a href="mailto:dkg2140@gmail.com" class="">dkg2140@gmail.com</a>><br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">Subject: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class=""><b class="">[petsc-users] MatCreateFFT in Fortran</b><br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">Date: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class="">August 18, 2016 at 5:15:16 PM CDT<br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">To: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class=""><<a href="mailto:petsc-users@mcs.anl.gov" class="">petsc-users@mcs.anl.gov</a>><br class=""></span></div><br class=""><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class="">Dear All,<br class=""><br class=""></div>I'm trying to implement Fast Fourier Transform (via FFTW interface in PETSC) to my program. Unfortunately, I have not found any example/tutorial for fortran. <br class=""></div>When I use MatCreateFFT it leads to a crash. My test code is:<br class=""><br class="">#include "petsc/finclude/petscsysdef.h"<br class="">#include "petsc/finclude/petscmatdef.h"<br class=""><br class="">program test<br class=""> use petscsys<br class=""> use petscmat<br class=""> implicit none<br class=""> Mat :: A<br class=""> integer :: ndim,vdim(3)<br class=""> integer :: ierr<br class=""><br class=""> ndim = 3<br class=""> vdim(1) = 30<br class=""> vdim(2) = 30<br class=""> vdim(3) = 30<br class=""><br class=""> call PetscInitialize(PETSC_NULL_CHARACTER,ierr)<br class=""><br class=""> call MatCreateFFT(PETSC_COMM_WORLD,ndim,vdim,MATFFTW,A,ierr)<br class=""><br class=""> call PetscFinalize(ierr)<br class=""><br class="">end program<br class=""><br class=""></div>When MatCreateFFT is executed program crushes and gives:<br class=""><br class="">"[Minas-Thirith:17525] *** An error occurred in MPI_Comm_size<br class="">[Minas-Thirith:17525] *** reported by process [3784507393,2]<br class="">[Minas-Thirith:17525] *** on communicator MPI_COMM_WORLD<br class="">[Minas-Thirith:17525] *** MPI_ERR_COMM: invalid communicator<br class="">[Minas-Thirith:17525] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort,<br class="">[Minas-Thirith:17525] *** and potentially your MPI job)<br class="">[Minas-Thirith:17520] 5 more processes have sent help message help-mpi-errors.txt / mpi_errors_are_fatal<br class="">[Minas-Thirith:17520] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages"<br class=""><br class=""></div>What am I doing wrong?<br class=""><br class=""></div>Best regards,<br class=""></div>K.G.<br class=""><div class=""><div class=""><div class=""><br class=""> </div></div></div></div>
</div></blockquote></div><br class=""></div></body></html>