aboutsummaryrefslogtreecommitdiff
path: root/semestr-3/anm/numerki/lista1
diff options
context:
space:
mode:
authorFranciszek Malinka <franciszek.malinka@gmail.com>2021-10-05 21:49:54 +0200
committerFranciszek Malinka <franciszek.malinka@gmail.com>2021-10-05 21:49:54 +0200
commitc5fcf7179a83ef65c86c6a4a390029149e518649 (patch)
treed29ffc5b86a0d257453cedcf87d91a13d8bf3b0d /semestr-3/anm/numerki/lista1
parentf8a88b6a4aba1f66d04711a9330eaba49a50c463 (diff)
Duzy commit ze smieciami
Diffstat (limited to 'semestr-3/anm/numerki/lista1')
-rw-r--r--semestr-3/anm/numerki/lista1/test.jl118
-rw-r--r--semestr-3/anm/numerki/lista1/test.jl.txt16
2 files changed, 134 insertions, 0 deletions
diff --git a/semestr-3/anm/numerki/lista1/test.jl b/semestr-3/anm/numerki/lista1/test.jl
new file mode 100644
index 0000000..235800a
--- /dev/null
+++ b/semestr-3/anm/numerki/lista1/test.jl
@@ -0,0 +1,118 @@
+using Random
+
+function prettyfloat64(f)
+ ss = bitstring(f)
+ s = ss[1]
+ c = ss[2:12]
+ m = ss[13:end]
+ return string(s, " ", c, " ", m)
+end
+
+function string_to_float64(st)
+ s = st[1]
+ c = st[2:12]
+ m = st[13:end]
+
+ if s == '0'
+ s = 1.0
+ else
+ s = -1.0
+ end
+
+ c = parse(Int, string("0b", c))
+ c = 2.0^(c - 1023)
+ m = Float64(parse(Int64, string("0b", m))) / Float64(2^52) + 1.0
+ return s * c * m
+end
+
+rng = MersenneTwister(123)
+
+function rand_f64_bitstring()
+ return randstring(rng, ['0', '1'], 64)
+end
+
+function test_string_to_float()
+ for i in 1:100000
+ # st = rand_f64_bitstring()
+ st = rand(rng)
+ res = string_to_float64(bitstring(st))
+ println(string(st, " -> ", prettyfloat64(st), " -> ", res))
+ if bitstring(res) == bitstring(st)
+ println("OK")
+ else
+ println("ZJEBALO SIE")
+ break
+ end
+ end
+end
+
+function find_felerny_x()
+ while true
+ x = rand(rng) + 1.0
+ res = x * (1.0 / x)
+ if res != 1.0
+ println(string("Znalazlem! ", x))
+ return x
+ end
+ end
+end
+
+function find_felerne_iksy(iterations)
+ bad_xs = 0
+ for i in 1:iterations
+ x = rand(rng) + 1.0
+ res = x * (1.0 / x)
+ if res != 1.0
+ bad_xs += 1
+ end
+ end
+ println("Znalazlem ", bad_xs, " felernych iksow, czyli ", bad_xs / iterations * 100.0, "% wszystkich losowanych wartosci.")
+end
+
+# find_felerne_iksy(10000000)
+
+
+function rel_error(x, rx)
+ x = Float64(x)
+ rx = Float64(rx)
+ relative = abs((x - rx) / x)
+ return relative
+end
+
+function test_poly(t, poly, x=4.71, exact=-14.636489)
+ println("Testing poly \"", poly, "\" with float type ", t, ".")
+ res = poly(x, t)
+ println(poly, "(", x, ") = ", res)
+ println("Relative error (exact: ", exact, "): ", rel_error(exact, res))
+end
+
+function w(x, t)
+ x = t(x)
+ return x^3 - (t(6)) * (x^2) + (t(3)) * x - t(0.149)
+end
+
+function w2(x, t)
+ x = t(x)
+ return ((x - t(6)) * x + t(3)) * x - t(0.149)
+end
+
+
+function task3()
+ println("Task3: ")
+
+ test_poly(Float16, w)
+ println()
+ test_poly(Float32, w)
+ println()
+ test_poly(Float64, w)
+ println()
+
+ test_poly(Float16, w2)
+ println()
+ test_poly(Float32, w2)
+ println()
+ test_poly(Float64, w2)
+ println()
+end
+
+
diff --git a/semestr-3/anm/numerki/lista1/test.jl.txt b/semestr-3/anm/numerki/lista1/test.jl.txt
new file mode 100644
index 0000000..81b4901
--- /dev/null
+++ b/semestr-3/anm/numerki/lista1/test.jl.txt
@@ -0,0 +1,16 @@
+function string_to_float64(st)
+ s = st[1]
+ c = st[2:12]
+ m = st[13:end]
+
+ if s == '0'
+ s = 1.0
+ else
+ s = -1.0
+ end
+
+ c = parse(Int, string("0b", c))
+ c = 2.0 ^ (c - 1023)
+ m = Float64(parse(Int64, string("0b", m))) / Float64(2^52) + 1.0
+ return s * c * m
+end