[closed] Determine the position of a intersection between a line and a triangle (raycasting)

I've found the answer!

The main problem with my code was that my cross product of 2 3D vectors was completely wrong.
The current math is the following with the value of P being stored in Vector3D I, if the method return 1 it was a successful operation, 0 means not within the triangle and 2 means that it's parallel to the triangle.

Given values:
Ray R = P0 (near point), P1 (far point)
Triangle T = V0, V1, V2 which are A B and C respectively 
All values are Vector3D objects

Vector3D Normal = T.V1 - T.V0;
Normal.cross(T.V2 - T.V0); //value stored in Normal
Normal.Normalize(); //divide all values by the length of the Vector (sqrt(x^2+y^2+z^2))

Vector3D d1 = R.P0 - T.V0;
Vector3D d2 = R.P1 - T.V0;

float dist1 = d1.dot(Normal);
float dist2 = d2.dot(Normal);

if(dist1 * dist2 >= 0.0f)
    return 0;//no cross

if(dist1==dist2)
    return 2;//parallel

*I = R.P0 + (R.P1-R.P0) * (-dist1/(dist2-dist1)); //Here the exact position is calculated

//check if the position in within the triangle and thus a valid value
Vector3D u = T.V1 - T.V0;
Vector3D v = T.V2 - T.V0;

// is I inside T?
float    uu, uv, vv, wu, wv, D;
uu = u.dot(u);
uv = u.dot(v);
vv = v.dot(v);
Vector3D w = *I - T.V0;
wu = w.dot(u);
wv = w.dot(v);
D = uv * uv - uu * vv;

// get and test parametric coords
float s, t;
s = (uv * wv - vv * wu) / D;
if (s < 0.0 || s > 1.0)         // I is outside T
    return 0;
t = (uv * wu - uu * wv) / D;
if (t < 0.0 || (s + t) > 1.0)  // I is outside T
    return 0;

return 1;                       // I is in T

The Vector3D math are calculated as such:

float Vector3D::dot(Vector3D v)
{
    return (x*v.x)+(y*v.y)+(z*v.z);
}    
void Vector3D::cross(Vector3D v)
{
    Vector3D o = Vector3D(x,y,z);

    x = o.y*v.z - o.z*v.y;
    y = o.z*v.x - o.x*v.z;
    z = o.x*v.y - o.y*v.x;
}    
void Vector3D::normalize()
{
    float l = length();

    x = x/l;
    y = y/l;
    z = z/l;
}    
float Vector3D::length() // |vector| = length
{
    return sqrtf(powf(x, 2)+powf(y,2)+powf(z,2));
}