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

628 lines
40 KiB
Fortran
Raw Normal View History

updated copyright header
2013-03-22 23:05:05 +05:30
! Copyright 2011-13 Max-Planck-Institut für Eisenforschung GmbH
!
! This file is part of DAMASK,
! the Düsseldorf Advanced MAterial Simulation Kit.
!
! DAMASK is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation, either version 3 of the License, or
! (at your option) any later version.
!
! DAMASK is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with DAMASK. If not, see <http://www.gnu.org/licenses/>.
!
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, &
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
mySpectralSolver, &
regridMode
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
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
#ifdef PETSc
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
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
#endif
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) :: &
some more documentation, mainly on the exit codes
2013-02-23 05:54:30 +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
totalIncsCounter = 0_pInt, & !< total No. of increments
convergedCounter = 0_pInt, & !< No. of converged increments
notConvergedCounter = 0_pInt, & !< No. of non-converged increments
resUnit = 0_pInt, & !< file unit for results writing
statUnit = 0_pInt, & !< file unit for statistics output
lastRestartWritten = 0_pInt !< total increment No. 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
new version of modular solver structure
2012-08-03 14:55:48 +05:30
select case (myspectralsolver)
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)
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
#ifdef PETSc
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
#endif
case default
call IO_error(error_ID = 891, ext_msg = trim(myspectralsolver))
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')
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) 'load', trim(loadCaseFile) ! ... and write header
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
write(resUnit) 'workingdir', trim(getSolverWorkingDirectoryName())
write(resUnit) 'geometry', trim(geometryFile)
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(resUnit) 'resolution', grid
write(resUnit) 'dimension', geomSize
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
write(resUnit) 'materialpoint_sizeResults', materialpoint_sizeResults
write(resUnit) 'loadcases', size(loadCases)
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
write(resUnit) 'frequencies', loadCases%outputfrequency ! one entry per currentLoadCase
write(resUnit) 'times', loadCases%time ! one entry per currentLoadCase
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
write(resUnit) 'logscales', loadCases%logscale
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 17:01:40 +05:30
write(resUnit) 'increments', loadCases%incs ! one entry per currentLoadCase
write(resUnit) 'startingIncrement', restartInc - 1_pInt ! start with writing out the previous inc
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
if (loadCases(currentLoadCase)%followFormerTrajectory) then
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.
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
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 = .false. ! 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
endif
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
removed old standalone AL solver and introduced new structure for solver zoo
2012-07-19 22:54:56 +05:30
totalIncsCounter = totalIncsCounter + 1_pInt
!--------------------------------------------------------------------------------------------------
! 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
select case(myspectralsolver)
case (DAMASK_spectral_SolverBasic_label)
#ifdef PETSc
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)
#endif
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
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
select case(myspectralsolver)
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)
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
#ifdef PETSc
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)
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
#endif
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
elseif (solres%termIll) then ! material point model cannot find a solution
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
if(regridMode > 0_pInt) call quit(-1_pInt*(lastRestartWritten+1_pInt)) ! regrid requested (mode 1 or 2)
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
call IO_error(850_pInt) ! no regrid (give up)
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
else
fixed bug in cutback scheme
2013-01-08 03:12:00 +05:30
if(regridMode == 2_pInt) call quit(-1_pInt*(lastRestartWritten+1_pInt)) ! regrid also if BVP solver do not converge
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
guess = .true. ! 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
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
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
next round in modularization
2012-07-25 19:31:39 +05:30
select case (myspectralsolver)
case (DAMASK_spectral_SolverBasic_label)
call basic_destroy()
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
#ifdef PETSc
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 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
#endif
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)
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
write(0,'(a,i6)') 'restart at ', stop_id*(-1_pInt)
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