def point_equal(pa, pb)
      # x1 / z1 == x2 / z2  <==>  x1 * z2 == x2 * z1
      return false if (pa[0] * pb[2] - pb[0] * pa[2]) % @@p != 0
      return false if (pa[1] * pb[2] - pb[1] * pa[2]) % @@p != 0
      return true
    end
def point_equal(pa, pb)
      # x1 / z1 == x2 / z2  <==>  x1 * z2 == x2 * z1
      return false if (pa[0] * pb[2] - pb[0] * pa[2]) % @@p != 0
      return false if (pa[1] * pb[2] - pb[1] * pa[2]) % @@p != 0
      return true
    end
return point A == point B