Fast InPolygon detection MEX

Fast detection points inside a polygonal region.
To my knowlegde, the code is faster than all the inpolygon functions found on MatLab File Exchange

[IN,ON,IN_strict] = InPolygon(px,py,cx,cy);

input: px - nPM x nPN (Matrix) X-coordinates of the points to be tested
py - nPM x nPN (Matrix) Y-coordinates of the points to be tested
cx - nC x 1 (Vector) X-coordinates of the polygon
cy - nC x 1 (Vector) Y-coordinates of the polygon
output:
IN - boolean matrix of 1 if px(i,j),py(i,j) is
in or on the polygon defined by cx,cy
ON - boolean matrix of 1 if px(i,j),py(i,j) is
on the polygon defined by cx,cy
IN_strict - boolean matrix of 1 if px(i,j),py(i,j) is
exclusively in the polygon defined by cx,cy