Skip to content
Snippets Groups Projects
Commit f46f3e1e authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Do not hand out names so freely

parent 62a23c9c
No related branches found
No related tags found
No related merge requests found
...@@ -151,16 +151,15 @@ int main() { ...@@ -151,16 +151,15 @@ int main() {
size_t const solver_maxIterations = 100; size_t const solver_maxIterations = 100;
double const solver_tolerance = 1e-4; double const solver_tolerance = 1e-4;
// Set up grid
typedef Dune::YaspGrid<dim> GridType; typedef Dune::YaspGrid<dim> GridType;
Dune::FieldVector<double, dim> const end_points( Dune::FieldVector<double, dim> const end_points(
1); // nth dimension (zero-indexed) goes from 0 to end_points[n] 1); // nth dimension (zero-indexed) goes from 0 to end_points[n]
Dune::FieldVector<int, dim> const elements( GridType grid(
2); // number of elements in each direction end_points,
Dune::FieldVector<bool, dim> const periodic(false); Dune::FieldVector<int, dim>(2), // number of elements in each direction
Dune::FieldVector<bool, dim>(false), // non-periodic in each direction
GridType grid(end_points, elements, periodic, 0); 0); // zero overlap (whatever that is)
grid.globalRefine(refinements); grid.globalRefine(refinements);
typedef P1NodalBasis<GridType::LeafGridView, double> P1Basis; typedef P1NodalBasis<GridType::LeafGridView, double> P1Basis;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment