public void EllipsoidAABBTest()
{
Point3d p = new Point3d(0, 0, 0);
Vector3d v1 = new Vector3d(3, 0, 0);
Vector3d v2 = new Vector3d(0, 2, 0);
Vector3d v3 = new Vector3d(0, 0, 4);
Ellipsoid e = new Ellipsoid(p, v1, v2, v3);
Rotation r = new Rotation(new Vector3d(1, 2, 3), PI / 4);
e = e.Rotate(r, p);
Box3d b = e.BoundingBox();
Plane3d s = new Plane3d(b.P1, b.P2, b.P3);
Assert.IsTrue(e.IntersectionWith(s).GetType() == typeof(Point3d));
s = new Plane3d(b.P8, b.P7, b.P6);
Assert.IsTrue(e.IntersectionWith(s).GetType() == typeof(Point3d));
s = new Plane3d(b.P1, b.P2, b.P5);
Assert.IsTrue(e.IntersectionWith(s).GetType() == typeof(Point3d));
s = new Plane3d(b.P2, b.P3, b.P6);
Assert.IsTrue(e.IntersectionWith(s).GetType() == typeof(Point3d));
s = new Plane3d(b.P3, b.P4, b.P7);
Assert.IsTrue(e.IntersectionWith(s).GetType() == typeof(Point3d));
s = new Plane3d(b.P1, b.P4, b.P8);
Assert.IsTrue(e.IntersectionWith(s).GetType() == typeof(Point3d));
}