straight geodesic

There is an issue with straight geodesics not being computed correctly

T = HyperbolicTiling(3,7,2)
T.generate()

fig, ax = plt.subplots(figsize=(8,8), dpi=100)

for poly in T: # loop over polygons
    for i in range(T.p): # loop over vertices
        z1 = poly.verticesP[i] # extract edges
        z2 = poly.verticesP[(i+1)%T.p]
        arc = geodesic_arc(z1,z2,lw=0.6) # compute arc
        plt.gca().add_patch(arc) # add to plot
        
plt.xlim(-bnd, bnd); plt.ylim(-bnd, bnd); 
plt.axis("off"); plt.gca().set_aspect('equal')
plt.tight_layout(); plt.show()