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

[Matthew Knepley]

On Thu, Aug 18, 2016 at 8:55 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>    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.
>

I am considering it. I have a student doing a summer project with it. If
that goes well, I will have
them do it.

  Matt


>    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.
>
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20160818/a0919174/attachment.html>