Add computation of MuGamma via the 2D Poisson-Type Equation

dune/microstructure/matrix_operations.hh

@@ -83,14 +83,21 @@ namespace MatrixOperations {
 	    return sum;
 	}
 
-	static double linearizedStVenantKirchhoffDensity(double mu, double lambda, MatrixRT E1, MatrixRT E2){
-		E1= sym(E1);
-		E2 = sym(E2);
-		double t1 = scalarProduct(E1,E2);
-		t1 *= 2* mu;
-		double t2 = trace(E1)*trace(E2); 
-		t2 *= lambda;
-		return t1 + t2;
+// 	static double linearizedStVenantKirchhoffDensity(double mu, double lambda, MatrixRT E1, MatrixRT E2){    // ?? Check with Robert
+// 		E1= sym(E1);
+// 		E2 = sym(E2);
+// 		double t1 = scalarProduct(E1,E2);
+// 		t1 *= 2* mu;
+// 		double t2 = trace(E1)*trace(E2); 
+// 		t2 *= lambda;
+// 		return t1 + t2;
+// 	}
+	
+    static double linearizedStVenantKirchhoffDensity(double mu, double lambda, MatrixRT E1, MatrixRT E2)  // CHANGED
+    {  
+        auto t1 = 2.0 * mu * sym(E1) + lambda * trace(sym(E1)) * Id();
+        return scalarProduct(t1,E2);
+
 	}
 
     // --- Generalization: Define Quadratic QuadraticForm

src/CMakeLists.txt

@@ -5,6 +5,7 @@
 
 set(programs Cell-Problem
 	     Cell-Problem_muGamma
+	     Compute_MuGamma
              )