diff --git a/code/config/numerics.config b/code/config/numerics.config index 39cfa858f..b64bdf2dd 100644 --- a/code/config/numerics.config +++ b/code/config/numerics.config @@ -52,11 +52,11 @@ discrepancyPower_RGC 5.0 fixed_seed 0 # put any number larger than zero, integer, if you want to have a pseudo random distribution ## spectral parameters ## -err_div_tol 1.0e-4 # 1.0e-4 proposed by Suquet -err_stress_tolrel 0.01 # relative tolerance for fullfillment of stress BC +err_div_tol 0.1 # Div(P)/avg(P)*meter +err_stress_tolrel 0.01 # relative tolerance for fulfillment of stress BC fftw_timelimit -1.0 # timelimit of plan creation for FFTW, see manual on www.fftw.org, Default -1.0: disable timelimit rotation_tol 1.0e-12 # tolerance of rotation specified in loadcase, Default 1.0e-12: first guess -fftw_plan_mode FFTW_PATIENT# reads the planing-rigor flag, see manual on www.fftw.org, Default FFTW_PATIENT: use patiant planner flag +fftw_plan_mode FFTW_PATIENT# reads the planing-rigor flag, see manual on www.fftw.org, Default FFTW_PATIENT: use patient planner flag itmax 20 # Maximum iteration number itmin 2 # Minimum iteration number memory_efficient 1 # Precalculate Gamma-operator (81 double per point) diff --git a/code/numerics.f90 b/code/numerics.f90 index 4d86be95f..2f00fb6d0 100644 --- a/code/numerics.f90 +++ b/code/numerics.f90 @@ -66,11 +66,11 @@ real(pReal) :: relevantStrain = 1.0e-7_pReal, & volDiscrMod_RGC = 1.0e+12_pReal, & ! stiffness of RGC volume discrepancy (zero = without volume discrepancy constraint) volDiscrPow_RGC = 5.0_pReal, & ! powerlaw penalty for volume discrepancy !* spectral parameters: - err_div_tol = 1.0e-4_pReal, & ! error of divergence in fourier space, Default 1.0e-4: proposed by Suquet + err_div_tol = 0.1_pReal, & ! Div(P)/avg(P)*meter err_stress_tolrel = 0.01_pReal , & ! relative tolerance for fullfillment of stress BC, Default: 0.01 allowing deviation of 1% of maximum stress fftw_timelimit = -1.0_pReal, & ! sets the timelimit of plan creation for FFTW, see manual on www.fftw.org, Default -1.0: disable timelimit rotation_tol = 1.0e-12_pReal ! tolerance of rotation specified in loadcase, Default 1.0e-12: first guess -character(len=64) :: fftw_plan_mode = 'FFTW_PATIENT' ! reads the planing-rigor flag, see manual on www.fftw.org, Default FFTW_PATIENT: use patiant planner flag +character(len=64) :: fftw_plan_mode = 'FFTW_PATIENT' ! reads the planing-rigor flag, see manual on www.fftw.org, Default FFTW_PATIENT: use patient planner flag integer(pInt) :: fftw_planner_flag = -1_pInt, & ! conversion of fftw_plan_mode to integer, basically what is usually done in the include file of fftw itmax = 20_pInt, & ! maximum number of iterations itmin = 2_pInt ! minimum number of iterations