diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index e91cb77891d3e3b86db62f6789a3c4b09969f10c..7e5a5c23b32920a1d78d0407b7c9ab2b480d97b5 100644 --- a/src/one-body-sample.cc +++ b/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);