Skip to content

Commit

Permalink
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Lighter construction of finite field elements from lists
Browse files Browse the repository at this point in the history
When doing intensive polynomial arithmetic with the NTL implementation
the constructor with lists is called a large number of times
and may spend a lot of time constructing the vector_space and
FreeModuleElement objects.

The very common call to vector_space(map=False) is optimized to be as cheap as
possible using the already cached object.

The common case of lists of length 0 and 1 is replaced by cheaper
shortcuts.
remyoudompheng committed Mar 27, 2023
1 parent 82e02a1 commit 7a6637a
Showing 2 changed files with 21 additions and 3 deletions.
14 changes: 13 additions & 1 deletion src/sage/rings/finite_rings/element_pari_ffelt.pyx
Original file line number Diff line number Diff line change
@@ -270,13 +270,19 @@ cdef class FiniteFieldElement_pari_ffelt(FinitePolyExtElement):
sage: k = FiniteField(3^11, 't', impl='pari_ffelt')
sage: k([ 0, 1/2 ])
2*t
sage: k([ 0, 1/2, 0, 0, 0, 0, 0, 0, 0, -1, 0 ])
2*t^9 + 2*t
sage: k([ k(0), k(1) ])
t
sage: k([ GF(3)(2), GF(3^5,'u')(1) ])
t + 2
sage: R.<x> = PolynomialRing(k)
sage: k([ x/x ])
1
sage: k([ R(-1), x/x ])
t + 2
sage: k([ R(-1), R(0), 0 ])
2
Check that zeros are created correctly (:trac:`11685`)::
@@ -496,7 +502,13 @@ cdef class FiniteFieldElement_pari_ffelt(FinitePolyExtElement):
self.construct_from(x.constant_coefficient())

elif isinstance(x, list):
if len(x) == self._parent.degree():
n = len(x)
if n == 0:
self.construct_from(None)
elif n == 1:
Fp = self._parent.base_ring()
self.construct_from(Fp(x[0]))
elif n == self._parent.degree():
self.construct_from(self._parent.vector_space(map=False)(x))
else:
Fp = self._parent.base_ring()
10 changes: 8 additions & 2 deletions src/sage/rings/finite_rings/finite_field_base.pyx
Original file line number Diff line number Diff line change
@@ -1229,9 +1229,8 @@ cdef class FiniteField(Field):
True
sage: all(to_V(h(c) * e) == c * to_V(e) for e in E for c in F)
True
"""
from sage.modules.all import VectorSpace
from sage.categories.morphism import is_Morphism
if subfield is not None:
if base is not None:
raise ValueError
@@ -1241,6 +1240,13 @@ cdef class FiniteField(Field):
deprecation(28481, "The default value for map will be changing to True. To keep the current behavior, explicitly pass map=False.")
map = False

if base is None and self.__vector_space is not None and not map:
# A very common case: return as early as possible.
return self.__vector_space

from sage.modules.all import VectorSpace
from sage.categories.morphism import is_Morphism

if base is None:
base = self.prime_subfield()
s = self.degree()

0 comments on commit 7a6637a

Please sign in to comment.