iam-git / WellMet (public) (License: MIT) (since 2021-08-31) (hash sha1)
WellMet is pure Python framework for spatial structural reliability analysis. Or, more specifically, for "failure probability estimation and detection of failure surfaces by adaptive sequential decomposition of the design domain".
Clone URLs: https://rocketgit.com/user/iam-git/WellMet ssh://rocketgit@ssh.rocketgit.com/user/iam-git/WellMet git://git.rocketgit.com/user/iam-git/WellMet
release trunk
List of commits:
Subject Hash Author Date (UTC)
dicebox.circumtri.CirQTri: switch to GaussCubatureIntegration 5b52dd25cb7a997af4953230116deb9efc577d56 I am 2023-02-11 04:32:48
simplex: implement GaussCubatureIntegration in the most memory-friendly way 689d253ae7e2a22242258fd5bef0a069caf7cf75 I am 2023-02-11 04:31:11
convex_hull.QHullCubature: implement memory-friendly outside (chi) integration ad8210a04b1e0903de7435cad16b1304707d0e6e I am 2023-02-09 22:22:05
qt_gui.qt_plot: require box recommendation to show next point 6f726047f7f08e884359020eaa1eac6f6cc125d2 I am 2023-02-09 11:51:44
dicebox.circumtri.CirQTri.get_circum_node: limit circumcenter by explore radius, not by just R 136ec73bb06da16c1f2bce64b3c349be4c8ba975 I am 2023-02-09 03:09:51
dicebox.circumtri: implement new brand CirQTri box 5879b8ad6317c65aa9481b59f76b6159f19e247a I am 2023-02-09 01:29:10
simplex.FullCubatureIntegration: store simplex probabilities in sorted dicts c0da90731ff3ede47d9b4eec0ad9b28a29027167 I am 2023-02-09 01:23:14
dicebox.circumtri: exploratory: even better idea 811ab11cd7172ff4a3807992f4928be2e8068ec0 I am 2023-02-08 15:31:23
dicebox.circumtri: exploratory, new idea 526d3f6185887ff48a341b0705d74dde4d15ca87 I am 2023-02-08 03:03:41
dicebox.circumtri: exploratory 806063d2e223c812280dc4845153450dd47faed3 I am 2023-02-06 17:15:15
dicebox.circumtri: exploratory efed2589f642d502f30e80f0e9b45dfeecd1c7c7 I am 2023-02-06 13:40:24
dicebox.circumtri: exploratory - again 34d3f4e47420e1c1e26b09570fb44d3270194751 I am 2023-02-06 12:50:45
qt_gui.qt_dicebox: change default q of circumtri classes 9fd5855e5d7cacf80d27fb383dd18a92d60e138b I am 2023-02-06 12:30:27
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
Commit 5b52dd25cb7a997af4953230116deb9efc577d56 - dicebox.circumtri.CirQTri: switch to GaussCubatureIntegration
Author: I am
Author date (UTC): 2023-02-11 04:32
Committer name: I am
Committer date (UTC): 2023-02-11 04:32
Parent(s): 689d253ae7e2a22242258fd5bef0a069caf7cf75
Signer:
Signing key:
Signing status: N
Tree: 0b9b262fe20c613f13cc13dad5f5a6e97141500e
File Lines added Lines deleted
wellmet/dicebox/circumtri.py 72 76
File wellmet/dicebox/circumtri.py changed (mode: 100644) (index f357382..8290e2f)
... ... class _CircumTri(_Exploration):
117 117 r = bx.sball.get_r(p_desired) r = bx.sball.get_r(p_desired)
118 118 return r return r
119 119
120
121
122
123
124
125
126
127
128
129
130
131 def _regen_inside(bx):
132 failsi = bx.failsi
133 # incremental triangulation require one more point
134 if (bx.nsim > bx.nvar + 1) and np.any(failsi) and not np.all(failsi):
135 #bx._logger(msg="triangulation started")
136 bx.__regen_inside()
137 else:
138 #č jíž není nutný
139 bx._logger(msg="triangulation skipped")
140
141 def __regen_inside(bx):
142 try:
143 if bx.holydays:
144 Integrator = sx.FullCubatureIntegration
145 else:
146 Integrator = sx.FastCubatureIntegration
147
148 bx.Tri = Integrator(bx.samplebox, bx.scheme,\
149 tri_space=bx.tri_space, incremental=True,\
150 on_safe_added=bx._on_safe_added,\
151 on_mixed_added=bx._on_mixed_added,\
152 on_delete_simplex=bx._invalidate_simplex)
153
154 bx.Tri.integrate() # nic nevrácí, všecko je přes kolbeky
155 #č tri - Deloneho triangulace
156 bx.tri = bx.Tri.tri #č všichní tam očekávajou QHull
157 bx._logger(msg="triangulation has been created")
158
159 except BaseException as e:
160 #č chcu zachytit spadnuti QHull na začatku,
161 #č kdy ještě není dostatek teček.
162 #č Jinak je třeba nechat QHull spadnout
163 if bx.nsim > 2*bx.nvar + 3:
164 #č no to teda ne!
165 raise
166 else:
167 #č lze přípustit chybu triangulace
168 bx._logger(msg='triangulation failed')
169
170 120
171 121 def is_mixed(bx, simplices=None): def is_mixed(bx, simplices=None):
172 122
 
