From 82cdc551a54979b49e5df231dc6b4d6d68190e69 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 18 Sep 2018 05:17:48 +0200 Subject: [PATCH] MPI Communicator only needed for PETSc 3.10. thx to Jaeyong for figuring this out --- src/FEM_mech.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/FEM_mech.f90 b/src/FEM_mech.f90 index e10071eb8..0fdaf6283 100644 --- a/src/FEM_mech.f90 +++ b/src/FEM_mech.f90 @@ -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)