test_that("constructor", {
  # Angle
  t <- degree(1)
  expect_equal(turn(turn = 1)@degree, 0)
  expect_equal(turn(turn = 2)@degree, 0)
  expect_equal(turn(turn = .5)@degree, 180)
  expect_equal(turn(turn = .5)@radian, pi)
  expect_equal(radian(radian = pi)@radian, pi)
  expect_equal(degree(degree = 180)@radian, pi)
  # Rectangle
  p_center <- ob_point(0, 0)
  p_northeast <- ob_point(2, 1)
  p_northwest <- ob_point(-2, 1)
  p_southwest <- ob_point(-2, -1)
  p_southeast <- ob_point(2, -1)
  width <- 4
  height <- 2
  r_center <- ob_rectangle(center = p_center,
                           width = width,
                           height = height)
  expect_identical(r_center,
                   ob_rectangle(center = p_center, northeast = p_northeast))
  expect_identical(r_center,
                   ob_rectangle(center = p_center, northwest = p_northwest))
  expect_identical(r_center,
                   ob_rectangle(center = p_center, southeast = p_southeast))
  expect_identical(r_center,
                   ob_rectangle(center = p_center, southwest = p_southwest))
  expect_identical(r_center,
                   ob_rectangle(
                     width = width,
                     height = height,
                     southwest = p_southwest
                   ))
  expect_identical(r_center,
                   ob_rectangle(
                     width = width,
                     height = height,
                     northwest = p_northwest
                   ))
  expect_identical(r_center,
                   ob_rectangle(
                     width = width,
                     height = height,
                     southeast = p_southeast
                   ))
  expect_identical(r_center,
                   ob_rectangle(
                     width = width,
                     height = height,
                     northeast = p_northeast
                   ))
  expect_identical(
    r_center,
    ob_rectangle(
      width = width,
      southeast = p_southeast,
      northeast = p_northeast
    )
  )
  expect_identical(
    r_center,
    ob_rectangle(
      width = width,
      southwest = p_southwest,
      northwest = p_northwest
    )
  )
  expect_identical(r_center,
                   ob_rectangle(southwest = p_southwest, northeast = p_northeast))
  expect_identical(r_center,
                   ob_rectangle(northwest = p_northwest, southeast = p_southeast))
  expect_identical(
    r_center,
    ob_rectangle(
      height = height,
      southwest = p_southwest,
      southeast = p_southeast
    )
  )
  expect_identical(
    r_center,
    ob_rectangle(
      height = height,
      northwest = p_northwest,
      northeast = p_northeast
    )
  )



})


test_that("adding", {
  p1 <- ob_point(1, 1)
  p2 <- ob_point(3, 4)
  p3 <- ob_point(4, 5)
  s1 <- ob_segment(p1, p2)
  s2 <- ob_segment(p2, p3)

  expect_identical(p1 + p2, p3)
  expect_identical(p3 - p2, p1)
  expect_identical(ob_segment(p1, p2) + p3, ob_segment(p1 + p3, p2 + p3))
  expect_identical(p3 + ob_segment(p1, p2), ob_segment(p1 + p3, p2 + p3))
  expect_equal(p3 - ob_segment(p1, p2), ob_segment(p3 - p1, p3 - p2))
  expect_identical(ob_circle(p1, 2) + p2, ob_circle(p1 + p2, 2))
  expect_identical(ob_circle(p1, 2) - p2, ob_circle(p1 - p2, 2))
  expect_identical(p2 + ob_circle(p1, 2), ob_circle(p1 + p2, 2))
  expect_identical(p2 - ob_circle(p1, 2), ob_circle(p2 - p1, 2))
  expect_identical(p2 + 2, ob_point(5, 6))
  expect_identical(2 + p2, ob_point(5, 6))
  expect_identical(p2 - 2, ob_point(1, 2))
  expect_identical(2 - p2, ob_point(-1, -2))
  expect_identical(s1 + s2, nudge(s1, s2))
})






test_that("distances", {
  p0 <- ob_point(0, 3)
  p1 <- ob_point(1, 1)
  p2 <- ob_point(4, 5)
  n1 <- ob_label(p1, label = "A")
  n2 <- ob_label(p2, label = "B")
  s1 <- ob_segment(p1, p2)
  c1 <- ob_circle(p1, radius = 1)
  c2 <- ob_circle(p2, radius = 2)
  expect_equal(distance(p0), 3)
  expect_equal(distance(p1, p2), 5)
  expect_equal(distance(s1), 5)
  expect_equal(distance(c1, c2), 2)
  expect_equal(distance(c1, c1), 0)
})






