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
|
/*
* lingot, a musical instrument tuner.
*
* Copyright (C) 2004-2018 Iban Cereijo.
* Copyright (C) 2004-2008 Jairo Chapela.
*
* This file is part of lingot.
*
* lingot is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* lingot is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with lingot; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*/
#include "lingot-complex.h"
void lingot_complex_add(const LingotComplex a, const LingotComplex b,
LingotComplex c) {
c[0] = a[0] + b[0];
c[1] = a[1] + b[1];
}
void lingot_complex_sub(const LingotComplex a, const LingotComplex b,
LingotComplex c) {
c[0] = a[0] - b[0];
c[1] = a[1] - b[1];
}
void lingot_complex_mul(const LingotComplex a, const LingotComplex b,
LingotComplex c) {
c[0] = a[0] * b[0] - a[1] * b[1];
c[1] = a[1] * b[0] + a[0] * b[1];
}
void lingot_complex_div(const LingotComplex a, const LingotComplex b,
LingotComplex c) {
FLT bm2 = b[0] * b[0] + b[1] * b[1];
c[0] = (a[0] * b[0] + a[1] * b[1]) / bm2;
c[1] = (a[1] * b[0] - a[0] * b[1]) / bm2;
}
void lingot_complex_mul_by(LingotComplex a, const LingotComplex b) {
double rr = a[0] * b[0] - a[1] * b[1];
a[1] = a[1] * b[0] + a[0] * b[1];
a[0] = rr;
}
void lingot_complex_div_by(LingotComplex a, const LingotComplex b) {
FLT bm2 = b[0] * b[0] + b[1] * b[1];
double rr = (a[0] * b[0] + a[1] * b[1]) / bm2;
a[1] = (a[1] * b[0] - a[0] * b[1]) / bm2;
a[0] = rr;
}
|