.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "_as_gen/tests_ci/run_GW_isolated_star.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr__as_gen_tests_ci_run_GW_isolated_star.py: Isolated star GW emmission ================================= Tests that the GW emission from an isolated star is indeed 0. .. GENERATED FROM PYTHON SOURCE LINES 7-25 .. code-block:: Python :lineno-start: 9 import os import shamrock if not shamrock.sys.is_initialized(): shamrock.change_loglevel(1) shamrock.sys.init("0:0") shamrock.enable_experimental_features() def is_in_sphere(pt): x, y, z = pt return (x**2 + y**2 + z**2) < 1 .. GENERATED FROM PYTHON SOURCE LINES 27-28 Use shamrock documentation style for matplotlib .. GENERATED FROM PYTHON SOURCE LINES 28-39 .. code-block:: Python :lineno-start: 28 shamrock.matplotlib.set_shamrock_mpl_style() si = shamrock.UnitSystem() sicte = shamrock.Constants(si) codeu = shamrock.UnitSystem( unit_time=1, unit_length=sicte.solar_radius(), unit_mass=sicte.sol_mass(), ) ucte = shamrock.Constants(codeu) .. GENERATED FROM PYTHON SOURCE LINES 40-41 Setup parameters .. GENERATED FROM PYTHON SOURCE LINES 41-58 .. code-block:: Python :lineno-start: 42 dr = 0.1 pmass = 1 C_cour = 0.3 C_force = 0.25 bsize = 4 render_gif = True dump_folder = "_to_trash" sim_name = "isolated_star" if shamrock.sys.world_rank() == 0: os.makedirs(dump_folder, exist_ok=True) .. GENERATED FROM PYTHON SOURCE LINES 59-60 Setup .. GENERATED FROM PYTHON SOURCE LINES 60-78 .. code-block:: Python :lineno-start: 60 ctx = shamrock.Context() ctx.pdata_layout_new() model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4") cfg = model.gen_default_config() cfg.set_units(codeu) # cfg.set_artif_viscosity_VaryingMM97(alpha_min = 0.1,alpha_max = 1,sigma_decay = 0.1, alpha_u = 1, beta_AV = 2) # cfg.set_artif_viscosity_ConstantDisc(alpha_u=alpha_u, alpha_AV=alpha_AV, beta_AV=beta_AV) cfg.set_eos_polytropic(1e-12, 5.0 / 3.0) cfg.set_artif_viscosity_VaryingCD10( alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2 ) cfg.set_self_gravity_fmm(order=1, opening_angle=0.5, reduction_level=3) cfg.set_softening_plummer(epsilon=1e-9) cfg.set_boundary_free() cfg.compute_GW(True) .. GENERATED FROM PYTHON SOURCE LINES 79-80 The important part to enable killing .. GENERATED FROM PYTHON SOURCE LINES 80-82 .. code-block:: Python :lineno-start: 80 cfg.add_kill_sphere(center=(0.0, 0.0, 0.0), radius=1.0) .. GENERATED FROM PYTHON SOURCE LINES 83-84 Rest of the setup .. GENERATED FROM PYTHON SOURCE LINES 84-111 .. code-block:: Python :lineno-start: 84 cfg.print_status() model.set_solver_config(cfg) model.init_scheduler(int(1e7), 1) bmin = (-bsize, -bsize, -bsize) bmax = (bsize, bsize, bsize) model.resize_simulation_box(bmin, bmax) model.set_particle_mass(pmass) setup = model.get_setup() lat = setup.make_generator_lattice_hcp(dr, (-bsize, -bsize, -bsize), (bsize, bsize, bsize)) thesphere = setup.make_modifier_filter(parent=lat, filter=is_in_sphere) setup.apply_setup(thesphere) model.set_value_in_a_box("uint", "f64", 1, bmin, bmax) model.set_cfl_cour(C_cour) model.set_cfl_force(C_force) model.change_htolerance(1.3) model.timestep() model.change_htolerance(1.1) .. rst-class:: sphx-glr-script-out .. code-block:: none ----- 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", "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", "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 .. GENERATED FROM PYTHON SOURCE LINES 112-113 Show the current solvergraph .. GENERATED FROM PYTHON SOURCE LINES 113-115 .. code-block:: Python :lineno-start: 113 shamrock.utils.plot.DotGraph(model.get_solver_dot_graph()) .. raw:: html
G cluster_1009074 time step cluster_1009052 sink accretion cluster_1009051 if_has_sinks cluster_1009073 point mass accretion cluster_1009072 if_has_ext_force_accretion cluster_1009071 ext force accretion cluster_1009063 sink predictor cluster_1009062 sink predictor body cluster_1009058 sink ext force cluster_1009057 sink ext force body cluster_1009027 attach fields cluster_1009035 leapfrog predictor cluster_1009030 half step cluster_1009033 half step cluster_1009039 part killing step n_1009044 SetEdge e_3583094 has_sinks n_1009044->e_3583094 n_1009018 SetEdge n_1009044->n_1009018 n_1009063 sink predictor e_3583094->n_1009063 n_1009058 sink ext force e_3583094->n_1009058 e_3583094->n_1009058 e_3583076 m n_1009018->e_3583076 n_1009018->e_3583076 n_1009019 SetEdge n_1009018->n_1009019 n_1009018->n_1009019 n_1009049 SinkParticlesAccreteQuantities e_3583076->n_1009049 e_3583073 patchdatas n_1009019->e_3583073 n_1009019->e_3583073 n_1009019->e_3583073 n_1009019->e_3583073 n_1009020 GetObjCntFromLayer n_1009019->n_1009020 n_1009019->n_1009020 n_1009019->n_1009020 n_1009019->n_1009020 e_3583073->n_1009020 e_3583073->n_1009020 e_3583073->n_1009020 e_3583073->n_1009020 n_1009021 GetFieldRefFromLayer e_3583073->n_1009021 e_3583073->n_1009021 e_3583073->n_1009021 e_3583073->n_1009021 n_1009022 GetFieldRefFromLayer e_3583073->n_1009022 e_3583073->n_1009022 e_3583073->n_1009022 n_1009023 GetFieldRefFromLayer e_3583073->n_1009023 e_3583073->n_1009023 e_3583073->n_1009023 n_1009024 GetFieldRefFromLayer e_3583073->n_1009024 e_3583073->n_1009024 n_1009025 GetFieldRefFromLayer e_3583073->n_1009025 e_3583073->n_1009025 n_1009026 GetFieldRefFromLayer e_3583073->n_1009026 e_3583073->n_1009026 e_3583074 Npart_patch n_1009020->e_3583074 n_1009020->e_3583074 n_1009020->e_3583074 n_1009020->e_3583074 n_1009020->n_1009021 n_1009020->n_1009021 n_1009020->n_1009021 n_1009020->n_1009021 n_1009048 SinkParticlesFlagAccreteHard e_3583074->n_1009048 e_3583074->n_1009049 n_1009050 SinkParticlesEvictAccretedParticles e_3583074->n_1009050 n_1009069 SinkParticlesFlagAccreteHard e_3583074->n_1009069 n_1009070 SinkParticlesEvictAccretedParticles e_3583074->n_1009070 n_1009028 ForwardEuler e_3583074->n_1009028 n_1009029 ForwardEuler e_3583074->n_1009029 n_1009034 ForwardEuler e_3583074->n_1009034 n_1009031 ForwardEuler e_3583074->n_1009031 n_1009032 ForwardEuler e_3583074->n_1009032 e_3583077 xyz n_1009021->e_3583077 n_1009021->e_3583077 n_1009021->e_3583077 n_1009021->e_3583077 n_1009021->n_1009022 n_1009021->n_1009022 n_1009021->n_1009022 n_1009021->n_1009069 e_3583077->n_1009048 e_3583077->n_1009049 e_3583077->n_1009069 n_1009037 GetParticlesOutsideSphere e_3583077->n_1009037 e_3583078 vxyz n_1009022->e_3583078 n_1009022->e_3583078 n_1009022->e_3583078 n_1009022->n_1009023 n_1009022->n_1009023 n_1009022->n_1009023 e_3583078->n_1009049 e_3583078->n_1009034 e_3583079 axyz n_1009023->e_3583079 n_1009023->e_3583079 n_1009023->e_3583079 n_1009023->n_1009048 n_1009023->n_1009024 n_1009023->n_1009024 e_3583079->n_1009049 e_3583079->n_1009028 e_3583079->n_1009031 e_3583095 sink_accretion_table n_1009048->e_3583095 n_1009048->n_1009049 e_3583066 sink_pos e_3583066->n_1009048 n_1009056 SinkSelfGravityHost e_3583066->n_1009056 e_3583066->n_1009056 e_3583072 sink_accretion_radius e_3583072->n_1009048 e_3583095->n_1009050 n_1009049->e_3583066 n_1009049->e_3583095 e_3583067 sink_vel n_1009049->e_3583067 e_3583068 sink_acc_sph n_1009049->e_3583068 e_3583071 sink_angular_momentum n_1009049->e_3583071 e_3583070 sink_mass n_1009049->e_3583070 n_1009049->n_1009050 e_3583064 dt e_3583064->n_1009049 n_1009059 MapEdge e_3583064->n_1009059 n_1009061 ForwardEulerHost e_3583064->n_1009061 n_1009040 MapEdge e_3583064->n_1009040 e_3583064->n_1009034 e_3583067->n_1009061 n_1009060 ForwardEulerHost2Deriv e_3583068->n_1009060 e_3583070->n_1009056 e_3583070->n_1009056 n_1009050->e_3583073 n_1009045 FreeAlloc n_1009050->n_1009045 n_1009045->e_3583077 n_1009045->e_3583077 n_1009046 FreeAlloc n_1009045->n_1009046 n_1009072_end n_1009045->n_1009072_end n_1009046->e_3583078 n_1009047 FreeAlloc n_1009046->n_1009047 n_1009047->e_3583079 n_1009066 SetEdge n_1009047->n_1009066 e_3583099 ext_force_accretion_pos n_1009066->e_3583099 n_1009067 SetEdge n_1009066->n_1009067 n_1009068 MapEdge e_3583099->n_1009068 e_3583099->n_1009069 e_3583100 ext_force_accretion_racc n_1009067->e_3583100 n_1009067->n_1009068 e_3583100->n_1009069 e_3583102 has_ext_force_accretion n_1009068->e_3583102 n_1009072 if_has_ext_force_accretion n_1009068->n_1009072 e_3583102->n_1009072 n_1009072->n_1009019 true n_1009072->n_1009072_end false e_3583101 ext_force_accretion_table n_1009069->e_3583101 n_1009069->n_1009070 e_3583101->n_1009070 n_1009070->e_3583073 n_1009070->n_1009045 n_1009072_end->n_1009063 n_1009063->n_1009058 true n_1009063_end n_1009063->n_1009063_end false n_1009053 SetEdge n_1009058->n_1009053 true n_1009058->n_1009053 true n_1009058_end n_1009058->n_1009058_end false n_1009058->n_1009058_end false e_3583096 G n_1009053->e_3583096 n_1009053->e_3583096 n_1009054 SetEdge n_1009053->n_1009054 n_1009053->n_1009054 e_3583096->n_1009056 e_3583096->n_1009056 e_3583097 epsilon_grav_sink n_1009054->e_3583097 n_1009054->e_3583097 n_1009055 SetEdge n_1009054->n_1009055 n_1009054->n_1009055 e_3583097->n_1009056 e_3583097->n_1009056 e_3583069 sink_acc_ext n_1009055->e_3583069 n_1009055->e_3583069 n_1009055->n_1009056 n_1009055->n_1009056 e_3583069->n_1009060 n_1009056->e_3583069 n_1009056->e_3583069 n_1009056->n_1009058_end n_1009056->n_1009058_end n_1009058_end->n_1009059 e_3583098 dt_half n_1009059->e_3583098 n_1009059->n_1009060 e_3583098->n_1009060 n_1009060->e_3583067 n_1009060->n_1009061 n_1009061->e_3583066 n_1009061->n_1009063_end n_1009063_end->n_1009040 n_1009040->n_1009018 e_3583075 dt_half n_1009040->e_3583075 e_3583075->n_1009028 e_3583075->n_1009029 e_3583075->n_1009031 e_3583075->n_1009032 e_3583080 uint n_1009024->e_3583080 n_1009024->e_3583080 n_1009024->n_1009025 n_1009024->n_1009025 e_3583081 duint n_1009025->e_3583081 n_1009025->e_3583081 n_1009025->n_1009026 n_1009025->n_1009026 e_3583081->n_1009029 e_3583081->n_1009032 n_1009026->n_1009058 e_3583082 hpart n_1009026->e_3583082 n_1009026->e_3583082 n_1009026->n_1009028 n_1009028->e_3583078 n_1009028->n_1009029 n_1009029->e_3583080 n_1009029->n_1009034 n_1009034->e_3583077 n_1009034->n_1009031 n_1009031->e_3583078 n_1009031->n_1009032 n_1009032->e_3583080 n_1009036 FreeAlloc n_1009032->n_1009036 e_3583083 part_to_remove n_1009036->e_3583083 n_1009036->n_1009037 n_1009038 KillParticles e_3583083->n_1009038 n_1009037->e_3583083 n_1009037->n_1009038 n_1009038->n_1009019 n_1009038->e_3583073


