X2660

TRILINEAR POLE OF PU(83)

< 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) { 410, 823 }: Central Line L(410, 823) { 4, 1248 }: Central Line L(4, 1248) { 416, 1936 }: Central Line L(416, 1936) { 1, 2656 }: Central Line L(1, 2656) A = Vertex A (0.0000, 0.0000) B = Vertex B (1.0000, 0.0000) C = Vertex C (0.6122, 0.5999) X_2659 = CROSSDIFFERENCE OF PU(83) (0.7496, 1.2245) X_410 = EULER X(29)-1st-SUBSTITUTION POINT (0.5800, 0.3115) X_823 = ISOGONAL CONJUGATE OF X(822) (0.5623, 0.6035) X_4 = ORTHOCENTER (0.6122, 0.3958) X_1248 = ISOGONAL CONJUGATE OF X(1047) (0.6139, 0.6190) X_416 = -HIRST INVERSE OF X(21) (0.0862, -0.9805) X_1936 = CROSSDIFFERENCE OF PU(15) (-1.1948, -3.0115) X_1 = INCENTER (0.5714, 0.2333) X_2656 = TRILINEAR POLE OF PU(82) (0.6337, -0.4192) X_2660 = TRILINEAR POLE OF PU(83) (0.6083, -0.1529) A B C X2659 X410 X823 X4 X1248 X416 X1936 X1 X2656 X2660 (a, b, c) = (5, 6, 7) | AB = 1

X2660 (TRILINEAR POLE OF PU(83)) constructed on acute reference triangle \((a, b, c) = (5, 6, 7)\) with unit baseline AB = \(1\).

