MPI Communicator only needed for PETSc 3.10.

thx to Jaeyong for figuring this out
This commit is contained in:
Martin Diehl 2018-09-18 05:17:48 +02:00
parent 245676728d
commit 82cdc551a5
1 changed files with 1 additions and 1 deletions

View File

@ -140,7 +140,7 @@ subroutine FEM_mech_init(fieldBC)
nc = dimPlex
call PetscQuadratureSetData(mechQuad,dimPlex,nc,nQuadrature,qPointsP,qWeightsP,ierr)
CHKERRQ(ierr)
call PetscFECreateDefault(PETSC_COMM_SELF,mech_mesh,dimPlex,nc,PETSC_TRUE,prefix, &
call PetscFECreateDefault(mech_mesh,dimPlex,nc,PETSC_TRUE,prefix, &
integrationOrder,mechFE,ierr); CHKERRQ(ierr)
call PetscFESetQuadrature(mechFE,mechQuad,ierr); CHKERRQ(ierr)
call PetscFEGetDimension(mechFE,nBasis,ierr); CHKERRQ(ierr)