lisyarus/psemek
1
0
Fork 1
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Add segment-box intersection

Browse source
This commit is contained in:
Nikita Lisitsa 2021-02-17 23:11:42 +03:00
parent dbeb613740
commit db2fb3095d
1 changed files with 24 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

24
libs/geom/include/psemek/geom/intersection.hpp
View file

@ -130,6 +130,30 @@ namespace psemek::geom
return !intersection(r, b).empty();
}
template <typename T, std::size_t N>
interval<T> intersection(segment<point<T, N>> const & s, box<T, N> const & b)
{
auto t = interval<T>{T(0), T(1)};
for (std::size_t i = 0; i < N; ++i)
{
T tmin = (b[i].min - s[0][i]) / (s[1][i] - s[0][i]);
T tmax = (b[i].max - s[0][i]) / (s[1][i] - s[0][i]);
if (tmin > tmax) std::swap(tmin, tmax);
t &= interval<T>{tmin, tmax};
}
return t;
}
template <typename T, std::size_t N>
bool intersect(segment<point<T, N>> const & s, box<T, N> const & b)
{
return !intersection(s, b).empty();
}
template <typename T, std::size_t N>
std::optional<T> intersection(ray<T, N> const & r, simplex<point<T, N>, N - 1> const & s)
{