1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586 |
- // David Eberly, Geometric Tools, Redmond WA 98052
- // Copyright (c) 1998-2020
- // Distributed under the Boost Software License, Version 1.0.
- // https://www.boost.org/LICENSE_1_0.txt
- // https://www.geometrictools.com/License/Boost/LICENSE_1_0.txt
- // Version: 4.0.2019.08.13
- #pragma once
- #include <Mathematics/Cone.h>
- #include <Mathematics/OrientedBox.h>
- #include <Mathematics/IntrAlignedBox3Cone3.h>
- // Test for intersection of a box and a cone. The cone can be infinite
- // 0 <= minHeight < maxHeight = std::numeric_limits<Real>::max()
- // or finite (cone frustum)
- // 0 <= minHeight < maxHeight < std::numeric_limits<Real>::max().
- // The algorithm is described in
- // https://www.geometrictools.com/Documentation/IntersectionBoxCone.pdf
- // and reports an intersection only when th intersection set has positive
- // volume. For example, let the box be outside the cone. If the box is
- // below the minHeight plane at the cone vertex and just touches the cone
- // vertex, no intersection is reported. If the box is above the maxHeight
- // plane and just touches the disk capping the cone, either at a single
- // point, a line segment of points or a polygon of points, no intersection
- // is reported.
- // TODO: These queries were designed when an infinite cone was defined
- // by choosing maxHeight of std::numeric_limits<Real>::max(). The Cone<N,Real>
- // class has been redesigned not to use std::numeric_limits to allow for
- // arithmetic systems that do not have representations for infinities
- // (such as BSNumber and BSRational). The intersection queries need to be
- // rewritten for the new class design. FOR NOW, the queries will work with
- // float/double when you create a cone using the cone-frustum constructor
- // Cone(ray, angle, minHeight, std::numeric_limits<Real>::max()).
- namespace WwiseGTE
- {
- template <typename Real>
- class TIQuery<Real, OrientedBox<3, Real>, Cone<3, Real>>
- :
- public TIQuery<Real, AlignedBox<3, Real>, Cone<3, Real>>
- {
- public:
- struct Result
- :
- public TIQuery<Real, AlignedBox<3, Real>, Cone<3, Real>>::Result
- {
- // No additional information to compute.
- };
- Result operator()(OrientedBox<3, Real> const& box, Cone<3, Real> const& cone)
- {
- // Transform the cone and box so that the cone vertex is at the
- // origin and the box is axis aligned. This allows us to call the
- // base class operator()(...).
- Vector<3, Real> diff = box.center - cone.ray.origin;
- Vector<3, Real> xfrmBoxCenter
- {
- Dot(box.axis[0], diff),
- Dot(box.axis[1], diff),
- Dot(box.axis[2], diff)
- };
- AlignedBox<3, Real> xfrmBox;
- xfrmBox.min = xfrmBoxCenter - box.extent;
- xfrmBox.max = xfrmBoxCenter + box.extent;
- Cone<3, Real> xfrmCone = cone;
- for (int i = 0; i < 3; ++i)
- {
- xfrmCone.ray.origin[i] = (Real)0;
- xfrmCone.ray.direction[i] = Dot(box.axis[i], cone.ray.direction);
- }
- // Test for intersection between the aligned box and the cone.
- auto bcResult = TIQuery<Real, AlignedBox<3, Real>, Cone<3, Real>>::operator()(xfrmBox, xfrmCone);
- Result result;
- result.intersect = bcResult.intersect;
- return result;
- }
- };
- // Template alias for convenience.
- template <typename Real>
- using TIOrientedBox3Cone3 = TIQuery<Real, OrientedBox<3, Real>, Cone<3, Real>>;
- }
|