Edinburgh Speech Tools 2.4-release
EST_Complex.cc
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
44EST_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
49EST_Complex operator + (const EST_Complex& z, float x)
50{
51 return EST_Complex(z.r + x, z.i);
52}
53
54EST_Complex operator + (float x, const EST_Complex &z)
55{
56 return EST_Complex(z.r + x, z.i);
57}
58
59// Subtraction
60
61EST_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
66EST_Complex operator - (const EST_Complex& z, float x)
67{
68 return EST_Complex(z.r -x, z.i);
69}
70
71EST_Complex operator - (float x, const EST_Complex &z)
72{
73 return EST_Complex(x - z.r, - z.i);
74}
75
76// Multiplication
77
78EST_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
83EST_Complex operator * (const EST_Complex& z, float x)
84{
85 return EST_Complex(z.r * x, z.i *x);
86}
87
88EST_Complex operator * (float x, const EST_Complex &z)
89{
90 return EST_Complex(z.r * x, z.i *x);
91}
92
93// Division
94
95EST_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
104EST_Complex operator / (const EST_Complex& z, float x)
105{
106 return EST_Complex(z.r / x, z.i / x);
107}
108
109EST_Complex operator / (float x, const EST_Complex &z)
110{
111 EST_Complex a(x, 0.0);
112 return (z / a);
113}
114
115
double mag() const
Polar magnitude, read only.
Definition: EST_Complex.h:73