<div dir="ltr"><div><div><div><div><div><div>Dear All,<br><br></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></div>When I use MatCreateFFT it leads to a crash. My test code is:<br><br>#include "petsc/finclude/petscsysdef.h"<br>#include "petsc/finclude/petscmatdef.h"<br><br>program test<br>    use petscsys<br>    use petscmat<br>    implicit none<br>    Mat :: A<br>    integer :: ndim,vdim(3)<br>    integer :: ierr<br><br>    ndim = 3<br>    vdim(1) = 30<br>    vdim(2) = 30<br>    vdim(3) = 30<br><br>    call PetscInitialize(PETSC_NULL_CHARACTER,ierr)<br><br>    call MatCreateFFT(PETSC_COMM_WORLD,ndim,vdim,MATFFTW,A,ierr)<br><br>    call PetscFinalize(ierr)<br><br>end program<br><br></div>When MatCreateFFT is executed program crushes and gives:<br><br>"[Minas-Thirith:17525] *** An error occurred in MPI_Comm_size<br>[Minas-Thirith:17525] *** reported by process [3784507393,2]<br>[Minas-Thirith:17525] *** on communicator MPI_COMM_WORLD<br>[Minas-Thirith:17525] *** MPI_ERR_COMM: invalid communicator<br>[Minas-Thirith:17525] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort,<br>[Minas-Thirith:17525] ***    and potentially your MPI job)<br>[Minas-Thirith:17520] 5 more processes have sent help message help-mpi-errors.txt / mpi_errors_are_fatal<br>[Minas-Thirith:17520] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages"<br><br></div>What am I doing wrong?<br><br></div>Best regards,<br></div>K.G.<br><div><div><div><br> </div></div></div></div>