From bcd6b732ed10d87afe738c33843de596d98317bd Mon Sep 17 00:00:00 2001
From: Klaus Boehnlein <klaus.boehnlein@tu-dresden.de>
Date: Wed, 16 Jun 2021 10:34:55 +0200
Subject: [PATCH] backup
---
src/dune-microstructure.cc | 19 ++++++++++++-------
1 file changed, 12 insertions(+), 7 deletions(-)
diff --git a/src/dune-microstructure.cc b/src/dune-microstructure.cc
index 1c8835b5..e8cbf42b 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}};
--
GitLab