aabb_intersect.glsl

void aabbIntersectFast(vec3 P, vec3 R, vec3 invR, vec3 aabb[2], out float tmin, out float tmax)                                                                                                                    
{                                                                               
    vec3 t1 = (aabb[0] - P) * invR;                                             
    vec3 t2 = (aabb[1] - P) * invR;                                             
    vec3 tmins = min(t1, t2);                                                   
    vec3 tmaxs = max(t1, t2);
    tmin = max(max(tmins.x, tmins.y), tmins.z);                                 
    tmax = min(min(tmaxs.x, tmaxs.y), tmaxs.z);                                 
}