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

438 lines
26 KiB
Fortran
Raw Normal View History

removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief Managing of parameters related to numerics
!--------------------------------------------------------------------------------------------------
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
module numerics
one implicit none is enough
2019-05-15 02:42:32 +05:30
use prec
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
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
implicit none
private
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pInt not needed
2019-06-11 18:09:51 +05:30
integer, protected, public :: &
iJacoStiffness = 1, & !< frequency of stiffness update
nMPstate = 10, & !< materialpoint state loop limit
randomSeed = 0, & !< fixed seeding for pseudo-random number generator, Default 0: use random seed
worldrank = 0, & !< MPI worldrank (/=0 for MPI simulations only)
worldsize = 1, & !< MPI worldsize (/=1 for MPI simulations only)
numerics_integrator = 1 !< method used for state integration Default 1: fix-point iteration
set OMP related variables to 4 byte integer as expected when compiling marc
2016-12-23 18:50:29 +05:30
integer(4), protected, public :: &
fixed type-casting problems (mainly with Lib_VTK_IO)
2013-09-26 22:51:46 +05:30
DAMASK_NumThreadsInt = 0 !< value stored in environment variable DAMASK_NUM_THREADS, set to zero if no OpenMP directive
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
real(pReal), protected, public :: &
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
defgradTolerance = 1.0e-7_pReal, & !< deviation of deformation gradient that is still allowed (used by CPFEM to determine outdated ffn1)
subStepMinHomog = 1.0e-3_pReal, & !< minimum (relative) size of sub-step allowed during cutback in homogenization
subStepSizeHomog = 0.25_pReal, & !< size of first substep when cutback in homogenization
stepIncreaseHomog = 1.5_pReal, & !< increase of next substep size when previous substep converged in homogenization
numerics_unitlength = 1.0_pReal, & !< determines the physical length of one computational length unit
absTol_RGC = 1.0e+4_pReal, & !< absolute tolerance of RGC residuum
relTol_RGC = 1.0e-3_pReal, & !< relative tolerance of RGC residuum
absMax_RGC = 1.0e+10_pReal, & !< absolute maximum of RGC residuum
relMax_RGC = 1.0e+2_pReal, & !< relative maximum of RGC residuum
pPert_RGC = 1.0e-7_pReal, & !< perturbation for computing RGC penalty tangent
xSmoo_RGC = 1.0e-5_pReal, & !< RGC penalty smoothing parameter (hyperbolic tangent)
viscPower_RGC = 1.0e+0_pReal, & !< power (sensitivity rate) of numerical viscosity in RGC scheme, Default 1.0e0: Newton viscosity (linear model)
viscModus_RGC = 0.0e+0_pReal, & !< stress modulus of RGC numerical viscosity, Default 0.0e0: No viscosity is applied
refRelaxRate_RGC = 1.0e-3_pReal, & !< reference relaxation rate in RGC viscosity
maxdRelax_RGC = 1.0e+0_pReal, & !< threshold of maximum relaxation vector increment (if exceed this then cutback)
maxVolDiscr_RGC = 1.0e-5_pReal, & !< threshold of maximum volume discrepancy allowed
volDiscrMod_RGC = 1.0e+12_pReal, & !< stiffness of RGC volume discrepancy (zero = without volume discrepancy constraint)
added mesh dependent characteristic length as a numerical parameter needed for gradient and phase field problems on unstructured meshes
2014-06-25 04:48:07 +05:30
volDiscrPow_RGC = 5.0_pReal, & !< powerlaw penalty for volume discrepancy
residualStiffness is not just an FEM parameter
2014-12-02 22:47:35 +05:30
charLength = 1.0_pReal, & !< characteristic length scale for gradient problems
residualStiffness = 1.0e-6_pReal !< non-zero residual damage
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
logical, protected, public :: &
unused variables
2019-04-11 10:16:31 +05:30
usePingPong = .true.
introduced preprocessor identifiers Marc, Abaqus, and Spectral to enable conditional compilation. This allows deleted dummy functions that are used by one solver only. Mainly affected modules are IO and mesh. Most of the changes in mesh result from reordering the functions when grouping them depending on their solver. Further advantage is that FE solver do not need FFTW and kdtree2 anymore. The include files for these two libraries moved to DAMASKROO/lib now as I figured out how to use a include path in the Makefile. Put all the files I got when testing compilation with abaqus in a folder which to become the abaqus compilation test.
2012-06-15 21:40:21 +05:30
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
!--------------------------------------------------------------------------------------------------
! field parameters:
real(pReal), protected, public :: &
err_struct_tolAbs = 1.0e-10_pReal, & !< absolute tolerance for mechanical equilibrium
err_struct_tolRel = 1.0e-4_pReal, & !< relative tolerance for mechanical equilibrium
err_thermal_tolAbs = 1.0e-2_pReal, & !< absolute tolerance for thermal equilibrium
err_thermal_tolRel = 1.0e-6_pReal, & !< relative tolerance for thermal equilibrium
err_damage_tolAbs = 1.0e-2_pReal, & !< absolute tolerance for damage evolution
cleaning
2019-03-07 11:39:00 +05:30
err_damage_tolRel = 1.0e-6_pReal !< relative tolerance for damage evolution
pInt not needed
2019-06-11 18:09:51 +05:30
integer, protected, public :: &
itmax = 250, & !< maximum number of iterations
itmin = 1, & !< minimum number of iterations
stagItMax = 10, & !< max number of field level staggered iterations
maxCutBack = 3 !< max number of cut backs
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
!--------------------------------------------------------------------------------------------------
! spectral parameters:
renaming "spectral" to "grid"
2019-03-23 13:57:58 +05:30
#ifdef Grid
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
real(pReal), protected, public :: &
absTol for divergence was way too tight
2018-08-21 02:04:43 +05:30
err_div_tolAbs = 1.0e-4_pReal, & !< absolute tolerance for equilibrium
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 22:50:05 +05:30
err_div_tolRel = 5.0e-4_pReal, & !< relative tolerance for equilibrium
err_curl_tolAbs = 1.0e-10_pReal, & !< absolute tolerance for compatibility
err_curl_tolRel = 5.0e-4_pReal, & !< relative tolerance for compatibility
err_stress_tolAbs = 1.0e3_pReal, & !< absolute tolerance for fullfillment of stress BC
err_stress_tolRel = 0.01_pReal, & !< relative tolerance for fullfillment of stress BC
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
rotation_tol = 1.0e-12_pReal, & !< tolerance of rotation specified in loadcase, Default 1.0e-12: first guess
polarAlpha = 1.0_pReal, & !< polarization scheme parameter 0.0 < alpha < 2.0. alpha = 1.0 ==> AL scheme, alpha = 2.0 ==> accelerated scheme
polarBeta = 1.0_pReal !< polarization scheme parameter 0.0 < beta < 2.0. beta = 1.0 ==> AL scheme, beta = 2.0 ==> accelerated scheme
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
character(len=1024), protected, public :: &
split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options
2015-06-09 18:58:50 +05:30
petsc_defaultOptions = '-mech_snes_type ngmres &
added spectral thermal and damage solvers
2015-06-03 23:00:31 +05:30
&-damage_snes_type ngmres &
split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options
2015-06-09 18:58:50 +05:30
&-thermal_snes_type ngmres ', &
petsc_options = ''
allowed arbitrary grid (formerly know as resolution), i.e. any integer>1 in each dimension is possible. removed square root correction for divergence (doesn't make sense) fixed FFT debug of spectral solver
2013-04-10 15:49:16 +05:30
logical, protected, public :: &
better define variables close to their usage
2019-03-26 12:06:55 +05:30
continueCalculation = .false. !< false:exit if BVP solver does not converge, true: continue calculation despite BVP solver not converging
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
mpie_spectral and mpie_interface: changed workingdir from pathToLoadFile to pathToGeomFile mpie_spectral and numerics: added switch to prevent pre calculation of gamma_hat. slower, but saves memory 3Dvisualize: started to add support for gmsh (not fully working yet) reconstruct: new version of f2py/Fortran subroutines for output of results from spectral method
2011-02-07 20:05:42 +05:30
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
!--------------------------------------------------------------------------------------------------
! FEM parameters:
#ifdef FEM
pInt not needed
2019-06-11 18:09:51 +05:30
integer, protected, public :: &
integrationOrder = 2, & !< order of quadrature rule required
structOrder = 2 !< order of displacement shape functions
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
logical, protected, public :: &
BBarStabilisation = .false.
increased size of petsc_options string
2014-10-21 21:55:53 +05:30
character(len=4096), protected, public :: &
split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options
2015-06-09 18:58:50 +05:30
petsc_defaultOptions = '-mech_snes_type newtonls &
numerics with respect to vacancy diffusion.
2014-12-08 21:07:44 +05:30
&-mech_snes_linesearch_type cp &
&-mech_snes_ksp_ew &
&-mech_snes_ksp_ew_rtol0 0.01 &
&-mech_snes_ksp_ew_rtolmax 0.01 &
&-mech_ksp_type fgmres &
&-mech_ksp_max_it 25 &
&-mech_pc_type ml &
&-mech_mg_levels_ksp_type chebyshev &
&-mech_mg_levels_pc_type sor &
cleaning, only having mech at the moment
2018-09-15 18:44:48 +05:30
&-mech_pc_ml_nullspace user ',&
split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options
2015-06-09 18:58:50 +05:30
petsc_options = ''
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
#endif
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
public :: numerics_init
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
contains
introduced a flexibility in cut-backing scheme in homogenization.f90 and in crystallite.f90: (1) subStepSizeHomog and subStepSizeCryst := size of substep when cut-back is applied (initially was hard-coded). (2) stepIncreaseHomog and stepIncreaseCryst := step increase when calculation for substep converge (was also hardcoded). introduced a possibility to choose different finite difference scheme, i.e., forward-, backward- and central-difference, for computing grain numerical tangent. note that central-difference scheme will slow down the computation significantly. please use it only if necessary. parameters to set these new features have been included in numerics.f90 and numerics.config, whereas corresponding error messages have been introduced in the IO.f90
2009-11-10 19:06:27 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief reads in parameters from numerics.config and sets openMP related parameters. Also does
! a sanity check
!--------------------------------------------------------------------------------------------------
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
subroutine numerics_init
made PETSc silent by removing -snes_view from default options. introduced PETSc option for debugging that introduces some debugging options into the petsc options and move PETSc initialization from numerics to DAMASK_spectral_utilities
2012-12-15 23:37:49 +05:30
use IO, only: &
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
IO_read_ASCII, &
made PETSc silent by removing -snes_view from default options. introduced PETSc option for debugging that introduces some debugging options into the petsc options and move PETSc initialization from numerics to DAMASK_spectral_utilities
2012-12-15 23:37:49 +05:30
IO_error, &
IO_isBlank, &
IO_stringPos, &
IO_stringValue, &
IO_lc, &
IO_floatValue, &
IO_intValue, &
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
IO_warning
employing use statements (almost no external functions anymore) adopted function calls to newer PETSc versions
2018-05-17 15:34:21 +05:30
#ifdef PETSc
#include <petsc/finclude/petscsys.h>
use petscsys
#endif
Marc now works also with the module reason, why it did NOT work earlier still not clear
2019-02-03 18:00:04 +05:30
!$ use OMP_LIB, only: omp_set_num_threads
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
!$ integer :: gotDAMASK_NUM_THREADS = 1
pInt not needed
2019-06-11 18:09:51 +05:30
integer :: i,j, ierr
integer, allocatable, dimension(:) :: chunkPos
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
character(len=pStringLen), dimension(:), allocatable :: fileContent
character(len=pStringLen) :: &
set tolerances to useful values and renamed them. allowed "keyword = value" in numerics.config (simply removing "=" from line)
2013-07-31 20:33:38 +05:30
tag ,&
line
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
logical :: fexist
fixed usage of OpenMP function library
2012-06-12 15:14:05 +05:30
!$ character(len=6) DAMASK_NumThreadsString ! environment variable DAMASK_NUM_THREADS
mpie_cpfem_marc.f90 >>> changing the compilation sequence of modules: numerics.f90 and IO.f90 move up, while math.f90 down. IO.f90 >>> adding "inRad" instead of using inRad from math.f90 module, since math.f90 module is now below IO.f90. numerics.f90 >>> now can read the parameter "fixed_seed" from numerics.config. this parameter is used to set a pseudo-random number generator/fixed seeding. default value is 0 (zero), which will give fully random number. math.f90 >> adding a possibility to fix the seeding (i.e., pseudo random number) numerics.config >>> add a new parameter: "fixed_seed"
2009-08-27 21:00:40 +05:30
simplified multi processor (MPI) reporting
2014-10-10 18:38:34 +05:30
#ifdef PETSc
only output from the root processor for parallel runs
2014-10-10 01:53:06 +05:30
call MPI_Comm_rank(PETSC_COMM_WORLD,worldrank,ierr);CHKERRQ(ierr)
call MPI_Comm_size(PETSC_COMM_WORLD,worldsize,ierr);CHKERRQ(ierr)
simplified multi processor (MPI) reporting
2014-10-10 18:38:34 +05:30
#endif
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(/,a)') ' <<<+- numerics init -+>>>'
simplified multi processor (MPI) reporting
2014-10-10 18:38:34 +05:30
merge declaration and initialization with standard values corrected handling of $DAMASK_NUM_THREADS
2012-01-31 01:46:19 +05:30
!$ call GET_ENVIRONMENT_VARIABLE(NAME='DAMASK_NUM_THREADS',VALUE=DAMASK_NumThreadsString,STATUS=gotDAMASK_NUM_THREADS) ! get environment variable DAMASK_NUM_THREADS...
improved error handling when getting number of threads
2014-02-28 16:00:07 +05:30
!$ if(gotDAMASK_NUM_THREADS /= 0) then ! could not get number of threads, set it to 1
pInt not needed
2019-06-11 18:09:51 +05:30
!$ call IO_warning(35,ext_msg='BEGIN:'//DAMASK_NumThreadsString//':END')
set OMP related variables to 4 byte integer as expected when compiling marc
2016-12-23 18:50:29 +05:30
!$ DAMASK_NumThreadsInt = 1_4
improved error handling when getting number of threads
2014-02-28 16:00:07 +05:30
!$ else
!$ read(DAMASK_NumThreadsString,'(i6)') DAMASK_NumThreadsInt ! read as integer
set OMP related variables to 4 byte integer as expected when compiling marc
2016-12-23 18:50:29 +05:30
!$ if (DAMASK_NumThreadsInt < 1_4) DAMASK_NumThreadsInt = 1_4 ! in case of string conversion fails, set it to one
improved error handling when getting number of threads
2014-02-28 16:00:07 +05:30
!$ endif
!$ call omp_set_num_threads(DAMASK_NumThreadsInt) ! set number of threads for parallel execution
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
inquire(file='numerics.config', exist=fexist)
fileExists: if (fexist) then
compiler complained
2017-04-25 16:04:14 +05:30
write(6,'(a,/)') ' using values from config file'
flush(6)
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
fileContent = IO_read_ASCII('numerics.config')
do j=1, size(fileContent)
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! read variables from config file and overwrite default parameters if keyword is present
use central functionality IO function reads numerics.config and debug.config
2019-03-09 04:37:57 +05:30
line = fileContent(j)
set tolerances to useful values and renamed them. allowed "keyword = value" in numerics.config (simply removing "=" from line)
2013-07-31 20:33:38 +05:30
do i=1,len(line)
if(line(i:i) == '=') line(i:i) = ' ' ! also allow keyword = value version
enddo
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name. functions in the core module are not used, including kdTree
2015-08-28 13:08:48 +05:30
chunkPos = IO_stringPos(line)
pInt not needed
2019-06-11 18:09:51 +05:30
tag = IO_lc(IO_stringValue(line,chunkPos,1)) ! extract key
behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name. functions in the core module are not used, including kdTree
2015-08-28 13:08:48 +05:30
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
select case(tag)
case ('defgradtolerance')
pInt not needed
2019-06-11 18:09:51 +05:30
defgradTolerance = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('ijacostiffness')
pInt not needed
2019-06-11 18:09:51 +05:30
iJacoStiffness = IO_intValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('nmpstate')
pInt not needed
2019-06-11 18:09:51 +05:30
nMPstate = IO_intValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('substepminhomog')
pInt not needed
2019-06-11 18:09:51 +05:30
subStepMinHomog = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('substepsizehomog')
pInt not needed
2019-06-11 18:09:51 +05:30
subStepSizeHomog = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('stepincreasehomog')
pInt not needed
2019-06-11 18:09:51 +05:30
stepIncreaseHomog = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('integrator')
pInt not needed
2019-06-11 18:09:51 +05:30
numerics_integrator = IO_intValue(line,chunkPos,2)
introduced plain mode (no ping pong) again and added test for MSC.MArc 2012 as prove that its working
2013-05-17 23:22:46 +05:30
case ('usepingpong')
pInt not needed
2019-06-11 18:09:51 +05:30
usepingpong = IO_intValue(line,chunkPos,2) > 0
forgot to check in; belongs to rev 1762
2012-09-24 21:52:25 +05:30
case ('unitlength')
pInt not needed
2019-06-11 18:09:51 +05:30
numerics_unitlength = IO_floatValue(line,chunkPos,2)
List of changes/modifications: * IO.f90 :: Adding error messages for RGC homogenization * crystallite.f90 :: Modifying convergent criteria in crystalline_updateState and crystalline_updateTemperature * material.f90 :: Adding IO_lc in homogenization_type(Name) * debug.f90 :: Adding debbugging statement and counter for material point loop * homogenization.f90 :: Adding homogenization_RGC blocks * homogenization_isostrain.f90 :: Modifying argument of homogenization_isostrain_stateInit * mpie_cpfem_marc.f90 :: Adding homogenization_RGC include * numerics.f90 :: Adding numerical parameters for RGC scheme * math.f90 :: Changing function name: math_permut to math_civita
2009-07-31 17:32:20 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! RGC parameters
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('atol_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
absTol_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('rtol_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
relTol_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('amax_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
absMax_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('rmax_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
relMax_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('perturbpenalty_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
pPert_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('relevantmismatch_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
xSmoo_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('viscositypower_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
viscPower_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('viscositymodulus_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
viscModus_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('refrelaxationrate_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
refRelaxRate_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('maxrelaxation_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
maxdRelax_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('maxvoldiscrepancy_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
maxVolDiscr_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('voldiscrepancymod_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
volDiscrMod_RGC = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('discrepancypower_rgc')
pInt not needed
2019-06-11 18:09:51 +05:30
volDiscrPow_RGC = IO_floatValue(line,chunkPos,2)
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
! random seeding parameter
changed fixed_seed to random_seed for clarity
2017-11-08 01:26:28 +05:30
case ('random_seed','fixed_seed')
pInt not needed
2019-06-11 18:09:51 +05:30
randomSeed = IO_intValue(line,chunkPos,2)
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
added mesh dependent characteristic length as a numerical parameter needed for gradient and phase field problems on unstructured meshes
2014-06-25 04:48:07 +05:30
!--------------------------------------------------------------------------------------------------
! gradient parameter
case ('charlength')
pInt not needed
2019-06-11 18:09:51 +05:30
charLength = IO_floatValue(line,chunkPos,2)
residualStiffness is not just an FEM parameter
2014-12-02 22:47:35 +05:30
case ('residualstiffness')
pInt not needed
2019-06-11 18:09:51 +05:30
residualStiffness = IO_floatValue(line,chunkPos,2)
added mesh dependent characteristic length as a numerical parameter needed for gradient and phase field problems on unstructured meshes
2014-06-25 04:48:07 +05:30
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
!--------------------------------------------------------------------------------------------------
! field parameters
case ('err_struct_tolabs')
pInt not needed
2019-06-11 18:09:51 +05:30
err_struct_tolAbs = IO_floatValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('err_struct_tolrel')
pInt not needed
2019-06-11 18:09:51 +05:30
err_struct_tolRel = IO_floatValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('err_thermal_tolabs')
pInt not needed
2019-06-11 18:09:51 +05:30
err_thermal_tolabs = IO_floatValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('err_thermal_tolrel')
pInt not needed
2019-06-11 18:09:51 +05:30
err_thermal_tolrel = IO_floatValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('err_damage_tolabs')
pInt not needed
2019-06-11 18:09:51 +05:30
err_damage_tolabs = IO_floatValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('err_damage_tolrel')
pInt not needed
2019-06-11 18:09:51 +05:30
err_damage_tolrel = IO_floatValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('itmax')
pInt not needed
2019-06-11 18:09:51 +05:30
itmax = IO_intValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('itmin')
pInt not needed
2019-06-11 18:09:51 +05:30
itmin = IO_intValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('maxcutback')
pInt not needed
2019-06-11 18:09:51 +05:30
maxCutBack = IO_intValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
case ('maxstaggerediter')
pInt not needed
2019-06-11 18:09:51 +05:30
stagItMax = IO_intValue(line,chunkPos,2)
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! spectral parameters
renaming "spectral" to "grid"
2019-03-23 13:57:58 +05:30
#ifdef Grid
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 22:50:05 +05:30
case ('err_div_tolabs')
pInt not needed
2019-06-11 18:09:51 +05:30
err_div_tolAbs = IO_floatValue(line,chunkPos,2)
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 22:50:05 +05:30
case ('err_div_tolrel')
pInt not needed
2019-06-11 18:09:51 +05:30
err_div_tolRel = IO_floatValue(line,chunkPos,2)
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
case ('err_stress_tolrel')
pInt not needed
2019-06-11 18:09:51 +05:30
err_stress_tolrel = IO_floatValue(line,chunkPos,2)
added missing ":" in mesh.f90, introduced absolute stress tolerance for spectral solver in numerics.config/numerics.f90
2012-04-11 18:27:25 +05:30
case ('err_stress_tolabs')
pInt not needed
2019-06-11 18:09:51 +05:30
err_stress_tolabs = IO_floatValue(line,chunkPos,2)
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
case ('continuecalculation')
pInt not needed
2019-06-11 18:09:51 +05:30
continueCalculation = IO_intValue(line,chunkPos,2) > 0
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 ('petsc_options')
behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name. functions in the core module are not used, including kdTree
2015-08-28 13:08:48 +05:30
petsc_options = trim(line(chunkPos(4):))
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 22:50:05 +05:30
case ('err_curl_tolabs')
pInt not needed
2019-06-11 18:09:51 +05:30
err_curl_tolAbs = IO_floatValue(line,chunkPos,2)
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 22:50:05 +05:30
case ('err_curl_tolrel')
pInt not needed
2019-06-11 18:09:51 +05:30
err_curl_tolRel = IO_floatValue(line,chunkPos,2)
restructured AL solver to resemble polarization scheme and introduced alpha and beta parameters in numerics to determine what kind of solution you want. alpha = beta = 1 results in the AL scheme. alpha = beta = 2 results in the accelerated scheme. safe values are anything in between. changed the solver to use stress calculations on the campatible field which seems to result in faster convergence.
2013-02-28 23:05:02 +05:30
case ('polaralpha')
pInt not needed
2019-06-11 18:09:51 +05:30
polarAlpha = IO_floatValue(line,chunkPos,2)
restructured AL solver to resemble polarization scheme and introduced alpha and beta parameters in numerics to determine what kind of solution you want. alpha = beta = 1 results in the AL scheme. alpha = beta = 2 results in the accelerated scheme. safe values are anything in between. changed the solver to use stress calculations on the campatible field which seems to result in faster convergence.
2013-02-28 23:05:02 +05:30
case ('polarbeta')
pInt not needed
2019-06-11 18:09:51 +05:30
polarBeta = IO_floatValue(line,chunkPos,2)
introduced preprocessor identifiers Marc, Abaqus, and Spectral to enable conditional compilation. This allows deleted dummy functions that are used by one solver only. Mainly affected modules are IO and mesh. Most of the changes in mesh result from reordering the functions when grouping them depending on their solver. Further advantage is that FE solver do not need FFTW and kdtree2 anymore. The include files for these two libraries moved to DAMASKROO/lib now as I figured out how to use a include path in the Makefile. Put all the files I got when testing compilation with abaqus in a folder which to become the abaqus compilation test.
2012-06-15 21:40:21 +05:30
#endif
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
!--------------------------------------------------------------------------------------------------
! FEM parameters
#ifdef FEM
case ('integrationorder')
pInt not needed
2019-06-11 18:09:51 +05:30
integrationorder = IO_intValue(line,chunkPos,2)
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
case ('structorder')
pInt not needed
2019-06-11 18:09:51 +05:30
structorder = IO_intValue(line,chunkPos,2)
split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options
2015-06-09 18:58:50 +05:30
case ('petsc_options')
behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name. functions in the core module are not used, including kdTree
2015-08-28 13:08:48 +05:30
petsc_options = trim(line(chunkPos(4):))
BBar stabilisation now optional (default is off)
2014-12-08 21:57:23 +05:30
case ('bbarstabilisation')
pInt not needed
2019-06-11 18:09:51 +05:30
BBarStabilisation = IO_intValue(line,chunkPos,2) > 0
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
#endif
foramtting of random_seed output to align with rest
2018-10-09 02:31:27 +05:30
end select
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
enddo
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
using config_numerics instead of global values
2019-03-25 14:06:59 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
else fileExists
write(6,'(a,/)') ' using standard values'
slightly polished makefile and numerics, recompile fftw with icc; still not running with ifort 14.01
2014-02-06 16:11:34 +05:30
flush(6)
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
endif fileExists
Introduced an alternative cutback scheme for the nonlocal model, that allows to keep the results of most of the ips that immediately converged, and only do cutbacks in regions where some ips are in big trouble. Only works for nonlocal model and state integrator 2.
2012-11-28 00:06:55 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
! writing parameters to output
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,es8.1)') ' defgradTolerance: ',defgradTolerance
write(6,'(a24,1x,i8)') ' iJacoStiffness: ',iJacoStiffness
WIP: cleaned no file reading getting rid of a number of obsolete dependencies
2019-02-23 01:07:41 +05:30
write(6,'(a24,1x,i8)') ' integrator: ',numerics_integrator
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,L8)') ' use ping pong scheme: ',usepingpong
write(6,'(a24,1x,es8.1,/)')' unitlength: ',numerics_unitlength
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 15:32:11 +05:30
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,es8.1)') ' subStepMinHomog: ',subStepMinHomog
write(6,'(a24,1x,es8.1)') ' subStepSizeHomog: ',subStepSizeHomog
write(6,'(a24,1x,es8.1)') ' stepIncreaseHomog: ',stepIncreaseHomog
write(6,'(a24,1x,i8,/)') ' nMPstate: ',nMPstate
List of changes/modifications: * IO.f90 :: Adding error messages for RGC homogenization * crystallite.f90 :: Modifying convergent criteria in crystalline_updateState and crystalline_updateTemperature * material.f90 :: Adding IO_lc in homogenization_type(Name) * debug.f90 :: Adding debbugging statement and counter for material point loop * homogenization.f90 :: Adding homogenization_RGC blocks * homogenization_isostrain.f90 :: Modifying argument of homogenization_isostrain_stateInit * mpie_cpfem_marc.f90 :: Adding homogenization_RGC include * numerics.f90 :: Adding numerical parameters for RGC scheme * math.f90 :: Changing function name: math_permut to math_civita
2009-07-31 17:32:20 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! RGC parameters
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,es8.1)') ' aTol_RGC: ',absTol_RGC
write(6,'(a24,1x,es8.1)') ' rTol_RGC: ',relTol_RGC
write(6,'(a24,1x,es8.1)') ' aMax_RGC: ',absMax_RGC
write(6,'(a24,1x,es8.1)') ' rMax_RGC: ',relMax_RGC
write(6,'(a24,1x,es8.1)') ' perturbPenalty_RGC: ',pPert_RGC
write(6,'(a24,1x,es8.1)') ' relevantMismatch_RGC: ',xSmoo_RGC
write(6,'(a24,1x,es8.1)') ' viscosityrate_RGC: ',viscPower_RGC
write(6,'(a24,1x,es8.1)') ' viscositymodulus_RGC: ',viscModus_RGC
write(6,'(a24,1x,es8.1)') ' maxrelaxation_RGC: ',maxdRelax_RGC
write(6,'(a24,1x,es8.1)') ' maxVolDiscrepancy_RGC: ',maxVolDiscr_RGC
write(6,'(a24,1x,es8.1)') ' volDiscrepancyMod_RGC: ',volDiscrMod_RGC
write(6,'(a24,1x,es8.1,/)') ' discrepancyPower_RGC: ',volDiscrPow_RGC
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! Random seeding parameter
foramtting of random_seed output to align with rest
2018-10-09 02:31:27 +05:30
write(6,'(a16,1x,i16,/)') ' random_seed: ',randomSeed
pInt not needed
2019-06-11 18:09:51 +05:30
if (randomSeed <= 0) &
changed fixed_seed to random_seed for clarity
2017-11-08 01:26:28 +05:30
write(6,'(a,/)') ' random seed will be generated!'
nicer FEM parameters
2014-11-03 16:06:07 +05:30
!--------------------------------------------------------------------------------------------------
! gradient parameter
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,es8.1)') ' charLength: ',charLength
write(6,'(a24,1x,es8.1)') ' residualStiffness: ',residualStiffness
nicer FEM parameters
2014-11-03 16:06:07 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! openMP parameter
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
!$ write(6,'(a24,1x,i8,/)') ' number of threads: ',DAMASK_NumThreadsInt
added matrix multiplication 3333x33 to math.f90 added some parameters for spectral method to numerics.f90 (tolerance) changed error message concerning spectral method in IO.f90 corrected calculation of stress BC in mpie_spectral.f90
2010-10-13 21:34:44 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
! field parameters
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,i8)') ' itmax: ',itmax
write(6,'(a24,1x,i8)') ' itmin: ',itmin
write(6,'(a24,1x,i8)') ' maxCutBack: ',maxCutBack
write(6,'(a24,1x,i8)') ' maxStaggeredIter: ',stagItMax
write(6,'(a24,1x,es8.1)') ' err_struct_tolAbs: ',err_struct_tolAbs
write(6,'(a24,1x,es8.1)') ' err_struct_tolRel: ',err_struct_tolRel
write(6,'(a24,1x,es8.1)') ' err_thermal_tolabs: ',err_thermal_tolabs
write(6,'(a24,1x,es8.1)') ' err_thermal_tolrel: ',err_thermal_tolrel
write(6,'(a24,1x,es8.1)') ' err_damage_tolabs: ',err_damage_tolabs
write(6,'(a24,1x,es8.1)') ' err_damage_tolrel: ',err_damage_tolrel
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
!--------------------------------------------------------------------------------------------------
! spectral parameters
renaming "spectral" to "grid"
2019-03-23 13:57:58 +05:30
#ifdef Grid
included changes to correct cutback issue of spectral solver
2018-02-16 20:06:18 +05:30
write(6,'(a24,1x,L8)') ' continueCalculation: ',continueCalculation
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,es8.1)') ' err_stress_tolAbs: ',err_stress_tolAbs
write(6,'(a24,1x,es8.1)') ' err_stress_tolRel: ',err_stress_tolRel
write(6,'(a24,1x,es8.1)') ' err_div_tolAbs: ',err_div_tolAbs
write(6,'(a24,1x,es8.1)') ' err_div_tolRel: ',err_div_tolRel
write(6,'(a24,1x,es8.1)') ' err_curl_tolAbs: ',err_curl_tolAbs
write(6,'(a24,1x,es8.1)') ' err_curl_tolRel: ',err_curl_tolRel
write(6,'(a24,1x,es8.1)') ' polarAlpha: ',polarAlpha
write(6,'(a24,1x,es8.1)') ' polarBeta: ',polarBeta
write(6,'(a24,1x,a)') ' PETSc_options: ',trim(petsc_defaultOptions)//' '//trim(petsc_options)
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
some polishing concerning openMP initialization number of threads is now written to *.out file
2010-12-02 16:34:29 +05:30
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
!--------------------------------------------------------------------------------------------------
! spectral parameters
#ifdef FEM
not needed for disabled output MPI rank>0
2016-06-29 20:05:49 +05:30
write(6,'(a24,1x,i8)') ' integrationOrder: ',integrationOrder
write(6,'(a24,1x,i8)') ' structOrder: ',structOrder
write(6,'(a24,1x,a)') ' PETSc_options: ',trim(petsc_defaultOptions)//' '//trim(petsc_options)
write(6,'(a24,1x,L8)') ' B-Bar stabilisation: ',BBarStabilisation
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
#endif
simplified multi processor (MPI) reporting
2014-10-10 18:38:34 +05:30
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
removed unused variables and declared external functions as external
2013-02-11 15:14:17 +05:30
!--------------------------------------------------------------------------------------------------
! sanity checks
pInt not needed
2019-06-11 18:09:51 +05:30
if (defgradTolerance <= 0.0_pReal) call IO_error(301,ext_msg='defgradTolerance')
if (iJacoStiffness < 1) call IO_error(301,ext_msg='iJacoStiffness')
if (nMPstate < 1) call IO_error(301,ext_msg='nMPstate')
if (subStepMinHomog <= 0.0_pReal) call IO_error(301,ext_msg='subStepMinHomog')
if (subStepSizeHomog <= 0.0_pReal) call IO_error(301,ext_msg='subStepSizeHomog')
if (stepIncreaseHomog <= 0.0_pReal) call IO_error(301,ext_msg='stepIncreaseHomog')
if (numerics_integrator <= 0 .or. numerics_integrator >= 6) &
call IO_error(301,ext_msg='integrator')
if (numerics_unitlength <= 0.0_pReal) call IO_error(301,ext_msg='unitlength')
if (absTol_RGC <= 0.0_pReal) call IO_error(301,ext_msg='absTol_RGC')
if (relTol_RGC <= 0.0_pReal) call IO_error(301,ext_msg='relTol_RGC')
if (absMax_RGC <= 0.0_pReal) call IO_error(301,ext_msg='absMax_RGC')
if (relMax_RGC <= 0.0_pReal) call IO_error(301,ext_msg='relMax_RGC')
if (pPert_RGC <= 0.0_pReal) call IO_error(301,ext_msg='pPert_RGC')
if (xSmoo_RGC <= 0.0_pReal) call IO_error(301,ext_msg='xSmoo_RGC')
if (viscPower_RGC < 0.0_pReal) call IO_error(301,ext_msg='viscPower_RGC')
if (viscModus_RGC < 0.0_pReal) call IO_error(301,ext_msg='viscModus_RGC')
if (refRelaxRate_RGC <= 0.0_pReal) call IO_error(301,ext_msg='refRelaxRate_RGC')
if (maxdRelax_RGC <= 0.0_pReal) call IO_error(301,ext_msg='maxdRelax_RGC')
if (maxVolDiscr_RGC <= 0.0_pReal) call IO_error(301,ext_msg='maxVolDiscr_RGC')
if (volDiscrMod_RGC < 0.0_pReal) call IO_error(301,ext_msg='volDiscrMod_RGC')
if (volDiscrPow_RGC <= 0.0_pReal) call IO_error(301,ext_msg='volDiscrPw_RGC')
if (residualStiffness < 0.0_pReal) call IO_error(301,ext_msg='residualStiffness')
if (itmax <= 1) call IO_error(301,ext_msg='itmax')
if (itmin > itmax .or. itmin < 1) call IO_error(301,ext_msg='itmin')
if (maxCutBack < 0) call IO_error(301,ext_msg='maxCutBack')
if (stagItMax < 0) call IO_error(301,ext_msg='maxStaggeredIter')
if (err_struct_tolRel <= 0.0_pReal) call IO_error(301,ext_msg='err_struct_tolRel')
if (err_struct_tolAbs <= 0.0_pReal) call IO_error(301,ext_msg='err_struct_tolAbs')
if (err_thermal_tolabs <= 0.0_pReal) call IO_error(301,ext_msg='err_thermal_tolabs')
if (err_thermal_tolrel <= 0.0_pReal) call IO_error(301,ext_msg='err_thermal_tolrel')
if (err_damage_tolabs <= 0.0_pReal) call IO_error(301,ext_msg='err_damage_tolabs')
if (err_damage_tolrel <= 0.0_pReal) call IO_error(301,ext_msg='err_damage_tolrel')
renaming "spectral" to "grid"
2019-03-23 13:57:58 +05:30
#ifdef Grid
pInt not needed
2019-06-11 18:09:51 +05:30
if (err_stress_tolrel <= 0.0_pReal) call IO_error(301,ext_msg='err_stress_tolRel')
if (err_stress_tolabs <= 0.0_pReal) call IO_error(301,ext_msg='err_stress_tolAbs')
if (err_div_tolRel < 0.0_pReal) call IO_error(301,ext_msg='err_div_tolRel')
if (err_div_tolAbs <= 0.0_pReal) call IO_error(301,ext_msg='err_div_tolAbs')
if (err_curl_tolRel < 0.0_pReal) call IO_error(301,ext_msg='err_curl_tolRel')
if (err_curl_tolAbs <= 0.0_pReal) call IO_error(301,ext_msg='err_curl_tolAbs')
restructured AL solver to resemble polarization scheme and introduced alpha and beta parameters in numerics to determine what kind of solution you want. alpha = beta = 1 results in the AL scheme. alpha = beta = 2 results in the accelerated scheme. safe values are anything in between. changed the solver to use stress calculations on the campatible field which seems to result in faster convergence.
2013-02-28 23:05:02 +05:30
if (polarAlpha <= 0.0_pReal .or. &
pInt not needed
2019-06-11 18:09:51 +05:30
polarAlpha > 2.0_pReal) call IO_error(301,ext_msg='polarAlpha')
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 22:50:05 +05:30
if (polarBeta < 0.0_pReal .or. &
pInt not needed
2019-06-11 18:09:51 +05:30
polarBeta > 2.0_pReal) call IO_error(301,ext_msg='polarBeta')
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
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
end subroutine numerics_init
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
removed (in IO.f90) a dangerous initialization statement. Please read http://www.cs.rpi.edu/~szymansk/OOF90/bugs.html#4 for more details. Other files are just a little bit polished
2012-03-07 15:37:29 +05:30
end module numerics