X2624

SUM OF PU(73)

< A B C >: Reference Triangle △ABC [ A B ]: Sideline AB (Base c = 1) [ B C ]: Sideline BC (Side a) [ C A ]: Sideline CA (Side b) { 44, 513 }: Central Line L(44, 513) { 163, 1983 }: Central Line L(163, 1983) { 1109, 2628 }: Central Line L(1109, 2628) A = Vertex A (0.0000, 0.0000) B = Vertex B (1.0000, 0.0000) C = Vertex C (0.6122, 0.5999) X_44 = -LINE CONJUGATE OF X(1) (2.1429, 2.0996) X_513 = ISOGONAL CONJUGATE OF X(100) (1.0529, -0.0806) X_163 = TRILINEAR PRODUCT X(6)*X(110) (0.3174, 0.5883) X_1983 = VEGA TRANSFORM OF X(523) (0.9661, 0.8719) X_1109 = TRILINEAR SQUARE OF X(523) (0.7246, 0.0705) X_2628 = MIDPOINT OF PU(73) (0.6720, 0.0304) X_2624 = SUM OF PU(73) (2.8696, 1.7039) A B C X44 X513 X163 X1983 X1109 X2628 X2624 (a, b, c) = (5, 6, 7) | AB = 1

X2624 (SUM OF PU(73)) constructed on acute reference triangle \((a, b, c) = (5, 6, 7)\) with unit baseline AB = \(1\).

"""GEOMETOR construction model for X_2624 = SUM OF PU(73)."""

        from __future__ import annotations

        from pathlib import Path

        import sympy.geometry as spg
        from geometor.model import Model

        from geometor.etc.ops import (
            central_line,
            construct_baseline_triangle,
            create_model,
            eval_trilinear_point,
            promote_point_to_top,
        )
        from geometor.etc.render import write_center_svg


        def construct_x_02624(
            model: Model, A: spg.Point, B: spg.Point, C: spg.Point
        ) -> spg.Point:
            """Construct X_2624 = SUM OF PU(73) on reference triangle ABC."""
            X_2624 = eval_trilinear_point(model, A, B, C, 'cot(A - B) + cot(A - C)', ID="X_2624", classes=["selected", "target"])
            X_44 = eval_trilinear_point(model, A, B, C, 'b + c - 2*a', ID="X_44", classes=["center", "incident"])
            X_513 = eval_trilinear_point(model, A, B, C, 'b - c', ID="X_513", classes=["center", "incident"])
            central_line(model, X_44, X_513, 44, 513, incident=True)
            X_163 = eval_trilinear_point(model, A, B, C, '(sin(2*A))/(sin(2*B) - sin(2*C))', ID="X_163", classes=["center", "incident"])
            X_1983 = eval_trilinear_point(model, A, B, C, '(sin(C - A) - sin(A - B))/sin(B - C)', ID="X_1983", classes=["center", "incident"])
            central_line(model, X_163, X_1983, 163, 1983, incident=True)
            X_1109 = eval_trilinear_point(model, A, B, C, '(sin(B - C))**(2)', ID="X_1109", classes=["center", "incident"])
            X_2628 = eval_trilinear_point(model, A, B, C, '((a*cot(B - A) + b*cot(C - B) + c*cot(A - C))*cot(C - A) + (a*cot(C - A) + b*cot(A - B) + c*cot(B - C))*cot(B - A))', ID="X_2628", classes=["center", "incident"])
            central_line(model, X_1109, X_2628, 1109, 2628, incident=True)
            return promote_point_to_top(model, X_2624)


        def build_model(
            sides: tuple[float, float, float] = (5.0, 6.0, 7.0),
            show_guides: bool = False,
        ) -> Model:
            """Build a `geometor.model.Model` for X_2624 on the given reference triangle."""
            model = create_model('X_2624 = SUM OF PU(73)')
            A, B, C = construct_baseline_triangle(model, *sides, show_guides=show_guides)
            construct_x_02624(model, A, B, C)
            return model


        if __name__ == "__main__":
            here = Path(__file__).resolve().parent
            model_acute = build_model((5.0, 6.0, 7.0))
            model_acute.save(str(here / "x-02624.json"))
            write_center_svg(
                model_acute,
                here / "x-02624.svg",
                subtitle="(a, b, c) = (5, 6, 7) | AB = 1",
            )