Area within a simple closed contour
Given a simple (i.e. not self-intersecting) closed contour with line/arc/elliptic arc entities, the enclosed area is found by Green's Theorem
This type of line integrals has closed forms for lines, arcs, and ellipses.
where (x0,y0) and (x1,y1) are the start and end point of the line, perspectively.
Ignoring the constant of integration, the antiderivative can be found as
where cx is the x-coordinate of arc center, and r, t are radius and arc angle, respectively.
As an ellipse can be scaled to a circle, the arc formula can be used using r as the major radius and t as the elliptic parameter angle. The resulted area by arc integral should be scaled back to elliptic by timing the ellipse ratio.
where cx is the x-coordinate of the ellipse center, and a, b are major and minor radius, respectively; y is the y-coordinate of the endpoint, A is direction of major axis, and t is the elliptic angle,
The nearest point on an ellipse to a given point
An ellipse in the coordinates orientated alone its major and minor axes is given as,
The squared distance from a point on ellipse to a given point(x,y),
The stationary points at the zero points of its first order derivative of t,
This stationary condition is a quartic equation of cos t. With variable change u = cos t,
γ2u4 − 2αγu3 + (α2 + β2 − γ2)u2 + 2αγu − α2 = 0
where α = 2ax, β = 2by, and γ = 2(a2 − b2).
For all solutions from the quartic equation, the minimum distance point is identified by the convex condition,
Given a circle, construct a tangent circle passing two given points
Given a circle, the path of center of all tangent circles passing a given point is either an ellipse or a hyperbola. A quadratic form of either ellipse or hyperbola is constructed for each of the given points, and the center of the circle to be constructed in on the intersections of these two quadratic forms. After locating of the circle center, the radius is found by the distance from either of the given points to the center.
LibreCAD has a powerful general quadratic form framework to ease the construction of such quadratic forms and location of intersections between two quadratic forms.
modular function in LibreCAD
the standard glibc fmod(x, a) is not convenient here, since we need a modular function work the same way for both positive and negative numbers, instead, we use,