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
|
#ifndef _COMPLEX_NUMBER_H_
#define _COMPLEX_NUMBER_H_
#include <math.h>
class ComplexNumber
{
public:
ComplexNumber() ;
ComplexNumber( double r, double i = 0.0 );
double real() const;
double imag() const;
friend ComplexNumber operator*(
const ComplexNumber &, const ComplexNumber & );
friend ComplexNumber operator+(
const ComplexNumber &, const ComplexNumber & );
friend ComplexNumber operator-(
const ComplexNumber &, const ComplexNumber & );
friend ComplexNumber operator/(
const ComplexNumber &, const ComplexNumber & );
private:
double d_real;
double d_imag;
};
inline ComplexNumber::ComplexNumber():
d_real( 0.0 ),
d_imag( -0.0 )
{
}
inline ComplexNumber::ComplexNumber( double re, double im ):
d_real( re ),
d_imag( im )
{
}
inline double ComplexNumber::real() const
{
return d_real;
}
inline double ComplexNumber::imag() const
{
return d_imag;
}
inline ComplexNumber operator+(
const ComplexNumber &x1, const ComplexNumber &x2 )
{
return ComplexNumber( x1.d_real + x2.d_real, x1.d_imag + x2.d_imag );
}
inline ComplexNumber operator-(
const ComplexNumber &x1, const ComplexNumber &x2 )
{
return ComplexNumber( x1.d_real - x2.d_real, x1.d_imag - x2.d_imag );
}
inline ComplexNumber operator*(
const ComplexNumber &x1, const ComplexNumber &x2 )
{
return ComplexNumber( x1.d_real * x2.d_real - x1.d_imag * x2.d_imag,
x1.d_real * x2.d_imag + x2.d_real * x1.d_imag );
}
inline ComplexNumber operator/(
const ComplexNumber &x1, const ComplexNumber &x2 )
{
double denom = x2.d_real * x2.d_real + x2.d_imag * x2.d_imag;
return ComplexNumber(
( x1.d_real * x2.d_real + x1.d_imag * x2.d_imag ) / denom,
( x1.d_imag * x2.d_real - x2.d_imag * x1.d_real ) / denom
);
}
#endif
|
