Note
Go to the end to download the full example code.
Isolated star GW emmission#
Tests that the GW emission from an isolated star is indeed 0.
9 import os
10
11 import shamrock
12
13 if not shamrock.sys.is_initialized():
14 shamrock.change_loglevel(1)
15 shamrock.sys.init("0:0")
16
17 shamrock.enable_experimental_features()
18
19
20 def is_in_sphere(pt):
21 x, y, z = pt
22 return (x**2 + y**2 + z**2) < 1
Use shamrock documentation style for matplotlib
28 shamrock.matplotlib.set_shamrock_mpl_style()
29
30 si = shamrock.UnitSystem()
31 sicte = shamrock.Constants(si)
32 codeu = shamrock.UnitSystem(
33 unit_time=1,
34 unit_length=sicte.solar_radius(),
35 unit_mass=sicte.sol_mass(),
36 )
37 ucte = shamrock.Constants(codeu)
Setup parameters
42 dr = 0.1
43 pmass = 1
44
45 C_cour = 0.3
46 C_force = 0.25
47
48 bsize = 4
49
50 render_gif = True
51
52 dump_folder = "_to_trash"
53 sim_name = "isolated_star"
54
55 if shamrock.sys.world_rank() == 0:
56 os.makedirs(dump_folder, exist_ok=True)
Setup
60 ctx = shamrock.Context()
61 ctx.pdata_layout_new()
62
63 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
64
65 cfg = model.gen_default_config()
66 cfg.set_units(codeu)
67 # cfg.set_artif_viscosity_VaryingMM97(alpha_min = 0.1,alpha_max = 1,sigma_decay = 0.1, alpha_u = 1, beta_AV = 2)
68 # cfg.set_artif_viscosity_ConstantDisc(alpha_u=alpha_u, alpha_AV=alpha_AV, beta_AV=beta_AV)
69 cfg.set_eos_polytropic(1e-12, 5.0 / 3.0)
70 cfg.set_artif_viscosity_VaryingCD10(
71 alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2
72 )
73 cfg.set_self_gravity_fmm(order=1, opening_angle=0.5, reduction_level=3)
74 cfg.set_softening_plummer(epsilon=1e-9)
75 cfg.set_boundary_free()
76 cfg.compute_GW(True)
The important part to enable killing
80 cfg.add_kill_sphere(center=(0.0, 0.0, 0.0), radius=1.0)
Rest of the setup
84 cfg.print_status()
85 model.set_solver_config(cfg)
86
87 model.init_scheduler(int(1e7), 1)
88
89 bmin = (-bsize, -bsize, -bsize)
90 bmax = (bsize, bsize, bsize)
91 model.resize_simulation_box(bmin, bmax)
92
93 model.set_particle_mass(pmass)
94
95 setup = model.get_setup()
96 lat = setup.make_generator_lattice_hcp(dr, (-bsize, -bsize, -bsize), (bsize, bsize, bsize))
97
98 thesphere = setup.make_modifier_filter(parent=lat, filter=is_in_sphere)
99
100 setup.apply_setup(thesphere)
101
102 model.set_value_in_a_box("uint", "f64", 1, bmin, bmax)
103
104 model.set_cfl_cour(C_cour)
105 model.set_cfl_force(C_force)
106
107 model.change_htolerance(1.3)
108 model.timestep()
109 model.change_htolerance(1.1)
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "free"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"enable_particle_reordering": false,
"eos_config": {
"K": 1e-12,
"Tvec": "f64_3",
"eos_type": "polytropic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"neigh_cache_strategy": "two_stage",
"particle_killing": [
{
"center": [
0.0,
0.0,
0.0
],
"radius": 1.0,
"type": "sphere"
}
],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"self_grav_config": {
"opening_angle": 0.5,
"order": 1,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "fmm"
},
"show_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"tree_reduction_level": 3,
"type_id": "sycl::vec<f64,3>",
"unit_sys": {
"unit_current": 1.0,
"unit_length": 695700000.0,
"unit_lumint": 1.0,
"unit_mass": 1.98847e+30,
"unit_qte": 1.0,
"unit_temperature": 1.0,
"unit_time": 1.0
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 735 ( 7.4e+02 ) Ntotal = 735 ( 7.4e+02 rank min = 8.8e+03 max = 7.4e+02) rate = 7.350000e+02 N.s^-1
SPH setup: the generation step took : 0.084416929 s
SPH setup: final particle count = 735 beginning injection ...
SPH setup: injected 735 / 735 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.007351364000000001 s
SPH setup: the setup took : 0.09320268000000001 s
Warning: .change_htolerance(val) is deprecated, [SPH][rank=0]
-> calling this is replaced internally by .change_htolerances(coarse=val, fine=min(val, 1.1))
see: https://shamrock-code.github.io/Shamrock/sphinx/user_guide/sph/smoothing_length_tolerance.html
---------------- t = 0, dt = 0 ----------------
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.13 unconverged cnt = 735
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.169 unconverged cnt = 735
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.21970000000000003 unconverged cnt = 334
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.28561000000000003 unconverged cnt = 2
################## hx = [0.00014041713145738393, 5.747004266219756e-05, -4.0876097653306834e-05, -0.00012826952061287196]
################## hp = [0.00042473407116827566, -5.339795094255786e-06, -0.0010508039109739116, -0.0016661941605910224]
################## ddq = [-0.07660472093148431, 7.020856572869196e-05, 6.413476030643598e-05, -0.07702945500265258, 0.0001069924638150102, -0.07493852677089329]
################## ddq_xy = [-0.07660472093148431, 7.020856572869196e-05, 6.413476030643598e-05, 7.020856572869196e-05, -0.07702945500265258, 0.0001069924638150102, 6.413476030643598e-05, 0.0001069924638150102, -0.07493852677089329]
Warning: .change_htolerance(val) is deprecated, [SPH][rank=0]
-> calling this is replaced internally by .change_htolerances(coarse=val, fine=min(val, 1.1))
see: https://shamrock-code.github.io/Shamrock/sphinx/user_guide/sph/smoothing_length_tolerance.html
Show the current solvergraph
Draw utilities#
118 import matplotlib.pyplot as plt
119 import numpy as np
120
121
122 def plot_state(iplot):
123 pos = ctx.collect_data()["xyz"]
124
125 if shamrock.sys.world_rank() == 0:
126 X = pos[:, 0]
127 Y = pos[:, 1]
128 Z = pos[:, 2]
129
130 plt.cla()
131
132 ax.set_xlim3d(bmin[0], bmax[0])
133 ax.set_ylim3d(bmin[1], bmax[1])
134 ax.set_zlim3d(bmin[2], bmax[2])
135
136 ax.scatter(X, Y, Z, s=1)
137
138 ax.minorticks_off()
139
140 ax.set_title(f"t = {model.get_time():.2f} ")
141
142 plt.savefig(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.png"))
Run the simulation#
148 nstop = 28 # To be increased when empty simulations will be fixed
149 dt_stop = 0.1
150
151 t_stop = [i * dt_stop for i in range(nstop + 1)]
152
153 # Init MPL
154 fig = plt.figure(dpi=120)
155 ax = fig.add_subplot(111, projection="3d")
156
157 iplot = 0
158 istop = 0
159 for ttarg in t_stop:
160 model.evolve_until(ttarg)
161 model.do_vtk_dump(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.vtk"), True)
162
163 # if do_plots:
164 plot_state(iplot)
165
166 iplot += 1
167 istop += 1
168
169 plt.close(fig)
---------------- t = 0, dt = 0.07535026072600358 ----------------
################## hx = [0.0001884035412103867, 4.589396996328958e-05, -0.00010891285347292917, -0.0002345365657757074]
################## hp = [0.00043231092361886714, -0.00017658053053626224, -0.0011875291883355575, -0.0015895863919797165]
################## ddq = [-0.07960475984393818, 9.420177060519335e-05, 0.0001172682828878537, -0.08003707076755705, -0.00011941581021014926, -0.07801517345195846]
################## ddq_xy = [-0.07960475984393818, 9.420177060519335e-05, 0.0001172682828878537, 9.420177060519335e-05, -0.08003707076755705, -0.00011941581021014926, 0.0001172682828878537, -0.00011941581021014926, -0.07801517345195846]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.1203472531970372
################## hx = [0.00018840340389168596, 4.589393704000073e-05, -0.00010891277317903747, -0.0002345363937793186]
################## hp = [0.0004323110858687751, -0.00017658002931334986, -0.0011875287696984827, -0.0015895863949014905]
################## ddq = [-0.07960475394310856, 9.420170194584298e-05, 0.0001172681968896593, -0.08003706502897734, -0.00011941537111687705, -0.07801516754820707]
################## ddq_xy = [-0.07960475394310856, 9.420170194584298e-05, 0.0001172681968896593, 9.420170194584298e-05, -0.08003706502897734, -0.00011941537111687705, 0.0001172681968896593, -0.00011941537111687705, -0.07801516754820707]
---------------- t = 0.07535026072600358, dt = 0.024649739273996427 ----------------
################## hx = [-9.971813891475513e-05, -4.4934814007552505e-05, 2.1888758025016527e-05, 8.284725502146213e-05]
################## hp = [0.000648210974184632, -6.942898831230159e-05, -0.0011215321164327555, -0.0014559952820562688]
################## ddq = [-0.07947231071746554, -4.9859069457377565e-05, -4.142362751073107e-05, -0.08012052169165017, -0.00022122722682209377, -0.07801631543540927]
################## ddq_xy = [-0.07947231071746554, -4.9859069457377565e-05, -4.142362751073107e-05, -4.9859069457377565e-05, -0.08012052169165017, -0.00022122722682209377, -4.142362751073107e-05, -0.00022122722682209377, -0.07801631543540927]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.019875162473671937
################## hx = [-9.971793854946274e-05, -4.4934778929956025e-05, 2.1888618415903438e-05, 8.284697813378859e-05]
################## hp = [0.0006482108213409082, -6.942906851617085e-05, -0.0011215321924822534, -0.001455995426591264]
################## ddq = [-0.0794723109424222, -4.985896927473137e-05, -4.1423489066894296e-05, -0.08012052176376311, -0.000221227145343334, -0.07801631551583094]
################## ddq_xy = [-0.0794723109424222, -4.985896927473137e-05, -4.1423489066894296e-05, -4.985896927473137e-05, -0.08012052176376311, -0.000221227145343334, -4.1423489066894296e-05, -0.000221227145343334, -0.07801631551583094]
---------------- t = 0.1, dt = 0.1 ----------------
################## hx = [2.3816193457783393e-05, 1.750415168534091e-05, 6.5018866046194586e-06, -6.242553741088516e-06]
################## hp = [0.0006612421565927451, 7.362426841882571e-05, -0.0011099519594206736, -0.0017059102990862468]
################## ddq = [-0.07934971983783422, 1.1908096728891697e-05, 3.1212768705442587e-06, -0.08001096199442696, 4.815361913863272e-06, -0.07764380953874797]
################## ddq_xy = [-0.07934971983783422, 1.1908096728891697e-05, 3.1212768705442587e-06, 1.1908096728891697e-05, -0.08001096199442696, 4.815361913863272e-06, 3.1212768705442587e-06, 4.815361913863272e-06, -0.07764380953874797]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.039235827187017354
################## hx = [2.3815518482377747e-05, 1.7503955837442577e-05, 6.502222361514634e-06, -6.241776345188746e-06]
################## hp = [0.0006612423464209954, 7.362364700049787e-05, -0.0011099517234065014, -0.0017059083943929965]
################## ddq = [-0.07934972074617508, 1.1907759241188874e-05, 3.120888172594374e-06, -0.08001096309259607, 4.813930128117996e-06, -0.07764381235178208]
################## ddq_xy = [-0.07934972074617508, 1.1907759241188874e-05, 3.120888172594374e-06, 1.1907759241188874e-05, -0.08001096309259607, 4.813930128117996e-06, 3.120888172594374e-06, 4.813930128117996e-06, -0.07764381235178208]
---------------- t = 0.2, dt = 0.10000000000000003 ----------------
################## hx = [-0.00013507822666621532, -0.0001679537037307761, -0.0001558261215148594, -0.00010194505587936212]
################## hp = [0.00038850325439905575, 3.501974132229623e-05, -0.0010012706843743817, -0.0016840775969943073]
################## ddq = [-0.07946546558318789, -6.753911333310766e-05, 5.097252793968105e-05, -0.07985396883758694, 0.00019013495337668354, -0.07778138798619358]
################## ddq_xy = [-0.07946546558318789, -6.753911333310766e-05, 5.097252793968105e-05, -6.753911333310766e-05, -0.07985396883758694, 0.00019013495337668354, 5.097252793968105e-05, 0.00019013495337668354, -0.07778138798619358]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.026797945974446614
################## hx = [-0.00013507642843802034, -0.00016795178629953702, -0.00015582459865472827, -0.00010194433563948125]
################## hp = [0.0003885060418110653, 3.502029384959063e-05, -0.0010012714449527636, -0.0016840774357936433]
################## ddq = [-0.07946546463649606, -6.753821421901017e-05, 5.097216781974062e-05, -0.07985397067830713, 0.0001901331308760134, -0.07778138720070242]
################## ddq_xy = [-0.07946546463649606, -6.753821421901017e-05, 5.097216781974062e-05, -6.753821421901017e-05, -0.07985397067830713, 0.0001901331308760134, 5.097216781974062e-05, 0.0001901331308760134, -0.07778138720070242]
---------------- t = 0.30000000000000004, dt = 0.09999999999999998 ----------------
################## hx = [4.521791631523122e-05, 0.00012671888793447204, 0.00017426563586590116, 0.00017511804739856596]
################## hp = [0.0009489372984436556, 0.0004626037071456524, -0.0005818544632754153, -0.001139979042398473]
################## ddq = [-0.0790423740324346, 2.260895815761561e-05, -8.755902369928298e-05, -0.07999131133087825, 4.144854614618229e-05, -0.07790239499003612]
################## ddq_xy = [-0.0790423740324346, 2.260895815761561e-05, -8.755902369928298e-05, 2.260895815761561e-05, -0.07999131133087825, 4.144854614618229e-05, -8.755902369928298e-05, 4.144854614618229e-05, -0.07790239499003612]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.021878829995611877
################## hx = [4.5215054458214935e-05, 0.00012671461845006017, 0.00017426110275899396, 0.00017511446531149843]
################## hp = [0.0009489293656237546, 0.0004625975537854724, -0.0005818601153397607, -0.0011399859726266975]
################## ddq = [-0.07904238154184981, 2.2607527229107468e-05, -8.755723265574921e-05, -0.07999131090747356, 4.145031146598902e-05, -0.07790239556922311]
################## ddq_xy = [-0.07904238154184981, 2.2607527229107468e-05, -8.755723265574921e-05, 2.2607527229107468e-05, -0.07999131090747356, 4.145031146598902e-05, -8.755723265574921e-05, 4.145031146598902e-05, -0.07790239556922311]
---------------- t = 0.4, dt = 0.09999999999999998 ----------------
################## hx = [-3.68665005525858e-05, 8.035032188593551e-05, 0.00017603734046353964, 0.00022455529582621574]
################## hp = [0.001041422375260198, 0.000505182429467653, -0.0006420653638331233, -0.0012530732113413406]
################## ddq = [-0.07896753443010207, -1.84332502762929e-05, -0.00011227764791310787, -0.08000895680536227, 4.316726818228683e-05, -0.07771446121876073]
################## ddq_xy = [-0.07896753443010207, -1.84332502762929e-05, -0.00011227764791310787, -1.84332502762929e-05, -0.08000895680536227, 4.316726818228683e-05, -0.00011227764791310787, 4.316726818228683e-05, -0.07771446121876073]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.015440448554529806
################## hx = [-3.6864790231982394e-05, 8.035121142387247e-05, 0.00017603717086783833, 0.00022455411253990734]
################## hp = [0.0010414220820839648, 0.0005051832572331305, -0.0006420629840441083, -0.001253070400470499]
################## ddq = [-0.07896753477416099, -1.8432395115991197e-05, -0.00011227705626995367, -0.08000895685624496, 4.316766647307664e-05, -0.07771446437369049]
################## ddq_xy = [-0.07896753477416099, -1.8432395115991197e-05, -0.00011227705626995367, -1.8432395115991197e-05, -0.08000895685624496, 4.316766647307664e-05, -0.00011227705626995367, 4.316766647307664e-05, -0.07771446437369049]
---------------- t = 0.5, dt = 0.10000000000000009 ----------------
################## hx = [-0.00012949902284358536, -0.00010878913557843436, -5.8929287289761965e-05, 6.720615938743784e-06]
################## hp = [0.0004321996951898305, -1.5263191490399295e-05, -0.0010620450709538354, -0.001661364063737028]
################## ddq = [-0.07960984155126449, -6.474951142179268e-05, -3.360307969371896e-06, -0.08004204124645432, 8.767416373662993e-05, -0.07794847748752746]
################## ddq_xy = [-0.07960984155126449, -6.474951142179268e-05, -3.360307969371896e-06, -6.474951142179268e-05, -0.08004204124645432, 8.767416373662993e-05, -3.360307969371896e-06, 8.767416373662993e-05, -0.07794847748752746]
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.013639617510270843
################## hx = [-0.00012949842578350473, -0.00010878590311546828, -5.8924285559751755e-05, 6.726046726281124e-06]
################## hp = [0.0004322110749959418, -1.5253789217287888e-05, -0.0010620374890561613, -0.001661356324681798]
################## ddq = [-0.07960982894831724, -6.474921289175236e-05, -3.363023363140566e-06, -0.08004204002331318, 8.767293127246096e-05, -0.07794847262363544]
################## ddq_xy = [-0.07960982894831724, -6.474921289175236e-05, -3.363023363140566e-06, -6.474921289175236e-05, -0.08004204002331318, 8.767293127246096e-05, -3.363023363140566e-06, 8.767293127246096e-05, -0.07794847262363544]
---------------- t = 0.6000000000000001, dt = 0.09999999999999998 ----------------
################## hx = [-0.00015037614456903938, -0.00018539721085640085, -0.0001707412442158071, -0.00011033529907290278]
################## hp = [0.0007362012959452602, 0.0002944241500056639, -0.0007100797435661883, -0.0012728064911984305]
################## ddq = [-0.07913080659351311, -7.518807228451969e-05, 5.516764953645138e-05, -0.07986700788945837, 6.983028509562311e-05, -0.07785800010231468]
################## ddq_xy = [-0.07913080659351311, -7.518807228451969e-05, 5.516764953645138e-05, -7.518807228451969e-05, -0.07986700788945837, 6.983028509562311e-05, 5.516764953645138e-05, 6.983028509562311e-05, -0.07785800010231468]
---------------- t = 0.7000000000000001, dt = 0.09999999999999998 ----------------
################## hx = [-0.00022868011900989887, -0.00021590726650605452, -0.00014528223630190165, -3.572894820606669e-05]
################## hp = [0.0006139377221974912, 0.0001741110088503026, -0.0008405989080294074, -0.001415482111561922]
################## ddq = [-0.0794412660608319, -0.00011434005950494944, 1.786447410303334e-05, -0.0800552037830294, 7.797490096429801e-05, -0.07802578394926998]
################## ddq_xy = [-0.0794412660608319, -0.00011434005950494944, 1.786447410303334e-05, -0.00011434005950494944, -0.0800552037830294, 7.797490096429801e-05, 1.786447410303334e-05, 7.797490096429801e-05, -0.07802578394926998]
---------------- t = 0.8, dt = 0.09999999999999998 ----------------
################## hx = [-0.00033097856347786344, -0.00023702906647741766, -7.956782253164489e-05, 9.921355520498501e-05]
################## hp = [0.000514388064492638, 0.00019244179201366527, -0.0008030147875855215, -0.0014765250947057218]
################## ddq = [-0.07963620951687218, -0.00016548928173893172, -4.9606777602492514e-05, -0.08015059758136482, 0.00020297559003749964, -0.07815968442216646]
################## ddq_xy = [-0.07963620951687218, -0.00016548928173893172, -4.9606777602492514e-05, -0.00016548928173893172, -0.08015059758136482, 0.00020297559003749964, -4.9606777602492514e-05, 0.00020297559003749964, -0.07815968442216646]
---------------- t = 0.9, dt = 0.09999999999999998 ----------------
################## hx = [-0.00026682764867283113, -0.00023684927387885066, -0.00014340732742113436, -1.1539503392219499e-05]
################## hp = [0.0005927575822671227, 0.00014141442888178092, -0.0009449004468863833, -0.001579872169269192]
################## ddq = [-0.07964373155181907, -0.00013341382433641557, 5.769751696109741e-06, -0.0802364891340862, 0.00010601800374158922, -0.07806385938254988]
################## ddq_xy = [-0.07964373155181907, -0.00013341382433641557, 5.769751696109741e-06, -0.00013341382433641557, -0.0802364891340862, 0.00010601800374158922, 5.769751696109741e-06, 0.00010601800374158922, -0.07806385938254988]
---------------- t = 1, dt = 0.10000000000000009 ----------------
################## hx = [-0.00016148748375513863, -0.0001626749170950917, -0.00012027373777061511, -4.5645307539829575e-05]
################## hp = [0.0007817373563132057, 0.0003186924850272652, -0.000777123461037385, -0.0014098945358160808]
################## ddq = [-0.0794966765928312, -8.074374187756931e-05, 2.2822653769914784e-05, -0.08027841394914441, 9.799146927507143e-05, -0.07808678205701512]
################## ddq_xy = [-0.0794966765928312, -8.074374187756931e-05, 2.2822653769914784e-05, -8.074374187756931e-05, -0.08027841394914441, 9.799146927507143e-05, 2.2822653769914784e-05, 9.799146927507143e-05, -0.07808678205701512]
---------------- t = 1.1, dt = 0.10000000000000009 ----------------
################## hx = [6.851167573411802e-05, 4.5822482695775544e-05, 1.0855192423910929e-05, -2.7020737891625073e-05]
################## hp = [0.0007398907017788781, 0.0001907506163443345, -0.0009646063487971423, -0.0015708232285040619]
################## ddq = [-0.07951963622432504, 3.425583786705901e-05, 1.3510368945812538e-05, -0.08025952692610391, 3.295330253763715e-05, -0.07794881299582097]
################## ddq_xy = [-0.07951963622432504, 3.425583786705901e-05, 1.3510368945812538e-05, 3.425583786705901e-05, -0.08025952692610391, 3.295330253763715e-05, 1.3510368945812538e-05, 3.295330253763715e-05, -0.07794881299582097]
---------------- t = 1.2000000000000002, dt = 0.09999999999999987 ----------------
################## hx = [9.220926467726594e-05, 5.338077133756762e-05, 2.4894342661764546e-07, -5.294958867445558e-05]
################## hp = [0.0005189500156768528, 1.5061272105693726e-05, -0.0009713576532705639, -0.0014538878350756623]
################## ddq = [-0.07934983931320001, 4.610463233863297e-05, 2.6474794337227794e-05, -0.07986878932887687, -1.2331371385142367e-05, -0.07789595147812435]
################## ddq_xy = [-0.07934983931320001, 4.610463233863297e-05, 2.6474794337227794e-05, 4.610463233863297e-05, -0.07986878932887687, -1.2331371385142367e-05, 2.6474794337227794e-05, -1.2331371385142367e-05, -0.07789595147812435]
---------------- t = 1.3, dt = 0.10000000000000009 ----------------
################## hx = [-0.00012626207996012026, -0.00013105077428567145, -0.00010072451947390365, -4.340921101109045e-05]
################## hp = [0.00020567839954772904, -0.00020015814808727855, -0.001159080477799072, -0.0017121662598758441]
################## ddq = [-0.07967730711166426, -6.313103998006013e-05, 2.170460550554522e-05, -0.07988298551121199, 8.501438514291502e-05, -0.07796514085178842]
################## ddq_xy = [-0.07967730711166426, -6.313103998006013e-05, 2.170460550554522e-05, -6.313103998006013e-05, -0.07988298551121199, 8.501438514291502e-05, 2.170460550554522e-05, 8.501438514291502e-05, -0.07796514085178842]
---------------- t = 1.4000000000000001, dt = 0.09999999999999987 ----------------
################## hx = [4.297643335451646e-07, 3.4092621658021084e-06, 5.475250953946398e-06, 6.074150670623015e-06]
################## hp = [0.00035791320154221484, -3.5229933220021046e-05, -0.0009948478356380359, -0.0015613226032938077]
################## ddq = [-0.07946051880634139, 2.148821667725823e-07, -3.0370753353115075e-06, -0.0798184320078836, 0.00010007306491015841, -0.07789919620304758]
################## ddq_xy = [-0.07946051880634139, 2.148821667725823e-07, -3.0370753353115075e-06, 2.148821667725823e-07, -0.0798184320078836, 0.00010007306491015841, -3.0370753353115075e-06, 0.00010007306491015841, -0.07789919620304758]
---------------- t = 1.5, dt = 0.10000000000000009 ----------------
################## hx = [-9.10514954989783e-06, -7.320351102170473e-05, -0.00011768705083212306, -0.00013063644041247353]
################## hp = [0.00046506191228014304, 7.753926739901159e-05, -0.0008654638027129545, -0.0014209442279437823]
################## ddq = [-0.07915317982422004, -4.552574774948915e-06, 6.531822020623676e-05, -0.07961824173650019, 9.697046969738319e-05, -0.07773223559627626]
################## ddq_xy = [-0.07915317982422004, -4.552574774948915e-06, 6.531822020623676e-05, -4.552574774948915e-06, -0.07961824173650019, 9.697046969738319e-05, 6.531822020623676e-05, 9.697046969738319e-05, -0.07773223559627626]
---------------- t = 1.6, dt = 0.10000000000000009 ----------------
################## hx = [-4.941088571330397e-06, 8.41452667311115e-05, 0.00015068496576605068, 0.0001768487499124652]
################## hp = [0.000800321571210294, 0.0003115332575472571, -0.0007007896802767484, -0.00122432430443771]
################## ddq = [-0.07899157795482269, -2.4705442856651986e-06, -8.84243749562326e-05, -0.07979189952603298, 2.0060791719317108e-05, -0.07776725365038498]
################## ddq_xy = [-0.07899157795482269, -2.4705442856651986e-06, -8.84243749562326e-05, -2.4705442856651986e-06, -0.07979189952603298, 2.0060791719317108e-05, -8.84243749562326e-05, 2.0060791719317108e-05, -0.07776725365038498]
---------------- t = 1.7000000000000002, dt = 0.09999999999999987 ----------------
################## hx = [-0.000118670964146953, -4.2990584858866334e-05, 4.420908692429523e-05, 0.00011956296956797459]
################## hp = [0.0006516451750623181, 0.00021440194906280447, -0.0007991695285332012, -0.0013754977801296864]
################## ddq = [-0.07901332637745356, -5.93354820734765e-05, -5.97814847839873e-05, -0.07966497155251588, 8.030077696865147e-05, -0.07763782859732388]
################## ddq_xy = [-0.07901332637745356, -5.93354820734765e-05, -5.97814847839873e-05, -5.93354820734765e-05, -0.07966497155251588, 8.030077696865147e-05, -5.97814847839873e-05, 8.030077696865147e-05, -0.07763782859732388]
---------------- t = 1.8, dt = 0.10000000000000009 ----------------
################## hx = [-0.0001010765463181838, -7.811227693590169e-05, -3.4217886029688516e-05, 1.8845159804879907e-05]
################## hp = [0.0006980114731442583, 0.00026479643241028545, -0.0007384671103844134, -0.0013085156124451325]
################## ddq = [-0.07880734803633724, -5.05382731590919e-05, -9.422579902439957e-06, -0.0795053595094815, 7.900083573113918e-05, -0.0774988324238921]
################## ddq_xy = [-0.07880734803633724, -5.05382731590919e-05, -9.422579902439957e-06, -5.05382731590919e-05, -0.0795053595094815, 7.900083573113918e-05, -9.422579902439957e-06, 7.900083573113918e-05, -0.0774988324238921]
---------------- t = 1.9000000000000001, dt = 0.09999999999999987 ----------------
################## hx = [-0.00016942714509535456, -0.00015880643869987688, -0.00010563367530190458, -2.415645391325552e-05]
################## hp = [0.0007022450181569245, 0.00024450709254585036, -0.0007530931897136375, -0.0012929555463620512]
################## ddq = [-0.07888403427860445, -8.471357254767728e-05, 1.2078226956627755e-05, -0.07958627929676138, 4.7414562366448836e-05, -0.0775910787322424]
################## ddq_xy = [-0.07888403427860445, -8.471357254767728e-05, 1.2078226956627755e-05, -8.471357254767728e-05, -0.07958627929676138, 4.7414562366448836e-05, 1.2078226956627755e-05, 4.7414562366448836e-05, -0.0775910787322424]
---------------- t = 2, dt = 0.10000000000000009 ----------------
################## hx = [-0.00010424104913757816, -5.7467749780564055e-05, 4.703986720986031e-06, 6.561529377944119e-05]
################## hp = [0.00043814164733857464, 1.0482476362432775e-05, -0.0009723055703920269, -0.001527434446170331]
################## ddq = [-0.07887435260023509, -5.212052456878908e-05, -3.2807646889720596e-05, -0.07931249424757367, 7.359466838105044e-05, -0.07734691815406476]
################## ddq_xy = [-0.07887435260023509, -5.212052456878908e-05, -3.2807646889720596e-05, -5.212052456878908e-05, -0.07931249424757367, 7.359466838105044e-05, -3.2807646889720596e-05, 7.359466838105044e-05, -0.07734691815406476]
---------------- t = 2.1, dt = 0.10000000000000009 ----------------
################## hx = [1.409347662295504e-06, 5.813188307882898e-05, 9.927802736988978e-05, 0.00011382270440103373]
################## hp = [0.0005971962996619978, 0.00015455769128264538, -0.0008629844207326716, -0.0014378879243686221]
################## ddq = [-0.07873840199256579, 7.04673831147752e-07, -5.6911352200516863e-05, -0.07933559829222779, 7.636317288073057e-05, -0.07730051406819717]
################## ddq_xy = [-0.07873840199256579, 7.04673831147752e-07, -5.6911352200516863e-05, 7.04673831147752e-07, -0.07933559829222779, 7.636317288073057e-05, -5.6911352200516863e-05, 7.636317288073057e-05, -0.07730051406819717]
---------------- t = 2.2, dt = 0.10000000000000009 ----------------
################## hx = [3.954696822569629e-05, -4.6680356135969545e-06, -4.7632243079987177e-05, -7.783342947941191e-05]
################## hp = [0.0008500684951114612, 0.00022710287205425376, -0.0008878299219095011, -0.0013797970928160347]
################## ddq = [-0.07809847505602396, 1.9773484112848146e-05, 3.8916714739705956e-05, -0.07894854355113542, -7.563199161261591e-05, -0.07671867796320793]
################## ddq_xy = [-0.07809847505602396, 1.9773484112848146e-05, 3.8916714739705956e-05, 1.9773484112848146e-05, -0.07894854355113542, -7.563199161261591e-05, 3.8916714739705956e-05, -7.563199161261591e-05, -0.07671867796320793]
---------------- t = 2.3000000000000003, dt = 0.10000000000000009 ----------------
################## hx = [4.4195939500588836e-05, -3.0463463396279475e-06, -4.9472366138275924e-05, -8.264230538251605e-05]
################## hp = [0.0007499921659523956, 0.00010508379306172737, -0.001046133518537595, -0.0015524424572462558]
################## ddq = [-0.07830166815020739, 2.2097969750294418e-05, 4.132115269125803e-05, -0.07905166031615979, -8.002042063451773e-05, -0.07674922569296114]
################## ddq_xy = [-0.07830166815020739, 2.2097969750294418e-05, 4.132115269125803e-05, 2.2097969750294418e-05, -0.07905166031615979, -8.002042063451773e-05, 4.132115269125803e-05, -8.002042063451773e-05, -0.07674922569296114]
---------------- t = 2.4000000000000004, dt = 0.09999999999999964 ----------------
################## hx = [-6.567881154431634e-05, -2.5254845714418796e-05, 2.193613562962987e-05, 6.324934714665965e-05]
################## hp = [0.0009122368274927239, 0.0002526835203691956, -0.0008451592641580374, -0.001283448741561749]
################## ddq = [-0.07849860905813685, -3.283940577215817e-05, -3.1624673573329824e-05, -0.07941084588562958, -0.00012774673165066972, -0.0772151603165751]
################## ddq_xy = [-0.07849860905813685, -3.283940577215817e-05, -3.1624673573329824e-05, -3.283940577215817e-05, -0.07941084588562958, -0.00012774673165066972, -3.1624673573329824e-05, -0.00012774673165066972, -0.0772151603165751]
---------------- t = 2.5, dt = 0.10000000000000009 ----------------
################## hx = [-6.61004034838681e-05, -5.131283992115934e-05, -2.2776042340228464e-05, 1.186357739254369e-05]
################## hp = [0.0009249693067361481, 0.0003133689730725458, -0.0007927841326539047, -0.001287336904716746]
################## ddq = [-0.0783376039321429, -3.305020174193405e-05, -5.931788696271847e-06, -0.07926257323887904, -6.757744119819467e-05, -0.07705026702742615]
################## ddq_xy = [-0.0783376039321429, -3.305020174193405e-05, -5.931788696271847e-06, -3.305020174193405e-05, -0.07926257323887904, -6.757744119819467e-05, -5.931788696271847e-06, -6.757744119819467e-05, -0.07705026702742615]
---------------- t = 2.6, dt = 0.10000000000000009 ----------------
################## hx = [-6.754929037397267e-05, -0.0001234848537498609, -0.00014633275028599863, -0.00012997090455677787]
################## hp = [0.0005986791259985774, 0.00014650829462050415, -0.0008728995764152497, -0.0014401366160729234]
################## ddq = [-0.07861223661650016, -3.3774645186986335e-05, 6.498545227838894e-05, -0.07921091574249874, 6.643350632485083e-05, -0.07717210000042724]
################## ddq_xy = [-0.07861223661650016, -3.3774645186986335e-05, 6.498545227838894e-05, -3.3774645186986335e-05, -0.07921091574249874, 6.643350632485083e-05, 6.498545227838894e-05, 6.643350632485083e-05, -0.07717210000042724]
---------------- t = 2.7, dt = 0.10000000000000009 ----------------
################## hx = [7.763694035831071e-06, -8.48879623347015e-05, -0.00015479395775052726, -0.00018322303719388192]
################## hp = [0.0005504352823825637, 3.745246190294973e-05, -0.0010337946748781687, -0.0015920589911796734]
################## ddq = [-0.07866373791737986, 3.8818470179155355e-06, 9.161151859694096e-05, -0.07921417319976243, 2.6143283802073486e-05, -0.07707167892620019]
################## ddq_xy = [-0.07866373791737986, 3.8818470179155355e-06, 9.161151859694096e-05, 3.8818470179155355e-06, -0.07921417319976243, 2.6143283802073486e-05, 9.161151859694096e-05, 2.6143283802073486e-05, -0.07707167892620019]
Total running time of the script: (0 minutes 13.139 seconds)
Estimated memory usage: 160 MB