Subject | Hash | Author | Date (UTC) |
---|---|---|---|
dicebox.circumtri: tune up exploratory rule one more time | bfaa8d65bd13956252a6a25382c621aca7a33e3f | I am | 2023-02-06 12:05:00 |
convex_hull.QHull.is_outside: mark everything as outside until convex hull is created | beecd924a63603868807a8a603cbeb04857e5cde | I am | 2023-02-06 12:03:35 |
dicebox.circumtri: get_exploratory_radius: turn it up some more | ecd638ca90bac88df515cf803a103e350b08041e | I am | 2023-02-03 10:58:37 |
mplot: ongoing ad-hoc changes | fd71b93abee959231ffd76324cf2fbc0b3bbbdfe | I am | 2023-02-02 17:59:27 |
dicebox.circumtri: change default q setting | 2a9bdc6604113462a935d409bca429b49ac52cbc | I am | 2023-02-02 11:55:41 |
dicebox.circumtri: minor change in explore radia | 9aa06dc11be7c935dd47715c0aa53b401be951d4 | I am | 2023-02-01 23:37:13 |
dicebox.circumtri: exploratory: another iteration | 78c849a2c1b07a55f87d32e3cd6570db222a50b6 | I am | 2023-02-01 01:02:00 |
dicebox: new experimental explore() | 1d570ef1dcc66b9b1015bd6c1f6877e1405e676c | I am | 2023-02-01 00:34:57 |
mplot: improve GRaph; draw vertex estimator only in convergence diagram | 5c22e808e4c6cabd8e0d077d5f295eea93957582 | I am | 2023-01-31 18:02:00 |
simplex.TopologyAnalysis: implement postprocessing | 5b8d6c87bd18e06e2e3e1aa651ebccc6d9ed2e69 | I am | 2023-01-30 20:59:59 |
testcases.testcases_nD_papers: add low pf nonlinear oscillator setting | 4a75d5ec830f88c35a91e53898e416e568e3913c | I am | 2023-01-30 16:13:17 |
qt_gui.qt_plot: replace dropped np.int by int | 197cbc3f6cf7ef92fc092a1e026cbace609eaabe | I am | 2023-01-28 18:07:43 |
dicebox.circumtri.increment: treat sample as an exploration sample if lies in a critical direction | 95158c4f2e450ddc9bf3e4d1c4151e33d901c4b3 | I am | 2023-01-28 12:13:54 |
simplex.TopologyAnalysis: integration fix | 5e520d41715c5260dd297b1356de9ea17ebef4ba | I am | 2023-01-27 21:22:43 |
simplex: preliminary work on topology analysis | 72aeb80bd3469fe8ed0ec56e1cb6bbb83ecfc396 | I am | 2023-01-26 21:52:08 |
dicebox.circumtri.QTri: use different cubature degrees for facets and simplices inside | aa77a812954891fae171db6c5a6b25ee0aac2275 | I am | 2023-01-25 12:08:20 |
convex_hull.QHull.get_chi_outside: count negative-valued facets | c3f645977e19d66f7a402001485254ba30fdb2ce | I am | 2023-01-25 10:22:00 |
convex_hull.QHull.get_chi_outside: check and mask negative valued integrals | 861f9bfdf5dab31c5c6b81566eb8db7d338307e9 | I am | 2023-01-24 22:27:43 |
dicebox.circumtri: estimation export hotfix | 4f86660297b643446bacd3e254ded9faca148542 | I am | 2023-01-24 15:02:34 |
dicebox.circumtri: assemble QTri box | 9e2114c0c6d5defcf3509a30fda3b6230eb9b9d1 | I am | 2023-01-24 14:48:48 |
File | Lines added | Lines deleted |
---|---|---|
wellmet/dicebox/circumtri.py | 7 | 4 |
File wellmet/dicebox/circumtri.py changed (mode: 100644) (index 66bbdcb..8011ccd) | |||
... | ... | class _CircumTri(_Exploration): | |
104 | 104 | bx.to_explore = 1 | bx.to_explore = 1 |
105 | 105 | bx.to_refine = 0 | bx.to_refine = 0 |
106 | 106 | if bx.nsim > 0: | if bx.nsim > 0: |
107 | bx.get_pf_estimation() #č updates p_mixed | ||
107 | bx.pf_estimation = bx.get_pf_estimation() #č updates p_mixed | ||
108 | 108 | bx.update_exploration_ratio() | bx.update_exploration_ratio() |
109 | 109 | ||
110 | 110 | bx._nsim = bx.nsim | bx._nsim = bx.nsim |
... | ... | class _CircumTri(_Exploration): | |
155 | 155 | ||
156 | 156 | def update_exploration_ratio(bx): | def update_exploration_ratio(bx): |
157 | 157 | bx.to_explore += np.sqrt(bx.shell_estimation.outside) | bx.to_explore += np.sqrt(bx.shell_estimation.outside) |
158 | bx.to_refine += np.sqrt(bx.p_mixed) | ||
158 | failure_share = bx.pf_estimation.vertex_estimation**(1/2 - 1/2 / bx.nvar) | ||
159 | mixed_share = bx.p_mixed**(1/2 / bx.nvar) | ||
160 | bx.to_refine += failure_share * mixed_share | ||
159 | 161 | ||
160 | 162 | def get_exploratory_radius(bx): | def get_exploratory_radius(bx): |
161 | 163 | R = bx.shell_stats.R | R = bx.shell_stats.R |
162 | 164 | # empirical rule to get desired behavior | # empirical rule to get desired behavior |
163 | sphere_area = bx.convex_hull.nsphere_surface_area * (R + (bx.q - 1)**2)**bx.nvar | ||
165 | sphere_area = bx.convex_hull.nsphere_surface_area * (R + bx.q)**bx.nvar | ||
164 | 166 | p_desired = bx.shell_estimation.outside / sphere_area | p_desired = bx.shell_estimation.outside / sphere_area |
165 | 167 | ||
166 | 168 | # get matematically clean radius of it | # get matematically clean radius of it |
... | ... | class _CircumTri(_Exploration): | |
365 | 367 | else: | else: |
366 | 368 | bx.to_refine = 0 | bx.to_refine = 0 |
367 | 369 | ||
368 | bx.box_estimations.append(bx.get_pf_estimation()) | ||
370 | bx.pf_estimation = bx.get_pf_estimation() | ||
371 | bx.box_estimations.append(bx.pf_estimation) | ||
369 | 372 | bx.update_exploration_ratio() | bx.update_exploration_ratio() |
370 | 373 | ||
371 | 374 |