diff --git a/dune/microstructure/EffectiveQuantitiesComputer.hh b/dune/microstructure/EffectiveQuantitiesComputer.hh
index 7dcb8d2bbb76bb4aa2735d3c7aebceea5f6721f3..d9f597b48437617f283d4960b95e811000b3e7f2 100644
--- a/dune/microstructure/EffectiveQuantitiesComputer.hh
+++ b/dune/microstructure/EffectiveQuantitiesComputer.hh
@@ -45,18 +45,11 @@ protected:
 
 public:
 
-	MatrixRT Q_;	    //effective moduli 
+	MatrixRT Qeff_;	    //effective moduli 
 	VectorRT Bhat_;		
 	VectorRT Beff_;	
 
 
-	// corrector parts
-	// VectorCT phi_1_;		
-	// VectorCT phi_2_;
-	// VectorCT phi_3_;
-	
-
- 
 	///////////////////////////////
 	// constructor
 	///////////////////////////////
@@ -74,16 +67,11 @@ public:
     // Getter
     ///////////////////////////////
 	CorrectorComputer<Basis,Material> getCorrectorComputer(){return correctorComputer_;}
+	const shared_ptr<Basis> getBasis() {return correctorComputer_.getBasis();}
 
-	const shared_ptr<Basis> getBasis()  
-	{
-		return correctorComputer_.getBasis();
-	}
-
-    auto getQeff(){return Q_;}
+    auto getQeff(){return Qeff_;}
     auto getBeff(){return Beff_;}
 
-
   // -----------------------------------------------------------------
   // --- Compute Effective Quantities
     void computeEffectiveQuantities()
@@ -112,7 +100,7 @@ public:
         auto localIndicatorFunction = material_.getLocalIndicatorFunction();
         
 
-        Q_ = 0 ;
+        Qeff_ = 0 ;
         Bhat_ = 0;
     
         for(size_t a=0; a < 3; a++)
@@ -197,13 +185,13 @@ public:
                 prestrain += elementPrestrain;
             
             }
-            Q_[a][b] = energy;    
+            Qeff_[a][b] = energy;    
             if (b==0)
                 Bhat_[a] = prestrain;
         }
         if(parameterSet.get<bool>("print_debug", false))
         {
-            printmatrix(std::cout, Q_, "Matrix Q_", "--");
+            printmatrix(std::cout, Qeff_, "Qeff_", "--");
             printvector(std::cout, Bhat_, "Bhat_", "--");
         }
 
@@ -212,15 +200,17 @@ public:
         //////////////////////////////
          
         // std::cout << "------- Information about Q matrix -----" << std::endl;        // TODO
-        // MatrixInfo<MatrixRT> matrixInfo(Q_,true,2,1);
+        // MatrixInfo<MatrixRT> matrixInfo(Qeff_,true,2,1);
         // std::cout << "----------------------------------------" << std::endl;
-        Q_.solve(Beff_,Bhat_);
+        Qeff_.solve(Beff_,Bhat_);
         if(parameterSet.get<bool>("print_debug", false))
             printvector(std::cout, Beff_, "Beff_", "--");
         
-    
         //LOG-Output
         auto& log = *(correctorComputer_.getLog());
+        log << "--- Effective moduli --- " << std::endl;
+        log << "Qeff_: \n" << Qeff_ << std::endl;
+        log << "------------------------ " << std::endl;
         log << "--- Prestrain Output --- " << std::endl;
         log << "Bhat_: " << Bhat_ << std::endl;
         log << "Beff_: " << Beff_ <<  " (Effective Prestrain)" << std::endl;
@@ -238,7 +228,7 @@ public:
     {
         std::cout << "write effective quantities to Matlab folder..." << std::endl;
         //writeMatrixToMatlab(Q, "../../Matlab-Programs/QMatrix.txt");
-        writeMatrixToMatlab(Q_, outputPath + "/QMatrix.txt");
+        writeMatrixToMatlab(Qeff_, outputPath + "/QMatrix.txt");
         // write effective Prestrain in Matrix for Output
         FieldMatrix<double,1,3> BeffMat;
         BeffMat[0] = Beff_;
@@ -318,7 +308,7 @@ public:
     //     auto prestrainGVF  = Dune::Functions::makeGridViewFunction(prestrain_, basis.gridView());
     //     auto prestrainFunctional = localFunction(prestrainGVF);   
 
-    //     Q_ = 0 ;
+    //     Qeff_ = 0 ;
     //     Bhat_ = 0;
     
     //     for(size_t a=0; a < 3; a++)
@@ -395,19 +385,19 @@ public:
     //             prestrain += elementPrestrain;
             
     //         }
-    //         Q_[a][b] = energy;    
+    //         Qeff_[a][b] = energy;    
     //         if (b==0)
     //             Bhat_[a] = prestrain;
     //     }
-    //     printmatrix(std::cout, Q_, "Matrix Q_", "--");
+    //     printmatrix(std::cout, Qeff_, "Matrix Qeff_", "--");
     //     printvector(std::cout, Bhat_, "Bhat_", "--");
     //     ///////////////////////////////
     //     // Compute effective Prestrain B_eff (by solving linear system)
     //     //////////////////////////////
     //     // std::cout << "------- Information about Q matrix -----" << std::endl;        // TODO
-    //     // MatrixInfo<MatrixRT> matrixInfo(Q_,true,2,1);
+    //     // MatrixInfo<MatrixRT> matrixInfo(Qeff_,true,2,1);
     //     // std::cout << "----------------------------------------" << std::endl;
-    //     Q_.solve(Beff_,Bhat_);
+    //     Qeff_.solve(Beff_,Bhat_);
     //     printvector(std::cout, Beff_, "Beff_", "--");
         
     //     //LOG-Output
diff --git a/src/Cell-Problem.cc b/src/Cell-Problem.cc
index b639fb8b407fbb817e50789a37c66031829ba225..ea11074be010547c982b4576054567a06d319a5f 100644
--- a/src/Cell-Problem.cc
+++ b/src/Cell-Problem.cc
@@ -58,7 +58,7 @@
 #include <dune/fufem/dunepython.hh>
 #include <python2.7/Python.h>
 
-#include <dune/fufem/functions/virtualgridfunction.hh> //TEST 
+// #include <dune/fufem/functions/virtualgridfunction.hh> //TEST 
 
 // #include <boost/multiprecision/cpp_dec_float.hpp>
 #include <any>