"""Independent deterministic reproduction for the Menkaure selective-aperture model.

The 28 x 35.5 in aperture is the principal same-state evidence. Quantities
using BUILT_UP_RUN are conditional because its overlap with that aperture
is not established by Perring.
"""
import math
L,W,H=96.0,37.0,35.0
RAW_H,OPEN_W,OPEN_H,BUILT_UP_RUN=57.0,35.5,28.0,201.0
ROUNDING=0.5
MIN_EDGE=min(L,W,H)
MIN_FACE=W*H
OPEN_AREA=OPEN_W*OPEN_H
SPACE_DIAGONAL=math.sqrt(L**2+W**2+H**2)
ADVERSE_EXCLUSION=(MIN_EDGE-ROUNDING)-(OPEN_H+ROUNDING)
CRITICAL_SYMMETRIC_ERROR=(MIN_EDGE-OPEN_H)/2
results={"aperture_width_clearance_in":OPEN_W-H,"aperture_height_clearance_in":OPEN_H-MIN_EDGE,"height_exclusion_cm":(MIN_EDGE-OPEN_H)*2.54,"height_exclusion_percent":100*(MIN_EDGE-OPEN_H)/MIN_EDGE,"adverse_rounding_exclusion_in":ADVERSE_EXCLUSION,"adverse_rounding_exclusion_cm":ADVERSE_EXCLUSION*2.54,"critical_symmetric_error_in":CRITICAL_SYMMETRIC_ERROR,"conditional_aperture_area_in2":OPEN_AREA,"conditional_minimum_projection_area_in2":MIN_FACE,"conditional_projection_area_deficit_in2":OPEN_AREA-MIN_FACE,"conditional_projection_area_deficit_percent":100*(OPEN_AREA-MIN_FACE)/MIN_FACE,"coffer_space_diagonal_in":SPACE_DIAGONAL,"conditional_run_to_space_diagonal_ratio":BUILT_UP_RUN/SPACE_DIAGONAL,"raw_height_clearance_in":RAW_H-W,"transverse_diagonal_in":math.hypot(W,H),"upper_socket_residual_in":40-W,"final_reach_residual_in":54-math.hypot(W,H)}
for k,v in results.items(): print(f"{k},{v:.9f}")
assert results["aperture_width_clearance_in"]==0.5
assert results["aperture_height_clearance_in"]==-7.0
assert results["height_exclusion_cm"]==17.78
assert results["adverse_rounding_exclusion_in"]==6.0
assert results["critical_symmetric_error_in"]==3.5
# The next two checks reproduce only the explicitly conditional sensitivity.
assert OPEN_AREA<MIN_FACE
assert BUILT_UP_RUN>SPACE_DIAGONAL