diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 4b56ab139f35aa23555ec3c9ddfcc39165025aeb..d31ab8aa6216ae1d922c87ef0a12738c26debea4 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -112,17 +112,16 @@ int main(int argc, char *argv[]) { // {{{ Set up grid using Grid = Dune::ALUGrid<dims, dims, Dune::simplex, Dune::nonconforming>; - Dune::FieldVector<typename Grid::ctype, dims> lowerLeft(0); - Dune::FieldVector<typename Grid::ctype, dims> upperRight(1); - upperRight[0] = parset.get<size_t>("body.width"); - upperRight[1] = parset.get<size_t>("body.height"); + Dune::FieldVector<size_t, dims> lowerLeft(0); + Dune::FieldVector<size_t, dims> upperRight(1); + upperRight[0] = 5; + upperRight[1] = 1; std::shared_ptr<Grid> grid; { Dune::array<unsigned int, dims> elements; - std::fill(elements.begin(), elements.end(), 1); - elements[0] = parset.get<size_t>("body.width"); - elements[1] = parset.get<size_t>("body.height"); + for (size_t i = 0; i < dims; ++i) + elements[i] = upperRight[i]; // 1x1 elements grid = Dune::StructuredGridFactory<Grid>::createSimplexGrid( lowerLeft, upperRight, elements); diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index 5f61cf30e857cfa8f21d6a0221c6d2b806e70c9d..e10ed132aa7c24f89fc1eb462361b7d3468e43dc 100644 --- a/src/one-body-sample.parset +++ b/src/one-body-sample.parset @@ -10,8 +10,6 @@ finalTime = 15 youngModulus = 5e7 poissonRatio = 0.3 # The closer we get to 0.5, the more wiggly everything gets density = 5000 -height = 1 -width = 5 shearViscosity = 0 bulkViscosity = 0