Function RayTriangle:Byte(ray:Ray, tri:vec3[], frontfacing_only:Byte=True)
   Local vec_dst:vec3=vec3.Create(ray.o.v[0]-tri[0].v[0], ray.o.v[1]-tri[0].v[1], ray.o.v[2]-tri[0].v[2])
   Local vec_r1:vec3=vec3.Create(tri[1].v[0]-tri[0].v[0], tri[1].v[1]-tri[0].v[1], tri[1].v[2]-tri[0].v[2])
   Local vec_r2:vec3=vec3.Create(tri[2].v[0]-tri[0].v[0], tri[2].v[1]-tri[0].v[1], tri[2].v[2]-tri[0].v[2])
   Local minus_d:vec3=vec3.Create(-ray.d.v[0], -ray.d.v[1], -ray.d.v[2])
   
   Local det:Float=vec3.Det(minus_d, vec_r1, vec_r2)
   If det=0 Or (frontfacing_only And det>0) Then Return False
   
   Local t:Float=vec3.Det(vec_dst, vec_r1, vec_r2)/det
   If t<ray.min_t Then Return False
   
   Local lambda:Float=vec3.Det(minus_d, vec_dst, vec_r2)/det
   If lambda<0 Or lambda>1 Then Return False
   
   Local my:Float=vec3.Det(minus_d, vec_r1, vec_dst)/det
   If my<0 Or my+lambda>1 Then Return False
   
   If t<ray.t Then
      ray.t=t
      Return 2
   EndIf
   Return 1
End Function

Type Ray
   
   Field o:vec3
   Field d:vec3
   Field t:Float
   Field min_t:Float
   
   Method New()
      t=1e10
   End Method
   
   Function Create:Ray(o:vec3, d:vec3, min_t:Float=0.0)
      Local r:Ray=New Ray
      r.o=o
      r.d=d
      r.d.Normalize()
      r.min_t=min_t
      Return r:Ray
   End Function
   
   
   
   
   Method Transform:Ray(mat:mat4)
      Local r:Ray=Create(vec3.Cast(mat.Transform(vec4.Cast(o, 1.0))), vec3.Cast(mat.Transform(vec4.Cast(d, 0.0))))
      r.t=t
      Return r:Ray
   End Method
   
   
   
   
End Type