1 /*************************************************************************/
2 /* */
3 /* Centre for Speech Technology Research */
4 /* University of Edinburgh, UK */
5 /* Copyright (c) 1996 */
6 /* All Rights Reserved. */
7 /* */
8 /* Permission is hereby granted, free of charge, to use and distribute */
9 /* this software and its documentation without restriction, including */
10 /* without limitation the rights to use, copy, modify, merge, publish, */
11 /* distribute, sublicense, and/or sell copies of this work, and to */
12 /* permit persons to whom this work is furnished to do so, subject to */
13 /* the following conditions: */
14 /* 1. The code must retain the above copyright notice, this list of */
15 /* conditions and the following disclaimer. */
16 /* 2. Any modifications must be clearly marked as such. */
17 /* 3. Original authors' names are not deleted. */
18 /* 4. The authors' names are not used to endorse or promote products */
19 /* derived from this software without specific prior written */
20 /* permission. */
21 /* */
22 /* THE UNIVERSITY OF EDINBURGH AND THE CONTRIBUTORS TO THIS WORK */
23 /* DISCLAIM ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING */
24 /* ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS, IN NO EVENT */
25 /* SHALL THE UNIVERSITY OF EDINBURGH NOR THE CONTRIBUTORS BE LIABLE */
26 /* FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES */
27 /* WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN */
28 /* AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, */
29 /* ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF */
30 /* THIS SOFTWARE. */
31 /* */
32 /*************************************************************************/
33 /* Author : Paul Taylor */
34 /* Date : December 1997 */
35 /*-----------------------------------------------------------------------*/
36 /* */
37 /*=======================================================================*/
38
39 #include "EST_Complex.h"
40
41
42 // Addition
43
operator +(const EST_Complex & z1,const EST_Complex & z2)44 EST_Complex operator + (const EST_Complex& z1, const EST_Complex &z2)
45 {
46 return EST_Complex(z1.r + z2.r, z1.i + z2.i);
47 }
48
operator +(const EST_Complex & z,float x)49 EST_Complex operator + (const EST_Complex& z, float x)
50 {
51 return EST_Complex(z.r + x, z.i);
52 }
53
operator +(float x,const EST_Complex & z)54 EST_Complex operator + (float x, const EST_Complex &z)
55 {
56 return EST_Complex(z.r + x, z.i);
57 }
58
59 // Subtraction
60
operator -(const EST_Complex & z1,const EST_Complex & z2)61 EST_Complex operator - (const EST_Complex& z1, const EST_Complex &z2)
62 {
63 return EST_Complex(z1.r - z2.r, z1.i - z2.i);
64 }
65
operator -(const EST_Complex & z,float x)66 EST_Complex operator - (const EST_Complex& z, float x)
67 {
68 return EST_Complex(z.r -x, z.i);
69 }
70
operator -(float x,const EST_Complex & z)71 EST_Complex operator - (float x, const EST_Complex &z)
72 {
73 return EST_Complex(x - z.r, - z.i);
74 }
75
76 // Multiplication
77
operator *(const EST_Complex & z1,const EST_Complex & z2)78 EST_Complex operator * (const EST_Complex& z1, const EST_Complex &z2)
79 {
80 return EST_Complex((z1.r * z2.r) - (z1.i * z2.i), (z1.r * z2.i) + (z1.i * z2.r));
81 }
82
operator *(const EST_Complex & z,float x)83 EST_Complex operator * (const EST_Complex& z, float x)
84 {
85 return EST_Complex(z.r * x, z.i *x);
86 }
87
operator *(float x,const EST_Complex & z)88 EST_Complex operator * (float x, const EST_Complex &z)
89 {
90 return EST_Complex(z.r * x, z.i *x);
91 }
92
93 // Division
94
operator /(const EST_Complex & z1,const EST_Complex & z2)95 EST_Complex operator / (const EST_Complex& z1, const EST_Complex &z2)
96 {
97 double m = z1.mag();
98
99 EST_Complex inv(z1.r / m, z1.i /m);
100
101 return inv * z2;
102 }
103
operator /(const EST_Complex & z,float x)104 EST_Complex operator / (const EST_Complex& z, float x)
105 {
106 return EST_Complex(z.r / x, z.i / x);
107 }
108
operator /(float x,const EST_Complex & z)109 EST_Complex operator / (float x, const EST_Complex &z)
110 {
111 EST_Complex a(x, 0.0);
112 return (z / a);
113 }
114
115
116