Added Input/Output folder

inputs/cellsolver.parset

+
+
+#path for logfile
+outputPath = "../../outputs/output.txt"
+
+#############################################
+#  Cell Domain
+#############################################
+# Domain 1: (-1/2, 1/2)^3  , Domain 2 : [0,1)^2 x (-1/2, 1/2)
+
+cellDomain = 1
+
+#############################################
+#  Grid parameters
+#############################################
+
+#nElements_Cell = 20 20 20	               # number elements in each direction (y1 y2 x3) 
+nElements_Cell = 15 15 15
+
+width = 1.0 	    
+
+
+gamma = 50.0
+
+#############################################
+#  Material parameters
+#############################################
+
+beta = 2.0    # ratio between material parameters mu1 & mu2 .... beta = 1.0 corresponds to homogeneous case 
+
+mu1 = 10.0
+lambda1 = 0.0 
+
+
+#############################################
+#  Prestrain parameters
+#############################################
+
+theta = 0.3    # volume fraction
+
+prestraintype = 2
+
+rho1 = 1.0
+
+alpha = 2.0    # ratio between prestrain parameters rho1 & rho2 
+
+
+
+# Prestrain Types: 
+
+#1 Isotropic Pressure
+    Func2Tensor B1_ = [this] (const Domain& x) {              // ISOTROPIC PRESSURE
+                if (abs(x[0]) > (theta/2)  && x[2] > 0)
+                    return MatrixRT{{p1, 0.0 , 0.0}, {0.0, p1, 0.0}, {0.0, 0.0, p1}};
+                if (abs(x[0]) < (theta/2)  && x[2] < 0)
+                    return MatrixRT{{p2, 0.0 , 0.0}, {0.0, p2, 0.0}, {0.0, 0.0, p2}};
+                else
+                    return MatrixRT{{0.0, 0.0 , 0.0}, {0.0, 0.0, 0.0}, {0.0, 0.0, 0.0}};
+
+            };
+
+
+
+#2 Isotropic Pressure
+
+Func2Tensor B2_ = [this] (const Domain& x) {              // ISOTROPIC PRESSURE
+
+                if (abs(x[0]) < (theta/2) && x[2] < 0 && x[2] > -(1.0/2.0) )
+                    return MatrixRT{{p1, 0.0 , 0.0}, {0.0, p1, 0.0}, {0.0, 0.0, p1}};
+                else
+                    return MatrixRT{{0.0, 0.0 , 0.0}, {0.0, 0.0, 0.0}, {0.0, 0.0, 0.0}};
+
+
+            };
+
+
+
+
+#############################################
+#  Assembly options
+#############################################
+
+set_IntegralZero = true
+
+
+#############################################
+#  Solver Type
+#############################################
+# 1: CG-Solver  # 2: GMRES   # 3: QR
+
+Solvertype = 1
+
+write_corrector_phi1 = false
+write_corrector_phi2 = false
+write_corrector_phi3 = false
+
+write_L2Error = true
+write_IntegralMean = true
+
+#############################################
+#  Define Analytic Solutions
+#############################################
+
+b1 = (-(theta/4.0)*mu1*mu2)/(theta*mu1+(1.0-theta)*mu2)
+
+
+
+
+
+
+
+

outputs/output.txt

+prestrain imp: 2
+rho1 = 1
+rho2 = 2
+alpha: 2
+gamma: 50
+Number of Elements in each direction: [15,15,15]
+beta: 2
+material parameters: 
+mu1: 10
+mu2: 20
+lambda1: 0
+lambda2:0
+size of FiniteElementBasis: 10800
+Corrector-Matrix M_1: 
+-8.73326e-10 3.3873e-09 0
+3.3873e-09 -9.35729e-22 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2: 
+-1.20023e-25 -9.12937e-18 0
+-9.12937e-18 1.96361e-17 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3: 
+-1.352e-11 -3.15777e-12 0
+-3.15777e-12 -3.55246e-18 0
+0 0 0
+
+
+Solution of Corrector problems:
+Computed Matrix Q: 
+1.97452 3.06664e-22 1.73925e-11
+-2.77279e-25 2.16049 1.34612e-18
+-8.52488e-12 1.3461e-18 1.97454
+
+Computed prestrain B_hat: 
+-0.922995 -1.48056 -1.03062e-07
+Computed prestrain B_eff: 
+-0.467452 -0.68529 -5.21975e-08
+computed q1: 1.97452
+computed q2: 2.16049
+computed q3: 1.97454
+computed b1: -0.467452
+computed b2: -0.68529
+computed b3: -5.21975e-08
+computed b1_hat: -0.922995
+computed b2_hat: -1.48056
+computed b3_hat: -1.03062e-07
+ --- analytic solutions: --- 
+b1 : -0.882353
+b2 : -1.5
+b3 : 0
+q1 : 1.96078
+q2 : 2.16667
+q3 should be between q1 and q2

outputs/output20x20x20.txt

+prestrain imp: 2
+rho1 = 1
+rho2 = 2
+alpha: 2
+gamma: 50
+Number of Elements in each direction: [20,20,20]
+beta: 2
+material parameters: 
+mu1: 10
+mu2: 20
+lambda1: 0
+lambda2:0
+size feBasis: 25200
+Corrector-Matrix M_1: 
+1.923e-10 8.39184e-10 0
+8.39184e-10 5.15403e-21 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2: 
+2.39727e-23 -4.18818e-17 0
+-4.18818e-17 4.08491e-18 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3: 
+2.71269e-11 -2.81909e-12 0
+-2.81909e-12 -1.6844e-17 0
+0 0 0
+
+
+Solution of Corrector problems:
+Computed Matrix Q: 
+1.96078 3.32328e-21 2.38035e-11
+1.20763e-25 2.16667 1.18138e-18
+1.23313e-11 1.18136e-18 1.9608
+
+Computed prestrain B_hat: 
+-0.882353 -1.5 -7.19027e-08
+Computed prestrain B_eff: 
+-0.45 -0.692308 -3.66673e-08
+computed q1: 1.96078
+computed q2: 2.16667
+computed q3: 1.9608
+computed b1: -0.45
+computed b2: -0.692308
+computed b3: -3.66673e-08
+computed b1_hat: -0.882353
+computed b2_hat: -1.5
+computed b3_hat: -7.19027e-08
+ --- analytic solutions: --- 
+b1 : -0.882353
+b2 : -1.5
+b3 : 0
+q1 : 1.96078
+q2 : 2.16667
+q3 should be between q1 and q2