barycentre_driver.hpp

// FEAT3: Finite Element Analysis Toolbox, Version 3
// Copyright (C) 2010 by Stefan Turek & the FEAT group
// FEAT3 is released under the GNU General Public License version 3,
// see the file 'copyright.txt' in the top level directory for details.
 
#pragma once
 
// includes, FEAT
#include <kernel/cubature/driver_base.hpp>
#include <kernel/util/meta_math.hpp>
 
namespace FEAT
{
  namespace Cubature
  {
    namespace Intern
    {
      template<typename Shape_>
      class BarycentreDriverBase :
        public DriverBase<Shape_>
      {
      public:
        static constexpr bool variadic = false;
        static constexpr int num_points = 1;
 
        static String name()
        {
          return "barycentre";
        }
 
        template<typename Functor_>
        static void alias(Functor_& functor)
        {
          functor.alias("midpoint");
        }
      };
    } // namespace Intern
 
    template<typename Shape_>
    class BarycentreDriver DOXY({});
 
    // Simplex specialization
    template<int dim_>
    class BarycentreDriver<Shape::Simplex<dim_> > :
      public Intern::BarycentreDriverBase<Shape::Simplex<dim_> >
    {
    public:
      template<
        typename Weight_,
        typename Coord_,
        typename Point_>
      static void fill(Rule<Shape::Simplex<dim_>, Weight_, Coord_, Point_>& rule)
      {
        rule.get_weight(0) = Weight_(1) / Weight_(MetaMath::Factorial<dim_>::value);
 
        // create coords of barycentre point
        for(int i(0); i < dim_; ++i)
        {
          rule.get_coord(0, i) = Coord_(1) / Coord_(dim_ + 1);
        }
      }
    }; // class BarycentreDriver<Simplex<...>>
 
    // Hypercube specialization
    template<int dim_>
    class BarycentreDriver<Shape::Hypercube<dim_> > :
      public Intern::BarycentreDriverBase<Shape::Hypercube<dim_> >
    {
    public:
      template<
        typename Weight_,
        typename Coord_,
        typename Point_>
      static void fill(Rule<Shape::Hypercube<dim_>, Weight_, Coord_, Point_>& rule)
      {
        rule.get_weight(0) = Weight_(1 << dim_);
 
        // create coords of barycentre point
        for(int i(0); i < dim_; ++i)
        {
          rule.get_coord(0, i) = Coord_(0);
        }
      }
    }; // class BarycentreDriver<Hypercube<...>>
  } // namespace Cubature
} // namespace FEAT