Intersection of Line and Generalized Ellipse

这个改进的函数用于计算直线与广义椭圆的交点。（This improved function returns the points of intersection between a line and a generalized ellipse.）

这是个改进的lineEllipse函数，用于计算直线与广义椭圆的交点。增加了异常参数处理，并且改进了交点有效性判断。

基于Paurakh的工作成果。具体参见：(Intersection of Line to Generalized Ellipse)。

This is an improved lineEllipse function for calculating the intersection of a line and a generalized ellipse. The processing of exception parameters is increased, and the validity judgment of intersection points is improved. Based on Paurakh's work. See: (Intersection of Line to Generalized Ellipse).

INPUT : a - major axis of Ellipse.

b - minor axis of Ellipse.

x0 - Abscissa of the center point of the ellipse.

y0 - Ordinate of the center point of the ellipse.

phi - Angle (in radian) between x-axis and the major axis,

counter-clockwise axis of rotation.

p, q - line is defined by (x/p + y/q = 1), with p and q as

x axis and y axis intercepts, p or q ∊ [-inf, inf].

OUTPUT : C1, C2 - respectively the two points of intersection of line and

ellipse are (C1, C2).

POST CONDITION:

If (line does not intersect) NaN is returned.

If a single point of intersection, (C1) = (C2) or one of the point

is NaN.

If two points of intersetion, two distinct (C1) and (C2) returned.

Copyright (c)2024, Aeron. Shanghai China.