.. GENERATED FROM PYTHON SOURCE LINES 116-118 Draw utilities ################################################### .. GENERATED FROM PYTHON SOURCE LINES 118-145 .. code-block:: Python :lineno-start: 118 import matplotlib.pyplot as plt import numpy as np def plot_state(iplot): pos = ctx.collect_data()["xyz"] if shamrock.sys.world_rank() == 0: X = pos[:, 0] Y = pos[:, 1] Z = pos[:, 2] plt.cla() ax.set_xlim3d(bmin[0], bmax[0]) ax.set_ylim3d(bmin[1], bmax[1]) ax.set_zlim3d(bmin[2], bmax[2]) ax.scatter(X, Y, Z, s=1) ax.minorticks_off() ax.set_title(f"t = {model.get_time():.2f} ") plt.savefig(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.png")) .. GENERATED FROM PYTHON SOURCE LINES 146-148 Run the simulation ################################################### .. GENERATED FROM PYTHON SOURCE LINES 148-170 .. code-block:: Python :lineno-start: 148 nstop = 28 # To be increased when empty simulations will be fixed dt_stop = 0.1 t_stop = [i * dt_stop for i in range(nstop + 1)] # Init MPL fig = plt.figure(dpi=120) ax = fig.add_subplot(111, projection="3d") iplot = 0 istop = 0 for ttarg in t_stop: model.evolve_until(ttarg) model.do_vtk_dump(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.vtk"), True) # if do_plots: plot_state(iplot) iplot += 1 istop += 1 plt.close(fig) .. rst-class:: sphx-glr-script-out .. code-block:: none ---------------- 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] .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 13.139 seconds) **Estimated memory usage:** 160 MB .. _sphx_glr_download__as_gen_tests_ci_run_GW_isolated_star.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: run_GW_isolated_star.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: run_GW_isolated_star.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: run_GW_isolated_star.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_