Subject | Hash | Author | Date (UTC) |
---|---|---|---|
simplex.CubatureIntegration: check for non-finite simplex volumes | 3a1473c7452c7dc3a9299ed88520412eb31924c1 | I am | 2023-01-02 22:51:14 |
dicebox.circumtri: return to psi potential | 9db1f5b0e412ae871ae5f5fbb5ff214a8ddd7f45 | I am | 2023-01-02 22:20:28 |
dicebox.circumtri: bughunting | 5f9d2b9e4d8077c9e4552d47bf1da789d8cfbcf8 | I am | 2022-12-29 15:02:14 |
dicebox.circumtri: process filtered integration nodes as candidates too | 315e2f9700d93f1379985dafa818b44df53422a2 | I am | 2022-12-29 13:51:52 |
dicebox.circumtri: return back weighted entropy | 0424243e29a600ecb1d046631badbc9797e11369 | I am | 2022-12-28 19:11:30 |
qt_gui.qt_graph_widgets: reflect renewed simplex estimators | 90c83c7c7e6e58d6b9044db1a5e55f9a8a7dd7cd | I am | 2022-12-28 16:15:40 |
dicebox: rework CircumTri | e360b2f092c11a42b5ac1081061f2891bff3d8de | I am | 2022-12-28 15:52:28 |
simplex: make _Triangulation more clean; rework CubatureIntegration | fb7dd9fef8592944fbd93a3f07b44c2355d45e20 | I am | 2022-12-28 15:48:20 |
introduce shell module, clean version of ghull | b13e9b2b15d109770e0c3fbdd5bb7b2a3b5bf5fc | I am | 2022-12-28 15:45:27 |
qt_gui.qt_plot.CandidatesWidget: regression fix | a75874b5a53e7b420672b582f3ff415b17a32e9b | I am | 2022-12-23 12:19:51 |
qt_gui.qt_pairwise.ContactWidget: add mask function | 6f88fec7db92e5cd21826a06f42a609a3e03cb6b | I am | 2022-12-21 15:20:01 |
qt_gui.qt_pairwise.ContactWidget: add option for mixed pairs only adjacency search. | 2f87e8a18908de607b8abc0779448ca3095c134c | I am | 2022-12-21 14:30:23 |
qt_gui.qt_dicebox: reflect changes in CircumTri class | 9821d389240a383f0b13b7dd92cd501e1c40c537 | I am | 2022-12-15 09:36:16 |
dicebox.circumtri: employ _Exploration class. Keep naive CircumTri class in renamed module | b8605fb24d55b3296ceac34afba541a5b7ee3070 | I am | 2022-12-15 09:33:59 |
convex_hull.QHull.get_exploration_vector: self -> hull fix | dae47ef37b5871a98d2ea274cc35aeb2a4eebe1c | I am | 2022-12-15 09:30:31 |
qt_gui/qt_dicebox: add DumpExplorationWidget | b687895f06c780418a12eea59bd052ad7cd348f3 | I am | 2022-12-15 01:31:14 |
dicebox: add helper _exploration module | 836ae6cd12f1a04f1e2fc3423e46777d16ebad7d | I am | 2022-12-15 01:30:20 |
convex_hull: add get_exploration_vector() method | 795259d91b64c732e60a83a1d9e3b89c2141b003 | I am | 2022-12-15 01:28:12 |
qt_gui.qt_gui: dot not create graph widgets if estimations are not present | f070ed140a135a980e4436ed3779455c0e1bf399 | I am | 2022-12-14 02:25:45 |
qt_gui.gl_plot.CandidatesWidget: apply changes for GLplot too | 03550040b60c67d247e26425e35fac565dd88886 | Aleksei Gerasimov | 2022-12-14 01:56:01 |
File | Lines added | Lines deleted |
---|---|---|
wellmet/simplex.py | 11 | 3 |
File wellmet/simplex.py changed (mode: 100644) (index 97a9d84..b46897b) | |||
... | ... | class CubatureIntegration(_Triangulation): | |
1237 | 1237 | x = quadpy.tn.transform(sx.tn_scheme.points, vertices_model.T) | x = quadpy.tn.transform(sx.tn_scheme.points, vertices_model.T) |
1238 | 1238 | vol = quadpy.tn.get_vol(vertices_model) | vol = quadpy.tn.get_vol(vertices_model) |
1239 | 1239 | ||
1240 | if not np.isfinite(vol): | ||
1241 | print("Incorrect volume has occured in simplex %s" % indices) | ||
1242 | vol = 0 | ||
1243 | |||
1240 | 1244 | nodes = sx.f.new_sample(x.T, sx.tri_space) | nodes = sx.f.new_sample(x.T, sx.tri_space) |
1241 | 1245 | fx = nodes.pdf(sx.tri_space) | fx = nodes.pdf(sx.tri_space) |
1242 | 1246 | ||
1243 | 1247 | p_mixed = vol * np.dot(fx, sx.tn_scheme.weights) | p_mixed = vol * np.dot(fx, sx.tn_scheme.weights) |
1244 | 1248 | ||
1245 | 1249 | ||
1246 | if p_mixed < 0: | ||
1247 | print("Negative measure has occured in simplex %s" % indices) | ||
1248 | print("Rozbíla se nám integrace, integráční schema je na houby.") | ||
1250 | if not (p_mixed > 0): | ||
1251 | print("Incorrect measure %s has occured in simplex %s" % (p_mixed, indices)) | ||
1252 | #print("Rozbíla se nám integrace, integráční schema je na houby.") | ||
1249 | 1253 | p_mixed = 0 | p_mixed = 0 |
1250 | 1254 | ||
1251 | 1255 | ||
... | ... | class CubatureIntegration(_Triangulation): | |
1272 | 1276 | x = quadpy.tn.transform(sx.tn_scheme.points, vertices_model.T) | x = quadpy.tn.transform(sx.tn_scheme.points, vertices_model.T) |
1273 | 1277 | vol = quadpy.tn.get_vol(vertices_model) | vol = quadpy.tn.get_vol(vertices_model) |
1274 | 1278 | ||
1279 | if not np.isfinite(vol): | ||
1280 | print("Incorrect volume has occured in simplex %s" % indices) | ||
1281 | vol = 0 | ||
1282 | |||
1275 | 1283 | nodes = sx.f.new_sample(x.T, sx.tri_space) | nodes = sx.f.new_sample(x.T, sx.tri_space) |
1276 | 1284 | fx = nodes.pdf(sx.tri_space) | fx = nodes.pdf(sx.tri_space) |
1277 | 1285 |