Skip to content
Snippets Groups Projects
Commit 9535048e authored by Dave Whipp's avatar Dave Whipp
Browse files

Modified use of weightel array to not always be deallocated

parent 30c0e1f6
No related branches found
No related tags found
No related merge requests found
......@@ -215,28 +215,10 @@ if (ndof==3) then
call mpi_allreduce (is_plastic,osolve%is_plastic,osolve%nleaves,mpi_logical,mpi_lor,mpi_comm_world,ierr)
call mpi_allreduce (vbounded,vbounded2,osolve%nleaves,mpi_logical,mpi_lor,mpi_comm_world,ierr)
if (iproc.eq.0) write (*,'(a,i10)') shift//'viscosity capped in', count(vbounded2)
!call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
!print *,'ndof: ',ndof
!print *,'weightel_glob min: ',minval(weightel_glob)
!print *,'weightel_glob max: ',maxval(weightel_glob)
!print *,'weightel min: ',minval(weightel)
!print *,'weightel max: ',maxval(weightel)
call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
end if
!weightel_glob=0.d0
!print *,'1. ndof: ',ndof
!print *,'1. weightel_glob min: ',minval(weightel_glob)
!print *,'1. weightel_glob max: ',maxval(weightel_glob)
!print *,'1. weightel min: ',minval(weightel)
!print *,'1. weightel max: ',maxval(weightel)
call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
!print *,'2. ndof: ',ndof
!print *,'2. weightel_glob min: ',minval(weightel_glob)
!print *,'2. weightel_glob max: ',maxval(weightel_glob)
!print *,'2. weightel min: ',minval(weightel)
!print *,'2. weightel max: ',maxval(weightel)
!call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
call heap (threadinfo,'is_plastic','build_system',size(is_plastic),'bool',-1)
deallocate (is_plastic)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment