ec_msm.sage 1.1 KB

1234567891011121314151617181920212223242526272829303132
  1. load('curve.sage')
  2. load('share.sage')
  3. load('ec_share.sage')
  4. load('beaver.sage')
  5. N = 10
  6. source = Source(p)
  7. points = [CurvePoint.random() for _ in range(0, N)]
  8. lhs_points = [pt - CurvePoint.random() for pt in points]
  9. rhs_points = [point - lhs for (point, lhs) in zip(points, lhs_points)]
  10. lhs_points_shares = [ECAuthenticatedShare(pt) for pt in lhs_points]
  11. rhs_points_shares = [ECAuthenticatedShare(pt) for pt in rhs_points]
  12. scalars = [random.randint(0,p) for i in range(0, N)]
  13. lhs_scalars = [s - random.randint(0,p) for s in scalars]
  14. rhs_scalars = [s - lhs for (s, lhs) in zip(scalars, lhs_scalars)]
  15. lhs_scalars_shares = [AuthenticatedShare(s) for s in lhs_scalars]
  16. rhs_scalars_shares = [AuthenticatedShare(s) for s in rhs_scalars]
  17. lhs_msm = MSM(lhs_points_shares, lhs_scalars_shares, source, 0)
  18. rhs_msm = MSM(rhs_points_shares, rhs_scalars_shares, source, 1)
  19. res = []
  20. for lhs, rhs in zip(lhs_msm.msm(), rhs_msm.msm()):
  21. first_share = lhs*rhs
  22. second_share = rhs*lhs
  23. res += [first_share.authenticated_open(second_share)]
  24. assert (sum(res) == sum([p*s for p, s in zip(points, scalars)]))