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
83
84
85
86
87
|
// now #include the definition of vcl_complex<>
// and define the doublify() function.
#include <vcl_complex.h>
vcl_complex<double> doublify(vcl_complex<float> const &z)
{
return vcl_complex<double>(vcl_real(z), vcl_imag(z));
}
#include <vcl_iostream.h>
int test_complex_main(int /*argc*/,char* /*argv*/[])
{
vcl_complex<double> dc1(1.1,1.2), dc2(2.1,2.2);
vcl_complex<float> fc1(1.1f,1.2f), fc2(2.1f,2.2f);
vcl_cout << dc1 << " + " << dc2 << " = " << (dc1+dc2) << vcl_endl
<< fc1 << " + " << fc2 << " = " << (fc1+fc2) << vcl_endl;
vcl_complex<double> dc3(vcl_real(dc1),vcl_imag(dc2));
vcl_complex<float> fc3(vcl_real(fc1),vcl_imag(fc2));
vcl_cout << dc3 << " / " << dc1 << " = " << dc3/dc1 << vcl_endl
<< fc3 << " / " << fc1 << " = " << fc3/fc1 << vcl_endl;
vcl_cout << "polar representation of " << dc3 << " is [" << vcl_abs(dc3) << ',' << vcl_arg(dc3) << "]\n"
<< "going back: " << dc3 << " must be = " << vcl_polar(vcl_abs(dc3), vcl_arg(dc3)) << vcl_endl;
vcl_complex<float> fcsr3 = vcl_sqrt(fc3);
vcl_cout << "sqrt(" << fc3 << ") is " << fcsr3 << ", so " << fcsr3 << " * " << fcsr3 << " = " << fcsr3*fcsr3 << vcl_endl;
// Should also test form of complex stream input and output. The standard says:
// [26.2.6.12] operator>> understands "u", "(u)" and "(u,v)" where u, v are real.
// [26.2.6.13] operator<< does f << '(' << x.real() << ',' << x.imag() << ')';
// In particular, complex numbers written with operator<< can be read again with
// operator>>.
// complex should have a type called value_type;
vcl_complex<float>::value_type tmp = 1.0f;
tmp += 2.0f; // to avoid unused variable warnings.
// Test the vcl_pow functions
bool success = true;
const vcl_complex<double> neg1(-1.0, 0.0);
const vcl_complex<double> i(0.0,1.0);
vcl_complex<double> sqrt_neg1 = vcl_pow(neg1, 0.5);
vcl_cout << "pow("<<neg1<<",0.5) = "<<sqrt_neg1<<
" and should be (0,1)"<<vcl_endl;
double error = vcl_abs(sqrt_neg1-i);
// need to be careful of quiet NANs
if ( error < 0.0 || 1e-6 < error)
{
vcl_cout << "** FAILURE **\n";
success = false;
}
const vcl_complex<double> half(0.5,0.0);
sqrt_neg1 = vcl_pow(neg1, half);
vcl_cout << "pow("<<neg1<<","<<half<<") = "<<sqrt_neg1<<
" and should be (0,1)"<<vcl_endl;
error = vcl_abs(sqrt_neg1-i);
if ( error < 0.0 || 1e-6 < error)
{
vcl_cout << "** FAILURE **\n";
success = false;
}
vcl_complex<double> zero(0.0,0.0);
vcl_cout << "Implementation defines vcl_pow((0,0),(0,0)) = "
<< vcl_pow(zero, zero) << vcl_endl;
{
vcl_complex<double> x(2, 3);
vcl_complex<double> xc = vcl_conj(x);
vcl_cout << "Conjugate " << x << " = " << xc << "\n";
if( xc != vcl_complex<double>(2,-3) ) {
success = false;
}
}
return success?0:1;
}
|
