Actual source code: ex196f90.F

petsc-3.7.3 2016-07-24
Report Typos and Errors
  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