diff --git a/src/dune-microstructure.cc b/src/dune-microstructure.cc
index 1c8835b5f36abc63da5c01075feedc9046d3c044..e8cbf42b2d07bbdb2d388c80c87eb75d035aed68 100644
--- a/src/dune-microstructure.cc
+++ b/src/dune-microstructure.cc
@@ -1111,7 +1111,7 @@ int main(int argc, char *argv[])
     FieldVector<double,dim> lower({0.0, 0.0, -1.0/2.0});
     FieldVector<double,dim> upper({1.0, 1.0, 1.0/2.0});
 
-    int nE = 3;
+    int nE = 1;
 
     std::array<int,dim> nElements={nE,nE,nE};  //#Elements in each direction
 
@@ -1237,17 +1237,22 @@ bool set_oneBasisFunction_Zero = true;
 // using Range = FieldVector<double,dim>;
 // auto sourceTerm = [](const FieldVector<double,dim>& x){return Range{0.0, -1.0};};
 
-Func2Tensor x3G_1 = [] (const Domain& z) {
-	  		 return MatrixRT{{1.0*z[2], 0.0, 0.0}, {0.0, 0.0, 0.0}, {0.0, 0.0, 0.0}}; };
+Func2Tensor x3G_1 = [] (const Domain& x) {
+	  		 return MatrixRT{{1.0*x[2], 0.0, 0.0}, {0.0, 0.0, 0.0}, {0.0, 0.0, 0.0}}; };
 
-Func2Tensor x3G_2 = [] (const Domain& z) {
-	  		 return MatrixRT{{0.0, 0.0, 0.0}, {0.0, 1.0*z[2], 0.0}, {0.0, 0.0, 0.0}}; };
+Func2Tensor x3G_2 = [] (const Domain& x) {
+	  		 return MatrixRT{{0.0, 0.0, 0.0}, {0.0, 1.0*x[2], 0.0}, {0.0, 0.0, 0.0}}; };
 
-Func2Tensor x3G_3 = [] (const Domain& z) {
-	  		 return MatrixRT{{0.0, 0.5*z[2], 0.0}, {0.5*z[2], 0.0, 0.0}, {0.0, 0.0, 0.0}}; };
+Func2Tensor x3G_3 = [] (const Domain& x) {
+	  		 return MatrixRT{{0.0, 0.5*x[2], 0.0}, {0.5*x[2], 0.0, 0.0}, {0.0, 0.0, 0.0}}; };
 
 
 
+
+
+///////////////////////////////////////////////
+// Basis for R_sym^{2x2}
+//////////////////////////////////////////////
 MatrixRT G_1 {{1.0, 0.0, 0.0}, {0.0, 0.0, 0.0}, {0.0, 0, 0.0}};
 MatrixRT G_2 {{0.0, 0.0, 0.0}, {0.0, 1.0, 0.0}, {0, 0.0, 0.0}};
 MatrixRT G_3 {{0.0, 0.5, 0.0}, {0.5, 0.0, 0.0}, {0.0, 0.0, 0.0}};