<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Aug 18, 2016 at 8:55 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div><br></div><div>   PETSc developers,</div><div><br></div>  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.</div></blockquote><div><br></div><div>I am considering it. I have a student doing a summer project with it. If that goes well, I will have</div><div>them do it.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>   Barry</div><div><br><div><br><blockquote type="cite"><div>Begin forwarded message:</div><br><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif;color:rgba(0,0,0,1.0)"><b>From: </b></span><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif">Krzysztof Gawarecki <<a href="mailto:dkg2140@gmail.com" target="_blank">dkg2140@gmail.com</a>><br></span></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif;color:rgba(0,0,0,1.0)"><b>Subject: </b></span><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif"><b>[petsc-users] MatCreateFFT in Fortran</b><br></span></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif;color:rgba(0,0,0,1.0)"><b>Date: </b></span><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif">August 18, 2016 at 5:15:16 PM CDT<br></span></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif;color:rgba(0,0,0,1.0)"><b>To: </b></span><span style="font-family:-webkit-system-font,Helvetica Neue,Helvetica,sans-serif"><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br></span></div><div><div class="h5"><br><div><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_<wbr>CHARACTER,ierr)<br><br>    call MatCreateFFT(PETSC_COMM_WORLD,<wbr>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>
</div></div></div></blockquote></div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>