"""GEOMETOR construction model for X_2660 = TRILINEAR POLE OF PU(83)."""

        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_barycentric_point,
            eval_trilinear_point,
            isogonal_conjugate,
            orthocenter_of,
            point_from_barycentrics,
            promote_point_to_top,
        )
        from geometor.etc.render import write_center_svg


        def construct_x_02660(
            model: Model, A: spg.Point, B: spg.Point, C: spg.Point
        ) -> spg.Point:
            """Construct X_2660 = TRILINEAR POLE OF PU(83) on reference triangle ABC."""
            X_2659 = eval_trilinear_point(model, A, B, C, '(1/((sec(A) + sec(B))*(sec(A) + sec(C))) - 1/(sec(B) + sec(C))**(2))', ID="X_2659", classes=["center"])
            X_2660 = isogonal_conjugate(model, A, B, C, X_2659, ID="X_2660", classes=["selected", "target"])
            X_410 = eval_barycentric_point(model, A, B, C, '(a + b) * (a + c) * (a**(2) + b**(2) - c**(2)) * (a**(2) - b**(2) + c**(2)) * (a**(8)*b**(2) - 3*a**(6)*b**(4) + 3*a**(4)*b**(6) - a**(2)*b**(8) + 3*a**(8)*b*c - a**(7)*b**(2)*c - 5*a**(6)*b**(3)*c + a**(5)*b**(4)*c + a**(4)*b**(5)*c + a**(3)*b**(6)*c + a**(2)*b**(7)*c - a*b**(8)*c + a**(8)*c**(2) - a**(7)*b*c**(2) - 3*a**(6)*b**(2)*c**(2) + a**(5)*b**(3)*c**(2) + 2*a**(4)*b**(4)*c**(2) + a**(3)*b**(5)*c**(2) + a**(2)*b**(6)*c**(2) - a*b**(7)*c**(2) - b**(8)*c**(2) - 5*a**(6)*b*c**(3) + a**(5)*b**(2)*c**(3) + 8*a**(4)*b**(3)*c**(3) - 2*a**(3)*b**(4)*c**(3) - a**(2)*b**(5)*c**(3) + a*b**(6)*c**(3) - 2*b**(7)*c**(3) - 3*a**(6)*c**(4) + a**(5)*b*c**(4) + 2*a**(4)*b**(2)*c**(4) - 2*a**(3)*b**(3)*c**(4) + a*b**(5)*c**(4) + b**(6)*c**(4) + a**(4)*b*c**(5) + a**(3)*b**(2)*c**(5) - a**(2)*b**(3)*c**(5) + a*b**(4)*c**(5) + 4*b**(5)*c**(5) + 3*a**(4)*c**(6) + a**(3)*b*c**(6) + a**(2)*b**(2)*c**(6) + a*b**(3)*c**(6) + b**(4)*c**(6) + a**(2)*b*c**(7) - a*b**(2)*c**(7) - 2*b**(3)*c**(7) - a**(2)*c**(8) - a*b*c**(8) - b**(2)*c**(8))', ID="X_410", classes=["center", "incident"])
            X_823 = eval_trilinear_point(model, A, B, C, '1/((sec(B))**(2) - (sec(C))**(2))', ID="X_823", classes=["center", "incident"])
            central_line(model, X_410, X_823, 410, 823, incident=True)
            X_4 = orthocenter_of(model, A, B, C, ID="X_4", classes=["center", "incident"])
            X_1248 = point_from_barycentrics(model, A, B, C, -0.013970382, -0.017899094, 1.0318695, ID="X_1248", classes=["center", "incident"])
            central_line(model, X_4, X_1248, 4, 1248, incident=True)
            X_416 = eval_barycentric_point(model, A, B, C, 'a*(a + b) * (a + c) * (a**(7)*b**(2) - a**(6)*b**(3) - 2*a**(5)*b**(4) + 2*a**(4)*b**(5) + a**(3)*b**(6) - a**(2)*b**(7) + a**(7)*b*c - a**(6)*b**(2)*c - a**(5)*b**(3)*c + 2*a**(4)*b**(4)*c - a**(3)*b**(5)*c - a**(2)*b**(6)*c + a*b**(7)*c + a**(7)*c**(2) - a**(6)*b*c**(2) + a**(5)*b**(2)*c**(2) - a**(4)*b**(3)*c**(2) - a**(3)*b**(4)*c**(2) + a**(2)*b**(5)*c**(2) - a*b**(6)*c**(2) + b**(7)*c**(2) - a**(6)*c**(3) - a**(5)*b*c**(3) - a**(4)*b**(2)*c**(3) + 2*a**(3)*b**(3)*c**(3) + a**(2)*b**(4)*c**(3) - a*b**(5)*c**(3) + b**(6)*c**(3) - 2*a**(5)*c**(4) + 2*a**(4)*b*c**(4) - a**(3)*b**(2)*c**(4) + a**(2)*b**(3)*c**(4) + 2*a*b**(4)*c**(4) - 2*b**(5)*c**(4) + 2*a**(4)*c**(5) - a**(3)*b*c**(5) + a**(2)*b**(2)*c**(5) - a*b**(3)*c**(5) - 2*b**(4)*c**(5) + a**(3)*c**(6) - a**(2)*b*c**(6) - a*b**(2)*c**(6) + b**(3)*c**(6) - a**(2)*c**(7) + a*b*c**(7) + b**(2)*c**(7))', ID="X_416", classes=["center", "incident"])
            X_1936 = eval_trilinear_point(model, A, B, C, '(cos(A))**(2) - cos(B)*cos(C)', ID="X_1936", classes=["center", "incident"])
            central_line(model, X_416, X_1936, 416, 1936, incident=True)
            return promote_point_to_top(model, X_2660)


        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_2660 on the given reference triangle."""
            model = create_model('X_2660 = TRILINEAR POLE OF PU(83)')
            A, B, C = construct_baseline_triangle(model, *sides, show_guides=show_guides)
            construct_x_02660(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-02660.json"))
            write_center_svg(
                model_acute,
                here / "x-02660.svg",
                subtitle="(a, b, c) = (5, 6, 7) | AB = 1",
            )