Subject | Hash | Author | Date (UTC) |
---|---|---|---|
voronoi: WIP | c3bfa60b6db2042df11528fe95a51d7cc88d3c6b | I am | 2022-10-19 11:18:13 |
voronoi: WIP | baff3c5ccaf857f60194917fe5d6e33c40b01265 | I am | 2022-10-19 09:28:33 |
voronoi.ConvexSpline: little fix of np.sign() for better sleep | 968f649e54d97abd596cd60749872fa3ef788c06 | I am | 2022-10-19 06:28:46 |
implement even more optimized Convex.. now ConvexSpline class | c08511aecb32b880e73c9f10d07feff4dc588e8b | I am | 2022-10-19 04:11:13 |
voronoi: WIP | cbb61d4544367a67c9a5c3fe6d64a9dc284eb364 | I am | 2022-10-18 05:05:44 |
voronoi: WIP | e9c247fef30cfb678bad9b955d6f6a94a0ff61e7 | I am | 2022-10-18 01:02:30 |
candynodex: hotfix to support numpy masking indexing | f359123b72a3d998d5eed4b9e611a5e402a18c8d | I am | 2022-10-17 09:02:33 |
add new lightweight candynodes module to replace old heavy ugly candybox | ca97509f47c100df90e1e06d82ed7c759627bfd0 | I am | 2022-10-17 06:31:33 |
qt_gui: move CandyBox creating to the individual DiceBox widgets | 29253625129219b5e550f82cae3da0e53ef5bd47 | I am | 2022-10-16 03:17:55 |
voronoi: WIP | 4855190869e06df740cc05c64e8b27b2c3c5a88d | I am | 2022-10-15 22:46:51 |
voronoi: WIP | 0298b3a0563587ea4fc24fa99b51d2e01d517203 | I am | 2022-10-15 20:00:25 |
voronoi: WIP | 4b80757fa26e3d40dceac888480ec3d814f6abc9 | I am | 2022-10-15 17:30:11 |
voronoi: WIP | 7cadd1d01d7f4bc3f51e785fe6ebf8cae6ff839a | I am | 2022-10-15 03:12:54 |
voronoi: clean up and comment out ContactSolver | 543f8482a925ff6cf8485b9be616b5bd3c714d13 | I am | 2022-10-14 04:24:40 |
qt_gui.qt_gui: use ConvexSolver instead of ContactSolver (which did not work correctly) | c1e49a9b2c3cc13886d317b460d7598f661cf216 | I am | 2022-10-14 04:01:06 |
voronoi: add fantastic ConvexSolver | 2853826e5ef1abc79d3ac2fb8289b13c45211a31 | I am | 2022-10-14 04:00:14 |
qt_gui.qt_plot: add (finally!) Numbers class | 4cdc658c0fcc857c23c0e39e91e7a1ce5e1b30a1 | I am | 2022-10-13 06:23:46 |
qt_gui.qt_gui: show contacts in distance matrix. Based on ContactSolver for now | 21bd6101888d9f06d7d6a7c6ba2732ff30fdd68d | I am | 2022-10-13 04:16:41 |
voronoi: ContactSolver is ready. Ale je to na nic. Pomalá, potvora | 28e2442b0101eac2407ed4a69c6f757ffd579cf1 | I am | 2022-10-13 04:14:59 |
voronoi: add some preliminary lines of ContactSolver. WIP | 6a203b278c9fa6de7da75c402c80f78d40164fdf | I am | 2022-10-12 23:38:59 |
File | Lines added | Lines deleted |
---|---|---|
wellmet/voronoi.py | 28 | 2 |
File wellmet/voronoi.py changed (mode: 100644) (index f642ead..310c8bc) | |||
... | ... | class ContactVoronoi: | |
426 | 426 | IS - kusy mají odlíšnou hmotnost | IS - kusy mají odlíšnou hmotnost |
427 | 427 | ||
428 | 428 | ||
429 | |||
429 | č base_r dává poloměr, do kterého/za který spadá jistý podíl bodíků | ||
430 | ' | ||
431 | r_ball/r_base = sigma | ||
432 | ' | ||
433 | č Když máme nějakou oblast s poloměrem alespoň r | ||
434 | č můžeme nastavit směrodatnou odchylku normálního rozdělení tak, | ||
435 | č aby jisté procento bodíků dopadalo dovnitř | ||
436 | č (integrujeme-li obyčejnou Voroneho buňku, tak r = mindist / 2) | ||
437 | č Co ale máme dělat, když máme výběrovou směrodatnou odchylku | ||
438 | č bodíků z cenzurovanného rozdělení? Survival bias. | ||
439 | č Můžu spočítat r_ball, ten nebude korrettní, protože vzorky | ||
440 | č nebyly z normálního rozdělení, ale z cenzurovanného | ||
441 | č A ani ten r_ball nepotřebuji. Chcu nastavit vzorkovací sigmu. | ||
442 | č Co jako? Mám dělat výzkum? Počítat třetí a čtverté výběrové momenty? | ||
443 | č Před tím se dělalo tohle: sigmas = sigmy / base_r | ||
444 | č Což, pokud rozumím, ve 2D ještě víc změnšovalo směrodatné odchylky | ||
445 | ' | ||
446 | č Mám výberovou "sigmu" (std) z cenzurovanného rozdělení, | ||
447 | č podle mé osy má odříznuté chvosty. | ||
448 | č Odříznuté chvosty podel jinejch os mě nezajimají | ||
449 | č std * r_base(ndim) = r_ball(ndim) | ||
450 | č sigma * r_base(ndim-1) = r_ball(ndim-1) | ||
451 | č Nabízím opravu: sigma = std * r_base(ndim) / r_base(ndim-1) | ||
452 | č Ve výsokých dimenzích rozdíl je minimální | ||
453 | č Ale kdo tuší, co a jak tam funguje a má fungovat? | ||
430 | 454 | ||
431 | 455 | ||
432 | 456 | č Score počítáme jako pravděpodobnostní obsah buňky dělený | č Score počítáme jako pravděpodobnostní obsah buňky dělený |
... | ... | class ContactVoronoi: | |
450 | 474 | self.on_delete_mixed = on_delete_mixed | self.on_delete_mixed = on_delete_mixed |
451 | 475 | ||
452 | 476 | #č chcu, aby behem dalších iterací pulka (p_base) tečiček dopadala do buňky | #č chcu, aby behem dalších iterací pulka (p_base) tečiček dopadala do buňky |
453 | self.base_r = sball.Sball(sample_box.nvar).get_r(1 - p_base) | ||
477 | self.base_r = sball.get_Radius_ps(sample_box.nvar, p_base) | ||
478 | #č pro funkci IS_stat.sample_like(), víz. diskuse v poznamkách | ||
479 | self.d = (sball.get_Radius_ps(sample_box.nvar - 1, p_base) / self.base_r)**2 | ||
454 | 480 | ||
455 | 481 | ||
456 | 482 | self._nsim = 0 | self._nsim = 0 |