From e7e959af471ef75ee89d4f2dc409f37c1d80b282 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Wed, 19 Sep 2018 20:06:35 +0200 Subject: [PATCH] using explicit interface from PETScDT --- src/FEM_mech.f90 | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/FEM_mech.f90 b/src/FEM_mech.f90 index 09696b9ab..d77874515 100644 --- a/src/FEM_mech.f90 +++ b/src/FEM_mech.f90 @@ -14,6 +14,7 @@ module FEM_mech use PETScsnes use PETScDM use PETScDMplex + use PETScDT use prec, only: & pInt, & pReal @@ -100,8 +101,8 @@ subroutine FEM_mech_init(fieldBC) IS, pointer :: pBcComps(:), pBcPoints(:) PetscSection :: section PetscInt :: field, faceSet, topologDim, nNodalPoints - PetscReal, pointer :: qPointsP(:), qWeightsP(:), & - nodalPointsP(:), nodalWeightsP(:) + PetscReal, dimension(:) , pointer :: qPointsP, qWeightsP, & + nodalPointsP, nodalWeightsP PetscReal, allocatable, target :: nodalPoints(:), nodalWeights(:) PetscScalar, pointer :: px_scal(:) PetscScalar, allocatable, target :: x_scal(:) @@ -111,7 +112,8 @@ subroutine FEM_mech_init(fieldBC) PetscInt :: cellStart, cellEnd, cell, basis character(len=7) :: prefix = 'mechFE_' PetscErrorCode :: ierr - + PetscReal, allocatable, target, dimension(:) :: qWeights + write(6,'(/,a)') ' <<<+- FEM_mech init -+>>>' write(6,'(a15,a)') ' Current time: ',IO_timeStamp() #include "compilation_info.f90" @@ -123,8 +125,6 @@ subroutine FEM_mech_init(fieldBC) !-------------------------------------------------------------------------------------------------- ! Setup FEM mech discretization - allocate(qPoints(dimPlex*FEM_Zoo_nQuadrature(dimPlex,integrationOrder))) - allocate(qWeights(FEM_Zoo_nQuadrature(dimPlex,integrationOrder))) qPoints = FEM_Zoo_QuadraturePoints(dimPlex,integrationOrder)%p qWeights = FEM_Zoo_QuadratureWeights(dimPlex,integrationOrder)%p nQuadrature = FEM_Zoo_nQuadrature(dimPlex,integrationOrder)