Actual source code: ex4f.F

petsc-3.6.4 2016-04-12
Report Typos and Errors
  1: !
  2: !    Test AO with on IS with 0 entries - Fortran version of ex4.c
  3: !
  4:       program main
  5:       implicit none
  6: #include <petsc/finclude/petscsys.h>
  7: #include <petsc/finclude/petscao.h>

  9:       PetscErrorCode ierr
 10:       AO             ao
 11:       PetscInt       localvert(4),nlocal
 12:       PetscMPIInt    rank
 13:       IS             is

 15:       call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
 16:       call MPI_Comm_rank(PETSC_COMM_WORLD,rank,ierr)

 18:       nlocal = 0
 19:       if (rank .eq. 0) then
 20:          nlocal = 4
 21:          localvert(1) = 0
 22:          localvert(2) = 1
 23:          localvert(3) = 2
 24:          localvert(4) = 3
 25:       endif

 27: !     Test AOCreateBasic()
 28:       call AOCreateBasic(PETSC_COMM_WORLD, nlocal, localvert,             &
 29:      &     PETSC_NULL_INTEGER,ao,ierr)
 30:       call AODestroy(ao,ierr)

 32: !     Test AOCreateMemoryScalable()
 33:       call AOCreateMemoryScalable(PETSC_COMM_WORLD, nlocal, localvert,    &
 34:      &     PETSC_NULL_INTEGER,ao,ierr)
 35:       call AODestroy(ao,ierr)

 37:       call AOCreate(PETSC_COMM_WORLD,ao,ierr)
 38:       call ISCreateStride(PETSC_COMM_WORLD,1,0,1,is,ierr)
 39:       call AOSetIS(ao,is,is,ierr)
 40:       call ISDestroy(is,ierr)
 41:       call AOSetType(ao,AOMEMORYSCALABLE,ierr)
 42:       call AODestroy(ao,ierr)

 44:       call PetscFinalize(ierr)
 45:       end