def _intersection_homogenous(homog_line_0, homog_line_1):
    """Find point of intersection of two lines in homogenous coordinates"""
    # NB: renamed from '_intersection'
    eps = 1e-13
    a,b,c=homog_line_0
    u,v,w=homog_line_1
    D=float(b*u-v*a)
    if abs(D)<eps:
        # parallel lines
        return None, None
    xp=-(w*b-c*v)/D
    yp= (w*a-c*u)/D

    return xp, yp