475 lines
22 KiB
Fortran
475 lines
22 KiB
Fortran
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
|
||
|
|
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
|
||
|
|
!> @author Christoph Kords, Max-Planck-Institut für Eisenforschung GmbH
|
||
|
|
!> @author Martin Diehl, Max-Planck-Institut für Eisenforschung GmbH
|
||
|
|
!> @brief Reading in and interpretating the debugging settings for the various modules
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
module debug
|
||
|
|
use prec, only: &
|
||
|
|
pInt, &
|
||
|
|
pReal, &
|
||
|
|
pLongInt
|
||
|
|
|
||
|
|
implicit none
|
||
|
|
private
|
||
|
|
integer(pInt), parameter, public :: &
|
||
|
|
debug_LEVELSELECTIVE = 2_pInt**0_pInt, &
|
||
|
|
debug_LEVELBASIC = 2_pInt**1_pInt, &
|
||
|
|
debug_LEVELEXTENSIVE = 2_pInt**2_pInt
|
||
|
|
integer(pInt), parameter, private :: &
|
||
|
|
debug_MAXGENERAL = debug_LEVELEXTENSIVE ! must be set to the last bitcode used by (potentially) all debug types
|
||
|
|
integer(pInt), parameter, public :: &
|
||
|
|
debug_SPECTRALRESTART = debug_MAXGENERAL*2_pInt**1_pInt, &
|
||
|
|
debug_SPECTRALFFTW = debug_MAXGENERAL*2_pInt**2_pInt, &
|
||
|
|
debug_SPECTRALDIVERGENCE = debug_MAXGENERAL*2_pInt**3_pInt, &
|
||
|
|
debug_SPECTRALROTATION = debug_MAXGENERAL*2_pInt**4_pInt, &
|
||
|
|
debug_SPECTRALPETSC = debug_MAXGENERAL*2_pInt**5_pInt
|
||
|
|
|
||
|
|
integer(pInt), parameter, public :: &
|
||
|
|
debug_DEBUG = 1_pInt, &
|
||
|
|
debug_MATH = 2_pInt, &
|
||
|
|
debug_FESOLVING = 3_pInt, &
|
||
|
|
debug_MESH = 4_pInt, & !< stores debug level for mesh part of DAMASK bitwise coded
|
||
|
|
debug_MATERIAL = 5_pInt, & !< stores debug level for material part of DAMASK bitwise coded
|
||
|
|
debug_LATTICE = 6_pInt, & !< stores debug level for lattice part of DAMASK bitwise coded
|
||
|
|
debug_CONSTITUTIVE = 7_pInt, & !< stores debug level for constitutive part of DAMASK bitwise coded
|
||
|
|
debug_CRYSTALLITE = 8_pInt, &
|
||
|
|
debug_HOMOGENIZATION = 9_pInt, &
|
||
|
|
debug_CPFEM = 10_pInt, &
|
||
|
|
debug_SPECTRAL = 11_pInt, &
|
||
|
|
debug_MARC = 12_pInt, &
|
||
|
|
debug_ABAQUS = 13_pInt
|
||
|
|
integer(pInt), parameter, private :: &
|
||
|
|
debug_MAXNTYPE = debug_ABAQUS !< must be set to the maximum defined debug type
|
||
|
|
|
||
|
|
integer(pInt),protected, dimension(debug_maxNtype+2_pInt), public :: & ! specific ones, and 2 for "all" and "other"
|
||
|
|
debug_level = 0_pInt
|
||
|
|
|
||
|
|
integer(pInt), public :: &
|
||
|
|
debug_cumLpCalls = 0_pInt, & !< total number of calls to LpAndItsTangent
|
||
|
|
debug_cumDeltaStateCalls = 0_pInt, & !< total number of calls to deltaState
|
||
|
|
debug_cumDotStateCalls = 0_pInt !< total number of calls to dotState
|
||
|
|
|
||
|
|
integer(pInt), protected, public :: &
|
||
|
|
debug_e = 1_pInt, &
|
||
|
|
debug_i = 1_pInt, &
|
||
|
|
debug_g = 1_pInt
|
||
|
|
|
||
|
|
integer(pLongInt), public :: &
|
||
|
|
debug_cumLpTicks = 0_pLongInt, & !< total cpu ticks spent in LpAndItsTangent
|
||
|
|
debug_cumDeltaStateTicks = 0_pLongInt, & !< total cpu ticks spent in deltaState
|
||
|
|
debug_cumDotStateTicks = 0_pLongInt !< total cpu ticks spent in dotState
|
||
|
|
|
||
|
|
integer(pInt), dimension(2), public :: &
|
||
|
|
debug_stressMaxLocation = 0_pInt, &
|
||
|
|
debug_stressMinLocation = 0_pInt, &
|
||
|
|
debug_jacobianMaxLocation = 0_pInt, &
|
||
|
|
debug_jacobianMinLocation = 0_pInt
|
||
|
|
|
||
|
|
integer(pInt), dimension(:), allocatable, public :: &
|
||
|
|
debug_CrystalliteLoopDistribution, & !< distribution of crystallite cutbacks
|
||
|
|
debug_MaterialpointStateLoopDistribution, &
|
||
|
|
debug_MaterialpointLoopDistribution
|
||
|
|
|
||
|
|
integer(pInt), dimension(:,:), allocatable, public :: &
|
||
|
|
debug_StressLoopLiDistribution, & !< distribution of stress iterations until convergence
|
||
|
|
debug_StressLoopLpDistribution, & !< distribution of stress iterations until convergence
|
||
|
|
debug_StateLoopDistribution !< distribution of state iterations until convergence
|
||
|
|
|
||
|
|
real(pReal), public :: &
|
||
|
|
debug_stressMax = -huge(1.0_pReal), &
|
||
|
|
debug_stressMin = huge(1.0_pReal), &
|
||
|
|
debug_jacobianMax = -huge(1.0_pReal), &
|
||
|
|
debug_jacobianMin = huge(1.0_pReal)
|
||
|
|
|
||
|
|
character(len=64), parameter, private :: &
|
||
|
|
debug_CONFIGFILE = 'debug.config' !< name of configuration file
|
||
|
|
|
||
|
|
#ifdef PETSc
|
||
|
|
character(len=1024), parameter, public :: &
|
||
|
|
PETSCDEBUG = ' -snes_view -snes_monitor '
|
||
|
|
#endif
|
||
|
|
public :: debug_init, &
|
||
|
|
debug_reset, &
|
||
|
|
debug_info
|
||
|
|
|
||
|
|
contains
|
||
|
|
|
||
|
|
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
!> @brief reads in parameters from debug.config and allocates arrays
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
subroutine debug_init
|
||
|
|
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
|
||
|
|
use numerics, only: &
|
||
|
|
worldrank, &
|
||
|
|
nStress, &
|
||
|
|
nState, &
|
||
|
|
nCryst, &
|
||
|
|
nMPstate, &
|
||
|
|
nHomog
|
||
|
|
use IO, only: &
|
||
|
|
IO_read, &
|
||
|
|
IO_error, &
|
||
|
|
IO_open_file_stat, &
|
||
|
|
IO_isBlank, &
|
||
|
|
IO_stringPos, &
|
||
|
|
IO_stringValue, &
|
||
|
|
IO_lc, &
|
||
|
|
IO_floatValue, &
|
||
|
|
IO_intValue, &
|
||
|
|
IO_timeStamp, &
|
||
|
|
IO_EOF
|
||
|
|
|
||
|
|
implicit none
|
||
|
|
integer(pInt), parameter :: FILEUNIT = 300_pInt
|
||
|
|
|
||
|
|
integer(pInt) :: i, what
|
||
|
|
integer(pInt), allocatable, dimension(:) :: chunkPos
|
||
|
|
character(len=65536) :: tag, line
|
||
|
|
|
||
|
|
mainProcess: if (worldrank == 0) then
|
||
|
|
write(6,'(/,a)') ' <<<+- debug init -+>>>'
|
||
|
|
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
|
||
|
|
#include "compilation_info.f90"
|
||
|
|
endif mainProcess
|
||
|
|
|
||
|
|
if (allocated(debug_StressLoopLpDistribution)) &
|
||
|
|
deallocate(debug_StressLoopLpDistribution)
|
||
|
|
allocate(debug_StressLoopLpDistribution(nStress+1,2))
|
||
|
|
debug_StressLoopLpDistribution = 0_pInt
|
||
|
|
if (allocated(debug_StressLoopLiDistribution)) &
|
||
|
|
deallocate(debug_StressLoopLiDistribution)
|
||
|
|
allocate(debug_StressLoopLiDistribution(nStress+1,2))
|
||
|
|
debug_StressLoopLiDistribution = 0_pInt
|
||
|
|
if (allocated(debug_StateLoopDistribution)) &
|
||
|
|
deallocate(debug_StateLoopDistribution)
|
||
|
|
allocate(debug_StateLoopDistribution(nState+1,2))
|
||
|
|
debug_StateLoopDistribution = 0_pInt
|
||
|
|
if (allocated(debug_CrystalliteLoopDistribution)) &
|
||
|
|
deallocate(debug_CrystalliteLoopDistribution)
|
||
|
|
allocate(debug_CrystalliteLoopDistribution(nCryst+1))
|
||
|
|
debug_CrystalliteLoopDistribution = 0_pInt
|
||
|
|
if (allocated(debug_MaterialpointStateLoopDistribution)) &
|
||
|
|
deallocate(debug_MaterialpointStateLoopDistribution)
|
||
|
|
allocate(debug_MaterialpointStateLoopDistribution(nMPstate))
|
||
|
|
debug_MaterialpointStateLoopDistribution = 0_pInt
|
||
|
|
if (allocated(debug_MaterialpointLoopDistribution)) &
|
||
|
|
deallocate(debug_MaterialpointLoopDistribution)
|
||
|
|
allocate(debug_MaterialpointLoopDistribution(nHomog+1))
|
||
|
|
debug_MaterialpointLoopDistribution = 0_pInt
|
||
|
|
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
! try to open the config file
|
||
|
|
|
||
|
|
line = ''
|
||
|
|
fileExists: if(IO_open_file_stat(FILEUNIT,debug_configFile)) then
|
||
|
|
do while (trim(line) /= IO_EOF) ! read thru sections of phase part
|
||
|
|
line = IO_read(FILEUNIT)
|
||
|
|
if (IO_isBlank(line)) cycle ! skip empty lines
|
||
|
|
chunkPos = IO_stringPos(line)
|
||
|
|
tag = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key
|
||
|
|
select case(tag)
|
||
|
|
case ('element','e','el')
|
||
|
|
debug_e = IO_intValue(line,chunkPos,2_pInt)
|
||
|
|
case ('integrationpoint','i','ip')
|
||
|
|
debug_i = IO_intValue(line,chunkPos,2_pInt)
|
||
|
|
case ('grain','g','gr')
|
||
|
|
debug_g = IO_intValue(line,chunkPos,2_pInt)
|
||
|
|
end select
|
||
|
|
|
||
|
|
what = 0_pInt
|
||
|
|
select case(tag)
|
||
|
|
case ('debug')
|
||
|
|
what = debug_DEBUG
|
||
|
|
case ('math')
|
||
|
|
what = debug_MATH
|
||
|
|
case ('fesolving', 'fe')
|
||
|
|
what = debug_FESOLVING
|
||
|
|
case ('mesh')
|
||
|
|
what = debug_MESH
|
||
|
|
case ('material')
|
||
|
|
what = debug_MATERIAL
|
||
|
|
case ('lattice')
|
||
|
|
what = debug_LATTICE
|
||
|
|
case ('constitutive')
|
||
|
|
what = debug_CONSTITUTIVE
|
||
|
|
case ('crystallite')
|
||
|
|
what = debug_CRYSTALLITE
|
||
|
|
case ('homogenization')
|
||
|
|
what = debug_HOMOGENIZATION
|
||
|
|
case ('cpfem')
|
||
|
|
what = debug_CPFEM
|
||
|
|
case ('spectral')
|
||
|
|
what = debug_SPECTRAL
|
||
|
|
case ('marc')
|
||
|
|
what = debug_MARC
|
||
|
|
case ('abaqus')
|
||
|
|
what = debug_ABAQUS
|
||
|
|
case ('all')
|
||
|
|
what = debug_MAXNTYPE + 1_pInt
|
||
|
|
case ('other')
|
||
|
|
what = debug_MAXNTYPE + 2_pInt
|
||
|
|
end select
|
||
|
|
if (what /= 0) then
|
||
|
|
do i = 2_pInt, chunkPos(1)
|
||
|
|
select case(IO_lc(IO_stringValue(line,chunkPos,i)))
|
||
|
|
case('basic')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_LEVELBASIC)
|
||
|
|
case('extensive')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_LEVELEXTENSIVE)
|
||
|
|
case('selective')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_LEVELSELECTIVE)
|
||
|
|
case('restart')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_SPECTRALRESTART)
|
||
|
|
case('fft','fftw')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_SPECTRALFFTW)
|
||
|
|
case('divergence')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_SPECTRALDIVERGENCE)
|
||
|
|
case('rotation')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_SPECTRALROTATION)
|
||
|
|
case('petsc')
|
||
|
|
debug_level(what) = ior(debug_level(what), debug_SPECTRALPETSC)
|
||
|
|
end select
|
||
|
|
enddo
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
close(FILEUNIT)
|
||
|
|
|
||
|
|
do i = 1_pInt, debug_maxNtype
|
||
|
|
if (debug_level(i) == 0) &
|
||
|
|
debug_level(i) = ior(debug_level(i), debug_level(debug_MAXNTYPE + 2_pInt)) ! fill undefined debug types with levels specified by "other"
|
||
|
|
|
||
|
|
debug_level(i) = ior(debug_level(i), debug_level(debug_MAXNTYPE + 1_pInt)) ! fill all debug types with levels specified by "all"
|
||
|
|
enddo
|
||
|
|
|
||
|
|
if (iand(debug_level(debug_debug),debug_LEVELBASIC) /= 0) &
|
||
|
|
write(6,'(a,/)') ' using values from config file'
|
||
|
|
else fileExists
|
||
|
|
if (iand(debug_level(debug_debug),debug_LEVELBASIC) /= 0) &
|
||
|
|
write(6,'(a,/)') ' using standard values'
|
||
|
|
endif fileExists
|
||
|
|
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
! output switched on (debug level for debug must be extensive)
|
||
|
|
if (iand(debug_level(debug_debug),debug_LEVELEXTENSIVE) /= 0) then
|
||
|
|
do i = 1_pInt, debug_MAXNTYPE
|
||
|
|
select case(i)
|
||
|
|
case (debug_DEBUG)
|
||
|
|
tag = ' Debug'
|
||
|
|
case (debug_MATH)
|
||
|
|
tag = ' Math'
|
||
|
|
case (debug_FESOLVING)
|
||
|
|
tag = ' FEsolving'
|
||
|
|
case (debug_MESH)
|
||
|
|
tag = ' Mesh'
|
||
|
|
case (debug_MATERIAL)
|
||
|
|
tag = ' Material'
|
||
|
|
case (debug_LATTICE)
|
||
|
|
tag = ' Lattice'
|
||
|
|
case (debug_CONSTITUTIVE)
|
||
|
|
tag = ' Constitutive'
|
||
|
|
case (debug_CRYSTALLITE)
|
||
|
|
tag = ' Crystallite'
|
||
|
|
case (debug_HOMOGENIZATION)
|
||
|
|
tag = ' Homogenizaiton'
|
||
|
|
case (debug_CPFEM)
|
||
|
|
tag = ' CPFEM'
|
||
|
|
case (debug_SPECTRAL)
|
||
|
|
tag = ' Spectral solver'
|
||
|
|
case (debug_MARC)
|
||
|
|
tag = ' MSC.MARC FEM solver'
|
||
|
|
case (debug_ABAQUS)
|
||
|
|
tag = ' ABAQUS FEM solver'
|
||
|
|
end select
|
||
|
|
|
||
|
|
if(debug_level(i) /= 0) then
|
||
|
|
write(6,'(3a)') ' debug level for ', trim(tag), ':'
|
||
|
|
if(iand(debug_level(i),debug_LEVELBASIC) /= 0) write(6,'(a)') ' basic'
|
||
|
|
if(iand(debug_level(i),debug_LEVELEXTENSIVE) /= 0) write(6,'(a)') ' extensive'
|
||
|
|
if(iand(debug_level(i),debug_LEVELSELECTIVE) /= 0) then
|
||
|
|
write(6,'(a)') ' selective on:'
|
||
|
|
write(6,'(a24,1x,i8)') ' element: ',debug_e
|
||
|
|
write(6,'(a24,1x,i8)') ' ip: ',debug_i
|
||
|
|
write(6,'(a24,1x,i8)') ' grain: ',debug_g
|
||
|
|
endif
|
||
|
|
if(iand(debug_level(i),debug_SPECTRALRESTART) /= 0) write(6,'(a)') ' restart'
|
||
|
|
if(iand(debug_level(i),debug_SPECTRALFFTW) /= 0) write(6,'(a)') ' FFTW'
|
||
|
|
if(iand(debug_level(i),debug_SPECTRALDIVERGENCE)/= 0) write(6,'(a)') ' divergence'
|
||
|
|
if(iand(debug_level(i),debug_SPECTRALROTATION) /= 0) write(6,'(a)') ' rotation'
|
||
|
|
if(iand(debug_level(i),debug_SPECTRALPETSC) /= 0) write(6,'(a)') ' PETSc'
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
endif
|
||
|
|
|
||
|
|
end subroutine debug_init
|
||
|
|
|
||
|
|
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
!> @brief resets all debug values
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
subroutine debug_reset
|
||
|
|
|
||
|
|
implicit none
|
||
|
|
|
||
|
|
debug_StressLoopLpDistribution = 0_pInt
|
||
|
|
debug_StressLoopLiDistribution = 0_pInt
|
||
|
|
debug_StateLoopDistribution = 0_pInt
|
||
|
|
debug_CrystalliteLoopDistribution = 0_pInt
|
||
|
|
debug_MaterialpointStateLoopDistribution = 0_pInt
|
||
|
|
debug_MaterialpointLoopDistribution = 0_pInt
|
||
|
|
debug_cumLpTicks = 0_pLongInt
|
||
|
|
debug_cumDeltaStateTicks = 0_pLongInt
|
||
|
|
debug_cumDotStateTicks = 0_pLongInt
|
||
|
|
debug_cumLpCalls = 0_pInt
|
||
|
|
debug_cumDeltaStateCalls = 0_pInt
|
||
|
|
debug_cumDotStateCalls = 0_pInt
|
||
|
|
debug_stressMaxLocation = 0_pInt
|
||
|
|
debug_stressMinLocation = 0_pInt
|
||
|
|
debug_jacobianMaxLocation = 0_pInt
|
||
|
|
debug_jacobianMinLocation = 0_pInt
|
||
|
|
debug_stressMax = -huge(1.0_pReal)
|
||
|
|
debug_stressMin = huge(1.0_pReal)
|
||
|
|
debug_jacobianMax = -huge(1.0_pReal)
|
||
|
|
debug_jacobianMin = huge(1.0_pReal)
|
||
|
|
|
||
|
|
end subroutine debug_reset
|
||
|
|
|
||
|
|
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
!> @brief writes debug statements to standard out
|
||
|
|
!--------------------------------------------------------------------------------------------------
|
||
|
|
subroutine debug_info
|
||
|
|
use numerics, only: &
|
||
|
|
nStress, &
|
||
|
|
nState, &
|
||
|
|
nCryst, &
|
||
|
|
nMPstate, &
|
||
|
|
nHomog
|
||
|
|
|
||
|
|
implicit none
|
||
|
|
integer(pInt) :: j,integral
|
||
|
|
integer(pLongInt) :: tickrate
|
||
|
|
character(len=1) :: exceed
|
||
|
|
|
||
|
|
call system_clock(count_rate=tickrate)
|
||
|
|
|
||
|
|
!$OMP CRITICAL (write2out)
|
||
|
|
debugOutputCryst: if (iand(debug_level(debug_CRYSTALLITE),debug_LEVELBASIC) /= 0) then
|
||
|
|
write(6,'(/,a,/)') ' DEBUG Info (from previous cycle)'
|
||
|
|
write(6,'(a33,1x,i12)') 'total calls to LpAndItsTangent :',debug_cumLpCalls
|
||
|
|
if (debug_cumLpCalls > 0_pInt) then
|
||
|
|
write(6,'(a33,1x,f12.3)') 'total CPU time/s :',&
|
||
|
|
real(debug_cumLpTicks,pReal)/real(tickrate,pReal)
|
||
|
|
write(6,'(a33,1x,f12.6)') 'avg CPU time/microsecs per call :',&
|
||
|
|
real(debug_cumLpTicks,pReal)*1.0e6_pReal/real(tickrate*debug_cumLpCalls,pReal)
|
||
|
|
endif
|
||
|
|
write(6,'(/,a33,1x,i12)') 'total calls to collectDotState :',debug_cumDotStateCalls
|
||
|
|
if (debug_cumdotStateCalls > 0_pInt) then
|
||
|
|
write(6,'(a33,1x,f12.3)') 'total CPU time/s :',&
|
||
|
|
real(debug_cumDotStateTicks,pReal)/real(tickrate,pReal)
|
||
|
|
write(6,'(a33,1x,f12.6)') 'avg CPU time/microsecs per call :',&
|
||
|
|
real(debug_cumDotStateTicks,pReal)*1.0e6_pReal/real(tickrate*debug_cumDotStateCalls,pReal)
|
||
|
|
endif
|
||
|
|
write(6,'(/,a33,1x,i12)') 'total calls to collectDeltaState:',debug_cumDeltaStateCalls
|
||
|
|
if (debug_cumDeltaStateCalls > 0_pInt) then
|
||
|
|
write(6,'(a33,1x,f12.3)') 'total CPU time/s :',&
|
||
|
|
real(debug_cumDeltaStateTicks,pReal)/real(tickrate,pReal)
|
||
|
|
write(6,'(a33,1x,f12.6)') 'avg CPU time/microsecs per call :',&
|
||
|
|
real(debug_cumDeltaStateTicks,pReal)*1.0e6_pReal/real(tickrate*debug_cumDeltaStateCalls,pReal)
|
||
|
|
endif
|
||
|
|
|
||
|
|
integral = 0_pInt
|
||
|
|
write(6,'(3/,a)') 'distribution_StressLoopLp : stress stiffness'
|
||
|
|
do j=1_pInt,nStress+1_pInt
|
||
|
|
if (any(debug_StressLoopLpDistribution(j,:) /= 0_pInt )) then
|
||
|
|
integral = integral + j*(debug_StressLoopLpDistribution(j,1) + debug_StressLoopLpDistribution(j,2))
|
||
|
|
exceed = ' '
|
||
|
|
if (j > nStress) exceed = '+' ! last entry gets "+"
|
||
|
|
write(6,'(i25,a1,i10,1x,i10)') min(nStress,j),exceed,debug_StressLoopLpDistribution(j,1),&
|
||
|
|
debug_StressLoopLpDistribution(j,2)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
write(6,'(a15,i10,2(1x,i10))') ' total',integral,sum(debug_StressLoopLpDistribution(:,1)), &
|
||
|
|
sum(debug_StressLoopLpDistribution(:,2))
|
||
|
|
|
||
|
|
integral = 0_pInt
|
||
|
|
write(6,'(3/,a)') 'distribution_StressLoopLi : stress stiffness'
|
||
|
|
do j=1_pInt,nStress+1_pInt
|
||
|
|
if (any(debug_StressLoopLiDistribution(j,:) /= 0_pInt )) then
|
||
|
|
integral = integral + j*(debug_StressLoopLiDistribution(j,1) + debug_StressLoopLiDistribution(j,2))
|
||
|
|
exceed = ' '
|
||
|
|
if (j > nStress) exceed = '+' ! last entry gets "+"
|
||
|
|
write(6,'(i25,a1,i10,1x,i10)') min(nStress,j),exceed,debug_StressLoopLiDistribution(j,1),&
|
||
|
|
debug_StressLoopLiDistribution(j,2)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
write(6,'(a15,i10,2(1x,i10))') ' total',integral,sum(debug_StressLoopLiDistribution(:,1)), &
|
||
|
|
sum(debug_StressLoopLiDistribution(:,2))
|
||
|
|
|
||
|
|
integral = 0_pInt
|
||
|
|
write(6,'(2/,a)') 'distribution_CrystalliteStateLoop :'
|
||
|
|
do j=1_pInt,nState+1_pInt
|
||
|
|
if (any(debug_StateLoopDistribution(j,:) /= 0)) then
|
||
|
|
integral = integral + j*(debug_StateLoopDistribution(j,1) + debug_StateLoopDistribution(j,2))
|
||
|
|
exceed = ' '
|
||
|
|
if (j > nState) exceed = '+' ! last entry gets "+"
|
||
|
|
write(6,'(i25,a1,i10,1x,i10)') min(nState,j),exceed,debug_StateLoopDistribution(j,1),&
|
||
|
|
debug_StateLoopDistribution(j,2)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
write(6,'(a15,i10,2(1x,i10))') ' total',integral,sum(debug_StateLoopDistribution(:,1)), &
|
||
|
|
sum(debug_StateLoopDistribution(:,2))
|
||
|
|
|
||
|
|
integral = 0_pInt
|
||
|
|
write(6,'(2/,a)') 'distribution_CrystalliteCutbackLoop :'
|
||
|
|
do j=1_pInt,nCryst+1_pInt
|
||
|
|
if (debug_CrystalliteLoopDistribution(j) /= 0) then
|
||
|
|
integral = integral + j*debug_CrystalliteLoopDistribution(j)
|
||
|
|
exceed = ' '
|
||
|
|
if (j > nCryst) exceed = '+'
|
||
|
|
write(6,'(i25,a1,i10)') min(nCryst,j),exceed,debug_CrystalliteLoopDistribution(j)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
write(6,'(a15,i10,1x,i10)') ' total',integral,sum(debug_CrystalliteLoopDistribution)
|
||
|
|
endif debugOutputCryst
|
||
|
|
|
||
|
|
debugOutputHomog: if (iand(debug_level(debug_HOMOGENIZATION),debug_LEVELBASIC) /= 0) then
|
||
|
|
integral = 0_pInt
|
||
|
|
write(6,'(2/,a)') 'distribution_MaterialpointStateLoop :'
|
||
|
|
do j=1_pInt,nMPstate
|
||
|
|
if (debug_MaterialpointStateLoopDistribution(j) /= 0) then
|
||
|
|
integral = integral + j*debug_MaterialpointStateLoopDistribution(j)
|
||
|
|
write(6,'(i25,1x,i10)') j,debug_MaterialpointStateLoopDistribution(j)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
write(6,'(a15,i10,1x,i10)') ' total',integral,sum(debug_MaterialpointStateLoopDistribution)
|
||
|
|
|
||
|
|
integral = 0_pInt
|
||
|
|
write(6,'(2/,a)') 'distribution_MaterialpointCutbackLoop :'
|
||
|
|
do j=1_pInt,nHomog+1_pInt
|
||
|
|
if (debug_MaterialpointLoopDistribution(j) /= 0) then
|
||
|
|
integral = integral + j*debug_MaterialpointLoopDistribution(j)
|
||
|
|
exceed = ' '
|
||
|
|
if (j > nHomog) exceed = '+'
|
||
|
|
write(6,'(i25,a1,i10)') min(nHomog,j),exceed,debug_MaterialpointLoopDistribution(j)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
write(6,'(a15,i10,1x,i10)') ' total',integral,sum(debug_MaterialpointLoopDistribution)
|
||
|
|
endif debugOutputHomog
|
||
|
|
|
||
|
|
debugOutputCPFEM: if (iand(debug_level(debug_CPFEM),debug_LEVELBASIC) /= 0) then
|
||
|
|
write(6,'(2/,a,/)') ' Extreme values of returned stress and jacobian'
|
||
|
|
write(6,'(a39)') ' value el ip'
|
||
|
|
write(6,'(a14,1x,e12.3,1x,i6,1x,i4)') ' stress min :', debug_stressMin, debug_stressMinLocation
|
||
|
|
write(6,'(a14,1x,e12.3,1x,i6,1x,i4)') ' max :', debug_stressMax, debug_stressMaxLocation
|
||
|
|
write(6,'(a14,1x,e12.3,1x,i6,1x,i4)') ' jacobian min :', debug_jacobianMin, debug_jacobianMinLocation
|
||
|
|
write(6,'(a14,1x,e12.3,1x,i6,1x,i4,/)') ' max :', debug_jacobianMax, debug_jacobianMaxLocation
|
||
|
|
endif debugOutputCPFEM
|
||
|
|
!$OMP END CRITICAL (write2out)
|
||
|
|
|
||
|
|
end subroutine debug_info
|
||
|
|
|
||
|
|
end module debug
|