Sacado_TemplateContainer.hpp

Go to the documentation of this file.

// @HEADER
// ***********************************************************************
//
//                           Sacado Package
//                 Copyright (2006) Sandia Corporation
//
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
// the U.S. Government retains certain rights in this software.
//
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301
// USA
// Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps
// (etphipp@sandia.gov).
//
// ***********************************************************************
// @HEADER
 
#ifndef SACADO_TEMPLATE_CONTAINER_HPP
#define SACADO_TEMPLATE_CONTAINER_HPP
 
#include <type_traits>
 
// While this code does not directly use C++11 features, it uses mpl::vector,
// which does
#include "Sacado_ConfigDefs.h"
 
#include "Sacado_mpl_size.hpp"
#include "Sacado_mpl_find.hpp"
#include "Sacado_mpl_for_each.hpp"
#include "Sacado_mpl_apply.hpp"
#include "Sacado_mpl_begin.hpp"
#include "Sacado_mpl_end.hpp"
#include "Sacado_mpl_deref.hpp"
#include "Sacado_mpl_next.hpp"
#include "Sacado_mpl_enable_if.hpp"
 
namespace Sacado {
 
  namespace Impl {
 
    // Forward declaration
    template <typename TypeSeq,
              typename ObjectT,
              typename Iter1 = typename mpl::begin<TypeSeq>::type,
              typename Iter2 =typename  mpl::end<TypeSeq>::type>
    struct TupleSeq;
 
    // Forward declaration
    template <typename T,
              typename TypeSeq,
              typename ObjectT,
              typename Iter1 = typename mpl::begin<TypeSeq>::type,
              typename Iter2 = typename mpl::end<TypeSeq>::type,
              typename Enabled = void>
    struct GetTupleSeq;
 
  } // namespace Impl
 


  template <typename TypeSeq, typename ObjectT>
  class TemplateContainer {

    typedef Impl::TupleSeq<TypeSeq, ObjectT> tuple_type;

    template <typename BuilderOpT>
    struct BuildObject {
      tuple_type& objects;
      const BuilderOpT& builder;
      BuildObject(tuple_type& objects_,
                  const BuilderOpT& builder_) :
        objects(objects_), builder(builder_) {}
      template <typename T> void operator()(T x) const {
        Impl::GetTupleSeq<T,TypeSeq,ObjectT>::apply(objects) = builder(x);
      }
    };
 
  public:

    typedef TypeSeq types;

    struct DefaultBuilderOp {

      template<class T>
      typename Sacado::mpl::apply<ObjectT,T>::type operator() (T) const {
        return typename Sacado::mpl::apply<ObjectT,T>::type();
      }
 
    };

    TemplateContainer() {}

    ~TemplateContainer() {}

    template<typename T>
    typename Sacado::mpl::apply<ObjectT,T>::type& get() {
      return Impl::GetTupleSeq<T,TypeSeq,ObjectT>::apply(objects);
    }

    template<typename T>
    const typename Sacado::mpl::apply<ObjectT,T>::type& get() const  {
      return Impl::GetTupleSeq<T,TypeSeq,ObjectT>::apply(objects);
    }

    template <typename BuilderOpT = DefaultBuilderOp>
    void build(const BuilderOpT& builder) {
      Sacado::mpl::for_each_no_kokkos<TypeSeq>(
        BuildObject<BuilderOpT>(objects,builder));
    }
 
  private:

    tuple_type objects;
 
  };
 
  // Wrap call to mpl::for_each so you don't have to specify the container
  // or type sequence
  template <typename TypeSeq, typename ObjectT, typename FunctorT>
  void container_for_each(TemplateContainer<TypeSeq,ObjectT>& container,
                          const FunctorT& op) {
    typedef TemplateContainer<TypeSeq,ObjectT> Container;
    Sacado::mpl::for_each<Container> f(op);
  }
 
  // Wrap call to mpl::for_each so you don't have to specify the container
  // or type sequence
  template <typename TypeSeq, typename ObjectT, typename FunctorT>
  void container_for_each_no_kokkos(TemplateContainer<TypeSeq,ObjectT>& container,
                                    const FunctorT& op) {
    typedef TemplateContainer<TypeSeq,ObjectT> Container;
    Sacado::mpl::for_each_no_kokkos<Container> f(op);
  }
 
  namespace mpl {
 
    // Give TemplateContainer begin<> and end<> iterators for for_each
 
    template <typename TypeSeq, typename ObjectT>
    struct begin< TemplateContainer<TypeSeq,ObjectT> > {
      typedef typename begin<TypeSeq>::type type;
    };
 
    template <typename TypeSeq, typename ObjectT>
    struct end< TemplateContainer<TypeSeq,ObjectT> > {
      typedef typename end<TypeSeq>::type type;
    };
 
  }
 
  namespace Impl {
 
    // Container class to store our tuple sequence
    template <typename TypeSeq,
              typename ObjectT,
              typename Iter1,
              typename Iter2>
    struct TupleSeq :
      TupleSeq<TypeSeq, ObjectT, typename mpl::next<Iter1>::type, Iter2>
    {
      typedef typename mpl::apply<ObjectT,
                                  typename mpl::deref<Iter1>::type>::type type;
      type tail;
    };
 
    template <typename TypeSeq,
              typename ObjectT,
              typename Iter1>
    struct TupleSeq<TypeSeq, ObjectT, Iter1, Iter1> {};
 
    // Helper class to get a value out of the tuple sequence from a given type
    template <typename T,
              typename TypeSeq,
              typename ObjectT,
              typename Iter1,
              typename Iter2,
              typename Enabled>
    struct GetTupleSeq {};
 
    template <typename T,
              typename TypeSeq,
              typename ObjectT,
              typename Iter1,
              typename Iter2>
    struct GetTupleSeq< T,
                        TypeSeq,
                        ObjectT,
                        Iter1,
                        Iter2,
                        typename mpl::enable_if_c<
                          std::is_same< T, typename mpl::deref<Iter1>::type
                                        >::value
                          >::type > {
      static typename TupleSeq<TypeSeq,ObjectT,Iter1,Iter2>::type&
      apply(TupleSeq<TypeSeq,ObjectT,Iter1,Iter2>& t) {
        return t.tail;
      }
 
      static const typename TupleSeq<TypeSeq,ObjectT,Iter1,Iter2>::type&
      apply(const TupleSeq<TypeSeq,ObjectT,Iter1,Iter2>& t) {
        return t.tail;
      }
    };
 
    template <typename T,
              typename TypeSeq,
              typename ObjectT,
              typename Iter1,
              typename Iter2>
    struct GetTupleSeq< T,
                        TypeSeq,
                        ObjectT,
                        Iter1,
                        Iter2,
                        typename mpl::enable_if_c<
                          !std::is_same< T, typename mpl::deref<Iter1>::type
                                         >::value
                          >::type > :
      GetTupleSeq< T, TypeSeq, ObjectT, typename mpl::next<Iter1>::type, Iter2>
    {};
 
    template <typename T,
              typename TypeSeq,
              typename ObjectT,
              typename Iter1>
    struct GetTupleSeq< T, TypeSeq, ObjectT, Iter1, Iter1, void> {};
 
 
  } // namespace Impl
 
}
 
#endif