complexnumber.h
1.79 KB
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