not used at the moment
This commit is contained in:
Martin Diehl
parent
cbecb71000
commit
5f97ae9080
|
|
@ -27,11 +27,7 @@ use PETScis
|
||||||
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
!--------------------------------------------------------------------------------------------------
|
||||||
! output data
|
! output data
|
||||||
PetscViewer, public :: resUnit
|
|
||||||
Vec, public :: coordinatesVec
|
Vec, public :: coordinatesVec
|
||||||
Vec, allocatable, public :: homogenizationResultsVec(:), &
|
|
||||||
crystalliteResultsVec(:,:), &
|
|
||||||
phaseResultsVec(:,:)
|
|
||||||
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
!--------------------------------------------------------------------------------------------------
|
||||||
! field labels information
|
! field labels information
|
||||||
|
|
@ -217,9 +213,6 @@ subroutine utilities_init()
|
||||||
|
|
||||||
wgt = 1.0/real(mesh_maxNips*mesh_NcpElemsGlobal,pReal)
|
wgt = 1.0/real(mesh_maxNips*mesh_NcpElemsGlobal,pReal)
|
||||||
|
|
||||||
call PetscViewerHDF5Open(PETSC_COMM_WORLD, trim(getSolverJobName())//'.h5', &
|
|
||||||
FILE_MODE_WRITE, resUnit, ierr); CHKERRQ(ierr)
|
|
||||||
call PetscViewerHDF5PushGroup(resUnit, '/', ierr); CHKERRQ(ierr)
|
|
||||||
call DMGetDimension(geomMesh,dimPlex,ierr); CHKERRQ(ierr)
|
call DMGetDimension(geomMesh,dimPlex,ierr); CHKERRQ(ierr)
|
||||||
allocate(nEntities(dimPlex+1), source=0)
|
allocate(nEntities(dimPlex+1), source=0)
|
||||||
allocate(nOutputNodes(worldsize), source = 0)
|
allocate(nOutputNodes(worldsize), source = 0)
|
||||||
|
|
@ -247,178 +240,6 @@ subroutine utilities_init()
|
||||||
write(headerID, '(a,i0)') 'number of cells : ', sum(nOutputCells)
|
write(headerID, '(a,i0)') 'number of cells : ', sum(nOutputCells)
|
||||||
endif
|
endif
|
||||||
|
|
||||||
allocate(connectivity(2**dimPlex,nOutputCells(worldrank+1)))
|
|
||||||
call DMPlexGetHeightStratum(geomMesh,0,cellStart,cellEnd,ierr)
|
|
||||||
CHKERRQ(ierr)
|
|
||||||
ctr = 0
|
|
||||||
select case (integrationOrder)
|
|
||||||
case(1_pInt)
|
|
||||||
do cell = cellStart, cellEnd-1 !< loop over all elements
|
|
||||||
call DMPlexGetTransitiveClosure(geomMesh,cell,PETSC_TRUE,points,ierr)
|
|
||||||
CHKERRQ(ierr)
|
|
||||||
if (dimPlex == 2) then
|
|
||||||
connectivity(:,ctr+1) = [points( 9), points(11), points(13), points(13)] - nEntities(dimPlex+1)
|
|
||||||
ctr = ctr + 1
|
|
||||||
else
|
|
||||||
connectivity(:,ctr+1) = [points(23), points(25), points(27), points(27), &
|
|
||||||
points(29), points(29), points(29), points(29)] - nEntities(dimPlex+1)
|
|
||||||
ctr = ctr + 1
|
|
||||||
endif
|
|
||||||
enddo
|
|
||||||
|
|
||||||
case(2_pInt)
|
|
||||||
do cell = cellStart, cellEnd-1 !< loop over all elements
|
|
||||||
call DMPlexGetTransitiveClosure(geomMesh,cell,PETSC_TRUE,points,ierr)
|
|
||||||
CHKERRQ(ierr)
|
|
||||||
if (dimPlex == 2) then
|
|
||||||
connectivity(:,ctr+1) = [points(9 ), points(3), points(1), points(7)]
|
|
||||||
connectivity(:,ctr+2) = [points(11), points(5), points(1), points(3)]
|
|
||||||
connectivity(:,ctr+3) = [points(13), points(7), points(1), points(5)]
|
|
||||||
ctr = ctr + 3
|
|
||||||
else
|
|
||||||
connectivity(:,ctr+1) = [points(23), points(11), points(3), points(15), points(17), points(5), points(1), points(7)]
|
|
||||||
connectivity(:,ctr+2) = [points(25), points(13), points(3), points(11), points(19), points(9), points(1), points(5)]
|
|
||||||
connectivity(:,ctr+3) = [points(27), points(15), points(3), points(13), points(21), points(7), points(1), points(9)]
|
|
||||||
connectivity(:,ctr+4) = [points(29), points(17), points(7), points(21), points(19), points(5), points(1), points(9)]
|
|
||||||
ctr = ctr + 4_pInt
|
|
||||||
endif
|
|
||||||
enddo
|
|
||||||
|
|
||||||
case default
|
|
||||||
do cell = cellStart, cellEnd-1; do ip = 0, mesh_maxNips-1
|
|
||||||
connectivity(:,ctr+1) = cell*mesh_maxNips + ip
|
|
||||||
ctr = ctr + 1
|
|
||||||
enddo; enddo
|
|
||||||
|
|
||||||
end select
|
|
||||||
connectivity = connectivity + sum(nOutputNodes(1:worldrank))
|
|
||||||
|
|
||||||
call VecCreateMPI(PETSC_COMM_WORLD,dimPlex*nOutputNodes(worldrank+1),dimPlex*sum(nOutputNodes), &
|
|
||||||
coordinatesVec,ierr);CHKERRQ(ierr)
|
|
||||||
call PetscObjectSetName(coordinatesVec, 'NodalCoordinates',ierr)
|
|
||||||
call VecSetFromOptions(coordinatesVec, ierr); CHKERRQ(ierr)
|
|
||||||
|
|
||||||
!allocate(mappingCells(worldsize), source = 0)
|
|
||||||
!do homog = 1, material_Nhomogenization
|
|
||||||
! mappingCells = 0_pInt; mappingCells(worldrank+1) = homogOutput(homog)%sizeIpCells
|
|
||||||
! call MPI_Allreduce(MPI_IN_PLACE,mappingCells,worldsize,MPI_INT,MPI_SUM,PETSC_COMM_WORLD,ierr)
|
|
||||||
! call VecCreateMPI(PETSC_COMM_WORLD,mappingCells(worldrank+1),sum(mappingCells), &
|
|
||||||
! homogenizationResultsVec(homog),ierr);CHKERRQ(ierr)
|
|
||||||
! if (sum(mappingCells) > 0) then
|
|
||||||
! call VecCreateMPI(PETSC_COMM_WORLD,mappingCells(worldrank+1)*2**dimPlex,sum(mappingCells)*2**dimPlex, &
|
|
||||||
! connectivityVec,ierr);CHKERRQ(ierr)
|
|
||||||
! call PetscObjectSetName(connectivityVec,'mapping_'//trim(homogenization_name(homog)),ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! call VecGetArrayF90(connectivityVec,results,ierr); CHKERRQ(ierr)
|
|
||||||
! results = 0.0_pReal; ctr = 1_pInt
|
|
||||||
! do cell = cellStart, cellEnd-1; do qPt = 1, mesh_maxNips
|
|
||||||
! if (material_homog(qPt,cell+1) == homog) then
|
|
||||||
! results(ctr:ctr+2**dimPlex-1) = real(reshape(connectivity(1:2**dimPlex,mesh_maxNips*cell+qPt), &
|
|
||||||
! shape=[2**dimPlex]))
|
|
||||||
! ctr = ctr + 2**dimPlex
|
|
||||||
! endif
|
|
||||||
! enddo; enddo
|
|
||||||
! call VecRestoreArrayF90(connectivityVec, results, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecAssemblyBegin(connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecAssemblyEnd (connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecView(connectivityVec, resUnit, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecDestroy(connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! endif
|
|
||||||
!enddo
|
|
||||||
!do cryst = 1, material_Ncrystallite; do grain = 1, homogenization_maxNgrains
|
|
||||||
! mappingCells = 0_pInt
|
|
||||||
! mappingCells(worldrank+1) = crystalliteOutput(cryst,grain)%sizeIpCells
|
|
||||||
! call MPI_Allreduce(MPI_IN_PLACE,mappingCells,worldsize,MPI_INT,MPI_SUM,PETSC_COMM_WORLD,ierr)
|
|
||||||
! call VecCreateMPI(PETSC_COMM_WORLD,mappingCells(worldrank+1),sum(mappingCells), &
|
|
||||||
! crystalliteResultsVec(cryst,grain),ierr);CHKERRQ(ierr)
|
|
||||||
! if (sum(mappingCells) > 0) then
|
|
||||||
! call VecCreateMPI(PETSC_COMM_WORLD,mappingCells(worldrank+1)*2**dimPlex,sum(mappingCells)*2**dimPlex, &
|
|
||||||
! connectivityVec,ierr);CHKERRQ(ierr)
|
|
||||||
! write(grainStr,'(a,i0)') 'Grain',grain
|
|
||||||
! call PetscObjectSetName(connectivityVec,'mapping_'// &
|
|
||||||
! trim(crystallite_name(cryst))//'_'// &
|
|
||||||
! trim(grainStr),ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! call VecGetArrayF90(connectivityVec, results, ierr); CHKERRQ(ierr)
|
|
||||||
! results = 0.0_pReal; ctr = 1_pInt
|
|
||||||
! do cell = cellStart, cellEnd-1; do qPt = 1, mesh_maxNips
|
|
||||||
! if (homogenization_Ngrains (mesh_element(3,cell+1)) >= grain .and. &
|
|
||||||
! microstructure_crystallite(mesh_element(4,cell+1)) == cryst) then
|
|
||||||
! results(ctr:ctr+2**dimPlex-1) = real(reshape(connectivity(1:2**dimPlex,mesh_maxNips*cell+qPt), &
|
|
||||||
! shape=[2**dimPlex]))
|
|
||||||
! ctr = ctr + 2**dimPlex
|
|
||||||
! endif
|
|
||||||
! enddo; enddo
|
|
||||||
! call VecRestoreArrayF90(connectivityVec, results, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecAssemblyBegin(connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecAssemblyEnd (connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecView(connectivityVec, resUnit, ierr); CHKERRQ(ierr)
|
|
||||||
! call VecDestroy(connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! endif
|
|
||||||
!enddo; enddo
|
|
||||||
!do phase = 1, material_Nphase; do grain = 1, homogenization_maxNgrains
|
|
||||||
! mappingCells = 0_pInt
|
|
||||||
! mappingCells(worldrank+1) = phaseOutput(phase,grain)%sizeIpCells
|
|
||||||
! call MPI_Allreduce(MPI_IN_PLACE,mappingCells,worldsize,MPI_INT,MPI_SUM,PETSC_COMM_WORLD,ierr)
|
|
||||||
! call VecCreateMPI(PETSC_COMM_WORLD,mappingCells(worldrank+1),sum(mappingCells), &
|
|
||||||
! phaseResultsVec(phase,grain),ierr);CHKERRQ(ierr)
|
|
||||||
! if (sum(mappingCells) > 0) then
|
|
||||||
! call VecCreateMPI(PETSC_COMM_WORLD,mappingCells(worldrank+1)*2**dimPlex,sum(mappingCells)*2**dimPlex, &
|
|
||||||
! connectivityVec,ierr);CHKERRQ(ierr)
|
|
||||||
! write(grainStr,'(a,i0)') 'Grain',grain
|
|
||||||
! call PetscObjectSetName(connectivityVec,&
|
|
||||||
! 'mapping_'//trim(phase_name(phase))//'_'// &
|
|
||||||
! trim(grainStr),ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! call VecGetArrayF90(connectivityVec, results, ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! results = 0.0_pReal; ctr = 1_pInt
|
|
||||||
! do cell = cellStart, cellEnd-1; do qPt = 1, mesh_maxNips
|
|
||||||
! if (material_phase(grain,qPt,cell+1) == phase) then
|
|
||||||
! results(ctr:ctr+2**dimPlex-1) = real(reshape(connectivity(1:2**dimPlex,mesh_maxNips*cell+qPt), &
|
|
||||||
! shape=[2**dimPlex]))
|
|
||||||
! ctr = ctr + 2**dimPlex
|
|
||||||
! endif
|
|
||||||
! enddo; enddo
|
|
||||||
! call VecRestoreArrayF90(connectivityVec, results, ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! call VecAssemblyBegin(connectivityVec, ierr);CHKERRQ(ierr)
|
|
||||||
! call VecAssemblyEnd (connectivityVec, ierr);CHKERRQ(ierr)
|
|
||||||
! call VecView(connectivityVec, resUnit, ierr);CHKERRQ(ierr)
|
|
||||||
! call VecDestroy(connectivityVec, ierr); CHKERRQ(ierr)
|
|
||||||
! endif
|
|
||||||
!enddo; enddo
|
|
||||||
!if (worldrank == 0_pInt) then
|
|
||||||
! do homog = 1, material_Nhomogenization
|
|
||||||
! call VecGetSize(homogenizationResultsVec(homog),mappingCells(1),ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! if (mappingCells(1) > 0) &
|
|
||||||
! write(headerID, '(a,i0)') 'number of homog_'// &
|
|
||||||
! trim(homogenization_name(homog))//'_'// &
|
|
||||||
! 'cells : ', mappingCells(1)
|
|
||||||
! enddo
|
|
||||||
! do cryst = 1, material_Ncrystallite; do grain = 1, homogenization_maxNgrains
|
|
||||||
! call VecGetSize(crystalliteResultsVec(cryst,grain),mappingCells(1),ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! write(grainStr,'(a,i0)') 'Grain',grain
|
|
||||||
! if (mappingCells(1) > 0) &
|
|
||||||
! write(headerID, '(a,i0)') 'number of cryst_'// &
|
|
||||||
! trim(crystallite_name(cryst))//'_'// &
|
|
||||||
! trim(grainStr)//'_'// &
|
|
||||||
! 'cells : ', mappingCells(1)
|
|
||||||
! enddo; enddo
|
|
||||||
! do phase = 1, material_Nphase; do grain = 1, homogenization_maxNgrains
|
|
||||||
! call VecGetSize(phaseResultsVec(phase,grain),mappingCells(1),ierr)
|
|
||||||
! CHKERRQ(ierr)
|
|
||||||
! write(grainStr,'(a,i0)') 'Grain',grain
|
|
||||||
! if (mappingCells(1) > 0) &
|
|
||||||
! write(headerID, '(a,i0)') 'number of phase_'// &
|
|
||||||
! trim(phase_name(phase))//'_'//trim(grainStr)//'_'// &
|
|
||||||
! 'cells : ', mappingCells(1)
|
|
||||||
! enddo; enddo
|
|
||||||
! close(headerID)
|
|
||||||
!endif
|
|
||||||
|
|
||||||
end subroutine utilities_init
|
end subroutine utilities_init
|
||||||
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue