integral.h

Go to the documentation of this file.

/*
 * integral.h - scalar integrals
 *
 * this file is part of PJFry library
 * Copyright 2011 Valery Yundin
 */
 
#ifndef QUL_INTEGRAL_H
#define QUL_INTEGRAL_H
 
#include "cache.h"
#include "common.h"
#include "kinem.h"
 
// #define USE_ONELOOP 1
// #define USE_QCDLOOP 1
 
class Initialize {
public:
  Initialize();
  ~Initialize();
};
 
ICache::Ival qlI1( const Kinem1& k );
ICache::Ival qlI2( const Kinem2& k );
ICache::Ival qlI3( const Kinem3& k );
ICache::Ival qlI4( const Kinem4& k );
 
#ifdef USE_QCDLOOP
#  define qlprec F77_FUNC( ffprec, FFPREC )
#  define qlflag F77_FUNC( ffflag, FFFLAG )
#  ifdef HAVE_QCDLOOP1_H
#    include <qcdloop1.h>
#  else
extern "C" {
void F77_FUNC( qlinit, QLINIT )();
void F77_FUNC( ffexi, FFEXI )();
 
#    ifdef USE_F2C
void F77_FUNC( qli1, QLI1 )( std::complex<double>* rslt, double* m1, double* mu2, int* ep );
void F77_FUNC( qli2, QLI2 )( std::complex<double>* rslt, double* p1, double* m1, double* m2,
                             double* mu2, int* ep );
void F77_FUNC( qli3, QLI3 )( std::complex<double>* rslt, double* p1, double* p2, double* p3,
                             double* m1, double* m2, double* m3, double* mu2, int* ep );
void F77_FUNC( qli4, QLI4 )( std::complex<double>* rslt, double* p1, double* p2, double* p3,
                             double* p4, double* s12, double* s23, double* m1, double* m2,
                             double* m3, double* m4, double* mu2, int* ep );
#    else
std::complex<double> F77_FUNC( qli1, QLI1 )( double* m1, double* mu2, int* ep );
std::complex<double> F77_FUNC( qli2, QLI2 )( double* p1, double* m1, double* m2, double* mu2,
                                             int* ep );
std::complex<double> F77_FUNC( qli3, QLI3 )( double* p1, double* p2, double* p3, double* m1,
                                             double* m2, double* m3, double* mu2, int* ep );
std::complex<double> F77_FUNC( qli4, QLI4 )( double* p1, double* p2, double* p3, double* p4,
                                             double* s12, double* s23, double* m1, double* m2,
                                             double* m3, double* m4, double* mu2, int* ep );
#    endif
 
extern struct {
  double xloss, precx, precc, xalogm, xclogm, xalog2, xclog2, reqprc;
} qlprec;
extern struct {
  int lwrite, ltest, l4also, ldc3c4, lmem, lwarn, ldot, nevent, ner, id, idsub, nwidth, nschem,
      onshel, idot;
} qlflag;
}
#  endif /* HAVE_QCDLOOP1_H */
#endif   /* USE_QCDLOOP */
 
#ifdef USE_ONELOOP
extern "C" {
void F77_FUNC_( avh_olo_mu_set, AVH_OLO_MU_SET )( double* mu );
void F77_FUNC_( avh_olo_onshell, AVH_OLO_ONSHELL )( double* thrs );
void F77_FUNC_( avh_olo_a0m, AVH_OLO_A0M )( std::complex<double>* rslt, double* m1 );
void F77_FUNC_( avh_olo_b0m, AVH_OLO_B0M )( std::complex<double>* rslt, double* p1, double* m1,
                                            double* m2 );
void F77_FUNC_( avh_olo_c0m, AVH_OLO_C0M )( std::complex<double>* rslt, double* p1, double* p2,
                                            double* p3, double* m1, double* m2, double* m3 );
void F77_FUNC_( avh_olo_d0m, AVH_OLO_D0M )( std::complex<double>* rslt, double* p1, double* p2,
                                            double* p3, double* p4, double* p12, double* p23,
                                            double* m1, double* m2, double* m3, double* m4 );
}
#endif /* USE_ONELOOP */
 
#endif /* QUL_INTEGRAL_H */