constant_vectorfield.h

Go to the documentation of this file.

#ifndef TATOOINE_CONSTANT_VECTORFIELD_H
#define TATOOINE_CONSTANT_VECTORFIELD_H
//==============================================================================
#include "field.h"
//==============================================================================
namespace tatooine {
//==============================================================================
template <typename Real, size_t N>
struct constant_vectorfield : field<constant_vectorfield<Real, N>, Real, N, N> {
  //============================================================================
  // typedefs
  //============================================================================
  using this_type   = constant_vectorfield<Real, N>;
  using parent_type = field<this_type, Real, N, N>;
  using typename parent_type::pos_type;
  using typename parent_type::tensor_type;
  //============================================================================
  // members
  //============================================================================
  const tensor_type m_vector;
  //============================================================================
  // ctors
  //============================================================================
  private:
   template <size_t... Is>
   constexpr constant_vectorfield(std::index_sequence<Is...> /*is*/)
       : m_vector{1, ((void)Is, 0)...} {}
 
  public:
   constexpr constant_vectorfield()
       : constant_vectorfield{std::make_index_sequence<N - 1>{}} {}
   //----------------------------------------------------------------------------
   constexpr constant_vectorfield(const vec<Real, N>& vector)
       : m_vector{vector} {}
   //----------------------------------------------------------------------------
   constexpr constant_vectorfield(vec<Real, N>&& vector)
       : m_vector{std::move(vector)} {}
   //============================================================================
   // methods
   //============================================================================
   constexpr tensor_type evaluate(const pos_type& /*x*/, Real /*t*/) const {
     return m_vector;
  }
  //----------------------------------------------------------------------------
  constexpr bool in_domain(const pos_type& /*x*/, Real /*t*/) const {
    return true;
  }
};
//==============================================================================
}  // namespace tatooine
//==============================================================================
#endif