bool Intersect::RayBounds (const Vector3f& origin, const Vector3f& dir, const Bounds& bounds) { const Vector3f& max = bounds.GetMax(); const Vector3f& center = bounds.GetCenter(); Vector3f ext (max - center); Vector3f diff (origin - center); if (Float::Abs(diff.x) > ext.x && Float::IsProductPositive(diff.x, dir.x)) return false; if (Float::Abs(diff.y) > ext.y && Float::IsProductPositive(diff.y, dir.y)) return false; if (Float::Abs(diff.z) > ext.z && Float::IsProductPositive(diff.z, dir.z)) return false; Vector3f adir (Float::Abs(dir.x), Float::Abs(dir.y), Float::Abs(dir.z)); float a = dir.y * diff.z - dir.z * diff.y; float b = ext.y * adir.z + ext.z * adir.y; if (Float::Abs(a) > b) return false; a = dir.z * diff.x - dir.x * diff.z; b = ext.x * adir.z + ext.z * adir.x; if (Float::Abs(a) > b) return false; a = dir.x * diff.y - dir.y * diff.x; b = ext.x * adir.y + ext.y * adir.x; return !(Float::Abs(a) > b); }