#include <ray.h>

Public Member Functions
void	Init (const SFVEC3F &o, const SFVEC3F &d)

bool	IntersectSphere (const SFVEC3F &aCenter, float aRadius, float &aOutT0, float &aOutT1) const

SFVEC3F	at (float t) const

SFVEC2F	at2D (float t) const

Public Attributes
SFVEC3F	m_Origin

unsigned int	rayID
	unique ray ID - not used - dummy More...

SFVEC3F	m_Dir

RAY_CLASSIFICATION	m_Classification

SFVEC3F	m_InvDir

float	ibyj

float	jbyi

float	kbyj

float	jbyk

float	ibyk

float	kbyi

float	c_xy

float	c_xz

float	c_yx

float	c_yz

float	c_zx

float	c_zy

unsigned int	m_dirIsNeg [3]

Detailed Description

Definition at line 62 of file ray.h.

Member Function Documentation

◆ at()

SFVEC3F RAY::at ( float t ) const

inline

Definition at line 84 of file ray.h.

84{ return m_Origin + m_Dir * t; }

RAY::m_Dir

SFVEC3F m_Dir

Definition: ray.h:67

RAY::m_Origin

SFVEC3F m_Origin

Definition: ray.h:64

References m_Dir, and m_Origin.

Referenced by NL_3D_VIEWER_PLUGIN_IMPL::GetHitLookAt(), CYLINDER::Intersect(), DUMMY_BLOCK::Intersect(), LAYER_ITEM::Intersect(), XY_PLANE::Intersect(), ROUND_SEGMENT::Intersect(), TRIANGLE::Intersect(), LAYER_ITEM::IntersectP(), ROUND_SEGMENT::IntersectP(), EDA_3D_CANVAS::move_pivot_based_on_cur_mouse_position(), RENDER_3D_RAYTRACE::renderBlockTracing(), and RENDER_3D_RAYTRACE::shadeHit().

◆ at2D()

SFVEC2F RAY::at2D ( float t ) const

inline

Definition at line 86 of file ray.h.

 {
 return SFVEC2F( m_Origin.x + m_Dir.x * t, m_Origin.y + m_Dir.y * t );
 }

References m_Dir, and m_Origin.

◆ Init()

void RAY::Init	(	const SFVEC3F &	o,
		const SFVEC3F &	d
	)

Definition at line 35 of file ray.cpp.

