From 42bc1b468a2bba010beb2549567f3aead1b1a4f4 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 22 Sep 2018 11:19:08 +0200 Subject: [PATCH] cleaning+adopting to petsc-3.10 --- src/FEM_mech.f90 | 7 +------ 1 file changed, 1 insertion(+), 6 deletions(-) diff --git a/src/FEM_mech.f90 b/src/FEM_mech.f90 index 630029ac7..620c17a79 100644 --- a/src/FEM_mech.f90 +++ b/src/FEM_mech.f90 @@ -243,17 +243,12 @@ subroutine FEM_mech_init(fieldBC) call PetscFEGetDualSpace(mechFE,mechDualSpace,ierr); CHKERRQ(ierr) call DMPlexGetHeightStratum(mech_mesh,0,cellStart,cellEnd,ierr) CHKERRQ(ierr) - write(6,*) 'cellDof', cellDof;flush(6) - write(6,*) 'cell start and end-1',cellStart,cellEnd-1;flush(6) do cell = cellStart, cellEnd-1 !< loop over all elements - write(6,*) 'cell',cell;flush(6) x_scal = 0.0 call DMPlexComputeCellGeometryAffineFEM(mech_mesh,cell,pV0,pCellJ,pInvcellJ,detJ,ierr) CHKERRQ(ierr) cellJMat = reshape(pCellJ,shape=[dimPlex,dimPlex]) - do basis = 0, nBasis-1 - write(6,*) 'nBasis-1',nBasis-1;flush(6) - write(6,*) 'basis',basis;flush(6) + do basis = 0, nBasis/dimPlex-1 call PetscDualSpaceGetFunctional(mechDualSpace,basis,functional,ierr) CHKERRQ(ierr) call PetscQuadratureGetData(functional,dimPlex,nc,nNodalPoints,nodalPointsP,nodalWeightsP,ierr)