compiled-knowledge 4.0.0a20__cp312-cp312-win32.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Potentially problematic release.


This version of compiled-knowledge might be problematic. Click here for more details.

Files changed (178) hide show
  1. ck/__init__.py +0 -0
  2. ck/circuit/__init__.py +17 -0
  3. ck/circuit/_circuit_cy.c +37523 -0
  4. ck/circuit/_circuit_cy.cp312-win32.pyd +0 -0
  5. ck/circuit/_circuit_cy.pxd +32 -0
  6. ck/circuit/_circuit_cy.pyx +768 -0
  7. ck/circuit/_circuit_py.py +836 -0
  8. ck/circuit/tmp_const.py +74 -0
  9. ck/circuit_compiler/__init__.py +2 -0
  10. ck/circuit_compiler/circuit_compiler.py +26 -0
  11. ck/circuit_compiler/cython_vm_compiler/__init__.py +1 -0
  12. ck/circuit_compiler/cython_vm_compiler/_compiler.c +19824 -0
  13. ck/circuit_compiler/cython_vm_compiler/_compiler.cp312-win32.pyd +0 -0
  14. ck/circuit_compiler/cython_vm_compiler/_compiler.pyx +380 -0
  15. ck/circuit_compiler/cython_vm_compiler/cython_vm_compiler.py +121 -0
  16. ck/circuit_compiler/interpret_compiler.py +223 -0
  17. ck/circuit_compiler/llvm_compiler.py +388 -0
  18. ck/circuit_compiler/llvm_vm_compiler.py +546 -0
  19. ck/circuit_compiler/named_circuit_compilers.py +57 -0
  20. ck/circuit_compiler/support/__init__.py +0 -0
  21. ck/circuit_compiler/support/circuit_analyser/__init__.py +13 -0
  22. ck/circuit_compiler/support/circuit_analyser/_circuit_analyser_cy.c +10618 -0
  23. ck/circuit_compiler/support/circuit_analyser/_circuit_analyser_cy.cp312-win32.pyd +0 -0
  24. ck/circuit_compiler/support/circuit_analyser/_circuit_analyser_cy.pyx +98 -0
  25. ck/circuit_compiler/support/circuit_analyser/_circuit_analyser_py.py +93 -0
  26. ck/circuit_compiler/support/input_vars.py +148 -0
  27. ck/circuit_compiler/support/llvm_ir_function.py +234 -0
  28. ck/example/__init__.py +53 -0
  29. ck/example/alarm.py +366 -0
  30. ck/example/asia.py +28 -0
  31. ck/example/binary_clique.py +32 -0
  32. ck/example/bow_tie.py +33 -0
  33. ck/example/cancer.py +37 -0
  34. ck/example/chain.py +38 -0
  35. ck/example/child.py +199 -0
  36. ck/example/clique.py +33 -0
  37. ck/example/cnf_pgm.py +39 -0
  38. ck/example/diamond_square.py +68 -0
  39. ck/example/earthquake.py +36 -0
  40. ck/example/empty.py +10 -0
  41. ck/example/hailfinder.py +539 -0
  42. ck/example/hepar2.py +628 -0
  43. ck/example/insurance.py +504 -0
  44. ck/example/loop.py +40 -0
  45. ck/example/mildew.py +38161 -0
  46. ck/example/munin.py +22982 -0
  47. ck/example/pathfinder.py +53747 -0
  48. ck/example/rain.py +39 -0
  49. ck/example/rectangle.py +161 -0
  50. ck/example/run.py +30 -0
  51. ck/example/sachs.py +129 -0
  52. ck/example/sprinkler.py +30 -0
  53. ck/example/star.py +44 -0
  54. ck/example/stress.py +64 -0
  55. ck/example/student.py +43 -0
  56. ck/example/survey.py +46 -0
  57. ck/example/triangle_square.py +54 -0
  58. ck/example/truss.py +49 -0
  59. ck/in_out/__init__.py +3 -0
  60. ck/in_out/parse_ace_lmap.py +216 -0
  61. ck/in_out/parse_ace_nnf.py +322 -0
  62. ck/in_out/parse_net.py +480 -0
  63. ck/in_out/parser_utils.py +185 -0
  64. ck/in_out/pgm_pickle.py +42 -0
  65. ck/in_out/pgm_python.py +268 -0
  66. ck/in_out/render_bugs.py +111 -0
  67. ck/in_out/render_net.py +177 -0
  68. ck/in_out/render_pomegranate.py +184 -0
  69. ck/pgm.py +3475 -0
  70. ck/pgm_circuit/__init__.py +1 -0
  71. ck/pgm_circuit/marginals_program.py +352 -0
  72. ck/pgm_circuit/mpe_program.py +237 -0
  73. ck/pgm_circuit/pgm_circuit.py +79 -0
  74. ck/pgm_circuit/program_with_slotmap.py +236 -0
  75. ck/pgm_circuit/slot_map.py +35 -0
  76. ck/pgm_circuit/support/__init__.py +0 -0
  77. ck/pgm_circuit/support/compile_circuit.py +83 -0
  78. ck/pgm_circuit/target_marginals_program.py +103 -0
  79. ck/pgm_circuit/wmc_program.py +323 -0
  80. ck/pgm_compiler/__init__.py +2 -0
  81. ck/pgm_compiler/ace/__init__.py +1 -0
  82. ck/pgm_compiler/ace/ace.py +299 -0
  83. ck/pgm_compiler/factor_elimination.py +395 -0
  84. ck/pgm_compiler/named_pgm_compilers.py +63 -0
  85. ck/pgm_compiler/pgm_compiler.py +19 -0
  86. ck/pgm_compiler/recursive_conditioning.py +231 -0
  87. ck/pgm_compiler/support/__init__.py +0 -0
  88. ck/pgm_compiler/support/circuit_table/__init__.py +17 -0
  89. ck/pgm_compiler/support/circuit_table/_circuit_table_cy.c +16396 -0
  90. ck/pgm_compiler/support/circuit_table/_circuit_table_cy.cp312-win32.pyd +0 -0
  91. ck/pgm_compiler/support/circuit_table/_circuit_table_cy.pyx +332 -0
  92. ck/pgm_compiler/support/circuit_table/_circuit_table_py.py +304 -0
  93. ck/pgm_compiler/support/clusters.py +568 -0
  94. ck/pgm_compiler/support/factor_tables.py +406 -0
  95. ck/pgm_compiler/support/join_tree.py +332 -0
  96. ck/pgm_compiler/support/named_compiler_maker.py +43 -0
  97. ck/pgm_compiler/variable_elimination.py +91 -0
  98. ck/probability/__init__.py +0 -0
  99. ck/probability/empirical_probability_space.py +50 -0
  100. ck/probability/pgm_probability_space.py +32 -0
  101. ck/probability/probability_space.py +622 -0
  102. ck/program/__init__.py +3 -0
  103. ck/program/program.py +137 -0
  104. ck/program/program_buffer.py +180 -0
  105. ck/program/raw_program.py +67 -0
  106. ck/sampling/__init__.py +0 -0
  107. ck/sampling/forward_sampler.py +211 -0
  108. ck/sampling/marginals_direct_sampler.py +113 -0
  109. ck/sampling/sampler.py +62 -0
  110. ck/sampling/sampler_support.py +232 -0
  111. ck/sampling/uniform_sampler.py +72 -0
  112. ck/sampling/wmc_direct_sampler.py +171 -0
  113. ck/sampling/wmc_gibbs_sampler.py +153 -0
  114. ck/sampling/wmc_metropolis_sampler.py +165 -0
  115. ck/sampling/wmc_rejection_sampler.py +115 -0
  116. ck/utils/__init__.py +0 -0
  117. ck/utils/iter_extras.py +163 -0
  118. ck/utils/local_config.py +270 -0
  119. ck/utils/map_list.py +128 -0
  120. ck/utils/map_set.py +128 -0
  121. ck/utils/np_extras.py +51 -0
  122. ck/utils/random_extras.py +64 -0
  123. ck/utils/tmp_dir.py +94 -0
  124. ck_demos/__init__.py +0 -0
  125. ck_demos/ace/__init__.py +0 -0
  126. ck_demos/ace/copy_ace_to_ck.py +15 -0
  127. ck_demos/ace/demo_ace.py +49 -0
  128. ck_demos/all_demos.py +88 -0
  129. ck_demos/circuit/__init__.py +0 -0
  130. ck_demos/circuit/demo_circuit_dump.py +22 -0
  131. ck_demos/circuit/demo_derivatives.py +43 -0
  132. ck_demos/circuit_compiler/__init__.py +0 -0
  133. ck_demos/circuit_compiler/compare_circuit_compilers.py +32 -0
  134. ck_demos/circuit_compiler/show_llvm_program.py +26 -0
  135. ck_demos/pgm/__init__.py +0 -0
  136. ck_demos/pgm/demo_pgm_dump.py +18 -0
  137. ck_demos/pgm/demo_pgm_dump_stress.py +18 -0
  138. ck_demos/pgm/demo_pgm_string_rendering.py +15 -0
  139. ck_demos/pgm/show_examples.py +25 -0
  140. ck_demos/pgm_compiler/__init__.py +0 -0
  141. ck_demos/pgm_compiler/compare_pgm_compilers.py +63 -0
  142. ck_demos/pgm_compiler/demo_compiler_dump.py +60 -0
  143. ck_demos/pgm_compiler/demo_factor_elimination.py +47 -0
  144. ck_demos/pgm_compiler/demo_join_tree.py +25 -0
  145. ck_demos/pgm_compiler/demo_marginals_program.py +53 -0
  146. ck_demos/pgm_compiler/demo_mpe_program.py +55 -0
  147. ck_demos/pgm_compiler/demo_pgm_compiler.py +38 -0
  148. ck_demos/pgm_compiler/demo_recursive_conditioning.py +33 -0
  149. ck_demos/pgm_compiler/demo_variable_elimination.py +33 -0
  150. ck_demos/pgm_compiler/demo_wmc_program.py +29 -0
  151. ck_demos/pgm_compiler/time_fe_compiler.py +93 -0
  152. ck_demos/pgm_inference/__init__.py +0 -0
  153. ck_demos/pgm_inference/demo_inferencing_basic.py +188 -0
  154. ck_demos/pgm_inference/demo_inferencing_mpe_cancer.py +45 -0
  155. ck_demos/pgm_inference/demo_inferencing_wmc_and_mpe_sprinkler.py +154 -0
  156. ck_demos/pgm_inference/demo_inferencing_wmc_student.py +110 -0
  157. ck_demos/programs/__init__.py +0 -0
  158. ck_demos/programs/demo_program_buffer.py +24 -0
  159. ck_demos/programs/demo_program_multi.py +24 -0
  160. ck_demos/programs/demo_program_none.py +19 -0
  161. ck_demos/programs/demo_program_single.py +23 -0
  162. ck_demos/programs/demo_raw_program_interpreted.py +21 -0
  163. ck_demos/programs/demo_raw_program_llvm.py +21 -0
  164. ck_demos/sampling/__init__.py +0 -0
  165. ck_demos/sampling/check_sampler.py +71 -0
  166. ck_demos/sampling/demo_marginal_direct_sampler.py +40 -0
  167. ck_demos/sampling/demo_uniform_sampler.py +38 -0
  168. ck_demos/sampling/demo_wmc_direct_sampler.py +40 -0
  169. ck_demos/utils/__init__.py +0 -0
  170. ck_demos/utils/compare.py +120 -0
  171. ck_demos/utils/convert_network.py +45 -0
  172. ck_demos/utils/sample_model.py +216 -0
  173. ck_demos/utils/stop_watch.py +384 -0
  174. compiled_knowledge-4.0.0a20.dist-info/METADATA +50 -0
  175. compiled_knowledge-4.0.0a20.dist-info/RECORD +178 -0
  176. compiled_knowledge-4.0.0a20.dist-info/WHEEL +5 -0
  177. compiled_knowledge-4.0.0a20.dist-info/licenses/LICENSE.txt +21 -0
  178. compiled_knowledge-4.0.0a20.dist-info/top_level.txt +2 -0
