Commit 3d14e53e authored by Matthieu Schaller's avatar Matthieu Schaller

Revert "debugging, fix to rounding error in calculating delta time from scale factors"

This reverts commit 388c4ea3.
parent 73bc7556
......@@ -116,7 +116,7 @@ SNIa_timescale_Gyr = 2.0
expected_energy_released_cgs = np.zeros(n_snapshots)
for i in range(n_snapshots):
age_Gyr = time[i] * unit_time_in_cgs / Gyr_in_cgs
total_sn = SNIa_efficiency * (1.0 - np.exp(-age_Gyr/SNIa_timescale_Gyr)) * star_initial_mass / const_solar_mass
total_sn = SNIa_efficiency * (SNIa_timescale_Gyr * Gyr_in_cgs / unit_time_in_cgs * (1 - exp(-age_Gyr/SNIa_timescale_Gyr))) * star_initial_mass * const_solar_mass
expected_energy_released_cgs[i] = total_sn * energy_per_sn * unit_energy_in_cgs
# Did we get it right?
......@@ -130,7 +130,6 @@ figure()
subplot(111)
plot(time*unit_time_in_cgs/Gyr_in_cgs, total_energy_cgs + total_kinetic_energy_cgs - total_energy_cgs[0] - total_kinetic_energy_cgs[0],color='k', linewidth=0.5, label="SWIFT")
plot(time*unit_time_in_cgs/Gyr_in_cgs, expected_energy_released_cgs,color = 'r', linewidth=0.5, label="expected")
#plot(np.arange(0,2,0.1),(1.0 - np.exp(-np.arange(0,2,0.1)/SNIa_timescale_Gyr)), color = 'k', linewidth = 0.5)
xlabel("Time (Gyr)")
ylabel("Total energy (erg)")
legend()
......
......@@ -375,7 +375,7 @@ void runner_do_stars_ghost(struct runner *r, struct cell *c, int timer) {
double star_age, current_time_begin = -1;
if (with_cosmology) {
star_age = cosmology_get_delta_time_from_scale_factors(
cosmo, sp->birth_scale_factor, (float) cosmo->a); // use float to avoid rounding errors
cosmo, sp->birth_scale_factor, cosmo->a);
} else {
current_time_begin =
get_integer_time_begin(e->ti_current - 1, sp->time_bin) *
......
......@@ -324,7 +324,7 @@ inline static void evolve_SNIa(float log10_min_mass, float log10_max_mass,
sp->mass_init;
sp->to_distribute.num_SNIa =
num_SNIa_per_msun / stars->feedback.const_solar_mass;
num_SNIa_per_msun * stars->feedback.const_solar_mass;
/* compute mass fractions of each metal */
for (int i = 0; i < chemistry_element_count; i++) {
......
......@@ -134,6 +134,8 @@ runner_iact_nonsym_stars_feedback(
const float current_mass = hydro_get_mass(pj);
float new_mass =
current_mass + si->to_distribute.mass * density_weighted_frac;
// for testing energy injection
//new_mass = current_mass;
hydro_set_mass(pj, new_mass);
/* Update total metallicity */
......@@ -224,13 +226,16 @@ runner_iact_nonsym_stars_feedback(
/* Energy feedback */
float u_init = hydro_get_physical_internal_energy(pj, xp, cosmo);
float heating_probability = -1.f, du = 0.f, d_energy = 0.f;
d_energy += si->to_distribute.d_energy * density_weighted_frac;
d_energy += si->to_distribute.d_energy;
if (stars_properties->feedback.continuous_heating) {
// We're doing ONLY continuous heating
d_energy += si->to_distribute.num_SNIa *
stars_properties->feedback.total_energy_SNe *
density_weighted_frac * si->mass_init;
du = d_energy / hydro_get_mass(pj);
hydro_set_physical_internal_energy(pj, xp, cosmo, u_init + du);
hydro_set_drifted_physical_internal_energy(pj, cosmo, u_init + du);
} else {
// We're doing stochastic heating
heating_probability = si->to_distribute.heating_probability;
......@@ -254,12 +259,6 @@ runner_iact_nonsym_stars_feedback(
hydro_set_drifted_physical_internal_energy(pj, cosmo, u_init + du);
}
}
/* Add contribution from thermal and kinetic energy of ejected material (and continuous SNIa feedback) */
u_init = hydro_get_physical_internal_energy(pj, xp, cosmo);
du = d_energy / hydro_get_mass(pj);
hydro_set_physical_internal_energy(pj, xp, cosmo, u_init + du);
hydro_set_drifted_physical_internal_energy(pj, cosmo, u_init + du);
}
#endif /* SWIFT_EAGLE_STARS_IACT_H */
# ID pos_x pos_y pos_z wcount wcount_dh
0 1e-6 1e-6 1e-6 4e-4 1.2e-2
0 1e-6 1e-6 1e-6 1e-4 2e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6
# ID pos_x pos_y pos_z wcount wcount_dh
0 1e-6 1e-6 1e-6 2e-4 1e-2
0 1e-6 1e-6 1e-6 1e-5 2.4e-3
0 1e-6 1e-6 1e-6 1e-6 1e-2
# ID pos_x pos_y pos_z wcount wcount_dh
0 1e-6 1e-6 1e-6 5e-4 1.4e-2
0 1e-6 1e-6 1e-6 1e-5 4e-3
0 1e-6 1e-6 1e-6 1e-6 1e0
# ID pos_x pos_y pos_z wcount wcount_dh
0 1e-6 1e-6 1e-6 5e-4 1.5e-2
0 1e-6 1e-6 1e-6 1e-5 5.86e-3
0 1e-6 1e-6 1e-6 1e-6 1e0
......@@ -76,7 +76,7 @@ int main(int argc, char *argv[]) {
stars_evolve_init(params, &stars_properties);
/* Init spart */
stars_first_init_spart(&sp, &stars_properties);
stars_first_init_spart(&sp);
/* Define an initial stellar mass. (for use when calling the feedback
* functions, the results are presented per initial stellar mass, so the
......
# ID pos_x pos_y pos_z v_x v_y v_z h rho div_v S u P c a_x a_y a_z h_dt v_sig dS/dt du/dt
0 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4
0 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-5 2e-5 2e-5 1e-5 1e-5 1e-5 1e-5
# ID pos_x pos_y pos_z v_x v_y v_z h rho div_v S u P c a_x a_y a_z h_dt v_sig dS/dt du/dt
0 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4
0 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 3.6e-3 2e-3 2e-3 1e-4 1e-4 1e-4 1e-4
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-4 2e-4 2e-4 1e-6 1e-6 1e-6 1e-6
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-6 6.2e-5 4e-4 1.2e-2 1e-5 6e-6 6e-6 8e-6
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-4 1e-4 2e-4 2e-4 1e-4 1e-4 1e-4
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-6 1e-4 2e-4 1e-2 1e-5 3e-6 3e-6 7e-6
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1.5e-3 1e-5 2e-3 6e-5 3e-3 2e-3 2e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-3 1e-6 1e0 1e-6 2e-6 2e-6 2e-6
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 3e-6 1e-4 5e-4 1.4e-2 1.1e-5 3e-6 3e-6 8e-6
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1.5e-6 1.4e-2 1e-5 2e-3 2.5e-4 3e-3 3e-3 3e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e0 1e-6 4e-6 4e-6 4e-6
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 3e-6 1e-4 5e-4 1.5e-2 1.4e-5 3e-6 3e-6 1e-5
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1.5e-6 2.5e-2 1e-5 5.86e-3 1.17e-3 3e-3 8e-3 3e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e0 1e-6 4e-6 4e-6 4e-6
# ID wcount h_dt
0 1e-2 1
0 1e-2 1
0 1e-2 1
# ID wcount h_dt
0 1 2.4e-3
0 1 2.4e-3
0 1 3e-5
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 4e-6 7e-5 1e-3 1.3e-2 2e-4 2e-4 2e-4 2e-4
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-6 2e-4 1e-4 2e-4 6e-4 2e-3 2e-3 2e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-4 1e-6 1e-4 5e-4 2e-4 2e-4 2e-4
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 3e-6 4e-5 1e-3 1e-2 2e-4 1e-4 1e-4 1e-4
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-6 7e-3 1e-4 3e-3 1e-2 6e-3 6e-3 6e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 2e-3 1e-6 1e0 5e-4 3e-3 3e-3 3e-3
# ID pos_x pos_y pos_z v_x v_y v_z h rho div_v S u P c a_x a_y a_z h_dt v_sig dS/dt du/dt
0 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4
0 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-5 1e-5 1e-5 1e-5 1e-5 1e-5 1e-5
# Number of cores: 1
# Iteration | Time | TimeStep | NumPartUpdates | ElapsedTicks (cycle.h)
1 0 0 1841127 2.49541e+11
2 1e-08 1e-08 1841127 2.41972e+11
3 2e-08 1e-08 1841127 2.46989e+11
4 3e-08 1e-08 1841127 2.47354e+11
5 4e-08 1e-08 1841127 2.48764e+11
6 5e-08 1e-08 1841127 2.39867e+11
7 6e-08 1e-08 1841127 2.41151e+11
8 7e-08 1e-08 1841127 2.47176e+11
9 8e-08 1e-08 1841127 2.40415e+11
10 9e-08 1e-08 1841127 2.44373e+11
11 1e-07 1e-08 1841127 2.43127e+11
12 1.1e-07 1e-08 1841127 2.44901e+11
13 1.2e-07 1e-08 1841127 2.41489e+11
14 1.3e-07 1e-08 1841127 2.43391e+11
15 1.4e-07 1e-08 1841127 2.4912e+11
16 1.5e-07 1e-08 1841127 2.48539e+11
17 1.6e-07 1e-08 1841127 2.45593e+11
18 1.7e-07 1e-08 1841127 2.46147e+11
19 1.8e-07 1e-08 1841127 2.38684e+11
20 1.9e-07 1e-08 1841127 2.4488e+11
21 2e-07 1e-08 1841127 2.37401e+11
22 2.1e-07 1e-08 1841127 2.4023e+11
23 2.2e-07 1e-08 1841127 2.38668e+11
24 2.3e-07 1e-08 1841127 2.47569e+11
25 2.4e-07 1e-08 1841127 2.38949e+11
26 2.5e-07 1e-08 1841127 2.42014e+11
27 2.6e-07 1e-08 1841127 2.38534e+11
28 2.7e-07 1e-08 1841127 2.38274e+11
29 2.8e-07 1e-08 1841127 2.49901e+11
30 2.9e-07 1e-08 1841127 2.39624e+11
31 3e-07 1e-08 1841127 2.43427e+11
32 3.1e-07 1e-08 1841127 2.4817e+11
33 3.2e-07 1e-08 1841127 2.44629e+11
34 3.3e-07 1e-08 1841127 2.40832e+11
35 3.4e-07 1e-08 1841127 2.45043e+11
36 3.5e-07 1e-08 1841127 2.43631e+11
37 3.6e-07 1e-08 1841127 2.43124e+11
38 3.7e-07 1e-08 1841127 2.49266e+11
39 3.8e-07 1e-08 1841127 2.40842e+11
40 3.9e-07 1e-08 1841127 2.46852e+11
41 4e-07 1e-08 1841127 2.44451e+11
42 4.1e-07 1e-08 1841127 2.4711e+11
43 4.2e-07 1e-08 1841127 2.42068e+11
44 4.3e-07 1e-08 1841127 2.48401e+11
45 4.4e-07 1e-08 1841127 2.40454e+11
46 4.5e-07 1e-08 1841127 2.41779e+11
47 4.6e-07 1e-08 1841127 2.47104e+11
48 4.7e-07 1e-08 1841127 2.48655e+11
49 4.8e-07 1e-08 1841127 2.37705e+11
50 4.9e-07 1e-08 1841127 2.50405e+11
51 5e-07 1e-08 1841127 2.43756e+11
52 5.1e-07 1e-08 1841127 2.38041e+11
53 5.2e-07 1e-08 1841127 2.39564e+11
54 5.3e-07 1e-08 1841127 2.45563e+11
55 5.4e-07 1e-08 1841127 2.48208e+11
56 5.5e-07 1e-08 1841127 2.41248e+11
57 5.6e-07 1e-08 1841127 2.37633e+11
58 5.7e-07 1e-08 1841127 2.37176e+11
59 5.8e-07 1e-08 1841127 2.42828e+11
60 5.9e-07 1e-08 1841127 2.37799e+11
61 6e-07 1e-08 1841127 2.40112e+11
62 6.1e-07 1e-08 1841127 2.49319e+11
63 6.2e-07 1e-08 1841127 2.49081e+11
64 6.3e-07 1e-08 1841127 2.477e+11
65 6.4e-07 1e-08 1841127 2.40199e+11
66 6.5e-07 1e-08 1841127 2.44252e+11
67 6.6e-07 1e-08 1841127 2.41898e+11
68 6.7e-07 1e-08 1841127 2.46553e+11
69 6.8e-07 1e-08 1841127 2.43482e+11
70 6.9e-07 1e-08 1841127 2.45474e+11
71 7e-07 1e-08 1841127 2.43577e+11
72 7.1e-07 1e-08 1841127 2.37344e+11
73 7.2e-07 1e-08 1841127 2.43476e+11
74 7.3e-07 1e-08 1841127 2.48881e+11
75 7.4e-07 1e-08 1841127 2.48735e+11
76 7.5e-07 1e-08 1841127 2.46232e+11
77 7.6e-07 1e-08 1841127 2.40849e+11
78 7.7e-07 1e-08 1841127 2.46472e+11
79 7.8e-07 1e-08 1841127 2.45576e+11
80 7.9e-07 1e-08 1841127 2.41463e+11
81 8e-07 1e-08 1841127 2.45939e+11
82 8.1e-07 1e-08 1841127 2.39002e+11
83 8.2e-07 1e-08 1841127 2.43038e+11
84 8.3e-07 1e-08 1841127 2.48318e+11
85 8.4e-07 1e-08 1841127 2.47645e+11
86 8.5e-07 1e-08 1841127 2.41054e+11
87 8.6e-07 1e-08 1841127 2.39986e+11
88 8.7e-07 1e-08 1841127 2.48552e+11
89 8.8e-07 1e-08 1841127 2.4166e+11
90 8.9e-07 1e-08 1841127 2.45632e+11
91 9e-07 1e-08 1841127 2.49184e+11
92 9.1e-07 1e-08 1841127 2.44541e+11
93 9.2e-07 1e-08 1841127 2.45425e+11
94 9.3e-07 1e-08 1841127 2.48154e+11
95 9.4e-07 1e-08 1841127 2.42583e+11
96 9.5e-07 1e-08 1841127 2.48707e+11
97 9.6e-07 1e-08 1841127 2.42242e+11
98 9.7e-07 1e-08 1841127 2.47416e+11
99 9.8e-07 1e-08 1841127 2.45948e+11
100 9.9e-07 1e-08 1841127 2.4869e+11
101 1e-06 1e-08 1841127 2.43131e+11
102 1.01e-06 1e-08 1841127 2.40918e+11
103 1.02e-06 1e-08 1841127 2.49165e+11
104 1.03e-06 1e-08 1841127 2.48708e+11
105 1.04e-06 1e-08 1841127 2.39101e+11
106 1.05e-06 1e-08 1841127 2.48381e+11
107 1.06e-06 1e-08 1841127 2.45242e+11
108 1.07e-06 1e-08 1841127 2.42491e+11
109 1.08e-06 1e-08 1841127 2.45098e+11
110 1.09e-06 1e-08 1841127 2.40874e+11
111 1.1e-06 1e-08 1841127 2.47036e+11
112 1.11e-06 1e-08 1841127 2.40896e+11
113 1.12e-06 1e-08 1841127 2.42822e+11
114 1.13e-06 1e-08 1841127 2.40389e+11
115 1.14e-06 1e-08 1841127 2.39214e+11
116 1.15e-06 1e-08 1841127 2.40481e+11
117 1.16e-06 1e-08 1841127 2.44091e+11
118 1.17e-06 1e-08 1841127 2.42227e+11
119 1.18e-06 1e-08 1841127 2.39134e+11
120 1.19e-06 1e-08 1841127 2.48641e+11
121 1.2e-06 1e-08 1841127 2.38446e+11
122 1.21e-06 1e-08 1841127 2.38798e+11
123 1.22e-06 1e-08 1841127 2.43307e+11
124 1.23e-06 1e-08 1841127 2.46562e+11
125 1.24e-06 1e-08 1841127 2.49537e+11
126 1.25e-06 1e-08 1841127 2.49218e+11
127 1.26e-06 1e-08 1841127 2.4558e+11
128 1.27e-06 1e-08 1841127 2.41956e+11
129 1.28e-06 1e-08 1841127 2.46577e+11
130 1.29e-06 1e-08 1841127 2.41805e+11
131 1.3e-06 1e-08 1841127 2.4071e+11
132 1.31e-06 1e-08 1841127 2.4009e+11
133 1.32e-06 1e-08 1841127 2.44535e+11
134 1.33e-06 1e-08 1841127 2.41269e+11
135 1.34e-06 1e-08 1841127 2.38852e+11
136 1.35e-06 1e-08 1841127 2.46536e+11
137 1.36e-06 1e-08 1841127 2.38541e+11
138 1.37e-06 1e-08 1841127 2.47088e+11
139 1.38e-06 1e-08 1841127 2.3928e+11
140 1.39e-06 1e-08 1841127 2.46866e+11
141 1.4e-06 1e-08 1841127 2.37898e+11
142 1.41e-06 1e-08 1841127 2.49495e+11
143 1.42e-06 1e-08 1841127 2.378e+11
144 1.43e-06 1e-08 1841127 2.43855e+11
145 1.44e-06 1e-08 1841127 2.39209e+11
146 1.45e-06 1e-08 1841127 2.40331e+11
147 1.46e-06 1e-08 1841127 2.47196e+11
148 1.47e-06 1e-08 1841127 2.46447e+11
149 1.48e-06 1e-08 1841127 2.45558e+11
150 1.49e-06 1e-08 1841127 2.49346e+11
151 1.5e-06 1e-08 1841127 2.45179e+11
152 1.51e-06 1e-08 1841127 2.46992e+11
153 1.52e-06 1e-08 1841127 2.38439e+11
154 1.53e-06 1e-08 1841127 2.38665e+11
155 1.54e-06 1e-08 1841127 2.44009e+11
156 1.55e-06 1e-08 1841127 2.37937e+11
157 1.56e-06 1e-08 1841127 2.37968e+11
158 1.57e-06 1e-08 1841127 2.39537e+11
159 1.58e-06 1e-08 1841127 2.43e+11
160 1.59e-06 1e-08 1841127 2.48704e+11
161 1.6e-06 1e-08 1841127 2.44625e+11
162 1.61e-06 1e-08 1841127 2.46884e+11
163 1.62e-06 1e-08 1841127 2.37603e+11
164 1.63e-06 1e-08 1841127 2.39332e+11
165 1.64e-06 1e-08 1841127 2.49885e+11
166 1.65e-06 1e-08 1841127 2.40072e+11
167 1.66e-06 1e-08 1841127 2.48453e+11
168 1.67e-06 1e-08 1841127 2.38501e+11
169 1.68e-06 1e-08 1841127 2.4996e+11
170 1.69e-06 1e-08 1841127 2.37827e+11
171 1.7e-06 1e-08 1841127 2.48082e+11
172 1.71e-06 1e-08 1841127 2.38273e+11
173 1.72e-06 1e-08 1841127 2.36946e+11
174 1.73e-06 1e-08 1841127 2.48783e+11
175 1.74e-06 1e-08 1841127 2.44543e+11
176 1.75e-06 1e-08 1841127 2.39674e+11
177 1.76e-06 1e-08 1841127 2.39833e+11
178 1.77e-06 1e-08 1841127 2.41799e+11
179 1.78e-06 1e-08 1841127 2.48682e+11
180 1.79e-06 1e-08 1841127 2.4746e+11
181 1.8e-06 1e-08 1841127 2.41678e+11
182 1.81e-06 1e-08 1841127 2.44008e+11
183 1.82e-06 1e-08 1841127 2.4564e+11
184 1.83e-06 1e-08 1841127 2.42708e+11
185 1.84e-06 1e-08 1841127 2.45847e+11
186 1.85e-06 1e-08 1841127 2.38661e+11
187 1.86e-06 1e-08 1841127 2.43879e+11
188 1.87e-06 1e-08 1841127 2.4662e+11
189 1.88e-06 1e-08 1841127 2.40857e+11
190 1.89e-06 1e-08 1841127 2.49643e+11
191 1.9e-06 1e-08 1841127 2.44895e+11
192 1.91e-06 1e-08 1841127 2.48156e+11
193 1.92e-06 1e-08 1841127 2.4673e+11
194 1.93e-06 1e-08 1841127 2.45526e+11
195 1.94e-06 1e-08 1841127 2.37373e+11
196 1.95e-06 1e-08 1841127 2.46674e+11
197 1.96e-06 1e-08 1841127 2.47827e+11
198 1.97e-06 1e-08 1841127 2.48903e+11
199 1.98e-06 1e-08 1841127 2.48118e+11
200 1.99e-06 1e-08 1841127 2.47462e+11
201 2e-06 1e-08 1841127 2.49679e+11
202 2.01e-06 1e-08 1841127 2.4648e+11
203 2.02e-06 1e-08 1841127 2.48482e+11
204 2.03e-06 1e-08 1841127 2.50963e+11
205 2.04e-06 1e-08 1841127 2.45634e+11
206 2.05e-06 1e-08 1841127 2.4314e+11
207 2.06e-06 1e-08 1841127 2.39173e+11
208 2.07e-06 1e-08 1841127 2.47692e+11
209 2.08e-06 1e-08 1841127 2.48322e+11
210 2.09e-06 1e-08 1841127 2.38007e+11
211 2.1e-06 1e-08 1841127 2.45285e+11
212 2.11e-06 1e-08 1841127 2.40321e+11
213 2.12e-06 1e-08 1841127 2.4532e+11
214 2.13e-06 1e-08 1841127 2.41105e+11
215 2.14e-06 1e-08 1841127 2.45989e+11
216 2.15e-06 1e-08 1841127 2.41174e+11
217 2.16e-06 1e-08 1841127 2.41287e+11
218 2.17e-06 1e-08 1841127 2.40051e+11
219 2.18e-06 1e-08 1841127 2.37963e+11
220 2.19e-06 1e-08 1841127 2.45039e+11
221 2.2e-06 1e-08 1841127 2.40084e+11
222 2.21e-06 1e-08 1841127 2.45689e+11
223 2.22e-06 1e-08 1841127 2.4369e+11
224 2.23e-06 1e-08 1841127 2.49472e+11
225 2.24e-06 1e-08 1841127 2.40499e+11
226 2.25e-06 1e-08 1841127 2.4389e+11
227 2.26e-06 1e-08 1841127 2.46203e+11
228 2.27e-06 1e-08 1841127 2.38354e+11
229 2.28e-06 1e-08 1841127 2.43015e+11
230 2.29e-06 1e-08 1841127 2.44692e+11
231 2.3e-06 1e-08 1841127 2.49186e+11
232 2.31e-06 1e-08 1841127 2.43954e+11
233 2.32e-06 1e-08 1841127 2.41146e+11
234 2.33e-06 1e-08 1841127 2.47836e+11
235 2.34e-06 1e-08 1841127 2.42775e+11
236 2.35e-06 1e-08 1841127 2.37132e+11
237 2.36e-06 1e-08 1841127 2.41883e+11
238 2.37e-06 1e-08 1841127 2.44703e+11
239 2.38e-06 1e-08 1841127 2.47572e+11
240 2.39e-06 1e-08 1841127 2.39929e+11
241 2.4e-06 1e-08 1841127 2.46093e+11
242 2.41e-06 1e-08 1841127 2.43104e+11
243 2.42e-06 1e-08 1841127 2.43483e+11
244 2.43e-06 1e-08 1841127 2.40992e+11
245 2.44e-06 1e-08 1841127 2.46368e+11
246 2.45e-06 1e-08 1841127 2.39375e+11
247 2.46e-06 1e-08 1841127 2.45623e+11
248 2.47e-06 1e-08 1841127 2.37358e+11
249 2.48e-06 1e-08 1841127 2.42598e+11
250 2.49e-06 1e-08 1841127 2.47466e+11
251 2.5e-06 1e-08 1841127 2.45522e+11
252 2.51e-06 1e-08 1841127 2.45115e+11
253 2.52e-06 1e-08 1841127 2.41498e+11
254 2.53e-06 1e-08 1841127 2.39537e+11
255 2.54e-06 1e-08 1841127 2.44713e+11
256 2.55e-06 1e-08 1841127 2.44787e+11
257 2.56e-06 1e-08 1841127 2.46263e+11
# Number of cores: 10
# Iteration | Time | TimeStep | NumPartUpdates | ElapsedTicks (cycle.h)
1 0 0 183963 3.45886e+10
2 1e-08 1e-08 183963 3.19859e+10
3 2e-08 1e-08 183963 3.27121e+10
4 3e-08 1e-08 183963 3.2731e+10
5 4e-08 1e-08 183963 3.29447e+10
6 5e-08 1e-08 183963 3.17616e+10
7 6e-08 1e-08 183963 3.23126e+10
8 7e-08 1e-08 183964 3.31991e+10
9 8e-08 1e-08 183963 3.23147e+10
10 9e-08 1e-08 183963 3.28677e+10
11 1e-07 1e-08 183964 3.26094e+10
12 1.1e-07 1e-08 183963 3.28295e+10
13 1.2e-07 1e-08 183963 3.23607e+10
14 1.3e-07 1e-08 183963 3.279e+10
15 1.4e-07 1e-08 186256 3.35803e+10
16 1.5e-07 1e-08 186256 3.31371e+10
17 1.6e-07 1e-08 186255 3.2841e+10
18 1.7e-07 1e-08 186255 3.27458e+10
19 1.8e-07 1e-08 186255 3.1756e+10
20 1.9e-07 1e-08 186255 3.26209e+10
21 2e-07 1e-08 186255 3.15661e+10
22 2.1e-07 1e-08 186255 3.22889e+10
23 2.2e-07 1e-08 186255 3.18118e+10
24 2.3e-07 1e-08 186255 3.29793e+10
25 2.4e-07 1e-08 186255 3.18468e+10
26 2.5e-07 1e-08 186255 3.21964e+10
27 2.6e-07 1e-08 186255 3.17404e+10
28 2.7e-07 1e-08 186254 3.17026e+10
29 2.8e-07 1e-08 186254 3.3258e+10
30 2.9e-07 1e-08 186254 3.19298e+10
31 3e-07 1e-08 186254 3.24714e+10
32 3.1e-07 1e-08 186254 3.30908e+10
33 3.2e-07 1e-08 186254 3.2542e+10
34 3.3e-07 1e-08 186254 3.20705e+10
35 3.4e-07 1e-08 186254 3.25962e+10
36 3.5e-07 1e-08 186254 3.26849e+10
37 3.6e-07 1e-08 186255 3.25801e+10
38 3.7e-07 1e-08 186255 3.3074e+10
39 3.8e-07 1e-08 186255 3.20638e+10
40 3.9e-07 1e-08 188306 3.31353e+10
41 4e-07 1e-08 188305 3.30021e+10
42 4.1e-07 1e-08 188306 3.3234e+10
43 4.2e-07 1e-08 188305 3.22633e+10
44 4.3e-07 1e-08 188306 3.32161e+10
45 4.4e-07 1e-08 188306 3.2068e+10
46 4.5e-07 1e-08 188305 3.23431e+10
47 4.6e-07 1e-08 188305 3.29472e+10
48 4.7e-07 1e-08 188305 3.31947e+10
49 4.8e-07 1e-08 188305 3.1646e+10
50 4.9e-07 1e-08 186254 3.38917e+10
51 5e-07 1e-08 186254 3.22862e+10
52 5.1e-07 1e-08 186254 3.14376e+10
53 5.2e-07 1e-08 186254 3.14109e+10
54 5.3e-07 1e-08 186253 3.22141e+10
55 5.4e-07 1e-08 186253 3.24937e+10
56 5.5e-07 1e-08 186253 3.15856e+10
57 5.6e-07 1e-08 186253 3.12077e+10
58 5.7e-07 1e-08 186253 3.11483e+10
59 5.8e-07 1e-08 186253 3.18011e+10
60 5.9e-07 1e-08 186254 3.14936e+10
61 6e-07 1e-08 186254 3.16083e+10
62 6.1e-07 1e-08 186254 3.27944e+10
63 6.2e-07 1e-08 186254 3.30691e+10
64 6.3e-07 1e-08 186254 3.25708e+10
65 6.4e-07 1e-08 186254 3.16659e+10
66 6.5e-07 1e-08 190790 3.3257e+10
67 6.6e-07 1e-08 188305 3.21765e+10
68 6.7e-07 1e-08 188305 3.24206e+10
69 6.8e-07 1e-08 188304 3.22709e+10
70 6.9e-07 1e-08 188304 3.23127e+10
71 7e-07 1e-08 188303 3.20304e+10
72 7.1e-07 1e-08 188303 3.12704e+10
73 7.2e-07 1e-08 188303 3.24073e+10
74 7.3e-07 1e-08 188304 3.27518e+10
75 7.4e-07 1e-08 188303 3.27234e+10
76 7.5e-07 1e-08 188303 3.24965e+10
77 7.6e-07 1e-08 188304 3.16855e+10
78 7.7e-07 1e-08 188304 3.26358e+10
79 7.8e-07 1e-08 188304 3.24595e+10
80 7.9e-07 1e-08 188304 3.18251e+10
81 8e-07 1e-08 188305 3.23904e+10
82 8.1e-07 1e-08 188304 3.14037e+10
83 8.2e-07 1e-08 188304 3.21084e+10
84 8.3e-07 1e-08 188303 3.27458e+10
85 8.4e-07 1e-08 188303 3.26365e+10
86 8.5e-07 1e-08 188303 3.16628e+10
87 8.6e-07 1e-08 188304 3.16257e+10
88 8.7e-07 1e-08 188304 3.30135e+10
89 8.8e-07 1e-08 188303 3.19482e+10
90 8.9e-07 1e-08 188303 3.22844e+10
91 9e-07 1e-08 188303 3.27983e+10
92 9.1e-07 1e-08 188303 3.29996e+10
93 9.2e-07 1e-08 190789 3.28601e+10
94 9.3e-07 1e-08 190789 3.32404e+10
95 9.4e-07 1e-08 190789 3.24515e+10
96 9.5e-07 1e-08 190789 3.32833e+10
97 9.6e-07 1e-08 190789 3.2485e+10
98 9.7e-07 1e-08 190789 3.33455e+10
99 9.8e-07 1e-08 190788 3.2905e+10
100 9.9e-07 1e-08 190788 3.32486e+10
101 1e-06 1e-08 190788 3.24984e+10
102 1.01e-06 1e-08 190788 3.28478e+10
103 1.02e-06 1e-08 190788 3.32879e+10
104 1.03e-06 1e-08 190788 3.32668e+10
105 1.04e-06 1e-08 190788 3.1933e+10
106 1.05e-06 1e-08 190787 3.32779e+10
107 1.06e-06 1e-08 190788 3.31057e+10
108 1.07e-06 1e-08 190788 3.27557e+10
109 1.08e-06 1e-08 190788 3.28775e+10
110 1.09e-06 1e-08 190788 3.23192e+10
111 1.1e-06 1e-08 190788 3.29314e+10
112 1.11e-06 1e-08 190788 3.22382e+10
113 1.12e-06 1e-08 190788 3.24646e+10
114 1.13e-06 1e-08 195902 3.25285e+10
115 1.14e-06 1e-08 190786 3.26959e+10
116 1.15e-06 1e-08 190786 3.24673e+10
117 1.16e-06 1e-08 190786 3.32737e+10
118 1.17e-06 1e-08 190786 3.27206e+10
119 1.18e-06 1e-08 190786 3.22661e+10
120 1.19e-06 1e-08 190786 3.32953e+10
121 1.2e-06 1e-08 190785 3.19402e+10
122 1.21e-06 1e-08 190785 3.22452e+10
123 1.22e-06 1e-08 190785 3.25929e+10
124 1.23e-06 1e-08 190785 3.29408e+10
125 1.24e-06 1e-08 190785 3.33972e+10
126 1.25e-06 1e-08 190785 3.33906e+10
127 1.26e-06 1e-08 190785 3.28745e+10
128 1.27e-06 1e-08 190785 3.25483e+10
129 1.28e-06 1e-08 190785 3.30151e+10
130 1.29e-06 1e-08 190784 3.25054e+10
131 1.3e-06 1e-08 195900 3.24412e+10
132 1.31e-06 1e-08 195899 3.22191e+10
133 1.32e-06 1e-08 195900 3.28022e+10
134 1.33e-06 1e-08 196352 3.29283e+10
135 1.34e-06 1e-08 196352 3.19157e+10
136 1.35e-06 1e-08 196352 3.31593e+10
137 1.36e-06 1e-08 190785 3.26543e+10
138 1.37e-06 1e-08 190785 3.35142e+10
139 1.38e-06 1e-08 190785 3.24383e+10
140 1.39e-06 1e-08 190786 3.34023e+10
141 1.4e-06 1e-08 196716 3.2094e+10
142 1.41e-06 1e-08 196716 3.33375e+10
143 1.42e-06 1e-08 196716 3.18309e+10
144 1.43e-06 1e-08 196716 3.26702e+10
145 1.44e-06 1e-08 196716 3.25279e+10
146 1.45e-06 1e-08 196716 3.24923e+10
147 1.46e-06 1e-08 197971 3.31385e+10
148 1.47e-06 1e-08 197972 3.28822e+10
149 1.48e-06 1e-08 197971 3.29654e+10
150 1.49e-06 1e-08 197971 3.3305e+10
151 1.5e-06 1e-08 197971 3.35048e+10
152 1.51e-06 1e-08 199044 3.36316e+10
153 1.52e-06 1e-08 197972 3.29329e+10
154 1.53e-06 1e-08 196354 3.2975e+10
155 1.54e-06 1e-08 197971 3.34298e+10
156 1.55e-06 1e-08 197971 3.26334e+10
157 1.56e-06 1e-08 196716 3.30714e+10
158 1.57e-06 1e-08 196716 3.31398e+10
159 1.58e-06 1e-08 196717 3.35842e+10
160 1.59e-06 1e-08 196717 3.51079e+10
161 1.6e-06 1e-08 196717 3.37391e+10
162 1.61e-06 1e-08 196718 3.42525e+10
163 1.62e-06 1e-08 196718 3.29824e+10
164 1.63e-06 1e-08 196719 3.32355e+10
165 1.64e-06 1e-08 198705 3.47048e+10
166 1.65e-06 1e-08 198706 3.31558e+10
167 1.66e-06 1e-08 197975 3.46702e+10
168 1.67e-06 1e-08 197975 3.49663e+10
169 1.68e-06 1e-08 197976 3.47679e+10
170 1.69e-06 1e-08 197977 3.27361e+10
171 1.7e-06 1e-08 197977 3.70371e+10
172 1.71e-06 1e-08 195905 3.4716e+10
173 1.72e-06 1e-08 195905 3.41017e+10
174 1.73e-06 1e-08 195905 3.58008e+10
175 1.74e-06 1e-08 195905 3.52594e+10
176 1.75e-06 1e-08 195906 3.46455e+10
177 1.76e-06 1e-08 195905 3.47926e+10
178 1.77e-06 1e-08 198706 3.39486e+10
179 1.78e-06 1e-08 198706 3.46261e+10
180 1.79e-06 1e-08 198706 3.45316e+10
181 1.8e-06 1e-08 198706 3.34237e+10
182 1.81e-06 1e-08 198705 3.39301e+10
183 1.82e-06 1e-08 198705 3.4729e+10
184 1.83e-06 1e-08 198705 3.36368e+10
185 1.84e-06 1e-08 198705 3.40141e+10
186 1.85e-06 1e-08 199045 3.25626e+10
187 1.86e-06 1e-08 199045 3.35973e+10
188 1.87e-06 1e-08 199045 3.39872e+10
189 1.88e-06 1e-08 199202 3.33767e+10
190 1.89e-06 1e-08 199202 3.43821e+10
191 1.9e-06 1e-08 199201 3.4279e+10
192 1.91e-06 1e-08 199202 3.40044e+10
193 1.92e-06 1e-08 199203 3.38974e+10
194 1.93e-06 1e-08 196718 3.36556e+10
195 1.94e-06 1e-08 196718 3.22672e+10
196 1.95e-06 1e-08 196718 3.35574e+10
197 1.96e-06 1e-08 199202 3.44555e+10
198 1.97e-06 1e-08 199202 3.44261e+10
199 1.98e-06 1e-08 199202 3.42634e+10
200 1.99e-06 1e-08 199202 3.40627e+10
201 2e-06 1e-08 199202 3.45602e+10
202 2.01e-06 1e-08 199202 3.39244e+10
203 2.02e-06 1e-08 199202 3.42394e+10
204 2.03e-06 1e-08 197973 3.49728e+10
205 2.04e-06 1e-08 197973 3.36462e+10
206 2.05e-06 1e-08 197973 3.31425e+10
207 2.06e-06 1e-08 197974 3.25651e+10
208 2.07e-06 1e-08 197974 3.39305e+10
209 2.08e-06 1e-08 197974 3.41107e+10
210 2.09e-06 1e-08 197974 3.35054e+10
211 2.1e-06 1e-08 199203 3.52838e+10
212 2.11e-06 1e-08 198706 3.35741e+10
213 2.12e-06 1e-08 198707 3.39315e+10
214 2.13e-06 1e-08 198707 3.34991e+10
215 2.14e-06 1e-08 198707 3.41083e+10
216 2.15e-06 1e-08 198706 3.32681e+10
217 2.16e-06 1e-08 198705 3.32514e+10
218 2.17e-06 1e-08 198705 3.30284e+10
219 2.18e-06 1e-08 199046 3.274e+10
220 2.19e-06 1e-08 199046 3.37776e+10
221 2.2e-06 1e-08 199046 3.30446e+10
222 2.21e-06 1e-08 199203 3.38451e+10
223 2.22e-06 1e-08 199203 3.36645e+10
224 2.23e-06 1e-08 199203 3.4392e+10
225 2.24e-06 1e-08 199203 3.28791e+10
226 2.25e-06 1e-08 199203 3.35251e+10
227 2.26e-06 1e-08 199203 3.38901e+10
228 2.27e-06 1e-08 199203 3.26563e+10
229 2.28e-06 1e-08 199203 3.34862e+10
230 2.29e-06 1e-08 199203 3.36889e+10
231 2.3e-06 1e-08 200910 3.59469e+10
232 2.31e-06 1e-08 200910 3.52715e+10
233 2.32e-06 1e-08 200910 3.44258e+10
234 2.33e-06 1e-08 200910 3.55451e+10
235 2.34e-06 1e-08 200909 3.47334e+10
236 2.35e-06 1e-08 200908 3.38269e+10
237 2.36e-06 1e-08 200907 3.46175e+10
238 2.37e-06 1e-08 200908 3.49418e+10
239 2.38e-06 1e-08 200908 3.54627e+10
240 2.39e-06 1e-08 196719 3.52117e+10
241 2.4e-06 1e-08 196718 3.56858e+10
242 2.41e-06 1e-08 196718 3.52443e+10
243 2.42e-06 1e-08 196719 3.52405e+10
244 2.43e-06 1e-08 196718 3.48564e+10
245 2.44e-06 1e-08 196718 3.5518e+10
246 2.45e-06 1e-08 200907 3.46646e+10
247 2.46e-06 1e-08 200907 3.58395e+10
248 2.47e-06 1e-08 200907 3.42763e+10
249 2.48e-06 1e-08 200907 3.50998e+10
250 2.49e-06 1e-08 200906 3.55626e+10
251 2.5e-06 1e-08 197972 3.58873e+10
252 2.51e-06 1e-08 200906 3.52845e+10
253 2.52e-06 1e-08 200906 3.45777e+10
254 2.53e-06 1e-08 200906 3.41315e+10
255 2.54e-06 1e-08 200906 3.51243e+10
256 2.55e-06 1e-08 200906 3.53302e+10
257 2.56e-06 1e-08 200906 3.56577e+10<