Skip to content

GitLab

Commit 388c4ea3 authored by Alexei Borissov's avatar Alexei Borissov
Browse files

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

parent 8818c220
Expand all Hide whitespace changes
Inline Side-by-side
Showing with 19 additions and 11481 deletions
examples/SubgridTests/StellarEvolution/check_continuous_heating.py
View file @ 388c4ea3
......@@ -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 * (SNIa_timescale_Gyr * Gyr_in_cgs / unit_time_in_cgs * (1 - exp(-age_Gyr/SNIa_timescale_Gyr))) * star_initial_mass * const_solar_mass
total_sn = SNIa_efficiency * (1.0 - np.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,6 +130,7 @@ 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()
......
src/runner.c
View file @ 388c4ea3
......@@ -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, cosmo->a);
cosmo, sp->birth_scale_factor, (float) cosmo->a); // use float to avoid rounding errors
} else {
current_time_begin =
get_integer_time_begin(e->ti_current - 1, sp->time_bin) *
......@@ -704,6 +704,7 @@ void runner_do_star_formation(struct runner *r, struct cell *c, int timer) {
/* Copy the properties of the gas particle to the star particle */
star_formation_copy_properties(p, xp, sp, e, sf_props, cosmo,
with_cosmology);
message("formed star %llu", sp->id);
}
}
......
src/stars/EAGLE/stars.h
View file @ 388c4ea3
......@@ -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++) {
......
src/stars/EAGLE/stars_iact.h
View file @ 388c4ea3
......@@ -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,16 +226,17 @@ 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;
d_energy += si->to_distribute.d_energy * density_weighted_frac;
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);
// for testing continuous energy injection
//d_energy = si->to_distribute.num_SNIa *
// stars_properties->feedback.total_energy_SNe *
// density_weighted_frac * si->mass_init;
} else {
// We're doing stochastic heating
heating_probability = si->to_distribute.heating_probability;
......@@ -257,6 +260,12 @@ 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 */
tests/star_tolerance_27_normal.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/star_tolerance_27_perturbed.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/star_tolerance_27_perturbed_h.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/star_tolerance_27_perturbed_h2.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/testFeedback.c
View file @ 388c4ea3
......@@ -76,7 +76,7 @@ int main(int argc, char *argv[]) {
stars_evolve_init(params, &stars_properties);
/* Init spart */
stars_first_init_spart(&sp);
stars_first_init_spart(&sp, &stars_properties);
/* Define an initial stellar mass. (for use when calling the feedback
* functions, the results are presented per initial stellar mass, so the
......
tests/tolerance_125_normal.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_125_perturbed.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_27_normal.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_27_perturbed.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_27_perturbed_h.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_27_perturbed_h2.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_pair_active.dat deleted 100644 → 0
View file @ 8818c220
# ID wcount h_dt
0 1e-2 1
0 1e-2 1
0 1e-2 1
tests/tolerance_pair_force_active.dat deleted 100644 → 0
View file @ 8818c220
# ID wcount h_dt
0 1 2.4e-3
0 1 2.4e-3
0 1 3e-5
tests/tolerance_periodic_BC_normal.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_periodic_BC_perturbed.dat deleted 100644 → 0
View file @ 8818c220
# 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
tests/tolerance_testInteractions.dat deleted 100644 → 0
View file @ 8818c220
# 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
theory/Papers/SPHERIC_2013/data/Cosmo_fixed/performance_1.dat deleted 100644 → 0
View file @ 8818c220
# 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
theory/Papers/SPHERIC_2013/data/Cosmo_fixed/performance_10.dat deleted 100644 → 0
View file @ 8818c220
# 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