Random values

src/one-body-sample.cc

@@ -194,16 +194,16 @@ int main() {
                           P1Basis>(leafView, p1Basis, frictionalNodes,
                                    nodalIntegrals);
 
-      // TODO: populate on S_F
+      // TODO: Random value
       std::vector<double> normalStress;
       normalStress.resize(grid.size(grid.maxLevel(), dim));
-      std::fill(normalStress.begin(), normalStress.end(), 0.0);
+      std::fill(normalStress.begin(), normalStress.end(), 2.0);
 
-      // TODO: populate on S_F
+      // TODO: Random value
       std::vector<double> coefficientOfFriction;
       coefficientOfFriction.resize(grid.size(grid.maxLevel(), dim));
       std::fill(coefficientOfFriction.begin(), coefficientOfFriction.end(),
-                0.0);
+                1.0);
 
       Dune::GlobalNonlinearity<dim, Dune::LinearFunction> myGlobalNonlinearity(
           coefficientOfFriction, normalStress, nodalIntegrals);