diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 75ec2f2f70553f980484cbe52ec8d780edf4b0b5..4842a445d6c09b89d4cf36f51f2665ea8704d57d 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -220,14 +220,17 @@ int main(int argc, char *argv[]) { // {{{ Set up grid typedef Dune::ALUGrid<dim, dim, Dune::simplex, Dune::nonconforming> GridType; + Dune::FieldVector<typename GridType::ctype, dim> lowerLeft(0); - Dune::FieldVector<typename GridType::ctype, dim> upperRight(1); // depth 1 + Dune::FieldVector<typename GridType::ctype, dim> upperRight(1); upperRight[0] = parset.get<size_t>("body.width"); upperRight[1] = parset.get<size_t>("body.height"); + Dune::array<unsigned int, dim> elements; std::fill(elements.begin(), elements.end(), 1); elements[0] = parset.get<size_t>("body.width"); elements[1] = parset.get<size_t>("body.height"); + auto grid = Dune::StructuredGridFactory<GridType>::createSimplexGrid( lowerLeft, upperRight, elements); @@ -436,11 +439,8 @@ int main(int argc, char *argv[]) { { // Write all kinds of data state_writer << alpha[specialNode][0] << " " << std::endl; - displacement_writer << u[specialNode][0] << " " << std::endl; - velocity_writer << ud[specialNode][0] << " " << std::endl; - coefficient_writer << mu + a *std::log(ud[specialNode].two_norm() / V0) + b * (alpha[specialNode] + std::log(V0 / L))