{
 m_Origin = o;
 m_Dir = d;
 m_InvDir = 1.0f / d;
 
 rayID = 0; // Not used, just set to 0
 //rayID = gs_next_rayID;
 //gs_next_rayID++;
 
 // An Efficient and Robust Ray–Box Intersection Algorithm
 // Amy Williams Steve Barrus R. Keith Morley Peter Shirley
 // University of Utah
 // http://people.csail.mit.edu/amy/papers/box-jgt.pdf
 m_dirIsNeg[0] = m_Dir.x < 0.0f;
 m_dirIsNeg[1] = m_Dir.y < 0.0f;
 m_dirIsNeg[2] = m_Dir.z < 0.0f;
 
 // ray slope
 
 // "Fast Ray / Axis-Aligned Bounding Box Overlap Tests using Ray Slopes"
 // by Martin Eisemann, Thorsten Grosch, Stefan Müller and Marcus Magnor
 // Computer Graphics Lab, TU Braunschweig, Germany and
 // University of Koblenz-Landau, Germany
 // Licence: "This source code is public domain, but please mention us if you use it."
 //
 // https://github.com/rjw57/mcvoxel/tree/master/third-party/rayslope
 // https://github.com/rjw57/mcvoxel/blob/master/third-party/rayslope/ray.cpp
 
 ibyj = m_Dir.x * m_InvDir.y;
 jbyi = m_Dir.y * m_InvDir.x;
 jbyk = m_Dir.y * m_InvDir.z;
 kbyj = m_Dir.z * m_InvDir.y;
 ibyk = m_Dir.x * m_InvDir.z;
 kbyi = m_Dir.z * m_InvDir.x;
 c_xy = m_Origin.y - jbyi * m_Origin.x;
 c_xz = m_Origin.z - kbyi * m_Origin.x;
 c_yx = m_Origin.x - ibyj * m_Origin.y;
 c_yz = m_Origin.z - kbyj * m_Origin.y;
 c_zx = m_Origin.x - ibyk * m_Origin.z;
 c_zy = m_Origin.y - jbyk * m_Origin.z;
 
 // ray slope classification
 if( m_Dir.x < 0 )
 {
 if( m_Dir.y < 0 )
 {
 if( m_Dir.z < 0 )
 {
 m_Classification = RAY_CLASSIFICATION::MMM;
 }
 else if( m_Dir.z > 0 )
 {
 m_Classification = RAY_CLASSIFICATION::MMP;
 }
 else
 {
 m_Classification = RAY_CLASSIFICATION::MMO;
 }
 }
 else
 {
 if( m_Dir.z < 0 )
 {
 m_Classification = RAY_CLASSIFICATION::MPM;
 
 if( m_Dir.y == 0 )
 m_Classification = RAY_CLASSIFICATION::MOM;
 }
 else
 {
 if( ( m_Dir.y == 0 ) && ( m_Dir.z == 0 ) )
 m_Classification = RAY_CLASSIFICATION::MOO;
 else if( m_Dir.z == 0 )
 m_Classification = RAY_CLASSIFICATION::MPO;
 else if( m_Dir.y == 0 )
 m_Classification = RAY_CLASSIFICATION::MOP;
 else
 m_Classification = RAY_CLASSIFICATION::MPP;
 }
 }
 }
 else
 {
 if( m_Dir.y < 0 )
 {
 if( m_Dir.z < 0 )
 {
 m_Classification = RAY_CLASSIFICATION::PMM;
 
 if( m_Dir.x == 0 )
 m_Classification = RAY_CLASSIFICATION::OMM;
 }
 else
 {
 if( ( m_Dir.x == 0 ) && ( m_Dir.z == 0 ) )
 m_Classification = RAY_CLASSIFICATION::OMO;
 else if( m_Dir.z == 0 )
 m_Classification = RAY_CLASSIFICATION::PMO;
 else if( m_Dir.x == 0 )
 m_Classification = RAY_CLASSIFICATION::OMP;
 else
 m_Classification = RAY_CLASSIFICATION::PMP;
 }
 }
  else
 {
 if( m_Dir.z < 0 )
 {
 if( ( m_Dir.x == 0 ) && ( m_Dir.y == 0 ) )
 m_Classification = RAY_CLASSIFICATION::OOM;
 else if( m_Dir.x == 0 )
 m_Classification = RAY_CLASSIFICATION::OPM;
 else if( m_Dir.y == 0 )
 m_Classification = RAY_CLASSIFICATION::POM;
 else
 m_Classification = RAY_CLASSIFICATION::PPM;
 }
 else
 {
 if( m_Dir.x == 0 )
 {
 if( m_Dir.y == 0 )
 m_Classification = RAY_CLASSIFICATION::OOP;
 else if( m_Dir.z == 0 )
 m_Classification = RAY_CLASSIFICATION::OPO;
 else
 m_Classification = RAY_CLASSIFICATION::OPP;
 }
 else
 {
 if( ( m_Dir.y == 0 ) && ( m_Dir.z == 0 ) )
 m_Classification = RAY_CLASSIFICATION::POO;
 else if( m_Dir.y == 0 )
 m_Classification = RAY_CLASSIFICATION::POP;
 else if( m_Dir.z == 0 )
 m_Classification = RAY_CLASSIFICATION::PPO;
 else
 m_Classification = RAY_CLASSIFICATION::PPP;
 }
 }
 }
 }
}

References c_xy, c_xz, c_yx, c_yz, c_zx, c_zy, ibyj, ibyk, jbyi, jbyk, kbyi, kbyj, m_Classification, m_Dir, m_dirIsNeg, m_InvDir, m_Origin, MMM, MMO, MMP, MOM, MOO, MOP, MPM, MPO, MPP, OMM, OMO, OMP, OOM, OOP, OPM, OPO, OPP, PMM, PMO, PMP, POM, POO, POP, PPM, PPO, PPP, and rayID.

Referenced by NL_3D_VIEWER_PLUGIN_IMPL::GetHitLookAt(), EDA_3D_CANVAS::getRayAtCurrentMousePosition(), RAYPACKET::RAYPACKET(), RAYPACKET_InitRays(), RAYPACKET_InitRays_with2DDisplacement(), RENDER_3D_RAYTRACE::renderPreview(), and RENDER_3D_RAYTRACE::shadeHit().

◆ IntersectSphere()

bool RAY::IntersectSphere	(	const SFVEC3F &	aCenter,
		float	aRadius,
		float &	aOutT0,
		float &	aOutT1
	)		const

Todo:: : not tested

Definition at line 211 of file ray.cpp.

{
 // Ray-sphere intersection: geometric
 SFVEC3F OC = aCenter - m_Origin;
 float p_dot_d = glm::dot( OC, m_Dir );
 
 if( p_dot_d < 0.0f )
 return 0.0f;
 
 float p_dot_p = glm::dot( OC, OC );
 float discriminant = p_dot_p - p_dot_d * p_dot_d;
 
 if( discriminant > aRadius*aRadius )
 return false;
 
 discriminant = sqrtf( aRadius*aRadius - discriminant );
 
 aOutT0 = p_dot_d - discriminant;
 aOutT1 = p_dot_d + discriminant;
 
 if( aOutT0 > aOutT1 )
 {
 float temp = aOutT0;
 aOutT0 = aOutT1;
 aOutT1 = temp;
 }
 
 return true;
}