Hide keyboard shortcuts

Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

""" 

Base class for polyhedra over ``RDF``. 

""" 

from __future__ import absolute_import 

 

from sage.rings.all import RDF 

from .base import Polyhedron_base 

 

 

 

class Polyhedron_RDF(Polyhedron_base): 

""" 

Base class for polyhedra over ``RDF``. 

 

TESTS:: 

 

sage: p = Polyhedron([(0,0)], base_ring=RDF); p 

A 0-dimensional polyhedron in RDF^2 defined as the convex hull of 1 vertex 

sage: TestSuite(p).run() 

""" 

# 1e-6 is the cddf+ default fuzzy zero cutoff 

 

def _is_zero(self, x): 

""" 

Test whether ``x`` is zero. 

 

INPUT: 

 

- ``x`` -- a number in the base ring. 

 

OUTPUT: 

 

Boolean. 

 

EXAMPLES:: 

 

sage: p = Polyhedron([(0,0)], base_ring=RDF) 

sage: p._is_zero(0) 

True 

sage: p._is_zero(1e-3) 

False 

 

This is a fuzzy zero for floating-point numbers:: 

 

sage: p._is_zero(1e-10) 

True 

""" 

return abs(x)<=1e-6 

 

def _is_nonneg(self, x): 

""" 

Test whether ``x`` is nonnegative. 

 

INPUT: 

 

- ``x`` -- a number in the base ring. 

 

OUTPUT: 

 

Boolean. 

 

EXAMPLES:: 

 

sage: p = Polyhedron([(0,0)], base_ring=RDF) 

sage: p._is_nonneg(1) 

True 

sage: p._is_nonneg(-1e-3) 

False 

 

This is a fuzzy zero for floating-point numbers:: 

 

sage: p._is_nonneg(-1e-10) 

True 

""" 

return x>=-1e-6 

 

def _is_positive(self, x): 

""" 

Test whether ``x`` is positive. 

 

INPUT: 

 

- ``x`` -- a number in the base ring. 

 

OUTPUT: 

 

Boolean. 

 

EXAMPLES:: 

 

sage: p = Polyhedron([(0,0)], base_ring=RDF) 

sage: p._is_positive(1) 

True 

sage: p._is_positive(0) 

True 

 

This is a fuzzy zero for floating-point numbers:: 

 

sage: p._is_positive(-1e-10) 

True 

""" 

return x>=-1e-6 

 

_base_ring = RDF