Skip to content
Snippets Groups Projects
Commit 143c279e authored by Timo Koch's avatar Timo Koch
Browse files

[test] Use extrusion in rotsym grid geometry test

parent 65bab142
No related branches found
No related tags found
1 merge request!2188New concept for extrusion of grid geometries
...@@ -31,10 +31,7 @@ ...@@ -31,10 +31,7 @@
#include <dune/grid/utility/structuredgridfactory.hh> #include <dune/grid/utility/structuredgridfactory.hh>
#include <dune/grid/yaspgrid.hh> #include <dune/grid/yaspgrid.hh>
#include <dumux/discretization/cellcentered/tpfa/fvgridgeometry.hh> #include <dumux/discretization/cellcentered/tpfa/fvgridgeometry.hh>
#include <dumux/discretization/rotationpolicy.hh> #include <dumux/discretization/extrusion.hh>
#include <dumux/discretization/rotationsymmetricscv.hh>
#include <dumux/discretization/rotationsymmetricscvf.hh>
#include <dumux/discretization/rotationsymmetricgridgeometrytraits.hh>
namespace Dumux { namespace Dumux {
...@@ -49,11 +46,11 @@ void runTest(const GG& gg, const double refVolume, const double refSurface) ...@@ -49,11 +46,11 @@ void runTest(const GG& gg, const double refVolume, const double refSurface)
fvGeometry.bind(element); fvGeometry.bind(element);
for (const auto& scv : scvs(fvGeometry)) for (const auto& scv : scvs(fvGeometry))
volume += scv.volume(); volume += GG::Extrusion::volume(scv);
for (const auto& scvf : scvfs(fvGeometry)) for (const auto& scvf : scvfs(fvGeometry))
if (scvf.boundary()) if (scvf.boundary())
surface += scvf.area(); surface += GG::Extrusion::area(scvf);
} }
// compare to reference // compare to reference
...@@ -75,7 +72,9 @@ int main (int argc, char *argv[]) try ...@@ -75,7 +72,9 @@ int main (int argc, char *argv[]) try
// test the disc policy // test the disc policy
{ {
using Grid = Dune::YaspGrid<1, Dune::EquidistantOffsetCoordinates<double, 1>>; using Grid = Dune::YaspGrid<1, Dune::EquidistantOffsetCoordinates<double, 1>>;
using GGTraits = RotationSymmetricGridGeometryTraits<CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>, RotationPolicy::disc>;
struct GGTraits : public CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>
{ using Extrusion = RotationalExtrusion<0>; };
using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>; using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>;
using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition; using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition;
...@@ -104,7 +103,8 @@ int main (int argc, char *argv[]) try ...@@ -104,7 +103,8 @@ int main (int argc, char *argv[]) try
// test the ball policy // test the ball policy
{ {
using Grid = Dune::YaspGrid<1, Dune::EquidistantOffsetCoordinates<double, 1>>; using Grid = Dune::YaspGrid<1, Dune::EquidistantOffsetCoordinates<double, 1>>;
using GGTraits = RotationSymmetricGridGeometryTraits<CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>, RotationPolicy::ball>; struct GGTraits : public CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>
{ using Extrusion = SphericalExtrusion; };
using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>; using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>;
using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition; using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition;
...@@ -133,7 +133,8 @@ int main (int argc, char *argv[]) try ...@@ -133,7 +133,8 @@ int main (int argc, char *argv[]) try
// test the toroid policy // test the toroid policy
{ {
using Grid = Dune::YaspGrid<2, Dune::EquidistantOffsetCoordinates<double, 2>>; using Grid = Dune::YaspGrid<2, Dune::EquidistantOffsetCoordinates<double, 2>>;
using GGTraits = RotationSymmetricGridGeometryTraits<CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>, RotationPolicy::toroid>; struct GGTraits : public CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>
{ using Extrusion = RotationalExtrusion<0>; };
using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>; using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>;
using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition; using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition;
...@@ -165,7 +166,8 @@ int main (int argc, char *argv[]) try ...@@ -165,7 +166,8 @@ int main (int argc, char *argv[]) try
// test the toroid policy for perfect cylinder // test the toroid policy for perfect cylinder
{ {
using Grid = Dune::YaspGrid<2, Dune::EquidistantOffsetCoordinates<double, 2>>; using Grid = Dune::YaspGrid<2, Dune::EquidistantOffsetCoordinates<double, 2>>;
using GGTraits = RotationSymmetricGridGeometryTraits<CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>, RotationPolicy::toroid>; struct GGTraits : public CCTpfaDefaultGridGeometryTraits<typename Grid::LeafGridView>
{ using Extrusion = RotationalExtrusion<0>; };
using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>; using GridGeometry = CCTpfaFVGridGeometry<typename Grid::LeafGridView, /*caching=*/false, GGTraits>;
using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition; using GlobalPosition = typename GridGeometry::SubControlVolume::GlobalPosition;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment