diff --git a/dune/microstructure/EffectiveQuantitiesComputer.hh b/dune/microstructure/EffectiveQuantitiesComputer.hh index f2035a55b127769f34892bfba65639494ca3b439..2a02574a9206e88782867b663ea9f7150ffb09bc 100644 --- a/dune/microstructure/EffectiveQuantitiesComputer.hh +++ b/dune/microstructure/EffectiveQuantitiesComputer.hh @@ -40,7 +40,7 @@ public: MatrixRT Qeff_; // Effective moduli VectorRT Bhat_; - VectorRT Beff_; // Effective prestrain + VectorRT Beff_; // Effective prestrain COEFFICIENTS /////////////////////////////// diff --git a/experiment/buckling_experiment/buckling_experiment.py b/experiment/buckling_experiment/buckling_experiment.py index 52da4e23bdb1f561ed9da835584f9dd99ebd734e..611ebaa33e56e745f7777b2f362041302dbe5198 100644 --- a/experiment/buckling_experiment/buckling_experiment.py +++ b/experiment/buckling_experiment/buckling_experiment.py @@ -60,26 +60,26 @@ parameterSet.instrumented = 0 -#TEST RIemannian TR - +# --- (optional) Riemannian Trust-region solver: # parameterSet.Solver = "RiemannianTR" -parameterSet.numIt = 200 -parameterSet.nu1 = 3 -# Number of postsmoothing steps -parameterSet.nu2 = 3 -# Number of coarse grid corrections -parameterSet.mu = 1 -# Number of base solver iterations -parameterSet.baseIt = 100 -# Tolerance of the multigrid solver -parameterSet.mgTolerance = 1e-10 -# Tolerance of the base grid solver -parameterSet.baseTolerance = 1e-8 -parameterSet.tolerance = 1e-12 -# Max number of steps of the trust region solver -parameterSet.maxTrustRegionSteps = 100 -# Initial trust-region radius -parameterSet.initialTrustRegionRadius = 1 +# parameterSet.numIt = 200 +# parameterSet.nu1 = 3 +# # Number of postsmoothing steps +# parameterSet.nu2 = 3 +# # Number of coarse grid corrections +# parameterSet.mu = 1 +# # Number of base solver iterations +# parameterSet.baseIt = 100 +# # Tolerance of the multigrid solver +# parameterSet.mgTolerance = 1e-10 +# # Tolerance of the base grid solver +# parameterSet.baseTolerance = 1e-8 +# parameterSet.tolerance = 1e-12 +# # Max number of steps of the trust region solver +# parameterSet.maxTrustRegionSteps = 100 +# # Initial trust-region radius +# parameterSet.initialTrustRegionRadius = 1 + ############################ @@ -431,3 +431,8 @@ parameterSet.cacheElementMatrices = 1 # --- check orthogonality (75) from paper: parameterSet.write_checkOrthogonality = 0 + +# --- write effective quantities (Qhom,Beff) to .txt-files +# Qhom is written as a matrix +# Beff is written as Coefficient-vector +parameterSet.write_EffectiveQuantitiesToTxt = True