... ... class CirQTri(_CircumTri):
404 354 return bx.shell_estimation[1:] return bx.shell_estimation[1:]
405 355
406 356
407 def get_circum_node(bx, simplex):
408 indices = np.array(simplex)
357 def _regen_inside(bx):
358 failsi = bx.failsi
359 # incremental triangulation require one more point
360 if (bx.nsim > bx.nvar + 1) and np.any(failsi) and not np.all(failsi):
361 try:
362 bx.Tri = sx.GaussCubatureIntegration(bx.samplebox, bx.scheme,
363 incremental=True, full=bx.holydays)
364
365 #č tri - Deloneho triangulace
366 bx.tri = bx.Tri.tri #č všichni tam očekávají QHull
367 bx._logger(msg="triangulation has been created")
368
369 except BaseException as e:
370 #č chcu zachytit spadnuti QHull na začatku,
371 #č kdy ještě není dostatek teček.
372 #č Jinak je třeba nechat QHull spadnout
373 if bx.nsim > 2*bx.nvar + 3:
374 #č no to teda ne!
375 raise
376 else:
377 #č lze přípustit chybu triangulace
378 bx._logger(msg='triangulation failed')
379 else:
380 #č jíž není nutný
381 bx._logger(msg="triangulation skipped")
382
383
384 def get_circum_node(bx, indices):
409 385 vertices_model = bx.Tri.tri.points[indices] vertices_model = bx.Tri.tri.points[indices]
410 386 try: try:
411 387 circum_center = bx.CC.get_circumcenter(vertices_model) circum_center = bx.CC.get_circumcenter(vertices_model)
 
... ... class CirQTri(_CircumTri):
431 407 def refine(bx): def refine(bx):
432 408 if "Tri" in dir(bx): if "Tri" in dir(bx):
433 409 # get the greatest mixed simplex # get the greatest mixed simplex
434 key, __value = bx.Tri.mixed_simplex_probabilities.peekitem()
435 return bx.get_circum_node(key)
410 return bx.get_circum_node(bx.Tri.max_mixed_indices)
436 411 else: else:
437 412 bx._logger(msg='refine called, but triangulation does not exist yet. Fallback to a random sample') bx._logger(msg='refine called, but triangulation does not exist yet. Fallback to a random sample')
438 413 return bx.f_model(1) return bx.f_model(1)
439 414
440 415
441 416 def holyday(bx): def holyday(bx):
442 if "Tri" in dir(bx) and len(bx.Tri.safe_simplex_probabilities):
443 # return node with the greatest potential
444 key, __value = bx.Tri.safe_simplex_probabilities.peekitem()
445 return bx.get_circum_node(key)
417 if "Tri" in dir(bx) and (bx.Tri.max_safe > 0):
418 return bx.get_circum_node(bx.Tri.max_safe_indices)
446 419 else: else:
447 420 bx._logger(msg='Fallback to random sample') bx._logger(msg='Fallback to random sample')
448 421 return bx.f_model(1) return bx.f_model(1)
 
