Commit da37bae0 authored by Lukas Riedel's avatar Lukas Riedel

added FEorder=3, removed FEorder=0

parent 3e456b86
......@@ -122,15 +122,15 @@ int main(int argc, char** argv)
auto grid = Dune::Dorie::build_grid_gmsh<Dune::UGGrid<2>>(inifile,helper);
auto gv = grid->leafGridView();
switch(FEorder){
case 0:
RichardsSolverSimplex<RF,0>(grid,gv,inifile,helper);
break;
case 1:
RichardsSolverSimplex<RF,1>(grid,gv,inifile,helper);
break;
case 2:
RichardsSolverSimplex<RF,2>(grid,gv,inifile,helper);
break;
case 3:
RichardsSolverSimplex<RF,3>(grid,gv,inifile,helper);
break;
default:
DUNE_THROW(Dune::NotImplemented,"Finite Element Order (grid.FEorder) not supported!");
}
......@@ -142,15 +142,15 @@ int main(int argc, char** argv)
// Avoid the creation of triangles upon refinement
grid->setClosureType(Dune::UGGrid<2>::ClosureType::NONE);
switch(FEorder){
case 0:
RichardsSolverRectangular<RF,0>(grid,gv,inifile,helper);
break;
case 1:
RichardsSolverRectangular<RF,1>(grid,gv,inifile,helper);
break;
case 2:
RichardsSolverRectangular<RF,2>(grid,gv,inifile,helper);
break;
case 3:
RichardsSolverRectangular<RF,3>(grid,gv,inifile,helper);
break;
default:
DUNE_THROW(Dune::NotImplemented,"Finite Element Order (grid.FEorder) not supported!");
}
......@@ -159,15 +159,15 @@ int main(int argc, char** argv)
auto grid = Dune::Dorie::build_grid_cube<Dune::YaspGrid<2>>(inifile,helper);
auto gv = grid->leafGridView();
switch(FEorder){
case 0:
RichardsSolverRectangular<RF,0>(grid,gv,inifile,helper);
break;
case 1:
RichardsSolverRectangular<RF,1>(grid,gv,inifile,helper);
break;
case 2:
RichardsSolverRectangular<RF,2>(grid,gv,inifile,helper);
break;
case 3:
RichardsSolverRectangular<RF,3>(grid,gv,inifile,helper);
break;
default:
DUNE_THROW(Dune::NotImplemented,"Finite Element Order (grid.FEorder) not supported!");
}
......@@ -183,15 +183,15 @@ int main(int argc, char** argv)
auto grid = Dune::Dorie::build_grid_gmsh<Dune::UGGrid<3>>(inifile,helper);
auto gv = grid->leafGridView();
switch(FEorder){
case 0:
RichardsSolverSimplex<RF,0>(grid,gv,inifile,helper);
break;
case 1:
RichardsSolverSimplex<RF,1>(grid,gv,inifile,helper);
break;
case 2:
RichardsSolverSimplex<RF,2>(grid,gv,inifile,helper);
break;
case 3:
RichardsSolverSimplex<RF,3>(grid,gv,inifile,helper);
break;
default:
DUNE_THROW(Dune::NotImplemented,"Finite Element Order (grid.FEorder) not supported!");
}
......@@ -203,15 +203,15 @@ int main(int argc, char** argv)
// Avoid the creation of triangles upon refinement
grid->setClosureType(Dune::UGGrid<3>::ClosureType::NONE);
switch(FEorder){
case 0:
RichardsSolverRectangular<RF,0>(grid,gv,inifile,helper);
break;
case 1:
RichardsSolverRectangular<RF,1>(grid,gv,inifile,helper);
break;
case 2:
RichardsSolverRectangular<RF,2>(grid,gv,inifile,helper);
break;
case 3:
RichardsSolverRectangular<RF,3>(grid,gv,inifile,helper);
break;
default:
DUNE_THROW(Dune::NotImplemented,"Finite Element Order (grid.FEorder) not supported!");
}
......@@ -220,15 +220,15 @@ int main(int argc, char** argv)
auto grid = Dune::Dorie::build_grid_cube<Dune::YaspGrid<3>>(inifile,helper);
auto gv = grid->leafGridView();
switch(FEorder){
case 0:
RichardsSolverRectangular<RF,0>(grid,gv,inifile,helper);
break;
case 1:
RichardsSolverRectangular<RF,1>(grid,gv,inifile,helper);
break;
case 2:
RichardsSolverRectangular<RF,2>(grid,gv,inifile,helper);
break;
case 3:
RichardsSolverRectangular<RF,3>(grid,gv,inifile,helper);
break;
default:
DUNE_THROW(Dune::NotImplemented,"Finite Element Order (grid.FEorder) not supported!");
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment