public void Test_S2Cap_GetRectBound()
{
// Empty and full caps.
Assert.True(S2Cap.Empty.GetRectBound().IsEmpty());
Assert.True(S2Cap.Full.GetRectBound().IsFull());
// Maximum allowable error for latitudes and longitudes measured in
// degrees. (Assert2.Near isn't sufficient.)
// Cap that includes the south pole.
S2LatLngRect rect = new S2Cap(GetLatLngPoint(-45, 57),
S1Angle.FromDegrees(50)).GetRectBound();
Assert2.Near(rect.LatLo().GetDegrees(), -90, kDegreeEps);
Assert2.Near(rect.LatHi().GetDegrees(), 5, kDegreeEps);
Assert.True(rect.Lng.IsFull());
// Cap that is tangent to the north pole.
rect = new S2Cap(new S2Point(1, 0, 1).Normalize(),
S1Angle.FromRadians(S2.M_PI_4 + 1e-16)).GetRectBound();
Assert2.Near(rect.Lat.Lo, 0, kEps);
Assert2.Near(rect.Lat.Hi, S2.M_PI_2, kEps);
Assert.True(rect.Lng.IsFull());
var p = new S2Point(1, 0, 1).Normalize();
var rb1 = S1Angle.FromDegrees(45 + 5e-15);
rect = new S2Cap(p, rb1).GetRectBound();
Assert2.Near(rect.LatLo().GetDegrees(), 0, kDegreeEps);
Assert2.Near(rect.LatHi().GetDegrees(), 90, kDegreeEps);
Assert.True(rect.Lng.IsFull());
// The eastern hemisphere.
var rb2 = S1Angle.FromRadians(S2.M_PI_2 + 2e-16);
rect = new S2Cap(new S2Point(0, 1, 0), rb2).GetRectBound();
Assert2.Near(rect.LatLo().GetDegrees(), -90, kDegreeEps);
Assert2.Near(rect.LatHi().GetDegrees(), 90, kDegreeEps);
Assert.True(rect.Lng.IsFull());
// A cap centered on the equator.
rect = new S2Cap(GetLatLngPoint(0, 50), S1Angle.FromDegrees(20)).GetRectBound();
Assert2.Near(rect.LatLo().GetDegrees(), -20, kDegreeEps);
Assert2.Near(rect.LatHi().GetDegrees(), 20, kDegreeEps);
Assert2.Near(rect.LngLo().GetDegrees(), 30, kDegreeEps);
Assert2.Near(rect.LngHi().GetDegrees(), 70, kDegreeEps);
// A cap centered on the north pole.
rect = new S2Cap(GetLatLngPoint(90, 123), S1Angle.FromDegrees(10)).GetRectBound();
Assert2.Near(rect.LatLo().GetDegrees(), 80, kDegreeEps);
Assert2.Near(rect.LatHi().GetDegrees(), 90, kDegreeEps);
Assert.True(rect.Lng.IsFull());
}
