bool IntrPlane3Ellipsoid3<Real>::EllipsoidIsCulled () const
{
    Matrix3<Real> MInverse;
    mEllipsoid->GetMInverse(MInverse);
    Real discr = MInverse.QForm(mPlane->Normal,mPlane->Normal);
    Real root = Math<Real>::Sqrt(Math<Real>::FAbs(discr));
    Real sDist = mPlane->DistanceTo(mEllipsoid->Center);
    return sDist <= -root;
}
bool IntrPlane3Ellipsoid3<Real>::Test ()
{
    Matrix3<Real> MInverse;
    mEllipsoid->GetMInverse(MInverse);
    Real discr = MInverse.QForm(mPlane->Normal,mPlane->Normal);
    Real root = Math<Real>::Sqrt(Math<Real>::FAbs(discr));
    Real sDist = mPlane->DistanceTo(mEllipsoid->Center);
    return Math<Real>::FAbs(sDist) <= root;
}