geof_planet

defmodule GEOF.Planet.Geometry.InterfieldCentroidsTest do
  use ExUnit.Case

  import GEOF.Planet.Geometry.InterfieldCentroids

  doctest GEOF.Planet.Geometry.InterfieldCentroids

  test "computes interfield centroids" do
    interfield_centroid_sphere = interfield_centroids(3)

    # Check that we can access centroids from indexes in any order of 3 mutually adjacent fields

    # Not mutually adjacent fields
    assert :error =
             Map.fetch(
               interfield_centroid_sphere,
               MapSet.new([
                 :north,
                 {:sxy, 0, 2, 0},
                 {:sxy, 1, 0, 0}
               ])
             )

    # Usual order
    assert {:ok, {:pos, c1_lat, c1_lon}} =
             Map.fetch(
               interfield_centroid_sphere,
               MapSet.new([
                 :north,
                 {:sxy, 0, 0, 0},
                 {:sxy, 1, 0, 0}
               ])
             )

    # Shuffled order
    assert {:ok, {:pos, c2_lat, c2_lon}} =
             Map.fetch(
               interfield_centroid_sphere,
               MapSet.new([
                 {:sxy, 0, 0, 0},
                 :north,
                 {:sxy, 1, 0, 0}
               ])
             )

    # Assert that the accessed centroids are identical
    assert c1_lat == c2_lat
    assert c1_lon == c2_lon
  end

  # TODO: is there a way to reasonably test the position values we're getting from this module?
end