在折线内的点的算法莫名其妙地失败(AutoCAD)

我正在使用Autocad的。net API，我有一个算法(我没有写)来确定一个点是否位于多边形内(仅限直线)。

我已经在相同的51个多边形上反复测试了我的命令。99%的效果都很好。每隔一段时间，它会在一个或多个多边形上失败，对于我在折线的边界框内创建的超过2000个点返回false。当折线是一个简单的矩形并且所有的点都分布在折线内的网格中时，我看到它失败了。在这种情况下，它应该在2000次以上返回true。它永远不会因为1分而不及格，它会因为所有的分而不及格。我已经确认了这些点是在我期望的地方正确创建的，多边形的顶点是在我期望的地方。当它失败时，最后一个点的最后一个角度变量正好是两倍PI。

我没有做任何多线程。我正在做的唯一可能"有趣"的事情是COM与Excel的互操作。这是在使用此算法的部分提交事务之后发生的，并且我确信我正在清理所有COM对象。我无法在没有COM互操作部分的情况下重现失败，但我认为我还没有测试到足够的证据。

你知道有什么问题吗?

bool IsInsidePolygon(Polyline polygon, Point3d pt)
    {
        int n = polygon.NumberOfVertices;
        double angle = 0;
        Point pt1, pt2;
        for (int i = 0; i < n; i++)
        {
            pt1.X = polygon.GetPoint2dAt(i).X - pt.X;
            pt1.Y = polygon.GetPoint2dAt(i).Y - pt.Y;
            pt2.X = polygon.GetPoint2dAt((i + 1) % n).X - pt.X;
            pt2.Y = polygon.GetPoint2dAt((i + 1) % n).Y - pt.Y;
            angle += Angle2D(pt1.X, pt1.Y, pt2.X, pt2.Y);
        }
        if (Math.Abs(angle) < Math.PI)
            return false;
        else
            return true;
    }
    public struct Point
    {
        public double X, Y;
    };
    public static double Angle2D(double x1, double y1, double x2, double y2)
    {
        double dtheta, theta1, theta2;
        theta1 = Math.Atan2(y1, x1);
        theta2 = Math.Atan2(y2, x2);
        dtheta = theta2 - theta1;
        while (dtheta > Math.PI)
            dtheta -= (Math.PI * 2);
        while (dtheta < -Math.PI)
            dtheta += (Math.PI * 2);
        return (dtheta);
    }