Order.hpp

#pragma once
 
#include <amdis/Output.hpp>
#include <amdis/common/Math.hpp>
#include <amdis/common/Apply.hpp>
 
#include <dune/common/typetree/nodeconcepts.hh>
 
namespace AMDiS
{
  template <class N>
    requires Dune::TypeTree::Concept::TreeNode<N>
  int order(N const& node)
  {
    if constexpr (Dune::TypeTree::Concept::LeafTreeNode<N>)
      return node.finiteElement().localBasis().order();
    else if constexpr (Dune::TypeTree::Concept::UniformInnerTreeNode<N>)
      return order(node.child(0u));
    else if constexpr (Dune::TypeTree::Concept::StaticDegreeInnerTreeNode<N>)
      return Ranges::applyIndices<std::size_t(N::degree())>([&](auto... ii) {
        return Math::max(order(node.child(ii))...);
      });
    else {
      warning("Unknown basis-node type. Assuming polynomial degree 1.");
      return 1;
    }
  }
 
} // end namespace AMDiS