[petsc-dev] TS ex11 crashing in parallel

Matthew Knepley knepley at gmail.com
Wed May 28 22:57:56 CDT 2014


On Wed, May 28, 2014 at 10:26 PM, Adrian Croucher <a.croucher at auckland.ac.nz
> wrote:

>
> On 29/05/14 14:48, Matthew Knepley wrote:
>
>
>  OK. There are a few other things missing or broken in the Fortran
>> interfaces that I will need at some point- e.g.
>> DMPlexConstructGhostCells(), PetscFVCreate(). At present I'm writing my
>> test code in C to avoid this issue, but will need to get it going in
>> fortran eventually.
>>
>
>  Okay,  I pushed those fixes to next.
>
>
> Oh great, thanks very much :-)
>
> Is there a trick to calling PetscFVCreate() from fortran though? The test
> code snippet below errors out when I run it:
>

You might have an old version of 'bfort' which does not know PetscFV is a
type and thus creates
incorrect bindings. I am sending you src/dm/dt/interface/ftn-auto/dtfvf.c
so you can test, but
you should

  rm -rf $PETSC_DIR/$PETSC_ARCH/externalpackages/sowing*
  rm -rf $PETSC_DIR/$PETSC_ARCH/conf/sowing

and reconfigure using

  ./$PETSC_ARCH/conf/reconfigure-$PETSC_ARCH.py

  Thanks,

     Matt


> program main
>
> ! testing PetscFVM from fortran
>
>   implicit none
> #include <finclude/petsc.h90>
>   PetscFV :: fvm
>   PetscErrorCode :: ierr
>
>   call PetscInitialize(PETSC_NULL_CHARACTER,ierr); CHKERRQ(ierr)
>   call PetscFVCreate(PETSC_COMM_WORLD, fvm, ierr); CHKERRQ(ierr)
>   call PetscFVDestroy(fvm, ierr); CHKERRQ(ierr)
>   call PetscFinalize(ierr); CHKERRQ(ierr)
>
> end program main
>
> The error is:
>
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [0]PETSC ERROR: Null argument, when expecting valid pointer
> [0]PETSC ERROR: Null Pointer: Parameter # 2
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.htmlfor trouble shooting.
> [0]PETSC ERROR: Petsc Development GIT revision: v3.4.4-5833-g396e404  GIT
> Date: 2014-05-28 21:50:40 -0500
> [0]PETSC ERROR: ./testfvm on a linux-gnu-c-opt named des108 by acro018 Thu
> May 29 15:01:42 2014
> [0]PETSC ERROR: Configure options --download-netcdf --download-exodusii
> --with-hdf5-dir=/usr/lib --with-numpy-dir=/usr/lib/pymodules/p
> ython2.7/numpy
> --with-scientificpython-dir=/usr/lib/python2.7/dist-packages/scipy
> --with-petsc4py-dir=/usr/local/lib/python2.7/dist-pa
> ckages/petsc4py --download-triangle --with-fortran-interfaces=1
> --download-ptscotch --download-chaco --download-ctetgen
> [0]PETSC ERROR: #1 PetscFVCreate() line 1356 in
> /home/acro018/software/PETSc/code/src/dm/dt/interface/dtfv.c
> --------------------------------------------------------------------------
>
> Cheers, Adrian
>
> --
> Dr Adrian Croucher
> Senior Research Fellow
> Department of Engineering Science
> University of Auckland, New Zealand
> email: a.croucher at auckland.ac.nz
> tel: +64 (0)9 923 4611
>
>


-- 
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/20140528/83745bf8/attachment.html>


More information about the petsc-dev mailing list