Actual source code: ex4f.F
petsc-3.6.4 2016-04-12
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