[petsc-dev] Fwd: [petsc-users] MatCreateFFT in Fortran

[Barry Smith]

   PETSc developers,

  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.

   Barry


> Begin forwarded message:
> 
> From: Krzysztof Gawarecki <dkg2140 at gmail.com>
> Subject: [petsc-users] MatCreateFFT in Fortran
> Date: August 18, 2016 at 5:15:16 PM CDT
> To: <petsc-users at mcs.anl.gov>
> 
> 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-dev/attachments/20160818/eb047cf6/attachment.html>