PHAL_FactoryTraits.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#ifndef PHAL_FACTORY_TRAITS_HPP
#define PHAL_FACTORY_TRAITS_HPP

// User Defined Evaluator Types

#ifdef ALBANY_LCM
#include "LCM/evaluators/bc/KfieldBC.hpp"
#include "LCM/evaluators/bc/TimeDepBC.hpp"
#include "LCM/evaluators/bc/TimeTracBC.hpp"
#include "LCM/evaluators/Time.hpp"
#include "LCM/evaluators/bc/TorsionBC.hpp"
#endif
#ifdef ALBANY_QCAD
#include "QCAD_PoissonDirichlet.hpp"
#include "QCAD_PoissonNeumann.hpp"
#endif
#include "PHAL_Dirichlet.hpp"
#include "PHAL_Neumann.hpp"
#include "PHAL_GatherCoordinateVector.hpp"
#include "PHAL_GatherSolution.hpp"
#include "PHAL_GatherAuxData.hpp"
#ifdef ALBANY_FELIX
  #include "PHAL_GatherBasalFriction.hpp"
  #include "PHAL_GatherThickness.hpp"
  #include "PHAL_GatherSHeight.hpp"
#endif
#include "PHAL_DirichletCoordinateFunction.hpp"


#include "boost/mpl/vector/vector50.hpp"
#include "boost/mpl/placeholders.hpp"

// \cond  Have doxygern ignore this namespace
using namespace boost::mpl::placeholders;
// \endcond

namespace PHAL {
  template<typename Traits>
  struct DirichletFactoryTraits {

    static const int id_dirichlet                      =  0;
    static const int id_dirichlet_aggregator           =  1;
    static const int id_dirichlet_coordinate_function  =  2;
    static const int id_qcad_poisson_dirichlet         =  3;
    static const int id_kfield_bc                      =  4; // Only for LCM probs
    static const int id_timedep_bc                     =  5; // Only for LCM probs
    static const int id_time                           =  6; // Only for LCM probs
    static const int id_torsion_bc                     =  7; // Only for LCM probs

#ifdef ALBANY_LCM
    typedef boost::mpl::vector8<
#else
      typedef boost::mpl::vector4<
#endif
        PHAL::Dirichlet<_,Traits>,                //  0
        PHAL::DirichletAggregator<_,Traits>,      //  1
        PHAL::DirichletCoordFunction<_,Traits>,    //  2
#ifdef ALBANY_QCAD
        QCAD::PoissonDirichlet<_,Traits>         //  3
#else
        PHAL::Dirichlet<_,Traits>                //  3 dummy
#endif
#ifdef ALBANY_LCM
        , LCM::KfieldBC<_,Traits>,                //  4
        LCM::TimeDepBC<_, Traits>,                //  5
        LCM::Time<_, Traits>,                     //  6
        LCM::TorsionBC<_, Traits>                 //  7
#endif
        > EvaluatorTypes;
};


  template<typename Traits>
  struct NeumannFactoryTraits {

    static const int id_neumann                   =  0;
    static const int id_neumann_aggregator        =  1;
    static const int id_qcad_poisson_neumann      =  2;
    static const int id_gather_coord_vector       =  3;
    static const int id_gather_solution           =  4;
    static const int id_timedep_bc                =  5; // Only for LCM probs
    static const int id_gather_basalFriction      =  6; // Only for FELIX probs
    static const int id_gather_thickness        =  7; // Only for FELIX probs
    static const int id_gather_surfaceHeight      =  8; // Only for FELIX probs

#ifdef ALBANY_FELIX
    typedef boost::mpl::vector9<
#else
#ifdef ALBANY_LCM
    typedef boost::mpl::vector6<
#else
    typedef boost::mpl::vector5<
#endif
#endif

       PHAL::Neumann<_,Traits>,                     //  0
       PHAL::NeumannAggregator<_,Traits>,           //  1
#ifdef ALBANY_QCAD
       QCAD::PoissonNeumann<_,Traits>,              //  2
#else
       PHAL::Neumann<_,Traits>,                     //  2 dummy
#endif
         PHAL::GatherCoordinateVector<_,Traits>,      //  3
         PHAL::GatherSolution<_,Traits>               //  4
#ifdef ALBANY_LCM
         , LCM::TimeTracBC<_, Traits>                 //  5
#else
#ifdef ALBANY_FELIX
         , PHAL::Neumann<_,Traits>            //  5 dummy
#endif
#endif
#ifdef ALBANY_FELIX
      , PHAL::GatherBasalFriction<_,Traits>         //  6
    , PHAL::GatherThickness<_,Traits>             //  7
        , PHAL::GatherSHeight<_,Traits>               //  8
#endif
    > EvaluatorTypes;
};

}

#endif