test_that("intersection", {
  s1 <- ob_segment(ob_point(0, 1), ob_point(1, 0))
  s2 <- ob_segment(ob_point(0, 0), ob_point(1, 1))
  s3 <- ob_segment(ob_point(0.1, .1), ob_point(.9, .9))
  s4 <- ob_segment(ob_point(1, 1), ob_point(2, 2))
  p1 <- ob_point(.5, .5)
  expect_equal(intersection(s1, s2), p1)
  expect_equal(intersection(s1@line, s2), p1)
  expect_equal(intersection(s1, s2@line), p1)
  expect_equal(intersection(s1, s3), p1)
  expect_equal(intersection(s1, s4), list())
  expect_equal(intersection(s1@line, s4@line), p1)
  c1 <- ob_circle(center = ob_point(1, 1), radius = 1)
  l1 <- ob_line(intercept = 1)
  p1 <- ob_point(0, 1)
  p2 <- ob_point(2, 1)
  expect_equal(intersection(l1, c1), bind(c(p2, p1)))
  expect_equal(intersection(c1, l1), bind(c(p2, p1)))
  expect_equal(intersection(ob_line(intercept = 2), c1), ob_point(1, 2))
  l1 <- ob_line(xintercept = 1)
  p1 <- ob_point(1, 0)
  p2 <- ob_point(1, 2)
  expect_equal(intersection(l1, c1), bind(c(p1, p2)))
  expect_equal(intersection(c1, l1), bind(c(p1, p2)))
  expect_equal(intersection(ob_line(xintercept = 2), c1), ob_point(2, 1))
  l1 <- ob_line(slope = 1, intercept = 2 * sin(degree(45)))
  c1 <- ob_circle(ob_point(0, 0), radius = 1)
  # intersect at tangent
  expect_equal(intersection(l1, c1), S7::convert(ob_polar(
    theta = radian(radian = pi * 3 / 4), r = 1
  ), ob_point))
  l1 <- ob_line(slope = .5, intercept = 0)
  c1 <- ob_circle(ob_point(0, 0), radius = 1)
  expect_equal(intersection(l1, c1), bind(c(
    ob_point(x = cos(atan(.5)), y = sin(atan(.5))), ob_point(x = -cos(atan(.5)), y = -sin(atan(.5)))
  )))
  e1 <- ob_ellipse(a = 1, b = 2)

  s5 <- ob_segment(ob_point(-2, -2), ob_point(2, 2))
  # intersection(s2, e1)
  # intersection(s5, e1)
  #   intersection(ob_segment(ob_point(0,0), ob_point(2,0)), ob_segment(ob_point(0,1), ob_point(2,1)))
  #
  #   intersection(ob_segment(ob_point(1,0), ob_point(1,2)), ob_segment(ob_point(0,1), ob_point(2,1)))
  #   intersection(ob_segment(ob_point(1,0), ob_point(1,2)), ob_segment(ob_point(1,1), ob_point(1,2)))
  #   intersection(ob_line(1,intercept = 0), ob_rectangle(ob_point(0,0), width = 2, height = 2))
  #   intersection(ob_segment(ob_point(0,-3), ob_point(0,10)), ob_rectangle(center = ob_point(0,0), width = 2, height = 2))
  # x <- ob_segment(ob_point(1,0), ob_point(1,2))
  # y <- ob_segment(ob_point(1,1), ob_point(1,2))
})




test_that("rotate", {
  # rotate a line with an angle
  expect_identical(rotate(ob_line(xintercept = 2), turn(turn = .5)), ob_line(xintercept = -2))
  # rotate a line with a numeric radian
  expect_equal(rotate(ob_line(xintercept = 2), turn(turn = .5)), rotate(ob_line(xintercept = 2), pi))

  # rotate a ob_point
  expect_equal(rotate(ob_point(1, 0), turn(turn = .5)), ob_point(-1, 0))

  # rotate a segment
  expect_equal(rotate(ob_segment(ob_point(0, 1), ob_point(1, 0)), theta = turn(.5)),
               ob_segment(ob_point(0, -1), ob_point(-1, 0)))

  # rotate a circle
  expect_equal(rotate(x = ob_circle(ob_point(1, 2)), theta = turn(.25)), ob_circle(ob_point(-2, 1)))

  expect_equal(rotate(x = ob_circle(ob_point(1, 2), n = 50), theta = turn(.25)),
               ob_circle(ob_point(-2, 1), n = 50))

  # rotate an ellipse
  expect_equal(rotate(
    x = ob_ellipse(
      center = ob_point(1, 2),
      a = 2,
      b = 1
    ),
    theta = turn(.25)
  ),
  ob_ellipse(
    ob_point(-2, 1),
    a = 2,
    b = 1,
    angle = turn(.25)
  ))
})

test_that("resect", {
  s1 <- ob_segment(ob_point(0, 0), ob_point(1, 0))
  expect_equal(resect(s1, distance(s1) * .05), ob_segment(ob_point(.05, 0), ob_point(.95, 0)))
})

test_that("inside", {
  o <- ob_point(0, 0)
  r1 <- ob_rectangle(center = o,
                     width = 2,
                     height = 4)
  expect_equal(inside(o, r1), 1)
  expect_equal(inside(ob_point(0, 2), r1), 0)
  expect_equal(inside(ob_point(0, 1.5), r1), 1)
  expect_equal(inside(ob_point(1.5, 1.5), r1), -1)
  expect_equal(inside(ob_point(.5, 2.5), r1), -1)

  c1 <- ob_circle(ob_point(1, 1), radius = sqrt(2))
  expect_equal(inside(ob_point(1, 1), c1), 1)
  expect_equal(inside(o, c1), 0)
  expect_equal(inside(ob_point(1, 0), c1), 1)
  expect_equal(inside(ob_point(1, 2), c1), 1)
  expect_equal(inside(ob_point(2, 1), c1), 1)
  expect_equal(inside(ob_point(2, 2), c1), 0)
  expect_equal(inside(ob_point(2, 3), c1), -1)

  e1 <- ob_ellipse(o, a = 1, b = 2)
  expect_equal(inside(o, e1), 1)
  re1 <- ob_ellipse(o,
                    a = 1,
                    b = 2,
                    angle = degree(45))

  p1 <- ob_point(0, 1.9)
  rp1 <- rotate(p1, degree(45))

  expect_equal(inside(p1, e1), 1)
  expect_equal(inside(p1, re1), -1)
  expect_equal(inside(rp1, e1), -1)
  expect_equal(inside(rp1, re1), 1)


})