inire/examples/04_sbends_and_radii.py

70 lines
2 KiB
Python
Raw Normal View History

2026-03-08 20:18:53 -07:00
from shapely.geometry import Polygon
from inire.geometry.collision import CollisionEngine
from inire.geometry.primitives import Port
from inire.router.astar import AStarRouter
from inire.router.config import CostConfig, RouterConfig
from inire.router.cost import CostEvaluator
from inire.router.danger_map import DangerMap
from inire.router.pathfinder import PathFinder
from inire.utils.visualization import plot_routing_results
def main() -> None:
print("Running Example 04: S-Bends and Multiple Radii...")
# 1. Setup Environment
2026-03-08 22:02:07 -07:00
bounds = (0, 0, 100, 100)
2026-03-08 20:18:53 -07:00
engine = CollisionEngine(clearance=2.0)
danger_map = DangerMap(bounds=bounds)
2026-03-08 22:02:07 -07:00
danger_map.precompute([])
# 2. Configure Router
2026-03-08 20:18:53 -07:00
evaluator = CostEvaluator(
2026-03-08 22:02:07 -07:00
engine,
2026-03-08 20:18:53 -07:00
danger_map,
unit_length_cost=1.0,
2026-03-08 22:02:07 -07:00
greedy_h_weight=1.5,
2026-03-08 20:18:53 -07:00
)
2026-03-08 22:02:07 -07:00
# We want a 45 degree switchover for S-bend.
# Offset O = 2 * R * (1 - cos(theta))
# If R = 10, O = 5.86
2026-03-08 20:18:53 -07:00
router = AStarRouter(
evaluator,
2026-03-08 22:02:07 -07:00
node_limit=50000,
bend_radii=[10.0, 30.0],
sbend_offsets=[5.0], # Use a simpler offset
sbend_radii=[10.0],
bend_penalty=10.0,
sbend_penalty=20.0,
snap_to_target_dist=50.0, # Large snap range
2026-03-08 20:18:53 -07:00
)
2026-03-08 22:02:07 -07:00
2026-03-08 20:18:53 -07:00
pf = PathFinder(router, evaluator)
2026-03-08 22:02:07 -07:00
2026-03-08 20:18:53 -07:00
# 3. Define Netlist
2026-03-08 22:02:07 -07:00
# start (10, 50), target (60, 55) -> 5um offset
2026-03-08 20:18:53 -07:00
netlist = {
2026-03-08 22:02:07 -07:00
"sbend_only": (Port(10, 50, 0), Port(60, 55, 0)),
"multi_radii": (Port(10, 10, 0), Port(90, 90, 0)),
2026-03-08 20:18:53 -07:00
}
2026-03-08 22:02:07 -07:00
net_widths = {"sbend_only": 2.0, "multi_radii": 2.0}
2026-03-08 20:18:53 -07:00
# 4. Route
results = pf.route_all(netlist, net_widths)
# 5. Check Results
for nid, res in results.items():
status = "Success" if res.is_valid else "Failed"
print(f"{nid}: {status}, collisions={res.collisions}")
# 6. Visualize
2026-03-08 22:13:10 -07:00
fig, ax = plot_routing_results(results, [], bounds, netlist=netlist)
2026-03-08 20:18:53 -07:00
fig.savefig("examples/sbends_radii.png")
print("Saved plot to examples/sbends_radii.png")
if __name__ == "__main__":
main()