FEApp_TemplateTypes.hpp

Go to the documentation of this file.

// $Id: FEApp_TemplateTypes.hpp,v 1.2 2008/05/16 19:56:56 etphipp Exp $ 
// $Source: /space/CVS/Trilinos/packages/sacado/example/FEApp/FEApp_TemplateTypes.hpp,v $ 
// @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., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps
// (etphipp@sandia.gov).
// 
// ***********************************************************************
// @HEADER

#ifndef FEAPP_TEMPLATETYPES_HPP
#define FEAPP_TEMPLATETYPES_HPP

// Includ the MPL vector
#include "Sacado_mpl_vector.hpp"

// Include all of our AD types
#include "Sacado_MathFunctions.hpp"
#include "Sacado_Fad_DFad.hpp"

// Typedef AD types to standard names
typedef double RealType;
typedef Sacado::Fad::DFad<double> FadType;

// Define which types we are using
#define REAL_ACTIVE 1
#define FAD_ACTIVE 1

// Conditionally enable Stokhos stochastic Galerkin expansions
#ifdef HAVE_SACADO_STOKHOS
#include "Sacado_PCE_OrthogPoly.hpp"
#include "Stokhos_HermiteEBasis.hpp"
typedef Sacado::PCE::OrthogPoly<double> SGType;
typedef Sacado::Fad::DFad< Sacado::PCE::OrthogPoly<double> > SGFadType;
#define SG_ACTIVE 1
#define SGFAD_ACTIVE 1
#else
#define SG_ACTIVE 0
#define SGFAD_ACTIVE 0
#endif

// Turn on/off explicit template instantiation
#define SACADO_ETI

// Build the MPL vector containing our valid types
typedef Sacado::mpl::vector<> ValidTypes0;
#if REAL_ACTIVE
typedef Sacado::mpl::push_back<ValidTypes0, RealType>::type ValidTypes1;
#else
typedef ValidTypes0 ValidTypes1;
#endif
#if FAD_ACTIVE
typedef Sacado::mpl::push_back<ValidTypes1, FadType>::type ValidTypes2;
#else
typedef ValidTypes1 ValidTypes2;
#endif
#if SG_ACTIVE
typedef Sacado::mpl::push_back<ValidTypes2, SGType>::type ValidTypes3;
#else
typedef ValidTypes2 ValidTypes3;
#endif
#if SGFAD_ACTIVE
typedef Sacado::mpl::push_back<ValidTypes3, SGFadType>::type ValidTypes4;
#else
typedef ValidTypes3 ValidTypes4;
#endif
typedef ValidTypes4 ValidTypes;

// Define macro for explicit template instantiation
#if REAL_ACTIVE
#define INSTANTIATE_TEMPLATE_CLASS_REAL(name) template class name<double>;
#else
#define INSTANTIATE_TEMPLATE_CLASS_REAL(name)
#endif

#if FAD_ACTIVE
#define INSTANTIATE_TEMPLATE_CLASS_FAD(name) template class name<FadType>;
#else
#define INSTANTIATE_TEMPLATE_CLASS_FAD(name)
#endif

#if SG_ACTIVE
#define INSTANTIATE_TEMPLATE_CLASS_SG(name) template class name<SGType>;
#else
#define INSTANTIATE_TEMPLATE_CLASS_SG(name)
#endif

#if SGFAD_ACTIVE
#define INSTANTIATE_TEMPLATE_CLASS_SGFAD(name) template class name<SGFadType>;
#else
#define INSTANTIATE_TEMPLATE_CLASS_SGFAD(name)
#endif

#define INSTANTIATE_TEMPLATE_CLASS(name) \
  INSTANTIATE_TEMPLATE_CLASS_REAL(name)  \
  INSTANTIATE_TEMPLATE_CLASS_FAD(name)   \
  INSTANTIATE_TEMPLATE_CLASS_SG(name)  \
  INSTANTIATE_TEMPLATE_CLASS_SGFAD(name)

#endif // FEAPP_TEMPLATETYPES_HPP

---