ck/example/child.py ADDED
@@ -0,0 +1,199 @@
1
+ from ck.pgm import PGM
2
+
3
+
4
+ class Child(PGM):
5
+
6
+ def __init__(self):
7
+ super().__init__(self.__class__.__name__)
8
+
9
+ pgm_rv0 = self.new_rv('BirthAsphyxia', ('yes', 'no'))
10
+ pgm_rv1 = self.new_rv('HypDistrib', ('Equal', 'Unequal'))
11
+ pgm_rv2 = self.new_rv('HypoxiaInO2', ('Mild', 'Moderate', 'Severe'))
12
+ pgm_rv3 = self.new_rv('CO2', ('Normal', 'Low', 'High'))
13
+ pgm_rv4 = self.new_rv('ChestXray', ('Normal', 'Oligaemic', 'Plethoric', 'Grd_Glass', 'Asy/Patch'))
14
+ pgm_rv5 = self.new_rv('Grunting', ('yes', 'no'))
15
+ pgm_rv6 = self.new_rv('LVHreport', ('yes', 'no'))
16
+ pgm_rv7 = self.new_rv('LowerBodyO2', ('<5', '5-12', '12+'))
17
+ pgm_rv8 = self.new_rv('RUQO2', ('<5', '5-12', '12+'))
18
+ pgm_rv9 = self.new_rv('CO2Report', ('<7.5', '>=7.5'))
19
+ pgm_rv10 = self.new_rv('XrayReport', ('Normal', 'Oligaemic', 'Plethoric', 'Grd_Glass', 'Asy/Patchy'))
20
+ pgm_rv11 = self.new_rv('Disease', ('PFC', 'TGA', 'Fallot', 'PAIVS', 'TAPVD', 'Lung'))
21
+ pgm_rv12 = self.new_rv('GruntingReport', ('yes', 'no'))
22
+ pgm_rv13 = self.new_rv('Age', ('0-3_days', '4-10_days', '11-30_days'))
23
+ pgm_rv14 = self.new_rv('LVH', ('yes', 'no'))
24
+ pgm_rv15 = self.new_rv('DuctFlow', ('Lt_to_Rt', 'None', 'Rt_to_Lt'))
25
+ pgm_rv16 = self.new_rv('CardiacMixing', ('None', 'Mild', 'Complete', 'Transp.'))
26
+ pgm_rv17 = self.new_rv('LungParench', ('Normal', 'Congested', 'Abnormal'))
27
+ pgm_rv18 = self.new_rv('LungFlow', ('Normal', 'Low', 'High'))
28
+ pgm_rv19 = self.new_rv('Sick', ('yes', 'no'))
29
+ pgm_factor0 = self.new_factor(pgm_rv0)
30
+ pgm_factor1 = self.new_factor(pgm_rv1, pgm_rv15, pgm_rv16)
31
+ pgm_factor2 = self.new_factor(pgm_rv2, pgm_rv16, pgm_rv17)
32
+ pgm_factor3 = self.new_factor(pgm_rv3, pgm_rv17)
33
+ pgm_factor4 = self.new_factor(pgm_rv4, pgm_rv17, pgm_rv18)
34
+ pgm_factor5 = self.new_factor(pgm_rv5, pgm_rv17, pgm_rv19)
35
+ pgm_factor6 = self.new_factor(pgm_rv6, pgm_rv14)
36
+ pgm_factor7 = self.new_factor(pgm_rv7, pgm_rv1, pgm_rv2)
37
+ pgm_factor8 = self.new_factor(pgm_rv8, pgm_rv2)
38
+ pgm_factor9 = self.new_factor(pgm_rv9, pgm_rv3)
39
+ pgm_factor10 = self.new_factor(pgm_rv10, pgm_rv4)
40
+ pgm_factor11 = self.new_factor(pgm_rv11, pgm_rv0)
41
+ pgm_factor12 = self.new_factor(pgm_rv12, pgm_rv5)
42
+ pgm_factor13 = self.new_factor(pgm_rv13, pgm_rv11, pgm_rv19)
43
+ pgm_factor14 = self.new_factor(pgm_rv14, pgm_rv11)
44
+ pgm_factor15 = self.new_factor(pgm_rv15, pgm_rv11)
45
+ pgm_factor16 = self.new_factor(pgm_rv16, pgm_rv11)
46
+ pgm_factor17 = self.new_factor(pgm_rv17, pgm_rv11)
47
+ pgm_factor18 = self.new_factor(pgm_rv18, pgm_rv11)
48
+ pgm_factor19 = self.new_factor(pgm_rv19, pgm_rv11)
49
+
50
+ pgm_function0 = pgm_factor0.set_dense()
51
+ pgm_function0.set_flat(0.1, 0.9)
52
+
53
+ pgm_function1 = pgm_factor1.set_dense()
54
+ pgm_function1.set_flat(
55
+ 0.95, 0.95, 0.95, 0.95, 0.95,
56
+ 0.95, 0.95, 0.95, 0.05, 0.5,
57
+ 0.95, 0.5, 0.05, 0.05, 0.05,
58
+ 0.05, 0.05, 0.05, 0.05, 0.05,
59
+ 0.95, 0.5, 0.05, 0.5
60
+ )
61
+
62
+ pgm_function2 = pgm_factor2.set_dense()
63
+ pgm_function2.set_flat(
64
+ 0.93, 0.15, 0.7, 0.1, 0.1,
65
+ 0.1, 0.1, 0.05, 0.1, 0.02,
66
+ 0.1, 0.02, 0.05, 0.8, 0.2,
67
+ 0.8, 0.75, 0.65, 0.7, 0.65,
68
+ 0.5, 0.18, 0.3, 0.18, 0.02,
69
+ 0.05, 0.1, 0.1, 0.15, 0.25,
70
+ 0.2, 0.3, 0.4, 0.8, 0.6,
71
+ 0.8
72
+ )
73
+
74
+ pgm_function3 = pgm_factor3.set_dense()
75
+ pgm_function3.set_flat(
76
+ 0.8, 0.65, 0.45, 0.1, 0.05,
77
+ 0.05, 0.1, 0.3, 0.5
78
+ )
79
+
80
+ pgm_function4 = pgm_factor4.set_dense()
81
+ pgm_function4.set_flat(
82
+ 0.9, 0.14, 0.15, 0.05, 0.05,
83
+ 0.05, 0.05, 0.05, 0.24, 0.03,
84
+ 0.8, 0.01, 0.02, 0.22, 0.02,
85
+ 0.05, 0.15, 0.33, 0.03, 0.02,
86
+ 0.79, 0.15, 0.08, 0.4, 0.05,
87
+ 0.05, 0.03, 0.01, 0.02, 0.04,
88
+ 0.7, 0.5, 0.4, 0.05, 0.05,
89
+ 0.34, 0.03, 0.02, 0.01, 0.08,
90
+ 0.15, 0.13, 0.8, 0.7, 0.06
91
+ )
92
+
93
+ pgm_function5 = pgm_factor5.set_dense()
94
+ pgm_function5.set_flat(
95
+ 0.2, 0.05, 0.4, 0.2, 0.8,
96
+ 0.6, 0.8, 0.95, 0.6, 0.8,
97
+ 0.2, 0.4
98
+ )
99
+
100
+ pgm_function6 = pgm_factor6.set_dense()
101
+ pgm_function6.set_flat(0.9, 0.05, 0.1, 0.95)
102
+
103
+ pgm_function7 = pgm_factor7.set_dense()
104
+ pgm_function7.set_flat(
105
+ 0.1, 0.3, 0.5, 0.4, 0.5,
106
+ 0.6, 0.3, 0.6, 0.4, 0.5,
107
+ 0.45, 0.35, 0.6, 0.1, 0.1,
108
+ 0.1, 0.05, 0.05
109
+ )
110
+
111
+ pgm_function8 = pgm_factor8.set_dense()
112
+ pgm_function8.set_flat(
113
+ 0.1, 0.3, 0.5, 0.3, 0.6,
114
+ 0.4, 0.6, 0.1, 0.1
115
+ )
116
+
117
+ pgm_function9 = pgm_factor9.set_dense()
118
+ pgm_function9.set_flat(
119
+ 0.9, 0.9, 0.1, 0.1, 0.1,
120
+ 0.9
121
+ )
122
+
123
+ pgm_function10 = pgm_factor10.set_dense()
124
+ pgm_function10.set_flat(
125
+ 0.8, 0.1, 0.1, 0.08, 0.08,
126
+ 0.06, 0.8, 0.02, 0.02, 0.02,
127
+ 0.06, 0.02, 0.8, 0.1, 0.1,
128
+ 0.02, 0.02, 0.02, 0.6, 0.1,
129
+ 0.06, 0.06, 0.06, 0.2, 0.7
130
+ )
131
+
132
+ pgm_function11 = pgm_factor11.set_dense()
133
+ pgm_function11.set_flat(
134
+ 0.2, 0.03061224, 0.3, 0.33673469, 0.25,
135
+ 0.29591837, 0.15, 0.23469388, 0.05, 0.05102041,
136
+ 0.05, 0.05102041
137
+ )
138
+
139
+ pgm_function12 = pgm_factor12.set_dense()
140
+ pgm_function12.set_flat(0.8, 0.1, 0.2, 0.9)
141
+
142
+ pgm_function13 = pgm_factor13.set_dense()
143
+ pgm_function13.set_flat(
144
+ 0.95, 0.85, 0.8, 0.7, 0.7,
145
+ 0.25, 0.8, 0.8, 0.8, 0.7,
146
+ 0.9, 0.8, 0.03, 0.1, 0.15,
147
+ 0.2, 0.15, 0.25, 0.15, 0.15,
148
+ 0.15, 0.2, 0.08, 0.15, 0.02,
149
+ 0.05, 0.05, 0.1, 0.15, 0.5,
150
+ 0.05, 0.05, 0.05, 0.1, 0.02,
151
+ 0.05
152
+ )
153
+
154
+ pgm_function14 = pgm_factor14.set_dense()
155
+ pgm_function14.set_flat(
156
+ 0.1, 0.1, 0.1, 0.9, 0.05,
157
+ 0.1, 0.9, 0.9, 0.9, 0.1,
158
+ 0.95, 0.9
159
+ )
160
+
161
+ pgm_function15 = pgm_factor15.set_dense()
162
+ pgm_function15.set_flat(
163
+ 0.15, 0.1, 0.8, 1.0, 0.33,
164
+ 0.2, 0.05, 0.8, 0.2, 0.0,
165
+ 0.33, 0.4, 0.8, 0.1, 0.0,
166
+ 0.0, 0.34, 0.4
167
+ )
168
+
169
+ pgm_function16 = pgm_factor16.set_dense()
170
+ pgm_function16.set_flat(
171
+ 0.4, 0.02, 0.02, 0.01, 0.01,
172
+ 0.4, 0.43, 0.09, 0.16, 0.02,
173
+ 0.03, 0.53, 0.15, 0.09, 0.8,
174
+ 0.95, 0.95, 0.05, 0.02, 0.8,
175
+ 0.02, 0.02, 0.01, 0.02
176
+ )
177
+
178
+ pgm_function17 = pgm_factor17.set_dense()
179
+ pgm_function17.set_flat(
180
+ 0.6, 0.8, 0.8, 0.8, 0.1,
181
+ 0.03, 0.1, 0.05, 0.05, 0.05,
182
+ 0.6, 0.25, 0.3, 0.15, 0.15,
183
+ 0.15, 0.3, 0.72
184
+ )
185
+
186
+ pgm_function18 = pgm_factor18.set_dense()
187
+ pgm_function18.set_flat(
188
+ 0.3, 0.2, 0.15, 0.1, 0.3,
189
+ 0.7, 0.65, 0.05, 0.8, 0.85,
190
+ 0.1, 0.1, 0.05, 0.75, 0.05,
191
+ 0.05, 0.6, 0.2
192
+ )
193
+
194
+ pgm_function19 = pgm_factor19.set_dense()
195
+ pgm_function19.set_flat(
196
+ 0.4, 0.3, 0.2, 0.3, 0.7,
197
+ 0.7, 0.6, 0.7, 0.8, 0.7,
198
+ 0.3, 0.3
199
+ )
ck/example/clique.py ADDED
@@ -0,0 +1,33 @@
1
+ import random as _random
2
+
3
+ from ck.pgm import PGM
4
+
5
+
6
+ class Clique(PGM):
7
+ """
8
+ This PGM is the 'Clique' factor graph.
9
+
10
+ The Clique factor graph has one factor that connects to all random variables.
11
+ Random variables are named x0, x1, x2, ... .
12
+ If include_unaries then, also includes one unary factor per random variable.
13
+ """
14
+
15
+ def __init__(
16
+ self,
17
+ vars_per_clique: int = 17,
18
+ states_per_var: int = 2,
19
+ include_unaries: bool = True,
20
+ random_seed: int = 123456,
21
+ ):
22
+ params = (vars_per_clique, states_per_var, include_unaries)
23
+ super().__init__(f'{self.__class__.__name__}({",".join(str(param) for param in params)})')
24
+
25
+ random_stream = _random.Random(random_seed).random
26
+
27
+ rvs = [self.new_rv(f'x{i}', states_per_var) for i in range(vars_per_clique)]
28
+
29
+ self.new_factor(*rvs).set_dense().set_stream(random_stream)
30
+
31
+ if include_unaries:
32
+ for rv in rvs:
33
+ self.new_factor(rv).set_dense().set_stream(random_stream)
ck/example/cnf_pgm.py ADDED
@@ -0,0 +1,39 @@
1
+ import random as _random
2
+
3
+ from ck.pgm import PGM
4
+ from ck.utils.random_extras import random_permutation as _random_permutation
5
+
6
+
7
+ class CNF_PGM(PGM):
8
+ """
9
+ This is a PGM that uses a factor graph to represent a
10
+ CNF propositional theory.
11
+
12
+ A CNF_PGM factor graph has random clause factors.
13
+ Random variables are named x0, x1, x2, ... .
14
+ """
15
+
16
+ def __init__(
17
+ self,
18
+ num_vars: int = 30,
19
+ num_clauses: int = 200,
20
+ min_clause_size: int = 3,
21
+ max_clause_size: int = 5,
22
+ random_seed: int = 129837697,
23
+ ):
24
+ params = (num_vars, num_clauses, min_clause_size, max_clause_size)
25
+ super().__init__(f'{self.__class__.__name__}({",".join(str(param) for param in params)})')
26
+
27
+ random = _random.Random(random_seed)
28
+
29
+ for i in range(num_vars):
30
+ self.new_rv(f'x{i}', 2)
31
+
32
+ pgm_rvs = self.rvs
33
+ for i in range(num_clauses):
34
+ size = random.randint(min_clause_size, max_clause_size)
35
+ perm = _random_permutation(num_vars, random)
36
+ rvs = tuple(pgm_rvs[j] for j in perm[:size])
37
+ states = [random.randint(0, 1) for _ in range(size)]
38
+
39
+ self.new_factor(*rvs).set_clause(*states)
@@ -0,0 +1,68 @@
1
+ import random as _random
2
+
3
+ from ck.pgm import PGM
4
+
5
+
6
+ class DiamondSquare(PGM):
7
+ r"""
8
+ This PGM is the 'DiamondSquare' factor graph.
9
+
10
+ The DiamondSquare is a factor graph with seven random variables (a, b, c, ..., h).
11
+ Binary factors are between pairs of random variables creating the pattern:
12
+ b
13
+ / \
14
+ / \
15
+ c -- d
16
+ / | | \
17
+ a | | h
18
+ \ | | /
19
+ e -- f
20
+ \ /
21
+ \ /
22
+ g
23
+ If include_unaries then, also includes one unary factor per random variable.
24
+ """
25
+
26
+ def __init__(
27
+ self,
28
+ states_per_var: int = 5,
29
+ include_unaries: bool = True,
30
+ random_seed: int = 123456,
31
+ ):
32
+ params = (states_per_var, include_unaries)
33
+ super().__init__(f'{self.__class__.__name__}({",".join(str(param) for param in params)})')
34
+
35
+ random_stream = _random.Random(random_seed).random
36
+ rand_iter = iter(random_stream, None)
37
+
38
+ a = self.new_rv('a', states_per_var)
39
+ b = self.new_rv('b', states_per_var)
40
+ c = self.new_rv('c', states_per_var)
41
+ d = self.new_rv('d', states_per_var)
42
+ e = self.new_rv('e', states_per_var)
43
+ f = self.new_rv('f', states_per_var)
44
+ g = self.new_rv('g', states_per_var)
45
+ h = self.new_rv('h', states_per_var)
46
+
47
+ self.new_factor(a, c).set_dense().set_iter(rand_iter)
48
+ self.new_factor(a, e).set_dense().set_iter(rand_iter)
49
+ self.new_factor(b, c).set_dense().set_iter(rand_iter)
50
+ self.new_factor(b, d).set_dense().set_iter(rand_iter)
51
+ self.new_factor(g, e).set_dense().set_iter(rand_iter)
52
+ self.new_factor(g, f).set_dense().set_iter(rand_iter)
53
+ self.new_factor(h, d).set_dense().set_iter(rand_iter)
54
+ self.new_factor(h, f).set_dense().set_iter(rand_iter)
55
+ self.new_factor(c, d).set_dense().set_iter(rand_iter)
56
+ self.new_factor(c, e).set_dense().set_iter(rand_iter)
57
+ self.new_factor(d, f).set_dense().set_iter(rand_iter)
58
+ self.new_factor(e, f).set_dense().set_iter(rand_iter)
59
+
60
+ if include_unaries:
61
+ self.new_factor(a).set_dense().set_iter(rand_iter)
62
+ self.new_factor(b).set_dense().set_iter(rand_iter)
63
+ self.new_factor(c).set_dense().set_iter(rand_iter)
64
+ self.new_factor(d).set_dense().set_iter(rand_iter)
65
+ self.new_factor(e).set_dense().set_iter(rand_iter)
66
+ self.new_factor(f).set_dense().set_iter(rand_iter)
67
+ self.new_factor(g).set_dense().set_iter(rand_iter)
68
+ self.new_factor(h).set_dense().set_iter(rand_iter)
@@ -0,0 +1,36 @@
1
+ from ck.pgm import PGM
2
+
3
+
4
+ class Earthquake(PGM):
5
+
6
+ def __init__(self):
7
+ super().__init__(self.__class__.__name__)
8
+
9
+ pgm_rv0 = self.new_rv('Burglary', ('True', 'False'))
10
+ pgm_rv1 = self.new_rv('Earthquake', ('True', 'False'))
11
+ pgm_rv2 = self.new_rv('Alarm', ('True', 'False'))
12
+ pgm_rv3 = self.new_rv('JohnCalls', ('True', 'False'))
13
+ pgm_rv4 = self.new_rv('MaryCalls', ('True', 'False'))
14
+ pgm_factor0 = self.new_factor(pgm_rv0)
15
+ pgm_factor1 = self.new_factor(pgm_rv1)
16
+ pgm_factor2 = self.new_factor(pgm_rv2, pgm_rv0, pgm_rv1)
17
+ pgm_factor3 = self.new_factor(pgm_rv3, pgm_rv2)
18
+ pgm_factor4 = self.new_factor(pgm_rv4, pgm_rv2)
19
+
20
+ pgm_function0 = pgm_factor0.set_dense()
21
+ pgm_function0.set_flat(0.01, 0.99)
22
+
23
+ pgm_function1 = pgm_factor1.set_dense()
24
+ pgm_function1.set_flat(0.02, 0.98)
25
+
26
+ pgm_function2 = pgm_factor2.set_dense()
27
+ pgm_function2.set_flat(
28
+ 0.95, 0.94, 0.29, 0.001, 0.05,
29
+ 0.06, 0.71, 0.999
30
+ )
31
+
32
+ pgm_function3 = pgm_factor3.set_dense()
33
+ pgm_function3.set_flat(0.9, 0.05, 0.1, 0.95)
34
+
35
+ pgm_function4 = pgm_factor4.set_dense()
36
+ pgm_function4.set_flat(0.7, 0.01, 0.3, 0.99)
ck/example/empty.py ADDED
@@ -0,0 +1,10 @@
1
+ from ck.pgm import PGM
2
+
3
+
4
+ class Empty(PGM):
5
+ """
6
+ This PGM has no random variables.
7
+ """
8
+
9
+ def __init__(self):
10
+ super().__init__(self.__class__.__name__)