[petsc-users] MatCreateFFT in Fortran

[Krzysztof Gawarecki]

The problem was related to fortran wrapper:
src/mat/impls/fft/ftn-custom/zfftf.c

After changing:
*ierr = MatCreateFFT(*comm,*ndim,dim,mattype,A);
to
*ierr = MatCreateFFT(MPI_Comm_f2c(*(MPI_Fint*)&*comm),*ndim,dim,mattype,A);

FFTW works fine in fortran.

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