1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
use crate::math::{Isometry, Real, Vector};
use crate::query::gjk::{self, CSOPoint, GJKResult, VoronoiSimplex};
use crate::query::ClosestPoints;
use crate::shape::SupportMap;
use na::Unit;
pub fn closest_points_support_map_support_map<G1: ?Sized, G2: ?Sized>(
pos12: &Isometry<Real>,
g1: &G1,
g2: &G2,
prediction: Real,
) -> ClosestPoints
where
G1: SupportMap,
G2: SupportMap,
{
match closest_points_support_map_support_map_with_params(
pos12,
g1,
g2,
prediction,
&mut VoronoiSimplex::new(),
None,
) {
GJKResult::ClosestPoints(pt1, pt2, _) => {
ClosestPoints::WithinMargin(pt1, pos12.inverse_transform_point(&pt2))
}
GJKResult::NoIntersection(_) => ClosestPoints::Disjoint,
GJKResult::Intersection => ClosestPoints::Intersecting,
GJKResult::Proximity(_) => unreachable!(),
}
}
pub fn closest_points_support_map_support_map_with_params<G1: ?Sized, G2: ?Sized>(
pos12: &Isometry<Real>,
g1: &G1,
g2: &G2,
prediction: Real,
simplex: &mut VoronoiSimplex,
init_dir: Option<Vector<Real>>,
) -> GJKResult
where
G1: SupportMap,
G2: SupportMap,
{
let dir = match init_dir {
None => -pos12.translation.vector,
Some(dir) => dir,
};
if let Some(dir) = Unit::try_new(dir, crate::math::DEFAULT_EPSILON) {
simplex.reset(CSOPoint::from_shapes(pos12, g1, g2, &dir));
} else {
simplex.reset(CSOPoint::from_shapes(
pos12,
g1,
g2,
&Vector::<Real>::x_axis(),
));
}
gjk::closest_points(pos12, g1, g2, prediction, true, simplex)
}