... ... class CirQTri(_CircumTri):
471 444
472 445 def get_max_potential(bx): def get_max_potential(bx):
473 446 if "Tri" in dir(bx): if "Tri" in dir(bx):
474 # get the greatest mixed simplex
475 __key, value = bx.Tri.mixed_simplex_probabilities.peekitem()
476 return value
447 # get probability of the greatest mixed simplex
448 return bx.Tri.max_mixed
477 449 else: else:
478 450 return -1 return -1
479 451
 
... ... class CirQTri(_CircumTri):
481 453 bx.to_explore = stats.norm.sf(bx.pf_estimation.r) bx.to_explore = stats.norm.sf(bx.pf_estimation.r)
482 454 bx.to_refine = bx.get_max_potential() bx.to_refine = bx.get_max_potential()
483 455
484
485 #č teď je to kolbek, který volá Triangulation
486 def _on_safe_added(bx, simplex, indices, vertices_model, nodes, vol):
487 pass
488
489
490 #č teď je to kolbek, který volá Triangulation
491 def _on_mixed_added(bx, simplex, indices, vertices_model, nodes, _vol, fr, wfr, _mfpdf):
492 pass
493
494 # callback
495 #č sx.on_delete_simplex(indices=indices)
496 def _invalidate_simplex(bx, simplex):
497 pass
498
499 456
500 457 def increment(bx, input_sample): def increment(bx, input_sample):
501 458 #č tri - Deloneho triangulace #č tri - Deloneho triangulace
 
... ... class _PsiTri:
560 517 bx._nsim = bx.nsim bx._nsim = bx.nsim
561 518
562 519
520 def _regen_inside(bx):
521 failsi = bx.failsi
522 # incremental triangulation require one more point
523 if (bx.nsim > bx.nvar + 1) and np.any(failsi) and not np.all(failsi):
524 #bx._logger(msg="triangulation started")
525 bx.__regen_inside()
526 else:
527 #č jíž není nutný
528 bx._logger(msg="triangulation skipped")
529
530 def __regen_inside(bx):
531 try:
532 if bx.holydays:
533 Integrator = sx.FullCubatureIntegration
534 else:
535 Integrator = sx.FastCubatureIntegration
536
537 bx.Tri = Integrator(bx.samplebox, bx.scheme,\
538 tri_space=bx.tri_space, incremental=True,\
539 on_safe_added=bx._on_safe_added,\
540 on_mixed_added=bx._on_mixed_added,\
541 on_delete_simplex=bx._invalidate_simplex)
542
543 bx.Tri.integrate() # nic nevrácí, všecko je přes kolbeky
544 #č tri - Deloneho triangulace
545 bx.tri = bx.Tri.tri #č všichní tam očekávajou QHull
546 bx._logger(msg="triangulation has been created")
547
548 except BaseException as e:
549 #č chcu zachytit spadnuti QHull na začatku,
550 #č kdy ještě není dostatek teček.
551 #č Jinak je třeba nechat QHull spadnout
552 if bx.nsim > 2*bx.nvar + 3:
553 #č no to teda ne!
554 raise
555 else:
556 #č lze přípustit chybu triangulace
557 bx._logger(msg='triangulation failed')
558
563 559 def get_max_potential(bx): def get_max_potential(bx):
564 560 if len(bx.candidates_index): if len(bx.candidates_index):
565 561 __key, max_potential = bx.potential_index.peekitem() __key, max_potential = bx.potential_index.peekitem()
Hints:
Before first commit, do not forget to setup your git environment:
git config --global user.name "your_name_here"
git config --global user.email "your@email_here"

Clone this repository using HTTP(S):
git clone https://rocketgit.com/user/iam-git/WellMet

Clone this repository using ssh (do not forget to upload a key first):
git clone ssh://rocketgit@ssh.rocketgit.com/user/iam-git/WellMet

Clone this repository using git:
git clone git://git.rocketgit.com/user/iam-git/WellMet

You are allowed to anonymously push to this repository.
This means that your pushed commits will automatically be transformed into a merge request:
... clone the repository ...
... make some changes and some commits ...
git push origin main