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

601 lines
39 KiB
Fortran
Raw Normal View History

further work on modularization, basic scheme is now calculating the same results as DAMASK_spectral.f90
2012-07-30 19:36:22 +05:30
!--------------------------------------------------------------------------------------------------
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
! $Id$
further work on modularization, basic scheme is now calculating the same results as DAMASK_spectral.f90
2012-07-30 19:36:22 +05:30
!--------------------------------------------------------------------------------------------------
!> @author Pratheek Shanthraj, Max-Planck-Institut für Eisenforschung GmbH
!> @author Martin Diehl, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief Driver controlling inner and outer load case looping of the various spectral solvers
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
!> @details doing cutbacking, forwarding in case of restart, reporting statistics, writing
!> results
further work on modularization, basic scheme is now calculating the same results as DAMASK_spectral.f90
2012-07-30 19:36:22 +05:30
!--------------------------------------------------------------------------------------------------
next round in modularizzation of solvers
2012-07-24 22:37:10 +05:30
program DAMASK_spectral_Driver
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
use, intrinsic :: &
iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran >4.6 at the moment)
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
use prec, only: &
pInt, &
some comments and minor improvements. removed tol_gravityNodePos from prec because it's not used tol_math_check is now used by spectral driver to check rotation
2013-09-14 16:29:35 +05:30
pReal, &
tol_math_check
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
use DAMASK_interface, only: &
DAMASK_interface_init, &
loadCaseFile, &
geometryFile, &
getSolverWorkingDirectoryName, &
getSolverJobName, &
appendToOutFile
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
use IO, only: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
IO_read, &
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
IO_isBlank, &
IO_open_file, &
IO_stringPos, &
IO_stringValue, &
IO_floatValue, &
IO_intValue, &
IO_error, &
IO_lc, &
added some warning and explicit size of arrays
2013-01-10 19:03:43 +05:30
IO_intOut, &
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 21:22:29 +05:30
IO_warning, &
continued with sourced allocation and EOL improvements
2013-12-16 17:28:03 +05:30
IO_timeStamp, &
IO_EOF
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 21:22:29 +05:30
use debug, only: &
debug_level, &
debug_spectral, &
debug_levelBasic
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
use math ! need to include the whole module for FFTW
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
use CPFEM, only: &
CPFEM_initAll
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
use FEsolving, only: &
restartWrite, &
restartInc
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
use numerics, only: &
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
maxCutBack, &
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
spectral_solver, &
changed handling of non-converged BVP solution: By default, exit if no solution is found to prevent "user errors" Removed ambiguous "regridMode" keyword which was used to trigger this behavior in favor of "continueCalculation", Set this to 1 to get back old behavior, e.g. report non-converged increments.
2014-03-31 15:34:11 +05:30
continueCalculation
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
use homogenization, only: &
materialpoint_sizeResults, &
materialpoint_results
new version of modular solver structure
2012-08-03 14:55:48 +05:30
use DAMASK_spectral_Utilities, only: &
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 21:22:29 +05:30
grid, &
geomSize, &
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
tBoundaryCondition, &
tSolutionState, &
rolled back phase field changes
2013-12-18 14:39:32 +05:30
cutBack
next round in modularizzation of solvers
2012-07-24 22:37:10 +05:30
use DAMASK_spectral_SolverBasic
added SolverBasicPETSC module to the new FFT solver (solves the basic scheme using PETSc when mysolver = basicpetsc in numerics.config) some bug fixes
2012-08-14 22:28:23 +05:30
use DAMASK_spectral_SolverBasicPETSC
updates for AL version of spectral solver
2012-08-06 14:23:12 +05:30
use DAMASK_spectral_SolverAL
added Polarisation scheme as solver
2013-07-24 18:36:16 +05:30
use DAMASK_spectral_SolverPolarisation
specified quit as external as it is not part of a module
2013-02-20 20:07:12 +05:30
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
implicit none
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
type tLoadCase
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
real(pReal), dimension (3,3) :: rotation = math_I3 !< rotation of BC
type(tBoundaryCondition) :: P, & !< stress BC
deformation !< deformation BC (Fdot or L)
real(pReal) :: time = 0.0_pReal, & !< length of increment
rolled back phase field changes
2013-12-18 14:39:32 +05:30
temperature = 300.0_pReal, & !< isothermal starting conditions
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms. dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!). not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 21:55:37 +05:30
density = 0.0_pReal !< density
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
integer(pInt) :: incs = 0_pInt, & !< number of increments
outputfrequency = 1_pInt, & !< frequency of result writes
restartfrequency = 0_pInt, & !< frequency of restart writes
minor changes on buggy latest commit (wrong use statement)
2013-01-09 03:24:25 +05:30
logscale = 0_pInt !< linear/logarithmic time inc flag
rolled back phase field changes
2013-12-18 14:39:32 +05:30
logical :: followFormerTrajectory = .true. !< follow trajectory of former loadcase
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
end type tLoadCase
further modularisation of spectral solver (only for basic scheme so far)
2012-07-23 15:42:31 +05:30
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
!--------------------------------------------------------------------------------------------------
! variables related to information from load case and geom file
some more finetuning, still not running on gfortran 4.5.0
2012-11-09 03:03:58 +05:30
real(pReal), dimension(9) :: temp_valueVector = 0.0_pReal !< temporarily from loadcase file when reading in tensors (initialize to 0.0)
logical, dimension(9) :: temp_maskVector = .false. !< temporarily from loadcase file when reading in tensors
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
integer(pInt), parameter :: maxNchunks = (1_pInt + 9_pInt)*3_pInt + & ! deformation, rotation, and stress
(1_pInt + 1_pInt)*5_pInt + & ! time, (log)incs, temp, restartfrequency, and outputfrequency
1_pInt, & ! dropguessing
continued with sourced allocation and EOL improvements
2013-12-16 17:28:03 +05:30
FILEUNIT = 234_pInt !< file unit, DAMASK IO does not support newunit feature
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
integer(pInt), dimension(1_pInt + maxNchunks*2_pInt) :: positions ! this is longer than needed for geometry parsing
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
integer(pInt) :: &
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
N_t = 0_pInt, & !< # of time indicators found in load case file
N_n = 0_pInt, & !< # of increment specifiers found in load case file
rolled back phase field changes
2013-12-18 14:39:32 +05:30
N_def = 0_pInt !< # of rate of deformation specifiers found in load case file
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
character(len=65536) :: &
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
line
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
!--------------------------------------------------------------------------------------------------
! loop variables, convergence etc.
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
real(pReal), dimension(3,3), parameter :: &
ones = 1.0_pReal, &
zeros = 0.0_pReal
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
integer(pInt), parameter :: &
subStepFactor = 2_pInt !< for each substep, divide the last time increment by 2.0
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
real(pReal) :: &
time = 0.0_pReal, & !< elapsed time
time0 = 0.0_pReal, & !< begin of interval
timeinc = 1.0_pReal, & !< current time interval
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
timeIncOld = 0.0_pReal, & !< previous time interval
remainingLoadCaseTime = 0.0_pReal !< remaining time of current load case
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
logical :: &
guess !< guess along former trajectory
integer(pInt) :: &
i, j, k, l, &
errorID, &
cutBackLevel = 0_pInt, & !< cut back level \f$ t = \frac{t_{inc}}{2^l} \f$
stepFraction = 0_pInt !< fraction of current time interval
integer(pInt) :: &
currentLoadcase = 0_pInt, & !< current load case
inc, & !< current increment in current load case
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
totalIncsCounter = 0_pInt, & !< total # of increments
convergedCounter = 0_pInt, & !< # of converged increments
notConvergedCounter = 0_pInt, & !< # of non-converged increments
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
resUnit = 0_pInt, & !< file unit for results writing
statUnit = 0_pInt, & !< file unit for statistics output
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
lastRestartWritten = 0_pInt !< total increment # at which last restart information was written
further modularisation of spectral solver (only for basic scheme so far)
2012-07-23 15:42:31 +05:30
character(len=6) :: loadcase_string
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
character(len=1024) :: incInfo !< string parsed to solution with information about current load case
type(tLoadCase), allocatable, dimension(:) :: loadCases !< array of all load cases
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
type(tSolutionState) solres
specified quit as external as it is not part of a module
2013-02-20 20:07:12 +05:30
external :: quit
new version of modular solver structure
2012-08-03 14:55:48 +05:30
!--------------------------------------------------------------------------------------------------
! init DAMASK (all modules)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
call CPFEM_initAll(temperature = 300.0_pReal, el = 1_pInt, ip = 1_pInt)
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 21:22:29 +05:30
write(6,'(/,a)') ' <<<+- DAMASK_spectral_driver init -+>>>'
write(6,'(a)') ' $Id$'
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
#include "compilation_info.f90"
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
!--------------------------------------------------------------------------------------------------
new version of modular solver structure
2012-08-03 14:55:48 +05:30
! reading basic information from load case file and allocate data structure containing load cases
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
call IO_open_file(FILEUNIT,trim(loadCaseFile))
rewind(FILEUNIT)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
do
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
line = IO_read(FILEUNIT)
continued with sourced allocation and EOL improvements
2013-12-16 17:28:03 +05:30
if (trim(line) == IO_EOF) exit
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
positions = IO_stringPos(line,maxNchunks)
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
do i = 1_pInt, positions(1) ! reading compulsory parameters for loadcase
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
select case (IO_lc(IO_stringValue(line,positions,i)))
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
case('l','velocitygrad','velgrad','velocitygradient','fdot','dotf','f')
N_def = N_def + 1_pInt
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
case('t','time','delta')
N_t = N_t + 1_pInt
case('n','incs','increments','steps','logincs','logincrements','logsteps')
N_n = N_n + 1_pInt
end select
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
enddo ! count all identifiers to allocate memory and do sanity check
enddo
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
if ((N_def /= N_n) .or. (N_n /= N_t)) & ! sanity check
call IO_error(error_ID=837_pInt,ext_msg = trim(loadCaseFile)) ! error message for incomplete loadcase
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
allocate (loadCases(N_n)) ! array of load cases
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases%P%myType='p'
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
!--------------------------------------------------------------------------------------------------
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
! reading the load case and assign values to the allocated data structure
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
rewind(FILEUNIT)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
do
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
line = IO_read(FILEUNIT)
continued with sourced allocation and EOL improvements
2013-12-16 17:28:03 +05:30
if (trim(line) == IO_EOF) exit
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
new version of modular solver structure
2012-08-03 14:55:48 +05:30
currentLoadCase = currentLoadCase + 1_pInt
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
positions = IO_stringPos(line,maxNchunks)
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 22:15:34 +05:30
do i = 1_pInt, positions(1)
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
select case (IO_lc(IO_stringValue(line,positions,i)))
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
case('fdot','dotf','l','velocitygrad','velgrad','velocitygradient','f') ! assign values for the deformation BC matrix
reseting value vector when reading in for loadcase file, otherwise NaNs might be in there
2012-11-09 02:05:31 +05:30
temp_valueVector = 0.0_pReal
simplified load case parsing
2012-11-29 18:56:17 +05:30
if (IO_lc(IO_stringValue(line,positions,i)) == 'fdot'.or. & ! in case of Fdot, set type to fdot
IO_lc(IO_stringValue(line,positions,i)) == 'dotf') then
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%deformation%myType = 'fdot'
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
else if (IO_lc(IO_stringValue(line,positions,i)) == 'f') then
loadCases(currentLoadCase)%deformation%myType = 'f'
simplified load case parsing
2012-11-29 18:56:17 +05:30
else
loadCases(currentLoadCase)%deformation%myType = 'l'
new version of modular solver structure
2012-08-03 14:55:48 +05:30
endif
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 22:15:34 +05:30
do j = 1_pInt, 9_pInt
temp_maskVector(j) = IO_stringValue(line,positions,i+j) /= '*' ! true if not a *
enddo
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
do j = 1_pInt,9_pInt
if (temp_maskVector(j)) temp_valueVector(j) = IO_floatValue(line,positions,i+j) ! read value where applicable
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
enddo
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
loadCases(currentLoadCase)%deformation%maskLogical = & ! logical mask in 3x3 notation
transpose(reshape(temp_maskVector,[ 3,3]))
loadCases(currentLoadCase)%deformation%maskFloat = & ! float (1.0/0.0) mask in 3x3 notation
merge(ones,zeros,loadCases(currentLoadCase)%deformation%maskLogical)
loadCases(currentLoadCase)%deformation%values = math_plain9to33(temp_valueVector) ! values in 3x3 notation
introduced error code for run_test.py
2013-01-24 00:03:46 +05:30
case('p','pk1','piolakirchhoff','stress', 's')
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
temp_valueVector = 0.0_pReal
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 22:15:34 +05:30
do j = 1_pInt, 9_pInt
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
temp_maskVector(j) = IO_stringValue(line,positions,i+j) /= '*' ! true if not an asterisk
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 22:15:34 +05:30
enddo
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
do j = 1_pInt,9_pInt
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
if (temp_maskVector(j)) temp_valueVector(j) = IO_floatValue(line,positions,i+j) ! read value where applicable
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
enddo
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%P%maskLogical = transpose(reshape(temp_maskVector,[ 3,3]))
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
loadCases(currentLoadCase)%P%maskFloat = merge(ones,zeros,&
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%P%maskLogical)
loadCases(currentLoadCase)%P%values = math_plain9to33(temp_valueVector)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
case('t','time','delta') ! increment time
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%time = IO_floatValue(line,positions,i+1_pInt)
rolled back phase field changes
2013-12-18 14:39:32 +05:30
case('temp','temperature') ! starting temperature
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%temperature = IO_floatValue(line,positions,i+1_pInt)
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms. dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!). not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 21:55:37 +05:30
case('den','density') ! starting density
loadCases(currentLoadCase)%density = IO_floatValue(line,positions,i+1_pInt)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
case('n','incs','increments','steps') ! number of increments
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%incs = IO_intValue(line,positions,i+1_pInt)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
case('logincs','logincrements','logsteps') ! number of increments (switch to log time scaling)
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%incs = IO_intValue(line,positions,i+1_pInt)
loadCases(currentLoadCase)%logscale = 1_pInt
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
case('freq','frequency','outputfreq') ! frequency of result writings
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%outputfrequency = IO_intValue(line,positions,i+1_pInt)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
case('r','restart','restartwrite') ! frequency of writing restart information
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
loadCases(currentLoadCase)%restartfrequency = &
max(0_pInt,IO_intValue(line,positions,i+1_pInt))
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
case('guessreset','dropguessing')
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
loadCases(currentLoadCase)%followFormerTrajectory = .false. ! do not continue to predict deformation along former trajectory
new version of modular solver structure
2012-08-03 14:55:48 +05:30
case('euler') ! rotation of currentLoadCase given in euler angles
some more finetuning, still not running on gfortran 4.5.0
2012-11-09 03:03:58 +05:30
temp_valueVector = 0.0_pReal
simplified the rate calculation interface in DAMASK_spectral_utilities.f90, corrected bug in rotation parsing and added debug option "rotation" to spectral solver options. This will show the current average stress and deformation additionally in lab coordinate system
2012-11-28 20:34:05 +05:30
l = 1_pInt ! assuming values given in degrees
k = 1_pInt ! assuming keyword indicating degree/radians present
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
select case (IO_lc(IO_stringValue(line,positions,i+1_pInt)))
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
case('deg','degree')
simplified the rate calculation interface in DAMASK_spectral_utilities.f90, corrected bug in rotation parsing and added debug option "rotation" to spectral solver options. This will show the current average stress and deformation additionally in lab coordinate system
2012-11-28 20:34:05 +05:30
case('rad','radian') ! don't convert from degree to radian
l = 0_pInt
case default
k = 0_pInt
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
end select
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 22:15:34 +05:30
do j = 1_pInt, 3_pInt
temp_valueVector(j) = IO_floatValue(line,positions,i+k+j)
enddo
simplified the rate calculation interface in DAMASK_spectral_utilities.f90, corrected bug in rotation parsing and added debug option "rotation" to spectral solver options. This will show the current average stress and deformation additionally in lab coordinate system
2012-11-28 20:34:05 +05:30
if (l == 1_pInt) temp_valueVector(1:3) = temp_valueVector(1:3) * inRad ! convert to rad
fixed bug in loadcase rotation: when not specifying "deg" the rotation was reset to 0,0,0..!
2012-09-06 19:35:28 +05:30
loadCases(currentLoadCase)%rotation = math_EulerToR(temp_valueVector(1:3)) ! convert rad Eulers to rotation matrix
new version of modular solver structure
2012-08-03 14:55:48 +05:30
case('rotation','rot') ! assign values for the rotation of currentLoadCase matrix
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
temp_valueVector = 0.0_pReal
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 22:15:34 +05:30
do j = 1_pInt, 9_pInt
temp_valueVector(j) = IO_floatValue(line,positions,i+j)
enddo
new version of modular solver structure
2012-08-03 14:55:48 +05:30
loadCases(currentLoadCase)%rotation = math_plain9to33(temp_valueVector)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
end select
enddo; enddo
rolled back phase field changes
2013-12-18 14:39:32 +05:30
close(FILEUNIT)
new version of modular solver structure
2012-08-03 14:55:48 +05:30
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
!--------------------------------------------------------------------------------------------------
! consistency checks and output of load case
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
loadCases(1)%followFormerTrajectory = .false. ! cannot guess along trajectory for first inc of first currentLoadCase
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
errorID = 0_pInt
new version of modular solver structure
2012-08-03 14:55:48 +05:30
checkLoadcases: do currentLoadCase = 1_pInt, size(loadCases)
write (loadcase_string, '(i6)' ) currentLoadCase
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
write(6,'(1x,a,i6)') 'load case: ', currentLoadCase
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
if (.not. loadCases(currentLoadCase)%followFormerTrajectory) &
write(6,'(2x,a)') 'drop guessing along trajectory'
new version of modular solver structure
2012-08-03 14:55:48 +05:30
if (loadCases(currentLoadCase)%deformation%myType=='l') then
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
do j = 1_pInt, 3_pInt
new version of modular solver structure
2012-08-03 14:55:48 +05:30
if (any(loadCases(currentLoadCase)%deformation%maskLogical(j,1:3) .eqv. .true.) .and. &
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
any(loadCases(currentLoadCase)%deformation%maskLogical(j,1:3) .eqv. .false.)) &
errorID = 832_pInt ! each row should be either fully or not at all defined
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
enddo
new version of modular solver structure
2012-08-03 14:55:48 +05:30
write(6,'(2x,a)') 'velocity gradient:'
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
else if (loadCases(currentLoadCase)%deformation%myType=='f') then
write(6,'(2x,a)') 'deformation gradient at end of load case:'
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
else
new version of modular solver structure
2012-08-03 14:55:48 +05:30
write(6,'(2x,a)') 'deformation gradient rate:'
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
endif
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
write(6,'(3(3(3x,f12.7,1x)/))',advance='no') &
merge(math_transpose33(loadCases(currentLoadCase)%deformation%values), &
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
reshape(spread(huge(1.0_pReal),1,9),[ 3,3]), & ! print *** (huge) for undefined
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
transpose(loadCases(currentLoadCase)%deformation%maskLogical))
if (any(loadCases(currentLoadCase)%P%maskLogical .eqv. &
loadCases(currentLoadCase)%deformation%maskLogical)) errorID = 831_pInt ! exclusive or masking only
if (any(loadCases(currentLoadCase)%P%maskLogical .and. &
transpose(loadCases(currentLoadCase)%P%maskLogical) .and. &
reshape([ .false.,.true.,.true.,.true.,.false.,.true.,.true.,.true.,.false.],[ 3,3]))) &
errorID = 838_pInt ! no rotation is allowed by stress BC
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
write(6,'(2x,a,/,3(3(3x,f12.7,1x)/))',advance='no') 'stress / GPa:',&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
1e-9_pReal*merge(math_transpose33(loadCases(currentLoadCase)%P%values),&
reshape(spread(huge(1.0_pReal),1,9),[ 3,3]),&
transpose(loadCases(currentLoadCase)%P%maskLogical))
if (any(abs(math_mul33x33(loadCases(currentLoadCase)%rotation, &
math_transpose33(loadCases(currentLoadCase)%rotation))-math_I3) >&
some comments and minor improvements. removed tol_gravityNodePos from prec because it's not used tol_math_check is now used by spectral driver to check rotation
2013-09-14 16:29:35 +05:30
reshape(spread(tol_math_check,1,9),[ 3,3]))&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
.or. abs(math_det33(loadCases(currentLoadCase)%rotation)) > &
some comments and minor improvements. removed tol_gravityNodePos from prec because it's not used tol_math_check is now used by spectral driver to check rotation
2013-09-14 16:29:35 +05:30
1.0_pReal + tol_math_check) errorID = 846_pInt ! given rotation matrix contains strain
new version of modular solver structure
2012-08-03 14:55:48 +05:30
if (any(loadCases(currentLoadCase)%rotation /= math_I3)) &
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
write(6,'(2x,a,/,3(3(3x,f12.7,1x)/))',advance='no') 'rotation of loadframe:',&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
math_transpose33(loadCases(currentLoadCase)%rotation)
new version of modular solver structure
2012-08-03 14:55:48 +05:30
write(6,'(2x,a,f12.6)') 'temperature:', loadCases(currentLoadCase)%temperature
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms. dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!). not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 21:55:37 +05:30
write(6,'(2x,a,f12.6)') 'density: ', loadCases(currentLoadCase)%density
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
if (loadCases(currentLoadCase)%time < 0.0_pReal) errorID = 834_pInt ! negative time increment
new version of modular solver structure
2012-08-03 14:55:48 +05:30
write(6,'(2x,a,f12.6)') 'time: ', loadCases(currentLoadCase)%time
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
if (loadCases(currentLoadCase)%incs < 1_pInt) errorID = 835_pInt ! non-positive incs count
new version of modular solver structure
2012-08-03 14:55:48 +05:30
write(6,'(2x,a,i5)') 'increments: ', loadCases(currentLoadCase)%incs
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
if (loadCases(currentLoadCase)%outputfrequency < 1_pInt) errorID = 836_pInt ! non-positive result frequency
write(6,'(2x,a,i5)') 'output frequency: ', &
loadCases(currentLoadCase)%outputfrequency
write(6,'(2x,a,i5,/)') 'restart frequency: ', &
loadCases(currentLoadCase)%restartfrequency
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
if (errorID > 0_pInt) call IO_error(error_ID = errorID, ext_msg = loadcase_string) ! exit with error message
next round in modularizzation of solvers
2012-07-24 22:37:10 +05:30
enddo checkLoadcases
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
!--------------------------------------------------------------------------------------------------
! doing initialization depending on selected solver
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
select case (spectral_solver)
new version of modular solver structure
2012-08-03 14:55:48 +05:30
case (DAMASK_spectral_SolverBasic_label)
fixed wrong temperature when using spectral solver
2013-01-09 23:38:08 +05:30
call basic_init(loadCases(1)%temperature)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
case (DAMASK_spectral_SolverBasicPETSc_label)
rolled back phase field changes
2013-12-18 14:39:32 +05:30
call basicPETSc_init(loadCases(1)%temperature)
updates for AL version of spectral solver
2012-08-06 14:23:12 +05:30
case (DAMASK_spectral_SolverAL_label)
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 21:22:29 +05:30
if(iand(debug_level(debug_spectral),debug_levelBasic)/= 0) &
call IO_warning(42_pInt, ext_msg='debug Divergence')
fixed wrong temperature when using spectral solver
2013-01-09 23:38:08 +05:30
call AL_init(loadCases(1)%temperature)
added Polarisation scheme as solver
2013-07-24 18:36:16 +05:30
case (DAMASK_spectral_SolverPolarisation_label)
if(iand(debug_level(debug_spectral),debug_levelBasic)/= 0) &
call IO_warning(42_pInt, ext_msg='debug Divergence')
call Polarisation_init(loadCases(1)%temperature)
new solver is now compiling without a PETSc installation, however only the plain basic solver is available then. also removed reporting of PETSc related variables in the case it is not installed/found
2012-08-29 00:49:47 +05:30
case default
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
call IO_error(error_ID = 891, ext_msg = trim(spectral_solver))
new version of modular solver structure
2012-08-03 14:55:48 +05:30
end select
new solver is now compiling without a PETSc installation, however only the plain basic solver is available then. also removed reporting of PETSc related variables in the case it is not installed/found
2012-08-29 00:49:47 +05:30
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
!--------------------------------------------------------------------------------------------------
! write header of output file
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
if (appendToOutFile) then ! after restart, append to existing results file
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
open(newunit=resUnit,file=trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
'.spectralOut',form='UNFORMATTED', position='APPEND', status='OLD')
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
open(newunit=statUnit,file=trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
'.sta',form='FORMATTED', position='APPEND', status='OLD')
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
else ! open new files ...
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
open(newunit=resUnit,file=trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
'.spectralOut',form='UNFORMATTED',status='REPLACE')
fixed interpretation 'spectralOut' header. Seem to have been buggy all the time, but manifested itself only after keyword change. Also added ':' after keywords to be more robust
2015-03-15 20:59:38 +05:30
write(resUnit) 'load:', trim(loadCaseFile) ! ... and write header
write(resUnit) 'workingdir:', trim(getSolverWorkingDirectoryName())
write(resUnit) 'geometry:', trim(geometryFile)
write(resUnit) 'grid:', grid
write(resUnit) 'size:', geomSize
write(resUnit) 'materialpoint_sizeResults:', materialpoint_sizeResults
write(resUnit) 'loadcases:', size(loadCases)
write(resUnit) 'frequencies:', loadCases%outputfrequency ! one entry per currentLoadCase
write(resUnit) 'times:', loadCases%time ! one entry per currentLoadCase
write(resUnit) 'logscales:', loadCases%logscale
write(resUnit) 'increments:', loadCases%incs ! one entry per currentLoadCase
write(resUnit) 'startingIncrement:', restartInc - 1_pInt ! start with writing out the previous inc
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
write(resUnit) 'eoh' ! end of header
fixed writing out of initial configuration for spectral solver in case of ngrains != 1, small formating improvements in other files
2013-07-01 12:10:09 +05:30
write(resUnit) materialpoint_results ! initial (non-deformed or read-in) results
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
open(newunit=statUnit,file=trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//&
'.sta',form='FORMATTED',status='REPLACE')
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
write(statUnit,'(a)') 'Increment Time CutbackLevel Converged IterationsNeeded' ! statistics file
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 21:22:29 +05:30
if (iand(debug_level(debug_spectral),debug_levelBasic) /= 0) &
write(6,'(/,a)') ' header of result file written out'
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
flush(6)
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
endif
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
new version of modular solver structure
2012-08-03 14:55:48 +05:30
!--------------------------------------------------------------------------------------------------
! loopping over loadcases
loadCaseLooping: do currentLoadCase = 1_pInt, size(loadCases)
time0 = time ! currentLoadCase start time
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
guess = loadCases(currentLoadCase)%followFormerTrajectory ! change of load case? homogeneous guess for the first inc
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
new version of modular solver structure
2012-08-03 14:55:48 +05:30
!--------------------------------------------------------------------------------------------------
! loop oper incs defined in input file for current currentLoadCase
incLooping: do inc = 1_pInt, loadCases(currentLoadCase)%incs
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
totalIncsCounter = totalIncsCounter + 1_pInt
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
!--------------------------------------------------------------------------------------------------
! forwarding time
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
timeIncOld = timeinc
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
if (loadCases(currentLoadCase)%logscale == 0_pInt) then ! linear scale
timeinc = loadCases(currentLoadCase)%time/loadCases(currentLoadCase)%incs ! only valid for given linear time scale. will be overwritten later in case loglinear scale is used
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
else
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
if (currentLoadCase == 1_pInt) then ! 1st currentLoadCase of logarithmic scale
new version of modular solver structure
2012-08-03 14:55:48 +05:30
if (inc == 1_pInt) then ! 1st inc of 1st currentLoadCase of logarithmic scale
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
timeinc = loadCases(1)%time*(2.0_pReal**real( 1_pInt-loadCases(1)%incs ,pReal)) ! assume 1st inc is equal to 2nd
new version of modular solver structure
2012-08-03 14:55:48 +05:30
else ! not-1st inc of 1st currentLoadCase of logarithmic scale
timeinc = loadCases(1)%time*(2.0_pReal**real(inc-1_pInt-loadCases(1)%incs ,pReal))
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
endif
new version of modular solver structure
2012-08-03 14:55:48 +05:30
else ! not-1st currentLoadCase of logarithmic scale
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
timeinc = time0 * &
( (1.0_pReal + loadCases(currentLoadCase)%time/time0 )**(real( inc,pReal)/&
real(loadCases(currentLoadCase)%incs ,pReal))&
-(1.0_pReal + loadCases(currentLoadCase)%time/time0 )**(real( (inc-1_pInt),pReal)/&
real(loadCases(currentLoadCase)%incs ,pReal)))
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
endif
endif
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
timeinc = timeinc / 2.0_pReal**real(cutBackLevel,pReal) ! depending on cut back level, decrease time step
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
doxygen comments for homogenization.f90, unified naming ip->i, el->e
2013-01-29 15:58:01 +05:30
forwarding: if(totalIncsCounter >= restartInc) then
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
stepFraction = 0_pInt
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
!--------------------------------------------------------------------------------------------------
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
! loop over sub incs
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
subIncLooping: do while (stepFraction/subStepFactor**cutBackLevel <1_pInt)
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
time = time + timeinc ! forward time
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
stepFraction = stepFraction + 1_pInt
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
remainingLoadCaseTime = time0 - time + loadCases(currentLoadCase)%time + timeInc
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
!--------------------------------------------------------------------------------------------------
separated forwarding and solution subroutines for better control at the load step looping level
2013-12-18 15:05:05 +05:30
! forward solution
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
select case(spectral_solver)
separated forwarding and solution subroutines for better control at the load step looping level
2013-12-18 15:05:05 +05:30
case (DAMASK_spectral_SolverBasic_label)
case (DAMASK_spectral_SolverBasicPETSC_label)
call BasicPETSC_forward (&
guess,timeinc,timeIncOld,remainingLoadCaseTime, &
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
rotation_BC = loadCases(currentLoadCase)%rotation)
case (DAMASK_spectral_SolverAL_label)
call AL_forward (&
guess,timeinc,timeIncOld,remainingLoadCaseTime, &
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
rotation_BC = loadCases(currentLoadCase)%rotation)
case (DAMASK_spectral_SolverPolarisation_label)
call Polarisation_forward (&
guess,timeinc,timeIncOld,remainingLoadCaseTime, &
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
rotation_BC = loadCases(currentLoadCase)%rotation)
fixed error in previous commit
2013-12-18 15:08:44 +05:30
end select
separated forwarding and solution subroutines for better control at the load step looping level
2013-12-18 15:05:05 +05:30
!--------------------------------------------------------------------------------------------------
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
! report begin of new increment
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
write(6,'(/,a)') ' ###########################################################################'
write(6,'(1x,a,es12.5'//&
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
',a,'//IO_intOut(inc)//',a,'//IO_intOut(loadCases(currentLoadCase)%incs)//&
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
',a,'//IO_intOut(stepFraction)//',a,'//IO_intOut(subStepFactor**cutBackLevel)//&
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
',a,'//IO_intOut(currentLoadCase)//',a,'//IO_intOut(size(loadCases))//')') &
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
'Time', time, &
's: Increment ', inc, '/', loadCases(currentLoadCase)%incs,&
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
'-', stepFraction, '/', subStepFactor**cutBackLevel,&
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
' of load case ', currentLoadCase,'/',size(loadCases)
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
flush(6)
added some warning and explicit size of arrays
2013-01-10 19:03:43 +05:30
write(incInfo,'(a,'//IO_intOut(totalIncsCounter)//',a,'//IO_intOut(sum(loadCases%incs))//&
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
',a,'//IO_intOut(stepFraction)//',a,'//IO_intOut(subStepFactor**cutBackLevel)//')') &
added some warning and explicit size of arrays
2013-01-10 19:03:43 +05:30
'Increment ',totalIncsCounter,'/',sum(loadCases%incs),&
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
'-',stepFraction, '/', subStepFactor**cutBackLevel
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
select case(spectral_solver)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
!--------------------------------------------------------------------------------------------------
! calculate solution
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
case (DAMASK_spectral_SolverBasic_label)
solres = basic_solution (&
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
incInfo,guess,timeinc,timeIncOld,remainingLoadCaseTime, &
fixed wrong temperature when using spectral solver
2013-01-09 23:38:08 +05:30
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
temperature_bc = loadCases(currentLoadCase)%temperature, &
rotation_BC = loadCases(currentLoadCase)%rotation)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
case (DAMASK_spectral_SolverBasicPETSC_label)
solres = BasicPETSC_solution (&
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
incInfo,guess,timeinc,timeIncOld,remainingLoadCaseTime, &
fixed wrong temperature when using spectral solver
2013-01-09 23:38:08 +05:30
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
temperature_bc = loadCases(currentLoadCase)%temperature, &
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms. dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!). not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 21:55:37 +05:30
rotation_BC = loadCases(currentLoadCase)%rotation, &
rolled back phase field changes
2013-12-18 14:39:32 +05:30
density = loadCases(currentLoadCase)%density)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
case (DAMASK_spectral_SolverAL_label)
solres = AL_solution (&
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
incInfo,guess,timeinc,timeIncOld,remainingLoadCaseTime, &
fixed wrong temperature when using spectral solver
2013-01-09 23:38:08 +05:30
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
temperature_bc = loadCases(currentLoadCase)%temperature, &
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms. dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!). not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 21:55:37 +05:30
rotation_BC = loadCases(currentLoadCase)%rotation, &
density = loadCases(currentLoadCase)%density)
added Polarisation scheme as solver
2013-07-24 18:36:16 +05:30
case (DAMASK_spectral_SolverPolarisation_label)
solres = Polarisation_solution (&
incInfo,guess,timeinc,timeIncOld,remainingLoadCaseTime, &
P_BC = loadCases(currentLoadCase)%P, &
F_BC = loadCases(currentLoadCase)%deformation, &
temperature_bc = loadCases(currentLoadCase)%temperature, &
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms. dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!). not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 21:55:37 +05:30
rotation_BC = loadCases(currentLoadCase)%rotation, &
density = loadCases(currentLoadCase)%density)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
end select
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
!--------------------------------------------------------------------------------------------------
! check solution
corrected statistics reporting in case of cutback
2013-04-26 14:43:39 +05:30
cutBack = .False.
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
if(solres%termIll .or. .not. solres%converged) then ! no solution found
if (cutBackLevel < maxCutBack) then ! do cut back
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
write(6,'(/,a)') ' cut back detected'
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
cutBack = .True.
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
stepFraction = (stepFraction - 1_pInt) * subStepFactor ! adjust to new denominator
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
cutBackLevel = cutBackLevel + 1_pInt
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
time = time - timeinc ! rewind time
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading. Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 15:14:23 +05:30
timeIncOld = timeinc
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
timeinc = timeinc/2.0_pReal
changed handling of non-converged BVP solution: By default, exit if no solution is found to prevent "user errors" Removed ambiguous "regridMode" keyword which was used to trigger this behavior in favor of "continueCalculation", Set this to 1 to get back old behavior, e.g. report non-converged increments.
2014-03-31 15:34:11 +05:30
elseif (solres%termIll) then ! material point model cannot find a solution, exit in any casy
call IO_warning(850_pInt)
call quit(-1_pInt*(lastRestartWritten+1_pInt)) ! quit and provide information about last restart inc written (e.g. for regridding)
elseif (continueCalculation == 1_pInt) then
guess = .true. ! accept non converged BVP solution
else ! default behavior, exit if spectral solver does not converge
call IO_warning(850_pInt)
call quit(-1_pInt*(lastRestartWritten+1_pInt)) ! quit and provide information about last restart inc written (e.g. for regridding) ! continue from non-converged solution and start guessing after accepted (sub)inc
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
endif
else
replaced real 'guessmode' with logical 'guess', caused trouble on gfortran 6.3 with ultra optimization
2012-11-09 01:02:00 +05:30
guess = .true. ! start guessing after first converged (sub)inc
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
endif
simplified logic for single run, now not checking for extra arguments but simply checking if start and end elem/IP are the same
2013-04-26 18:53:36 +05:30
if (.not. cutBack) &
write(statUnit,*) totalIncsCounter, time, cutBackLevel, &
solres%converged, solres%iterationsNeeded ! write statistics about accepted solution
small changes regarding restart and HDF5
2014-03-26 14:11:45 +05:30
flush(statUnit)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
enddo subIncLooping
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
cutBackLevel = max(0_pInt, cutBackLevel - 1_pInt) ! try half number of subincs next inc
if(solres%converged) then ! report converged inc
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
convergedCounter = convergedCounter + 1_pInt
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
write(6,'(/,a,'//IO_intOut(totalIncsCounter)//',a)') &
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
' increment ', totalIncsCounter, ' converged'
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
else
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
write(6,'(/,a,'//IO_intOut(totalIncsCounter)//',a)') & ! report non-converged inc
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
' increment ', totalIncsCounter, ' NOT converged'
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
notConvergedCounter = notConvergedCounter + 1_pInt
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
endif; flush(6)
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters. spectral solver got cut back facilities + improved output to screen
2012-10-02 20:56:56 +05:30
if (mod(inc,loadCases(currentLoadCase)%outputFrequency) == 0_pInt) then ! at output frequency
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
write(6,'(1/,a)') ' ... writing results to file ......................................'
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
write(resUnit) materialpoint_results ! write result to file
small changes regarding restart and HDF5
2014-03-26 14:11:45 +05:30
flush(resUnit)
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
endif
made doxygen working for MSC.Marc again, small changes on the todo-statements for doxygen, they don't work on single lines of codes but only on module/variables
2013-06-11 22:05:04 +05:30
if( loadCases(currentLoadCase)%restartFrequency > 0_pInt .and. & ! at frequency of writing restart information set restart parameter for FEsolving
mod(inc,loadCases(currentLoadCase)%restartFrequency) == 0_pInt) then ! ToDo first call to CPFEM_general will write?
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
restartWrite = .true.
introduced option for regridding to numerics.f90, working as follows: first cut back is tried if material point model oder BVP solver does not converge. If no regridding is enabled after max cut back, in case of non-converged material point the simulation stops and in case of non-converged BVP solver it continues. set regridMode to 2 to enable regridding if BVP solver OR materialPoint model do not converge, set regridMode to 1 to enable regridding if materialPoint model do not converge, non-converged BVP solver will be ignored as in the standard case. For regridding, the load case need to have a restart freq set. enabled restarting for Basic PETSc variant
2012-12-14 20:48:04 +05:30
lastRestartWritten = inc
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
endif
doxygen comments for homogenization.f90, unified naming ip->i, el->e
2013-01-29 15:58:01 +05:30
else forwarding
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
time = time + timeinc
replaced real 'guessmode' with logical 'guess', caused trouble on gfortran 6.3 with ultra optimization
2012-11-09 01:02:00 +05:30
guess = .true.
doxygen comments for homogenization.f90, unified naming ip->i, el->e
2013-01-29 15:58:01 +05:30
endif forwarding
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
next round in modularizzation of solvers
2012-07-24 22:37:10 +05:30
enddo incLooping
enddo loadCaseLooping
new version of modular solver structure
2012-08-03 14:55:48 +05:30
changed naming convention (keeping backward compatibility in numerics.config): myspectralsolver --> spectral_solver (called "spectralsolver") myfilter --> spectral_filter (called "spectralfilter")
2014-09-04 01:29:47 +05:30
select case (spectral_solver)
next round in modularization
2012-07-25 19:31:39 +05:30
case (DAMASK_spectral_SolverBasic_label)
call basic_destroy()
added SolverBasicPETSC module to the new FFT solver (solves the basic scheme using PETSc when mysolver = basicpetsc in numerics.config) some bug fixes
2012-08-14 22:28:23 +05:30
case (DAMASK_spectral_SolverBasicPETSC_label)
call BasicPETSC_destroy()
updates for AL version of spectral solver
2012-08-06 14:23:12 +05:30
case (DAMASK_spectral_SolverAL_label)
call AL_destroy()
added Polarisation scheme as solver
2013-07-24 18:36:16 +05:30
case (DAMASK_spectral_SolverPolarisation_label)
call Polarisation_destroy()
new version of modular solver structure
2012-08-03 14:55:48 +05:30
end select
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
!--------------------------------------------------------------------------------------------------
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
! report summary of whole calculation
fixed some buggy last minute changes in the last commit
2013-08-08 14:43:29 +05:30
write(6,'(/,a)') ' ###########################################################################'
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-10 03:49:32 +05:30
write(6,'(1x,i6.6,a,i6.6,a,f5.1,a)') convergedCounter, ' out of ', &
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
notConvergedCounter + convergedCounter, ' (', &
real(convergedCounter, pReal)/&
real(notConvergedCounter + convergedCounter,pReal)*100.0_pReal, &
' %) increments converged!'
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 14:14:21 +05:30
close(resUnit)
changed handling of non-converged BVP solution: By default, exit if no solution is found to prevent "user errors" Removed ambiguous "regridMode" keyword which was used to trigger this behavior in favor of "continueCalculation", Set this to 1 to get back old behavior, e.g. report non-converged increments.
2014-03-31 15:34:11 +05:30
close(statUnit)
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6 added comments
2012-12-15 21:51:10 +05:30
if (notConvergedCounter > 0_pInt) call quit(3_pInt) ! error if some are not converged
call quit(0_pInt) ! no complains ;)
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
next round in modularizzation of solvers
2012-07-24 22:37:10 +05:30
end program DAMASK_spectral_Driver
intermediate, not working state of new solver structure
2012-07-20 21:03:13 +05:30
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
!--------------------------------------------------------------------------------------------------
!> @author Martin Diehl, Max-Planck-Institut für Eisenforschung GmbH
!> @brief quit subroutine to mimic behavior of FEM solvers
!> @details exits the Spectral solver and reports time and duration. Exit code 0 signals
!> everything went fine. Exit code 1 signals an error, message according to IO_error. Exit code
!> 2 signals request for regridding, increment of last saved restart information is written to
!> stderr. Exit code 3 signals no severe problems, but some increments did not converge
!--------------------------------------------------------------------------------------------------
relocated spectral_quit() into driver code. added (dummy) core_quit() for Python damask module.
2013-02-19 20:24:34 +05:30
subroutine quit(stop_id)
use prec, only: &
pInt
implicit none
integer(pInt), intent(in) :: stop_id
integer, dimension(8) :: dateAndTime ! type default integer
call date_and_time(values = dateAndTime)
write(6,'(/,a)') 'DAMASK terminated on:'
write(6,'(a,2(i2.2,a),i4.4)') 'Date: ',dateAndTime(3),'/',&
dateAndTime(2),'/',&
dateAndTime(1)
write(6,'(a,2(i2.2,a),i2.2)') 'Time: ',dateAndTime(5),':',&
dateAndTime(6),':',&
dateAndTime(7)
if (stop_id == 0_pInt) stop 0 ! normal termination
if (stop_id < 0_pInt) then ! trigger regridding
changed handling of non-converged BVP solution: By default, exit if no solution is found to prevent "user errors" Removed ambiguous "regridMode" keyword which was used to trigger this behavior in favor of "continueCalculation", Set this to 1 to get back old behavior, e.g. report non-converged increments.
2014-03-31 15:34:11 +05:30
write(0,'(a,i6)') 'restart information available at ', stop_id*(-1_pInt)
relocated spectral_quit() into driver code. added (dummy) core_quit() for Python damask module.
2013-02-19 20:24:34 +05:30
stop 2
endif
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
if (stop_id == 3_pInt) stop 3 ! not all incs converged
relocated spectral_quit() into driver code. added (dummy) core_quit() for Python damask module.
2013-02-19 20:24:34 +05:30
stop 1 ! error (message from IO_error)
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +05:30
end subroutine quit