summaryrefslogtreecommitdiff
path: root/apps/plugins/fft/math.c
diff options
context:
space:
mode:
Diffstat (limited to 'apps/plugins/fft/math.c')
-rw-r--r--apps/plugins/fft/math.c13
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
3int64_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}