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
|
// Copyright 2009-2020 Intel Corporation
// SPDX-License-Identifier: Apache-2.0
#pragma once
#include "../common/ray.h"
namespace embree
{
namespace isa
{
struct HalfPlane
{
const Vec3fa P; //!< plane origin
const Vec3fa N; //!< plane normal
__forceinline HalfPlane(const Vec3fa& P, const Vec3fa& N)
: P(P), N(N) {}
__forceinline BBox1f intersect(const Vec3fa& ray_org, const Vec3fa& ray_dir) const
{
Vec3fa O = Vec3fa(ray_org) - P;
Vec3fa D = Vec3fa(ray_dir);
float ON = dot(O,N);
float DN = dot(D,N);
bool eps = abs(DN) < min_rcp_input;
float t = -ON*rcp(DN);
float lower = select(eps || DN < 0.0f, float(neg_inf), t);
float upper = select(eps || DN > 0.0f, float(pos_inf), t);
return BBox1f(lower,upper);
}
};
template<int M>
struct HalfPlaneN
{
const Vec3vf<M> P; //!< plane origin
const Vec3vf<M> N; //!< plane normal
__forceinline HalfPlaneN(const Vec3vf<M>& P, const Vec3vf<M>& N)
: P(P), N(N) {}
__forceinline BBox<vfloat<M>> intersect(const Vec3fa& ray_org, const Vec3fa& ray_dir) const
{
Vec3vf<M> O = Vec3vf<M>((Vec3fa)ray_org) - P;
Vec3vf<M> D = Vec3vf<M>((Vec3fa)ray_dir);
vfloat<M> ON = dot(O,N);
vfloat<M> DN = dot(D,N);
vbool<M> eps = abs(DN) < min_rcp_input;
vfloat<M> t = -ON*rcp(DN);
vfloat<M> lower = select(eps | DN < 0.0f, vfloat<M>(neg_inf), t);
vfloat<M> upper = select(eps | DN > 0.0f, vfloat<M>(pos_inf), t);
return BBox<vfloat<M>>(lower,upper);
}
};
}
}
|