corrected nbasis

This commit is contained in:
Pratheek Shanthraj 2018-09-22 12:53:02 +02:00
parent becd959913
commit 6c21c60c70
1 changed files with 1 additions and 0 deletions

View File

@ -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)