diff options
Diffstat (limited to 'apps/plugins/fft/math.c')
-rw-r--r-- | apps/plugins/fft/math.c | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/apps/plugins/fft/math.c b/apps/plugins/fft/math.c new file mode 100644 index 0000000000..c00804fcc8 --- /dev/null +++ b/apps/plugins/fft/math.c | |||
@@ -0,0 +1,13 @@ | |||
1 | #include "math.h" | ||
2 | |||
3 | int64_t fsqrt64(int64_t a, unsigned int fracbits) | ||
4 | { | ||
5 | int64_t b = a/2 + (1 << fracbits); /* initial approximation */ | ||
6 | unsigned int n; | ||
7 | const unsigned int iterations = 3; /* very rough approximation */ | ||
8 | |||
9 | for (n = 0; n < iterations; ++n) | ||
10 | b = (b + (((int64_t)(a) << fracbits)/b))/2; | ||
11 | |||
12 | return b; | ||
13 | } | ||