rayTriangleIntersect (Ray * r, Vec3 * v0, Vec3 * v1, Vec3 * v2){
  Vec3 v0v1 = sub (*v1, *v0);
  Vec3 v0v2 = sub (*v2, *v0);

  Vec3 pvec = cross (r->dir, v0v2);

  float det = dot (v0v1, pvec);
  if (det < 0.000001)
    return -INFINITY;

  float invDet = 1.0 / det;

  Vec3 tvec = sub (r->orig, *v0);

  float u = dot (tvec, pvec) * invDet;

  if (u < -0.0001 || u > 1.0001)
    return -INFINITY;

  Vec3 qvec = cross (tvec, v0v1);

  float v = dot (r->dir, qvec) * invDet;

  if (v < -0.0001 || u + v > 1.0001)
    return -INFINITY;

  return dot (v0v2, qvec) * invDet;   }
    

v0 = (0.0,0.0,0.1);
v1 = (0.01,0.0,0.1);
v2 = (0.002,0.01,0.1);
for int 2DY = 0 to 500; ++2DY {
for int 2DX = 0 to 500; ++2DX {
RayDir.xyz=(float 2DX,float 2DY ,-1.0); 
float r = rayTriangleIntersect(Ray,v0,v1,v2);
if (r != -INFINITY) drawone_pixel(2DX,2DY) }
}

    

vec2 f = float(2DX,2DY);
f = (2.0 * f / 500.0 - 1.0);
RayDir.xyz=(f.x,f.y ,-1.0);