The method for number fields iterates over the points of bounded height for each coordinate. This includes all the appropriate points, but includes some other points that are larger than the specified bound.

B=3
K.<v>=QuadraticField(3)
P.<x,y,z>=ProjectiveSpace(K,2)
for Q in P.points_of_bounded_height(bound=B):
    if exp(Q.global_height()) > B+0.001:
        print(Q,exp(Q.global_height()))

See also #31400

Depends on #34212

CC: @pfili @EnderWannabe @bhutz

Component: algebraic geometry

Keywords: gsoc2022

Author: Jing Guo

Branch/Commit: c053c2a

Reviewer: Alexander Galarraga, Ben Hutz

Issue created by migration from https://trac.sagemath.org/ticket/32686