diff --git a/inputs/cellsolver.parset b/inputs/cellsolver.parset
index afc14b794b27b365fc97f110cbd79c85f00aa63c..c9a46000116651bf2f3f4ebd1bc041446e6d84f8 100644
--- a/inputs/cellsolver.parset
+++ b/inputs/cellsolver.parset
@@ -9,8 +9,8 @@
 #  Choose Output-path for logfile
 #############################################
 ### Remove/Comment this when running via Python-Script:
-#outputPath=/home/stefan/DUNE/dune-microstructure/outputs
-outputPath=/home/klaus/Desktop/Dune-Testing/dune-microstructure/outputs
+outputPath=/home/stefan/DUNE/dune-microstructure/outputs
+#outputPath=/home/klaus/Desktop/Dune-Testing/dune-microstructure/outputs
 
 
 #############################################
@@ -83,8 +83,8 @@ theta=0.25
 
 
 #--- choose composite-Implementation:
-#geometryFunctionPath = /home/stefan/DUNE/dune-microstructure/geometries
-geometryFunctionPath = /home/klaus/Desktop/Dune-Testing/dune-microstructure/geometries
+geometryFunctionPath = /home/stefan/DUNE/dune-microstructure/geometries
+#geometryFunctionPath = /home/klaus/Desktop/Dune-Testing/dune-microstructure/geometries
 material_prestrain_imp= "material_neukamm"   #(Python-indicator-function with same name determines material phases)
 #material_prestrain_imp= "two_phase_material_2"  #(Python-indicator-function with same name determines material phases)
 #material_prestrain_imp= "two_phase_material_3"  #(Python-indicator-function with same name determines material phases)
@@ -138,6 +138,6 @@ Solver_verbosity = 0  #(default = 2)  degree of information for solver output
 #write_corrector_phi1 = false
 #write_corrector_phi2 = false
 #write_corrector_phi3 = false
-write_corrector_phi1 = true
-write_corrector_phi2 = true
-write_corrector_phi3 = true
+#write_corrector_phi1 = true
+#write_corrector_phi2 = true
+#write_corrector_phi3 = true
diff --git a/src/Cell-Problem.cc b/src/Cell-Problem.cc
index 787eb9d05db12a1728bf333c6b47f3e9278a7760..682244ca658f50a2f5766735179b42ed66bc5ee5 100644
--- a/src/Cell-Problem.cc
+++ b/src/Cell-Problem.cc
@@ -2217,8 +2217,7 @@ int main(int argc, char *argv[])
             Q_2[a][b] = computeFullQ(Basis_CE, muLocal, lambdaLocal,gamma,mContainer[a],mContainer[b],*phiDerContainer[a],*phiDerContainer[b],x3MatrixBasis[a],x3MatrixBasis[b]);
         }
     printmatrix(std::cout, Q_2, "Matrix Q_2", "--");
-        
-    Q = Q_2; 
+//     Q = Q_2; 
     
     //--- VARIANT 3
     // Compute effective elastic law Q