36 |
|
|
37 |
#ifdef ALLOW_PETSC |
#ifdef ALLOW_PETSC |
38 |
#include "finclude/petsc.h" |
#include "finclude/petsc.h" |
39 |
#include "finclude/petscvec.h" |
! UNCOMMENT IF V3.0 |
40 |
#include "finclude/petscmat.h" |
!#include "finclude/petscvec.h" |
41 |
#include "finclude/petscksp.h" |
!#include "finclude/petscmat.h" |
42 |
#include "finclude/petscpc.h" |
!#include "finclude/petscksp.h" |
43 |
|
!#include "finclude/petscpc.h" |
44 |
#endif |
#endif |
45 |
C === Global variables === |
C === Global variables === |
46 |
|
|
192 |
call VecAssemblyBegin(solution, ierr) |
call VecAssemblyBegin(solution, ierr) |
193 |
call VecAssemblyEnd(solution, ierr) |
call VecAssemblyEnd(solution, ierr) |
194 |
|
|
195 |
|
! IF USING v3.0 THEN |
196 |
call MatCreateMPIAIJ (PETSC_COMM_WORLD, |
! call MatCreateMPIAIJ (PETSC_COMM_WORLD, |
197 |
|
call MatCreateAIJ (PETSC_COMM_WORLD, |
198 |
& local_dofs, local_dofs, |
& local_dofs, local_dofs, |
199 |
& global_dofs, global_dofs, |
& global_dofs, global_dofs, |
200 |
& 18, PETSC_NULL_INTEGER, |
& 18, PETSC_NULL_INTEGER, |