shape_flags.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_OPERATORS_SHAPE_FLAGS_HPP_
#define MEDUSA_BITS_OPERATORS_SHAPE_FLAGS_HPP_
 
#include <string>
#include <tuple>
#include <medusa/bits/approximations/Operators.hpp>
 
namespace mm {
 
namespace sh {
 
typedef unsigned int shape_flags;
static const shape_flags d1 = 1;  
static const shape_flags lap = 2;  
static const shape_flags d2 = 4;  
static const shape_flags div = d1;  
static const shape_flags grad = d1;  
static const shape_flags graddiv = d2;  
static const shape_flags all = d1 | d2 | lap;  
 
inline std::string str(shape_flags f) {
    std::string s = "";
    if (f & sh::d1) s += " d1";
    if (f & sh::lap) s += " lap";
    if (f & sh::d2) s += " d2";
    return s;
}
 
template <shape_flags mask, int dim>
struct operator_tuple {
    static_assert(mask < 0, "Type for this mask not defined");
    typedef void type;   
};
template <int dim> struct operator_tuple<0, dim> {typedef std::tuple<> type; };  // NOLINT
template <int dim> struct operator_tuple<lap, dim> {typedef std::tuple<Lap<dim>> type; };  // NOLINT
template <int dim> struct operator_tuple<lap|d1, dim> { typedef std::tuple<Lap<dim>, Der1s<dim>> type; };  // NOLINT
template <int dim> struct operator_tuple<lap|d2, dim> { typedef std::tuple<Lap<dim>, Der2s<dim>> type; };  // NOLINT
template <int dim> struct operator_tuple<lap|d1|d2, dim> { typedef std::tuple<Lap<dim>, Der1s<dim>, Der2s<dim>> type; };  // NOLINT
template <int dim> struct operator_tuple<d1, dim> { typedef std::tuple<Der1s<dim>> type; };  // NOLINT
template <int dim> struct operator_tuple<d1|d2, dim> { typedef std::tuple<Der1s<dim>, Der2s<dim>> type; };  // NOLINT
template <int dim> struct operator_tuple<d2, dim> { typedef std::tuple<Der2s<dim>> type; };  // NOLINT
 
}  // namespace sh
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_OPERATORS_SHAPE_FLAGS_HPP_