Actual source code: ex196f90.F
petsc-3.7.2 2016-06-05
1: !
2: !
3: ! This program demonstrates use of MatSeqAIJGetArrayF90()
4: !
5: program main
6: implicit none
7: #include <petsc/finclude/petscsys.h>
8: #include <petsc/finclude/petscvec.h>
9: #include <petsc/finclude/petscmat.h>
10: #include <petsc/finclude/petscviewer.h>
11: #include <petsc/finclude/petscmat.h90>
13: Mat A
14: PetscErrorCode ierr
15: PetscViewer v
16: PetscScalar, pointer :: aa(:)
17: character*(256) f
18: PetscBool flg
20: call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
22: call PetscOptionsGetString(PETSC_NULL_OBJECT,PETSC_NULL_CHARACTER, &
23: & '-f',f,flg,ierr)
24: call PetscViewerBinaryOpen(PETSC_COMM_WORLD,f,FILE_MODE_READ,v, &
25: & ierr)
27: call MatCreate(PETSC_COMM_WORLD,A,ierr)
28: call MatSetType(A, MATSEQAIJ,ierr)
29: call MatLoad(A,v,ierr)
31: call MatView(A,PETSC_VIEWER_STDOUT_WORLD,ierr)
33: call MatSeqAIJGetArrayF90(A,aa,ierr)
34: print*,aa(3)
36: call MatDestroy(A,ierr)
37: call PetscViewerDestroy(v,ierr)
39: call PetscFinalize(ierr)
40: end