From 9477dbe667f250ecd23f8fc0d56b942191526421 Mon Sep 17 00:00:00 2001
From: Franciszek Malinka <franciszek.malinka@gmail.com>
Date: Thu, 25 Feb 2021 14:42:55 +0100
Subject: Stare semestry, niepoukladane

---
 Semestr 3/anm/cordic/taylor.jl | 103 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 103 insertions(+)
 create mode 100644 Semestr 3/anm/cordic/taylor.jl

(limited to 'Semestr 3/anm/cordic/taylor.jl')

diff --git a/Semestr 3/anm/cordic/taylor.jl b/Semestr 3/anm/cordic/taylor.jl
new file mode 100644
index 0000000..5d2d26d
--- /dev/null
+++ b/Semestr 3/anm/cordic/taylor.jl	
@@ -0,0 +1,103 @@
+module taylor
+using Base
+using Printf
+ 
+ITER = 12
+HYPERBOLIC_MAX = 0.2
+ 
+function series(x, parity, change_sign, iterations)
+    elements = ones(Float64, 2 * iterations)
+    res = 0.0
+    i = 2
+    while i <= 2 * iterations + parity - 2
+        elements[i + 1] = elements[i] / i
+        if change_sign && (i % 2 == parity)
+            elements[i + 1] = -elements[i + 1]
+        end -
+        i += 1
+    end
+    i = 2 * iterations + parity - 2
+    while i >= parity
+        res *= x * x
+        res += elements[i + 1]
+        i -= 2
+    end
+    if parity == 1
+        res *= x
+    end
+    return res
+end
+ 
+function real_sin(r, iterations)
+    r = r - floor(r / (2 * pi)) * 2 * pi
+    if r > pi
+        return -real_sin(r - pi, iterations)
+    end
+    if r > pi / 2
+        return real_cos(r - pi / 2, iterations)
+    end
+    if r > pi / 4
+        return real_cos(pi / 2 - r, iterations)
+    end
+ 
+    return series(r, 1, true, iterations)
+end
+ 
+function real_cos(r, iterations)
+    r = r - floor(r / (2 * pi)) * 2 * pi
+    if r > pi
+        return -real_cos(r - pi, iterations)
+    end
+    if r > pi / 2
+        return -real_sin(r - pi / 2, iterations)
+    end
+    if r > pi / 4
+        return real_sin(pi / 2 - r, iterations)
+    end
+ 
+    return series(r, 0, true, iterations)
+end
+ 
+function real_sinh(r, iterations)
+    if abs(r) > HYPERBOLIC_MAX
+        return 2 * real_sinh(r / 2, iterations) * real_cosh(r / 2, iterations)
+    end
+    return series(r, 1, false, iterations)
+end
+ 
+function real_cosh(r, iterations)
+    if abs(r) > HYPERBOLIC_MAX
+        s = real_sinh(r / 2, iterations)
+        c = real_cosh(r / 2, iterations)
+        return s * s + c * c
+    end
+    return series(r, 0, false, iterations)
+end
+ 
+function complex_sin(a, b, iterations)
+    return (real_sin(a, iterations) * real_cosh(b, iterations), 
+            real_cos(a, iterations) * real_sinh(b, iterations)) 
+end
+ 
+function complex_cos(a, b, iterations)
+    return (real_cos(a, iterations) * real_cosh(b, iterations),
+            -real_sin(a, iterations) * real_sinh(b, iterations))
+end
+ 
+# c = a + bi
+function csin(a, b)
+    return complex_sin(a, b, ITER)
+end
+ 
+function ccos(a, b)
+    return complex_cos(a, b, ITER)
+end
+ 
+function rsinh(r)
+    return real_sinh(r, ITER)
+end
+ 
+function rcosh(r)
+    return real_cosh(r, ITER)
+end
+end
\ No newline at end of file
-- 
cgit v1.2.3