loop_library_interface.hpp

Go to the documentation of this file.

// ====================================================================
// This file is part of FlexibleSUSY.
//
// FlexibleSUSY is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published
// by the Free Software Foundation, either version 3 of the License,
// or (at your option) any later version.
//
// FlexibleSUSY 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
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with FlexibleSUSY.  If not, see
// <http://www.gnu.org/licenses/>.
// ====================================================================
 
#ifndef LOOP_LIBRARY_INTERFACE_H
#define LOOP_LIBRARY_INTERFACE_H
 
#include <boost/preprocessor/seq/enum.hpp>
#include <boost/preprocessor/seq/for_each.hpp>
#include <boost/preprocessor/seq/size.hpp>
#include <boost/preprocessor/seq/transform.hpp>
#include <array>
#include <complex>
 
#define ARGS_TYPE(R, COMMA, ELEM) std::complex<double> ELEM,
#define VIRTUAL(R, ARGS, NAME) virtual std::complex<double> NAME ARGS = 0;
 
#define A_ARGS_SEQ (m02_in)
#define B_ARGS_SEQ (p10_in)(m02_in)(m12_in)
#define C_ARGS_SEQ (p10_in)(p21_in)(p20_in)(m02_in)(m12_in)(m22_in)
#define D_ARGS_SEQ                                                             \
   (p10_in)(p21_in)(p32_in)(p30_in)(p20_in)(p31_in)(m02_in)(m12_in)(m22_in)(   \
      m32_in)
 
#define A_ARGS BOOST_PP_SEQ_FOR_EACH(ARGS_TYPE, , A_ARGS_SEQ) double scl2_in
#define B_ARGS BOOST_PP_SEQ_FOR_EACH(ARGS_TYPE, , B_ARGS_SEQ) double scl2_in
#define C_ARGS BOOST_PP_SEQ_FOR_EACH(ARGS_TYPE, , C_ARGS_SEQ) double scl2_in
#define D_ARGS BOOST_PP_SEQ_FOR_EACH(ARGS_TYPE, , D_ARGS_SEQ) double scl2_in
 
#define A_CSEQ (0)
#define B_CSEQ (0)(1)(00)
#define C_CSEQ (0)(1)(2)(00)(11)(12)(22)
#define D_CSEQ (0)(1)(2)(3)(00)(11)(12)(13)(22)(23)(33)
 
#define A_N BOOST_PP_SEQ_SIZE(A_CSEQ)
#define B_N BOOST_PP_SEQ_SIZE(B_CSEQ)
#define C_N BOOST_PP_SEQ_SIZE(C_CSEQ)
#define D_N BOOST_PP_SEQ_SIZE(D_CSEQ)
 
#define APPEND(s, data, elem) data##elem
#define A_SEQ BOOST_PP_SEQ_TRANSFORM(APPEND, A, A_CSEQ)
#define B_SEQ BOOST_PP_SEQ_TRANSFORM(APPEND, B, B_CSEQ)
#define C_SEQ BOOST_PP_SEQ_TRANSFORM(APPEND, C, C_CSEQ)
#define D_SEQ BOOST_PP_SEQ_TRANSFORM(APPEND, D, D_CSEQ)
 
namespace flexiblesusy
{
namespace looplibrary
{
 
using Acoeff_t = std::array<std::complex<double>, A_N>;
using Bcoeff_t = std::array<std::complex<double>, B_N>;
using Ccoeff_t = std::array<std::complex<double>, C_N>;
using Dcoeff_t = std::array<std::complex<double>, D_N>;
 
constexpr int aa0  = 0;
 
constexpr int bb0  = 0;
constexpr int bb1  = 1;
constexpr int bb00 = 2;
 
constexpr int cc0  = 0;
constexpr int cc1  = 1;
constexpr int cc2  = 2;
constexpr int cc00 = 3;
constexpr int cc11 = 4;
constexpr int cc12 = 5;
constexpr int cc22 = 6;
 
constexpr int dd0  = 0;
constexpr int dd1  = 1;
constexpr int dd2  = 2;
constexpr int dd3  = 3;
constexpr int dd00 = 4;
constexpr int dd11 = 5;
constexpr int dd12 = 6;
constexpr int dd13 = 7;
constexpr int dd22 = 8;
constexpr int dd23 = 9;
constexpr int dd33 = 10;
 
enum Acoeffs : int {
   BOOST_PP_SEQ_ENUM(BOOST_PP_SEQ_TRANSFORM(APPEND, a, A_CSEQ))
};
enum Bcoeffs : int {
   BOOST_PP_SEQ_ENUM(BOOST_PP_SEQ_TRANSFORM(APPEND, b, B_CSEQ))
};
enum Ccoeffs : int {
   BOOST_PP_SEQ_ENUM(BOOST_PP_SEQ_TRANSFORM(APPEND, c, C_CSEQ))
};
enum Dcoeffs : int {
   BOOST_PP_SEQ_ENUM(BOOST_PP_SEQ_TRANSFORM(APPEND, d, D_CSEQ))
};
 
class Loop_library_interface
{
public:
   BOOST_PP_SEQ_FOR_EACH(VIRTUAL, (A_ARGS), A_SEQ)
   BOOST_PP_SEQ_FOR_EACH(VIRTUAL, (B_ARGS), B_SEQ)
   BOOST_PP_SEQ_FOR_EACH(VIRTUAL, (C_ARGS), C_SEQ)
   BOOST_PP_SEQ_FOR_EACH(VIRTUAL, (D_ARGS), D_SEQ)
   virtual void A(Acoeff_t&, A_ARGS) = 0;
   virtual void B(Bcoeff_t&, B_ARGS) = 0;
   virtual void C(Ccoeff_t&, C_ARGS) = 0;
   virtual void D(Dcoeff_t&, D_ARGS) = 0;
   virtual ~Loop_library_interface() {}
};
} // namespace looplibrary
} // namespace flexiblesusy
 
#endif // LOOP_LIBRARY_INTERFACE_H