diff --git a/inputs/cellsolver.parset b/inputs/cellsolver.parset
new file mode 100644
index 0000000000000000000000000000000000000000..3e4f5f20bd47c6008534aa430fee9ac7f6e89274
--- /dev/null
+++ b/inputs/cellsolver.parset
@@ -0,0 +1,112 @@
+
+
+#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)
+
+
+
+
+
+
+
+
diff --git a/outputs/debugLog.txt b/outputs/debugLog.txt
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/outputs/output.txt b/outputs/output.txt
new file mode 100644
index 0000000000000000000000000000000000000000..5800dd61c928f56f450a2aab60b5de00506d71d7
--- /dev/null
+++ b/outputs/output.txt
@@ -0,0 +1,57 @@
+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
diff --git a/outputs/output20x20x20.txt b/outputs/output20x20x20.txt
new file mode 100644
index 0000000000000000000000000000000000000000..24404acda473b1d76a8abc20d6f9f503c038fa53
--- /dev/null
+++ b/outputs/output20x20x20.txt
@@ -0,0 +1,57 @@
+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