| 1 | /* Copyright (C) 2010-2020 The RetroArch team |
| 2 | * |
| 3 | * --------------------------------------------------------------------------------------- |
| 4 | * The following license statement only applies to this file (complex.h). |
| 5 | * --------------------------------------------------------------------------------------- |
| 6 | * |
| 7 | * Permission is hereby granted, free of charge, |
| 8 | * to any person obtaining a copy of this software and associated documentation files (the "Software"), |
| 9 | * to deal in the Software without restriction, including without limitation the rights to |
| 10 | * use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, |
| 11 | * and to permit persons to whom the Software is furnished to do so, subject to the following conditions: |
| 12 | * |
| 13 | * The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software. |
| 14 | * |
| 15 | * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, |
| 16 | * INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, |
| 17 | * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. |
| 18 | * IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, |
| 19 | * WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, |
| 20 | * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. |
| 21 | */ |
| 22 | #ifndef __LIBRETRO_SDK_MATH_COMPLEX_H__ |
| 23 | #define __LIBRETRO_SDK_MATH_COMPLEX_H__ |
| 24 | |
| 25 | #include <stdint.h> |
| 26 | |
| 27 | #include <retro_inline.h> |
| 28 | |
| 29 | typedef struct |
| 30 | { |
| 31 | float real; |
| 32 | float imag; |
| 33 | } fft_complex_t; |
| 34 | |
| 35 | static INLINE fft_complex_t fft_complex_mul(fft_complex_t a, |
| 36 | fft_complex_t b) |
| 37 | { |
| 38 | fft_complex_t out; |
| 39 | out.real = a.real * b.real - a.imag * b.imag; |
| 40 | out.imag = a.imag * b.real + a.real * b.imag; |
| 41 | |
| 42 | return out; |
| 43 | |
| 44 | } |
| 45 | |
| 46 | static INLINE fft_complex_t fft_complex_add(fft_complex_t a, |
| 47 | fft_complex_t b) |
| 48 | { |
| 49 | fft_complex_t out; |
| 50 | out.real = a.real + b.real; |
| 51 | out.imag = a.imag + b.imag; |
| 52 | |
| 53 | return out; |
| 54 | |
| 55 | } |
| 56 | |
| 57 | static INLINE fft_complex_t fft_complex_sub(fft_complex_t a, |
| 58 | fft_complex_t b) |
| 59 | { |
| 60 | fft_complex_t out; |
| 61 | out.real = a.real - b.real; |
| 62 | out.imag = a.imag - b.imag; |
| 63 | |
| 64 | return out; |
| 65 | |
| 66 | } |
| 67 | |
| 68 | static INLINE fft_complex_t fft_complex_conj(fft_complex_t a) |
| 69 | { |
| 70 | fft_complex_t out; |
| 71 | out.real = a.real; |
| 72 | out.imag = -a.imag; |
| 73 | |
| 74 | return out; |
| 75 | } |
| 76 | |
| 77 | #endif |