From 6c21c60c70bbcb83fe1895d8d35c720ea8fbb5e3 Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Sat, 22 Sep 2018 12:53:02 +0200 Subject: [PATCH] corrected nbasis --- src/FEM_mech.f90 | 1 + 1 file changed, 1 insertion(+) diff --git a/src/FEM_mech.f90 b/src/FEM_mech.f90 index 3dd6cfb1e..6192a1e49 100644 --- a/src/FEM_mech.f90 +++ b/src/FEM_mech.f90 @@ -140,6 +140,7 @@ subroutine FEM_mech_init(fieldBC) ! Polar decomposition failed in run time call PetscFESetQuadrature(mechFE,mechQuad,ierr); CHKERRQ(ierr) call PetscFEGetDimension(mechFE,nBasis,ierr); CHKERRQ(ierr) + nBasis = nBasis/nc call DMGetDS(mech_mesh,mechDS,ierr); CHKERRQ(ierr) call PetscDSAddDiscretization(mechDS,mechFE,ierr); CHKERRQ(ierr) call PetscDSGetTotalDimension(mechDS,cellDof,ierr); CHKERRQ(ierr)