223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

WIP: separate reading in and parsing

This commit is contained in:
Martin Diehl 2018-05-31 18:31:44 +02:00
parent 72b69959de
commit 804febe7f9
2 changed files with 71 additions and 45 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/math.f90
View File

@ -160,8 +160,7 @@ module math
math_rotate_forward33, & math_rotate_forward33, &
math_rotate_backward33, & math_rotate_backward33, &
math_rotate_forward3333, & math_rotate_forward3333, &
math_limit, & math_limit
math_expand
private :: & private :: &
math_check, & math_check, &
halton halton

113
src/plastic_phenopowerlaw.f90
View File

@ -26,14 +26,11 @@ module plastic_phenopowerlaw
totalNslip, & !< no. of slip system used in simulation totalNslip, & !< no. of slip system used in simulation
totalNtwin !< no. of twin system used in simulation totalNtwin !< no. of twin system used in simulation
real(pReal), dimension(:,:,:), allocatable, private :: & real(pReal), dimension(:,:,:), allocatable, private :: &
interaction_SlipSlip, & !< interaction factors slip - slip (input parameter)
interaction_SlipSlip, & !< interaction factors slip - slip (input parameter)
interaction_SlipTwin, & !< interaction factors slip - twin (input parameter) interaction_SlipTwin, & !< interaction factors slip - twin (input parameter)
interaction_TwinSlip, & !< interaction factors twin - slip (input parameter) interaction_TwinSlip, & !< interaction factors twin - slip (input parameter)
interaction_TwinTwin !< interaction factors twin - twin (input parameter) interaction_TwinTwin !< interaction factors twin - twin (input parameter)
enum, bind(c) enum, bind(c)
@ -52,7 +49,16 @@ module plastic_phenopowerlaw
integer(kind(undefined_ID)), dimension(:,:), allocatable, private :: & integer(kind(undefined_ID)), dimension(:,:), allocatable, private :: &
plastic_phenopowerlaw_outputID !< ID of each post result output plastic_phenopowerlaw_outputID !< ID of each post result output
type :: tKeyValues
character(len=64) :: &
key = ''
character(len=65536) :: &
rawValues = ''
end type
type, private :: tParameters !< container type for internal constitutive parameters type, private :: tParameters !< container type for internal constitutive parameters
type(tKeyValues) :: &
keyValues
real(pReal) :: & real(pReal) :: &
gdot0_slip, & !< reference shear strain rate for slip gdot0_slip, & !< reference shear strain rate for slip
gdot0_twin, & !< reference shear strain rate for twin gdot0_twin, & !< reference shear strain rate for twin
@ -172,6 +178,9 @@ subroutine plastic_phenopowerlaw_init(fileUnit)
offset_slip, index_myFamily, index_otherFamily, & offset_slip, index_myFamily, index_otherFamily, &
mySize=0_pInt,sizeState,sizeDotState, sizeDeltaState, & mySize=0_pInt,sizeState,sizeDotState, sizeDeltaState, &
startIndex, endIndex startIndex, endIndex
type(tKeyValues) :: keyValuesTemp
character(len=65536) :: & character(len=65536) :: &
tag = '', & tag = '', &
line = '', & line = '', &
@ -184,67 +193,85 @@ subroutine plastic_phenopowerlaw_init(fileUnit)
write(6,'(a15,a)') ' Current time: ',IO_timeStamp() write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
#include "compilation_info.f90" #include "compilation_info.f90"
maxNinstance = int(count(phase_plasticity == PLASTICITY_PHENOPOWERLAW_ID),pInt) maxNinstance = int(count(phase_plasticity == PLASTICITY_PHENOPOWERLAW_ID),pInt) ! ToDo: this does not happen
if (maxNinstance == 0_pInt) return if (maxNinstance == 0_pInt) return
if (iand(debug_level(debug_constitutive),debug_levelBasic) /= 0_pInt) & if (iand(debug_level(debug_constitutive),debug_levelBasic) /= 0_pInt) &
write(6,'(a16,1x,i5,/)') '# instances:',maxNinstance write(6,'(a16,1x,i5,/)') '# instances:',maxNinstance
allocate(plastic_phenopowerlaw_sizePostResults(maxNinstance), source=0_pInt) allocate(plastic_phenopowerlaw_sizePostResults(maxNinstance), source=0_pInt)
allocate(plastic_phenopowerlaw_sizePostResult(maxval(phase_Noutput),maxNinstance), & allocate(plastic_phenopowerlaw_sizePostResult(maxval(phase_Noutput),maxNinstance),source=0_pInt)
source=0_pInt) allocate(plastic_phenopowerlaw_Noutput(maxNinstance), source=0_pInt)
allocate(plastic_phenopowerlaw_output(maxval(phase_Noutput),maxNinstance)) allocate(plastic_phenopowerlaw_output(maxval(phase_Noutput),maxNinstance))
plastic_phenopowerlaw_output = '' plastic_phenopowerlaw_output = ''
allocate(plastic_phenopowerlaw_outputID(maxval(phase_Noutput),maxNinstance),source=undefined_ID)
allocate(plastic_phenopowerlaw_Noutput(maxNinstance), source=0_pInt)
allocate(plastic_phenopowerlaw_outputID(maxval(phase_Noutput),maxNinstance),source=undefined_ID)
allocate(totalNslip(maxNinstance), source=0_pInt) allocate(totalNslip(maxNinstance), source=0_pInt)
allocate(totalNtwin(maxNinstance), source=0_pInt) allocate(totalNtwin(maxNinstance), source=0_pInt)
allocate(param(maxNinstance)) ! one container of parameters per instance allocate(param(maxNinstance)) ! one container of parameters per instance
rewind(fileUnit) rewind(fileUnit)
phase = 0_pInt phase = 0_pInt
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= material_partPhase) ! wind forward to <phase> windForward: do while (IO_lc(IO_getTag(line,'<','>')) /= material_partPhase)
line = IO_read(fileUnit) line = IO_read(fileUnit)
enddo enddo windForward
getKeys: do while (trim(line) /= IO_EOF) ! read through sections of phase part
parsingFile: do while (trim(line) /= IO_EOF) ! read through sections of phase part
line = IO_read(fileUnit) line = IO_read(fileUnit)
if (IO_isBlank(line)) cycle ! skip empty lines if (IO_isBlank(line) .or. phase == 0_pInt) cycle ! skip empty lines
if (IO_getTag(line,'[',']') /= '') phase = phase + 1_pInt ! next phase
phase = phase + 1_pInt ! advance phase section counter
instance = phase_plasticityInstance(phase) ! instance of present phase
cycle
endif
if (phase_plasticity(phase) /= PLASTICITY_PHENOPOWERLAW_ID) cycle
if (IO_getTag(line,'<','>') /= '') then ! stop at next part if (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read line = IO_read(fileUnit, .true.) ! reset IO_read
exit exit
endif endif
if (IO_getTag(line,'[',']') /= '') then ! next phase chunkPos = IO_stringPos(line)
phase = phase + 1_pInt ! advance phase section counter keyValuesTemp%key = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key
if (phase_plasticity(phase) == PLASTICITY_PHENOPOWERLAW_ID) then if(chunkPos(1) > 1) keyValuesTemp%rawValues = IO_lc(line(chunkPos(4),:))
instance = phase_plasticityInstance(phase) ! which instance of my plasticity is present phase param(instance)%keyValues = [(instance)%keyValues,keyValuesTemp]
Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt) ! maximum number of slip families according to lattice type of current phase enddo getKeys
Nchunks_TwinFamilies = count(lattice_NtwinSystem(:,phase) > 0_pInt) ! maximum number of twin families according to lattice type of current phase
Nchunks_SlipSlip = maxval(lattice_interactionSlipSlip(:,:,phase))
Nchunks_SlipTwin = maxval(lattice_interactionSlipTwin(:,:,phase))
Nchunks_TwinSlip = maxval(lattice_interactionTwinSlip(:,:,phase))
Nchunks_TwinTwin = maxval(lattice_interactionTwinTwin(:,:,phase))
Nchunks_nonSchmid = lattice_NnonSchmid(phase)
if(allocated(tempPerSlip)) deallocate(tempPerSlip)
!allocate(param(instance)%H_int,source=tempPerSlip) gfortran 5 does not support this
allocate(param(instance)%H_int(Nchunks_SlipFamilies),source=0.0_pReal)
allocate(param(instance)%interaction_SlipSlip(Nchunks_SlipSlip),source=0.0_pReal)
allocate(param(instance)%interaction_SlipTwin(Nchunks_SlipTwin),source=0.0_pReal)
allocate(param(instance)%interaction_TwinSlip(Nchunks_TwinSlip),source=0.0_pReal)
allocate(param(instance)%interaction_TwinTwin(Nchunks_TwinTwin),source=0.0_pReal)
allocate(param(instance)%nonSchmidCoeff(Nchunks_nonSchmid),source=0.0_pReal)
allocate(tempPerSlip(Nchunks_SlipFamilies)) parseString: do instance = 1_pInt, maxNinstance
endif do i = 1_pInt, size(param(instance)%keyValues); key = param(instance)%keyValues(i)
cycle ! skip to next line enddo
enddo parseStrings
myPhase: if (phase_plasticity(phase) == PLASTICITY_phenopowerlaw_ID) then
instance = phase_plasticityInstance(phase)
! if (phase_plasticity(phase) == PLASTICITY_PHENOPOWERLAW_ID) then
! instance = phase_plasticityInstance(phase) ! which instance of my plasticity is present phase
! chunkPos = IO_stringPos(line)
! configTemp%key = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key
! if(chunkPos(1) > 1) configTemp%rawValues = IO_lc(line(chunkPos(4),:))
! config = [config,configTemp]
! Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt) ! maximum number of slip families according to lattice type of current phase
! Nchunks_TwinFamilies = count(lattice_NtwinSystem(:,phase) > 0_pInt) ! maximum number of twin families according to lattice type of current phase
! Nchunks_SlipSlip = maxval(lattice_interactionSlipSlip(:,:,phase))
! Nchunks_SlipTwin = maxval(lattice_interactionSlipTwin(:,:,phase))
! Nchunks_TwinSlip = maxval(lattice_interactionTwinSlip(:,:,phase))
! Nchunks_TwinTwin = maxval(lattice_interactionTwinTwin(:,:,phase))
! Nchunks_nonSchmid = lattice_NnonSchmid(phase)
! if(allocated(tempPerSlip)) deallocate(tempPerSlip)
! !allocate(param(instance)%H_int,source=tempPerSlip) gfortran 5 does not support this
! allocate(param(instance)%H_int(Nchunks_SlipFamilies),source=0.0_pReal)
! allocate(param(instance)%interaction_SlipSlip(Nchunks_SlipSlip),source=0.0_pReal)
! allocate(param(instance)%interaction_SlipTwin(Nchunks_SlipTwin),source=0.0_pReal)
! allocate(param(instance)%interaction_TwinSlip(Nchunks_TwinSlip),source=0.0_pReal)
! allocate(param(instance)%interaction_TwinTwin(Nchunks_TwinTwin),source=0.0_pReal)
! allocate(param(instance)%nonSchmidCoeff(Nchunks_nonSchmid),source=0.0_pReal)
! allocate(tempPerSlip(Nchunks_SlipFamilies))
! endif
! cycle ! skip to next line
endif endif
if (phase > 0_pInt ) then; if (phase_plasticity(phase) == PLASTICITY_PHENOPOWERLAW_ID) then ! one of my phases. Do not short-circuit here (.and. between if-statements), it's not safe in Fortran if (phase > 0_pInt ) then; if (phase_plasticity(phase) == PLASTICITY_PHENOPOWERLAW_ID) then ! one of my phases. Do not short-circuit here (.and. between if-statements), it's not safe in Fortran
tag = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key
chunkPos = IO_stringPos(line)
tag = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key
select case(tag) select case(tag)
case ('(output)') case ('(output)')