Cosserat-Continuum-Nonplanar

@osander: Can you have a look here?

This MR cleans up the usage of MIXED_SPACE in cosserat-continuum.

Now all combinations of dim = dimworld or dim != dimworld and MIXED_SPACE = 0 or MIXED_SPACE = 1 at least compile.

Also, this MR points out the one combination that still can't be solved: Problems with dim != dimworld do not work for MIXED_SPACE = 1!

As soon as I need this for locking tests I will add this; I will have to change Nonplanarcosseratshellenergy so it can be assembled with the GFEMixedAssembler.

dune/gfe/mixedriemanniantrsolver.cc

@@ -414,8 +414,10 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis,Basis0,TargetSpace0,Basis1,
             double oldEnergy = 0;
             Dune::Timer solutionTimer;
             int ii = 0;
+            CorrectionType diff{corr_global};
             try {
               for (; ii<innerIterations_; ii++) {
+                diff=corr_global;
                 residual[_0] = rhs_global[_0];
                 stiffnessMatrix[_0][_1].mmv(corr_global[_1], residual[_0]);
                 mmgStep0->setRhs(residual[_0]);
@@ -436,7 +438,8 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis,Basis0,TargetSpace0,Basis1,
                   std::cout << "Warning: energy increase!" << std::endl;
   
                 oldEnergy = energy;
-                if (corr_global.infinity_norm() < innerTolerance_)
+                diff -= corr_global;
+                if (diff.infinity_norm() < innerTolerance_)
                   break;
               }
             } catch (Dune::Exception &e) {