test_powerset.py 4.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141
  1. from sympy.core.expr import unchanged
  2. from sympy.core.singleton import S
  3. from sympy.core.symbol import Symbol
  4. from sympy.sets.contains import Contains
  5. from sympy.sets.fancysets import Interval
  6. from sympy.sets.powerset import PowerSet
  7. from sympy.sets.sets import FiniteSet
  8. from sympy.testing.pytest import raises, XFAIL
  9. def test_powerset_creation():
  10. assert unchanged(PowerSet, FiniteSet(1, 2))
  11. assert unchanged(PowerSet, S.EmptySet)
  12. raises(ValueError, lambda: PowerSet(123))
  13. assert unchanged(PowerSet, S.Reals)
  14. assert unchanged(PowerSet, S.Integers)
  15. def test_powerset_rewrite_FiniteSet():
  16. assert PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet) == \
  17. FiniteSet(S.EmptySet, FiniteSet(1), FiniteSet(2), FiniteSet(1, 2))
  18. assert PowerSet(S.EmptySet).rewrite(FiniteSet) == FiniteSet(S.EmptySet)
  19. assert PowerSet(S.Naturals).rewrite(FiniteSet) == PowerSet(S.Naturals)
  20. def test_finiteset_rewrite_powerset():
  21. assert FiniteSet(S.EmptySet).rewrite(PowerSet) == PowerSet(S.EmptySet)
  22. assert FiniteSet(
  23. S.EmptySet, FiniteSet(1),
  24. FiniteSet(2), FiniteSet(1, 2)).rewrite(PowerSet) == \
  25. PowerSet(FiniteSet(1, 2))
  26. assert FiniteSet(1, 2, 3).rewrite(PowerSet) == FiniteSet(1, 2, 3)
  27. def test_powerset__contains__():
  28. subset_series = [
  29. S.EmptySet,
  30. FiniteSet(1, 2),
  31. S.Naturals,
  32. S.Naturals0,
  33. S.Integers,
  34. S.Rationals,
  35. S.Reals,
  36. S.Complexes]
  37. l = len(subset_series)
  38. for i in range(l):
  39. for j in range(l):
  40. if i <= j:
  41. assert subset_series[i] in \
  42. PowerSet(subset_series[j], evaluate=False)
  43. else:
  44. assert subset_series[i] not in \
  45. PowerSet(subset_series[j], evaluate=False)
  46. @XFAIL
  47. def test_failing_powerset__contains__():
  48. # XXX These are failing when evaluate=True,
  49. # but using unevaluated PowerSet works fine.
  50. assert FiniteSet(1, 2) not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  51. assert S.Naturals not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  52. assert S.Naturals not in PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet)
  53. assert S.Naturals0 not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  54. assert S.Naturals0 not in PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet)
  55. assert S.Integers not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  56. assert S.Integers not in PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet)
  57. assert S.Rationals not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  58. assert S.Rationals not in PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet)
  59. assert S.Reals not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  60. assert S.Reals not in PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet)
  61. assert S.Complexes not in PowerSet(S.EmptySet).rewrite(FiniteSet)
  62. assert S.Complexes not in PowerSet(FiniteSet(1, 2)).rewrite(FiniteSet)
  63. def test_powerset__len__():
  64. A = PowerSet(S.EmptySet, evaluate=False)
  65. assert len(A) == 1
  66. A = PowerSet(A, evaluate=False)
  67. assert len(A) == 2
  68. A = PowerSet(A, evaluate=False)
  69. assert len(A) == 4
  70. A = PowerSet(A, evaluate=False)
  71. assert len(A) == 16
  72. def test_powerset__iter__():
  73. a = PowerSet(FiniteSet(1, 2)).__iter__()
  74. assert next(a) == S.EmptySet
  75. assert next(a) == FiniteSet(1)
  76. assert next(a) == FiniteSet(2)
  77. assert next(a) == FiniteSet(1, 2)
  78. a = PowerSet(S.Naturals).__iter__()
  79. assert next(a) == S.EmptySet
  80. assert next(a) == FiniteSet(1)
  81. assert next(a) == FiniteSet(2)
  82. assert next(a) == FiniteSet(1, 2)
  83. assert next(a) == FiniteSet(3)
  84. assert next(a) == FiniteSet(1, 3)
  85. assert next(a) == FiniteSet(2, 3)
  86. assert next(a) == FiniteSet(1, 2, 3)
  87. def test_powerset_contains():
  88. A = PowerSet(FiniteSet(1), evaluate=False)
  89. assert A.contains(2) == Contains(2, A)
  90. x = Symbol('x')
  91. A = PowerSet(FiniteSet(x), evaluate=False)
  92. assert A.contains(FiniteSet(1)) == Contains(FiniteSet(1), A)
  93. def test_powerset_method():
  94. # EmptySet
  95. A = FiniteSet()
  96. pset = A.powerset()
  97. assert len(pset) == 1
  98. assert pset == FiniteSet(S.EmptySet)
  99. # FiniteSets
  100. A = FiniteSet(1, 2)
  101. pset = A.powerset()
  102. assert len(pset) == 2**len(A)
  103. assert pset == FiniteSet(FiniteSet(), FiniteSet(1),
  104. FiniteSet(2), A)
  105. # Not finite sets
  106. A = Interval(0, 1)
  107. assert A.powerset() == PowerSet(A)
  108. def test_is_subset():
  109. # covers line 101-102
  110. # initialize powerset(1), which is a subset of powerset(1,2)
  111. subset = PowerSet(FiniteSet(1))
  112. pset = PowerSet(FiniteSet(1, 2))
  113. bad_set = PowerSet(FiniteSet(2, 3))
  114. # assert "subset" is subset of pset == True
  115. assert subset.is_subset(pset)
  116. # assert "bad_set" is subset of pset == False
  117. assert not pset.is_subset(bad_set)