Skip to content
Snippets Groups Projects
Commit 103962a8 authored by oliver.sander_at_tu-dresden.de's avatar oliver.sander_at_tu-dresden.de
Browse files

The PeriodicBasis code is in the Experimental namespace now

parent edff98e0
No related branches found
No related tags found
No related merge requests found
......@@ -50,7 +50,7 @@ const int order = 1;
using namespace Dune;
template <class GridView>
Functions::BasisFactory::PeriodicIndexSet getPeriodicIndices(const GridView& gridView)
Functions::BasisFactory::Experimental::PeriodicIndexSet getPeriodicIndices(const GridView& gridView)
{
// Check whether two points are equal on R/Z x R/Z
auto equivalent = [](const FieldVector<double,2>& x, const FieldVector<double,2>& y)
......@@ -69,7 +69,7 @@ Functions::BasisFactory::PeriodicIndexSet getPeriodicIndices(const GridView& gri
std::cout << "Grid has " << boundaryVertices.size() << " boundary vertices" << std::endl;
Functions::BasisFactory::PeriodicIndexSet periodicIndices;
Functions::BasisFactory::Experimental::PeriodicIndexSet periodicIndices;
if (order!=1)
DUNE_THROW(NotImplemented, "Periodicity detection only implemented for order==1");
......@@ -280,7 +280,7 @@ int main (int argc, char *argv[]) try
};
using namespace Dune::Functions::BasisFactory;
Functions::BasisFactory::PeriodicIndexSet periodicIndices;
Functions::BasisFactory::Experimental::PeriodicIndexSet periodicIndices;
if (periodic)
periodicIndices = getPeriodicIndices(gridView);
......@@ -288,7 +288,7 @@ int main (int argc, char *argv[]) try
auto feBasis = makeBasis(
gridView,
power<dim>(
Functions::BasisFactory::periodic(lagrange<order>(), periodicIndices),
Functions::BasisFactory::Experimental::periodic(lagrange<order>(), periodicIndices),
blockedInterleaved()
));
using FEBasis = decltype(feBasis);
......@@ -573,11 +573,11 @@ int main (int argc, char *argv[]) try
// For the restriction operators: FE bases on all levels
// This particular object is only needed to get the type of the level basis
PeriodicIndexSet dummyPeriodicIndices;
Experimental::PeriodicIndexSet dummyPeriodicIndices;
auto dummyLevelBasis = makeBasis(
grid->levelGridView(0),
power<dim>(
Functions::BasisFactory::periodic(lagrange<order>(), dummyPeriodicIndices),
Functions::BasisFactory::Experimental::periodic(lagrange<order>(), dummyPeriodicIndices),
blockedInterleaved()
));
......@@ -587,13 +587,13 @@ int main (int argc, char *argv[]) try
std::vector<std::shared_ptr<LevelBasis> > levelBases(numLevels);
for (int j=0; j<numLevels; j++)
{
PeriodicIndexSet periodicIndices;
Experimental::PeriodicIndexSet periodicIndices;
if (periodic)
periodicIndices = getPeriodicIndices(grid->levelGridView(j));
levelBases[j] = std::make_shared<LevelBasis>(makeBasis(
grid->levelGridView(j),
power<dim>(
Functions::BasisFactory::periodic(lagrange<order>(), periodicIndices),
Functions::BasisFactory::Experimental::periodic(lagrange<order>(), periodicIndices),
blockedInterleaved()
)));
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment