From d40de259097c5e8d8fd35539560ca7c3d47523e7 Mon Sep 17 00:00:00 2001 From: Ryan Rueger Date: Sat, 1 Mar 2025 20:25:41 +0100 Subject: Initial Commit MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Co-Authored-By: Damien Robert Co-Authored-By: Frederik Vercauteren Co-Authored-By: Jonathan Komada Eriksen Co-Authored-By: Pierrick Dartois Co-Authored-By: Riccardo Invernizzi Co-Authored-By: Ryan Rueger [0.01s] Co-Authored-By: Benjamin Wesolowski Co-Authored-By: Arthur Herlédan Le Merdy Co-Authored-By: Boris Fouotsa --- theta_lib/theta_structures/Tuple_point.py | 106 ++++++++++++++++++++++++++++++ 1 file changed, 106 insertions(+) create mode 100644 theta_lib/theta_structures/Tuple_point.py (limited to 'theta_lib/theta_structures/Tuple_point.py') diff --git a/theta_lib/theta_structures/Tuple_point.py b/theta_lib/theta_structures/Tuple_point.py new file mode 100644 index 0000000..aac248b --- /dev/null +++ b/theta_lib/theta_structures/Tuple_point.py @@ -0,0 +1,106 @@ +from sage.all import * +from ..utilities.discrete_log import weil_pairing_pari + +class TuplePoint: + def __init__(self,*args): + if len(args)==1: + self._points=list(args[0]) + else: + self._points=list(args) + + def points(self): + return self._points + + def parent_curves(self): + return [x.curve() for x in self._points] + + def parent_curve(self,i): + return self._points[i].curve() + + def n_points(self): + return len(self._points) + + def is_zero(self): + return all([self._points[i]==0 for i in range(self.n_points())]) + + def __repr__(self): + return str(self._points) + + def __getitem__(self,i): + return self._points[i] + + def __setitem__(self,i,P): + self._points[i]=P + + def __eq__(self,other): + n_self=self.n_points() + n_other=self.n_points() + return n_self==n_other and all([self._points[i]==other._points[i] for i in range(n_self)]) + + def __add__(self,other): + n_self=self.n_points() + n_other=self.n_points() + + if n_self!=n_other: + raise ValueError("Cannot add TuplePoint of distinct lengths {} and {}.".format(n_self,n_other)) + + points=[] + for i in range(n_self): + points.append(self._points[i]+other._points[i]) + return self.__class__(points) + + def __sub__(self,other): + n_self=self.n_points() + n_other=self.n_points() + + if n_self!=n_other: + raise ValueError("Cannot substract TuplePoint of distinct lengths {} and {}.".format(n_self,n_other)) + + points=[] + for i in range(n_self): + points.append(self._points[i]-other._points[i]) + return self.__class__(points) + + def __neg__(self): + n_self=self.n_points() + points=[] + for i in range(n_self): + points.append(-self._points[i]) + return self.__class__(points) + + def __mul__(self,m): + n_self=self.n_points() + points=[] + for i in range(n_self): + points.append(m*self._points[i]) + return self.__class__(points) + + def __rmul__(self,m): + return self*m + + def double_iter(self,n): + result=self + for i in range(n): + result=2*result + return result + + def weil_pairing(self,other,n): + n_self=self.n_points() + n_other=self.n_points() + + if n_self!=n_other: + raise ValueError("Cannot compute the Weil pairing of TuplePoint of distinct lengths {} and {}.".format(n_self,n_other)) + + zeta=1 + for i in range(n_self): + zeta*=weil_pairing_pari(self._points[i],other._points[i],n) + + return zeta + + + + + + + + -- cgit v1.2.3-70-g09d2