X3018¶
(ORTHIC AXIS)∩(FERMAT AXIS)
X3018 ((ORTHIC AXIS)∩(FERMAT AXIS)) constructed on acute reference triangle \((a, b, c) = (5, 6, 7)\) with unit baseline AB = \(1\).¶
"""GEOMETOR construction model for X_3018 = (ORTHIC AXIS)∩(FERMAT AXIS)."""
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,
centroid_of,
construct_baseline_triangle,
create_model,
crosspoint,
eval_trilinear_point,
promote_point_to_top,
symmedian_point_of,
)
from geometor.etc.render import write_center_svg
def construct_x_03018(
model: Model, A: spg.Point, B: spg.Point, C: spg.Point
) -> spg.Point:
"""Construct X_3018 = (ORTHIC AXIS)∩(FERMAT AXIS) on reference triangle ABC."""
X_2 = centroid_of(model, A, B, C, ID="X_2", classes=["center"])
X_477 = eval_trilinear_point(model, A, B, C, '1/(4*cos(A) + cos(A - 2*B) + cos(A - 2*C) - 3*cos(B - C))', ID="X_477", classes=["center"])
X_3018 = crosspoint(model, A, B, C, X_2, X_477, ID="X_3018", classes=["selected", "target"])
X_6 = symmedian_point_of(model, A, B, C, ID="X_6", classes=["center", "incident"])
X_13 = eval_trilinear_point(model, A, B, C, 'csc(A + pi /3)', ID="X_13", classes=["center", "incident"])
central_line(model, X_6, X_13, 6, 13, incident=True)
X_230 = eval_trilinear_point(model, A, B, C, 'b*c*(a**(2)*(2*a**(2) - b**(2) - c**(2)) + (b**(2) - c**(2))**(2))', ID="X_230", classes=["center", "incident"])
X_231 = eval_trilinear_point(model, A, B, C, '(cos(B))*(4*(sin(C))**(2) - 1)*(cos(A)*sin(B) - sin(A)*cos(B)) - (cos(C))*(4*(sin(B))**(2) - 1)*(cos(C)*sin(A) - sin(C)*cos(A))', ID="X_231", classes=["center", "incident"])
central_line(model, X_230, X_231, 230, 231, incident=True)
X_112 = eval_trilinear_point(model, A, B, C, 'a/(sin(2*B) - sin(2*C))', ID="X_112", classes=["center", "incident"])
X_393 = eval_trilinear_point(model, A, B, C, 'b*c*(tan(A))**(2)', ID="X_393", classes=["center", "incident"])
central_line(model, X_112, X_393, 112, 393, incident=True)
return promote_point_to_top(model, X_3018)
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_3018 on the given reference triangle."""
model = create_model('X_3018 = (ORTHIC AXIS)∩(FERMAT AXIS)')
A, B, C = construct_baseline_triangle(model, *sides, show_guides=show_guides)
construct_x_03018(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-03018.json"))
write_center_svg(
model_acute,
here / "x-03018.svg",
subtitle="(a, b, c) = (5, 6, 7) | AB = 1",
)