1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
// -*- C++ -*-
//
// Complex.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2011 Leif Lonnblad
//
// ThePEG is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_Complex_H
#define ThePEG_Complex_H
//
// This is file wraps the standard complex header and makes some
// convenient typedefs in the ThePEG namespace.
//
#include <complex>
namespace ThePEG {
using std::complex;
/** ThePEG code should use Complex for all complex scalars */
typedef std::complex<double> Complex;
/** @cond TRAITSPECIALIZATIONS */
template <typename T, typename U>
struct BinaryOpTraits;
template <typename T, typename U>
struct BinaryOpTraits<complex<T>, U> {
/** The type resulting from multiplication of the template type with
itself. */
typedef complex<typename BinaryOpTraits<T,U>::MulT> MulT;
/** The type resulting from division of one template type with
another. */
typedef complex<typename BinaryOpTraits<T,U>::DivT> DivT;
};
template <typename T, typename U>
struct BinaryOpTraits<T, complex<U> > {
/** The type resulting from multiplication of the template type with
itself. */
typedef complex<typename BinaryOpTraits<T,U>::MulT> MulT;
/** The type resulting from division of one template type with
another. */
typedef complex<typename BinaryOpTraits<T,U>::DivT> DivT;
};
template <typename T, typename U>
struct BinaryOpTraits<complex<T>, complex<U> > {
/** The type resulting from multiplication of the template type with
itself. */
typedef complex<typename BinaryOpTraits<T,U>::MulT> MulT;
/** The type resulting from division of one template type with
another. */
typedef complex<typename BinaryOpTraits<T,U>::DivT> DivT;
};
template <typename T>
struct BinaryOpTraits<complex<T>, complex<T> > {
/** The type resulting from multiplication of the template type with
itself. */
typedef complex<typename BinaryOpTraits<T,T>::MulT> MulT;
/** The type resulting from division of one template type with
another. */
typedef complex<typename BinaryOpTraits<T,T>::DivT> DivT;
/** @endcond */
};
}
#endif /* ThePEG_Complex_H */
|
