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
|
// Copyright 2012 John Maddock. Distributed under the Boost
// Software License, Version 1.0. (See accompanying file
// LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
#include <type_traits>
#include "test.hpp"
#include <boost/multiprecision/mpc.hpp>
#include <boost/math/constants/constants.hpp>
using boost::multiprecision::mpc_complex_100;
template <class Complex>
void test_overloads()
{
typedef typename Complex::value_type Real;
Complex ya = {5.2, 7.4};
Complex yb = {8.2, 7.3};
Real h = 0.0001;
auto I0 = (ya + yb) * h;
Complex I1 = I0 / 2 + yb * h;
//I1 = I0; // not supposed to work.
Complex z{2, 3};
typename Complex::value_type theta = 0.2;
int n = 2;
using std::sin;
Complex arg = z * sin(theta) - n * theta;
using std::exp;
Real v = 0.2;
Real cotv = 7.8;
Real cscv = 8.2;
Complex den = z + v * cscv * exp(-v * cotv);
boost::multiprecision::number<boost::multiprecision::backends::mpc_complex_backend<100> > a = 2;
boost::multiprecision::number<boost::multiprecision::backends::mpc_complex_backend<100> > b = 3;
/*
if (a <= b) {
b = a;
}*/
}
template <class F, class Real>
typename std::result_of_t<F(Real)> some_functional(F f, Real a, Real b)
{
if (a <= -boost::math::tools::max_value<Real>())
{
return f(a);
}
return f(b);
}
template <class Complex>
void test_functional()
{
typedef typename Complex::value_type Real;
auto f = [](Real x) -> Complex { Complex z(x, 3); return z; };
Real a = 0;
Real b = 1;
Complex result = some_functional(f, a, b);
}
int main()
{
test_overloads<mpc_complex_100>();
test_functional<mpc_complex_100>();
}
|
