euler/test.R at main · tobiste/euler · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
library(euler)
library(dplyr)

r1 <- c(33, 15, -1) %>% to_euler()
r2 <- c(65, -100, 2) %>% to_euler()
# expected results: c(42.35960, 10.60413, 10.85098)

euler_perf <- microbenchmark::microbenchmark(
  relative_euler_schaeben(r1, r2),
  relative_euler_schaeben2(r1, r2),
  relative_euler_schaeben3(r1, r2),
  relative_euler_py_schaeben(r1, r2),
  relative_euler_py2_schaeben(r1, r2),
  relative_euler_greiner(r1, r2),
  relative_euler_lepichon(r1, r2),
  times = 200
)
microbenchmark:::boxplot.microbenchmark(euler_perf)
summary(euler_perf) %>% as_tibble() %>% arrange(median) #%>% write.csv(file = "performance.csv")

lepichon <- relative_euler_lepichon(r1, r2)
greiner <- relative_euler_greiner(r1, r2)
schaeben2 <- relative_euler_schaeben2(r1, r2)
schaeben3 <- relative_euler_schaeben3(r1, r2)
python <- relative_euler_py_schaeben(r1, r2)
python2 <- relative_euler_py2_schaeben(r1, r2)
schaeben <- relative_euler_schaeben(r1, r2)

schaeben$angle == schaeben2$angle
schaeben$angle == python$angle
schaeben$angle == python2$angle
schaeben$angle == greiner$angle
schaeben$angle == lepichon$angle

sprintf("%.12f", schaeben$angle)
sprintf("%.12f", greiner$angle)

sprintf("%.1f", 1000 * schaeben$angle * 2 * pi * 63781e6 / 360) # mm
sprintf("%.1f", 1000 * greiner$angle * 2 * pi * 63781e6 / 360) # mm
# 1 mm  difference!