mast_data_types.h

Go to the documentation of this file.

/*
 * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
 * Copyright (C) 2013-2020  Manav Bhatia and MAST authors
 *
 * 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 Street, Fifth Floor, Boston, MA  02110-1301  USA
 */

#ifndef __mast__data_types__
#define __mast__data_types__

// MAST configuration
#include <mast/base/mast_config.h>

// Eigen includes
#include <Eigen/Dense>
#include <Eigen/SparseCore>

#if MAST_ENABLE_ADOLC == 1
// AdolC includes
#include <adolc/adtl.h>
#include <adolc/adouble.h>
// Eigen support for Adol-C
#include <unsupported/Eigen/AdolcForward>
// adouble for trace mode
typedef adouble              adouble_t;
// traceless adouble
typedef adtl::adouble        adouble_tl_t;
#endif

typedef double               real_t;
typedef unsigned int         uint_t;
typedef int                  int_t;
typedef std::complex<real_t> complex_t;
#define ComplexStepDelta     1.e-12

namespace MAST {

template <typename NodalScalarType, typename SolScalarType>
struct DeducedScalarType { };

template <>
struct DeducedScalarType<real_t, real_t> { using type = real_t;};

template <>
struct DeducedScalarType<std::complex<real_t>, real_t> { using type = std::complex<real_t>;};

template <>
struct DeducedScalarType<real_t, std::complex<real_t>> { using type = std::complex<real_t>;};

template <>
struct DeducedScalarType<std::complex<real_t>, std::complex<real_t>> { using type = std::complex<real_t>;};

#if MAST_ENABLE_ADOLC == 1
template <>
struct DeducedScalarType<adouble_tl_t, real_t> { using type = adouble_tl_t;};

template <>
struct DeducedScalarType<real_t, adouble_tl_t> { using type = adouble_tl_t;};

template <>
struct DeducedScalarType<adouble_tl_t, adouble_tl_t> { using type = adouble_tl_t;};
#endif
}

#endif // __mast__data_types__