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

Merge branch 'minor-modernizations' into 'master'

Several minor modernization changes

See merge request !94
parents c7939b2d 114dc5b6
Branches
No related tags found
1 merge request!94Several minor modernization changes
Pipeline #39617 passed
...@@ -220,11 +220,11 @@ class MassAssembler : public LocalOperatorAssembler < GridType, TrialLocalFE, An ...@@ -220,11 +220,11 @@ class MassAssembler : public LocalOperatorAssembler < GridType, TrialLocalFE, An
double z = quad[pt].weight() * integrationElement; double z = quad[pt].weight() * integrationElement;
for(int i=0; i<rows; ++i) for(int i=0; i<rows; ++i)
{ {
double zi = values[i]*z; auto zi = values[i]*z;
for (int j=i+1; j<cols; ++j) for (int j=i+1; j<cols; ++j)
{ {
double zij = values[j] * zi; auto zij = values[j] * zi;
Dune::MatrixVector::addToDiagonal(localMatrix[i][j],zij); Dune::MatrixVector::addToDiagonal(localMatrix[i][j],zij);
Dune::MatrixVector::addToDiagonal(localMatrix[j][i],zij); Dune::MatrixVector::addToDiagonal(localMatrix[j][i],zij);
} }
......
...@@ -81,7 +81,7 @@ mark(const std::vector<RefinementIndicator<GridType>*>& refinementIndicators, ...@@ -81,7 +81,7 @@ mark(const std::vector<RefinementIndicator<GridType>*>& refinementIndicators,
const std::vector<GridType*>& grids, const std::vector<GridType*>& grids,
double fraction) double fraction)
{ {
typedef typename Dune::LeafMultipleCodimMultipleGeomTypeMapper<GridType> AllCodimMapper; typedef typename Dune::MultipleCodimMultipleGeomTypeMapper<typename GridType::LeafGridView> AllCodimMapper;
typedef typename std::multimap<double, std::pair<int,int>, std::greater<double> > ErrorMap; typedef typename std::multimap<double, std::pair<int,int>, std::greater<double> > ErrorMap;
if (grids.size() != refinementIndicators.size()) if (grids.size() != refinementIndicators.size())
......
...@@ -109,7 +109,7 @@ mark(const std::vector<std::shared_ptr<RefinementIndicator<GridType> > >& refine ...@@ -109,7 +109,7 @@ mark(const std::vector<std::shared_ptr<RefinementIndicator<GridType> > >& refine
minGridError[i] = std::numeric_limits<field_type>::max(); minGridError[i] = std::numeric_limits<field_type>::max();
maxGridError[i] = -std::numeric_limits<field_type>::max(); maxGridError[i] = -std::numeric_limits<field_type>::max();
Dune::LeafMultipleCodimMultipleGeomTypeMapper<GridType> p0Mapper(*grids[i], Dune::mcmgElementLayout()); Dune::MultipleCodimMultipleGeomTypeMapper<typename GridType::LeafGridView> p0Mapper(grids[i]->leafGridView(), Dune::mcmgElementLayout());
const auto& leafView = grids[i]->leafGridView(); const auto& leafView = grids[i]->leafGridView();
...@@ -169,7 +169,7 @@ mark(const std::vector<std::shared_ptr<RefinementIndicator<GridType> > >& refine ...@@ -169,7 +169,7 @@ mark(const std::vector<std::shared_ptr<RefinementIndicator<GridType> > >& refine
for (size_t i=0; i<grids.size(); i++) { for (size_t i=0; i<grids.size(); i++) {
Dune::LeafMultipleCodimMultipleGeomTypeMapper<GridType> p0Mapper(*grids[i], Dune::mcmgElementLayout()); Dune::MultipleCodimMultipleGeomTypeMapper<typename GridType::LeafGridView> p0Mapper(grids[i]->leafGridView(), Dune::mcmgElementLayout());
for (const auto& e : elements(grids[i]->leafGridView())) for (const auto& e : elements(grids[i]->leafGridView()))
if (refinedElements[i][p0Mapper.index(e)][0]) if (refinedElements[i][p0Mapper.index(e)][0])
......
...@@ -23,9 +23,9 @@ class CompositeQuadratureRule: ...@@ -23,9 +23,9 @@ class CompositeQuadratureRule:
CompositeQuadratureRule(const Dune::QuadratureRule<ct,dim>& quad, int refinement) CompositeQuadratureRule(const Dune::QuadratureRule<ct,dim>& quad, int refinement)
{ {
typedef Dune::StaticRefinement<Dune::Impl::SimplexTopology<dim>::type::id, typedef Dune::StaticRefinement<Dune::GeometryTypes::simplex(dim).id(),
ct, ct,
Dune::Impl::SimplexTopology<dim>::type::id, Dune::GeometryTypes::simplex(dim).id(),
dim> Refinement; dim> Refinement;
int numberOfSubelements = (1<<(dim*refinement)); int numberOfSubelements = (1<<(dim*refinement));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment