update examples

Author: JianqiangDing

example/benchmark_boundary_analysis_vs_naive_partition_cmp.py

@@ -0,0 +1,54 @@
+import numpy as np
+from pybdr.geometry import Interval, Zonotope, Geometry
+from pybdr.geometry.operation import cvt2, partition, boundary
+from pybdr.algorithm import ASB2008CDC
+from pybdr.dynamic_system import NonLinSys
+from pybdr.model import *
+from pybdr.util.visualization import plot, plot_cmp
+from pybdr.util.functional import performance_counter, performance_counter_start
+
+if __name__ == '__main__':
+    time_start = performance_counter_start()
+    # init dynamic system
+    system = NonLinSys(Model(lotka_volterra_2d, [2, 1]))
+
+    # settings for the computation
+    options = ASB2008CDC.Options()
+    options.t_end = 2.2
+    options.step = 0.005
+    options.tensor_order = 3
+    options.taylor_terms = 4
+
+    options.u = Zonotope([0], np.diag([0]))
+    options.u_trans = options.u.c
+
+    # settings for the using geometry
+    Zonotope.REDUCE_METHOD = Zonotope.REDUCE_METHOD.GIRARD
+    Zonotope.ORDER = 50
+
+    z = Interval.identity(2) * 0.5 + 3
+
+    options.r0 = [cvt2(z, Geometry.TYPE.ZONOTOPE)]
+    _, tp_whole, _, _ = ASB2008CDC.reach(system, options)
+
+    # --------------------------------------------------------
+    # options.r0 = partition(z, 1, Geometry.TYPE.ZONOTOPE)
+    # 4
+    # ASB2008CDC cost: 43.344963666000005s
+    # --------------------------------------------------------
+    # options.r0 = partition(z, 0.5, Geometry.TYPE.ZONOTOPE)
+    # 9
+    # ASB2008CDC cost: 3868912500001s
+    # --------------------------------------------------------
+    options.r0 = partition(z, 0.2, Geometry.TYPE.ZONOTOPE)
+    # 36
+    # ASB2008CDC cost: 317.59988937500003s
+    # --------------------------------------------------------
+    # options.r0 = partition(z, 0.1, Geometry.TYPE.ZONOTOPE)
+
+    print(len(options.r0))
+
+    _, tp_part_00, _, _ = ASB2008CDC.reach(system, options)
+    performance_counter(time_start, "ASB2008CDC")
+
+    plot_cmp([tp_whole, tp_part_00], [0, 1], cs=["#FF5722", "#303F9F"])

example/benchmark_brusselator_cmp.py

@@ -28,7 +28,7 @@
     Zonotope.REDUCE_METHOD = Zonotope.REDUCE_METHOD.GIRARD
     Zonotope.ORDER = 50
 
-    z = Interval.identity(2) * 0.1 + 0.2;
+    z = Interval.identity(2) * 0.1 + 0.2
 
     no_boundary_analysis = False
 

example/benchmark_lotka_volterra_2d_cmp.py

@@ -8,8 +8,10 @@
 from pybdr.geometry.operation import boundary, cvt2
 from pybdr.model import *
 from pybdr.util.visualization import plot, plot_cmp
+from pybdr.util.functional import performance_counter_start, performance_counter
 
 if __name__ == '__main__':
+    time_start = performance_counter_start()
     # init dynamic system
     system = NonLinSys(Model(lotka_volterra_2d, [2, 1]))
 
@@ -34,8 +36,24 @@
     ti_whole, tp_whole, _, _ = ASB2008CDC.reach(system, options)
 
     # reachable sets computation with boundary analysis
+    # --------------------------------------------------------
+    # options.r0 = boundary(z, 2, Geometry.TYPE.ZONOTOPE)
+    # 4
+    # ASB2008CDC cost: 43.888805708s
+    # --------------------------------------------------------
     options.r0 = boundary(z, 1, Geometry.TYPE.ZONOTOPE)
+    # 8
+    # ASB2008CDC cost: 77.14907875s
+    # --------------------------------------------------------
     ti_bound, tp_bound, _, _ = ASB2008CDC.reach(system, options)
+    performance_counter(time_start, 'ASB2008CDC')
 
     # visualize the results
-    plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#303F9F"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#303F9F"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#D6D6D6"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#808000"])
+    plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#228B22"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#663399"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#008080"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#FFD700"])
+    # plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#2E8B57"])