| 3911 | | pts = self.torsion_points() |
|---|
| 3912 | | pts = [P for P in pts if P[0].is_integral()] |
|---|
| 3913 | | pts.sort() |
|---|
| 3914 | | return pts |
|---|
| | 3911 | int_points = self.torsion_points() |
|---|
| | 3912 | int_points = [P for P in int_points if not P.is_zero()] |
|---|
| | 3913 | int_points = [P for P in int_points if P[0].is_integral()] |
|---|
| | 3914 | if both_signs: |
|---|
| | 3915 | int_points = list(set(int_points)) # remove duplicates |
|---|
| | 3916 | else: |
|---|
| | 3917 | xlist = set([P[0] for P in int_points]) |
|---|
| | 3918 | int_points = [self.lift_x(x) for x in xlist] |
|---|
| | 3919 | int_points.sort() |
|---|
| | 3920 | if verbose: |
|---|
| | 3921 | print 'Total number of integral points:',len(int_points) |
|---|
| | 3922 | return int_points |
|---|
| | 3923 | |
|---|