point_in_cone()

static int point_in_cone ( const POINT3D *  A1, const POINT3D *  A2, const POINT3D *  P  )

static

Utility function for checking if P is within the cone defined by A1/A2.

Definition at line 3057 of file lwgeodetic.c.

References dot_product(), LW_FALSE, LW_TRUE, normalize(), and vector_sum().

Referenced by edge_intersects().

{
        POINT3D AC; /* Center point of A1/A2 */
        double min_similarity, similarity;

        /* The normalized sum bisects the angle between start and end. */
        vector_sum(A1, A2, &AC);
        normalize(&AC);

        /* The projection of start onto the center defines the minimum similarity */
        min_similarity = dot_product(A1, &AC);

        /* The projection of candidate p onto the center */
        similarity = dot_product(P, &AC);

        /* If the point is more similar than the end, the point is in the cone */
        if ( similarity > min_similarity || fabs(similarity - min_similarity) < 2e-16 )
        {
                return LW_TRUE;
        }
        return LW_FALSE;
}

Here is the call graph for this function:

Here is the caller graph for this function: