[petsc-users] MatCreateFFT in Fortran

[Krzysztof Gawarecki]

Dear All,

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.
When I use MatCreateFFT it leads to a crash. My test code is:

#include "petsc/finclude/petscsysdef.h"
#include "petsc/finclude/petscmatdef.h"

program test
    use petscsys
    use petscmat
    implicit none
    Mat :: A
    integer :: ndim,vdim(3)
    integer :: ierr

    ndim = 3
    vdim(1) = 30
    vdim(2) = 30
    vdim(3) = 30

    call PetscInitialize(PETSC_NULL_CHARACTER,ierr)

    call MatCreateFFT(PETSC_COMM_WORLD,ndim,vdim,MATFFTW,A,ierr)

    call PetscFinalize(ierr)

end program

When MatCreateFFT is executed program crushes and gives:

"[Minas-Thirith:17525] *** An error occurred in MPI_Comm_size
[Minas-Thirith:17525] *** reported by process [3784507393,2]
[Minas-Thirith:17525] *** on communicator MPI_COMM_WORLD
[Minas-Thirith:17525] *** MPI_ERR_COMM: invalid communicator
[Minas-Thirith:17525] *** MPI_ERRORS_ARE_FATAL (processes in this
communicator will now abort,
[Minas-Thirith:17525] ***    and potentially your MPI job)
[Minas-Thirith:17520] 5 more processes have sent help message
help-mpi-errors.txt / mpi_errors_are_fatal
[Minas-Thirith:17520] Set MCA parameter "orte_base_help_aggregate" to 0 to
see all help / error messages"

What am I doing wrong?

Best regards,
K.G.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160819/d88a6d97/attachment.html>