判断二维坐标系中任意一个点是否在三角形内

给定平面上一点p(x0,y0)，判断该点是否在三角形ABC中，三角形顶点坐标分别为A(xa,xb),B(xb,yb),C(xc,yc)。可以使用面积法来判断，方法如下：其中S(A,B,C)表示三角形ABC的面积。
1、 若abs( S(A,B,C) ) = abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) ，则P在三角形ABC的内部或边上；如果还有abs( S(P,B,C) )、abs( S(A,P,C) ) 和abs( S(A,B,P) )全都大于0，则说明P在三角形ABC的内部，否则P在三角形ABC的边上，具体为：S(P,B,C)为0，则说明P在BC边上，S(A,P,C)为0，则说明P在AC边上，S(A,B,P)为0，则说明P在AB边上；
2、 若abs( S(A,B,C) ) < abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) ，则P在三角形ABC的外部；
3、 对abs( S(A,B,C) ) > abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) 情况在理论上是不存在的
此处又引出另一个问题，如何求平面中三角形的面积？这个可以使用叉乘法来实现，即S(A,B,C)为向量AB叉乘AC所得向量模的1/2，再对该值求绝对值就是三角形ABC的面积。

#include <stdio.h>
#include <math.h>

#define ABS_FLOAT_0  0.0001

struct point_float
{
    float x;
    float y;
};

/*
 * 计算三角形面积
 */
float GetTriangleSquar(point_float pf0, point_float pf1, point_float pf2)
{
    point_float AB, BC;
    AB.x = pf1.x - pf0.x;
    AB.y = pf1.y - pf0.y;
    BC.x = pf2.x - pf1.x;
    BC.y = pf2.y - pf1.y;
    return fabs((AB.x * BC.y - AB.y * BC.x)) / 2.0f;
}

/*
 * 判断给定一点是否在三角形内或边上
*/
int IsInTriangle(point_float A, point_float B, point_float C, point_float D)
{
    float SABC, SADB, SBDC, SADC, Sum;
    SABC = GetTriangleSquar(A, B, C);
    SADB = GetTriangleSquar(A, D, B);
    SBDC = GetTriangleSquar(B, D, C);
    SADC = GetTriangleSquar(A, D, C);

    Sum = SADB + SBDC + SADC;

    if ((-ABS_FLOAT_0 < (SABC - Sum)) && ((SABC - Sum) < ABS_FLOAT_0))
    {
        return 1;
    }
    else
    {
        return 0;
    }
}

int main(void)
{
    point_float A, B, C, P;
    A.x = 6;
    A.y = 10;
    B.x = 3;
    B.y = 4;
    C.x = 10;
    C.y = 2;
    P.x = 6;
    P.y = 10;

    if (IsInTriangle(A, B, C, P))
    {
        printf("P is in ABC");
    }
    else
    {
        printf("P is not in ABC");
    }
    return 0;
}