jubjub.py 764 B

1234567891011121314151617181920212223242526272829303132
  1. from finite_fields.modp import IntegersModP
  2. q = 0x73eda753299d7d483339d80809a1d80553bda402fffe5bfeffffffff00000001
  3. modq = IntegersModP(q)
  4. a = modq(-1)
  5. d = -(modq(10240)/modq(10241))
  6. params = (a, d)
  7. def is_jubjub(params, x, y):
  8. a, d = params
  9. return a * x**2 + y**2 == 1 + d * x**2 * y**2
  10. def add(params, point_1, point_2):
  11. # From here: https://z.cash/technology/jubjub/
  12. a, d = params
  13. x1, y1 = point_1
  14. x2, y2 = point_2
  15. x3 = (x1 * y2 + y1 * x2) / (1 + d * x1 * x2 * y1 * y2)
  16. y3 = (y1 * y2 + x1 * x2) / (1 - d * x1 * x2 * y1 * y2)
  17. return (x3, y3)
  18. x = 0x15a36d1f0f390d8852a35a8c1908dd87a361ee3fd48fdf77b9819dc82d90607e
  19. y = 0x015d8c7f5b43fe33f7891142c001d9251f3abeeb98fad3e87b0dc53c4ebf1891
  20. print(is_jubjub(params, x, y))