From 2d01a121bbbd7c7122dd2e216c23628c8857c1db Mon Sep 17 00:00:00 2001
From: Klaus <klaus.boehnlein@tu-dresden.de>
Date: Mon, 13 May 2024 21:54:05 +0200
Subject: [PATCH] add comment on Beff coefficients as output in .txt files

---
 .../EffectiveQuantitiesComputer.hh            |  2 +-
 .../buckling_experiment.py                    | 43 +++++++++++--------
 2 files changed, 25 insertions(+), 20 deletions(-)

diff --git a/dune/microstructure/EffectiveQuantitiesComputer.hh b/dune/microstructure/EffectiveQuantitiesComputer.hh
index f2035a55..2a02574a 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 52da4e23..611ebaa3 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
-- 
GitLab