diff --git a/latex/images/particles_left_narrow_sweep.pdf b/latex/images/particles_left_narrow_sweep.pdf index f049bdb..d5c2bbc 100644 Binary files a/latex/images/particles_left_narrow_sweep.pdf and b/latex/images/particles_left_narrow_sweep.pdf differ diff --git a/latex/images/particles_left_narrow_sweep_fine.pdf b/latex/images/particles_left_narrow_sweep_fine.pdf new file mode 100644 index 0000000..31da0eb Binary files /dev/null and b/latex/images/particles_left_narrow_sweep_fine.pdf differ diff --git a/latex/images/particles_left_narrow_sweep_interactions.pdf b/latex/images/particles_left_narrow_sweep_interactions.pdf index cb06e28..f68274c 100644 Binary files a/latex/images/particles_left_narrow_sweep_interactions.pdf and b/latex/images/particles_left_narrow_sweep_interactions.pdf differ diff --git a/latex/images/particles_left_narrow_sweep_interactions_fine.pdf b/latex/images/particles_left_narrow_sweep_interactions_fine.pdf new file mode 100644 index 0000000..40d09e3 Binary files /dev/null and b/latex/images/particles_left_narrow_sweep_interactions_fine.pdf differ diff --git a/latex/images/particles_left_wide_sweep.pdf b/latex/images/particles_left_wide_sweep.pdf index 5cf0bff..c757fc8 100644 Binary files a/latex/images/particles_left_wide_sweep.pdf and b/latex/images/particles_left_wide_sweep.pdf differ diff --git a/src/main.cpp b/src/main.cpp index 0537470..e7e982b 100644 --- a/src/main.cpp +++ b/src/main.cpp @@ -42,6 +42,7 @@ vec3 analytical_solution_particle_1(double t) double w_n = (w_0 - std::sqrt(w_0 * w_0 - 2. * w_z2)) / 2.; double A_p = (25. + w_n * 20.) / (w_n - w_p); double A_n = -(25. + w_p * 20.) / (w_n - w_p); + std::cout << A_p << "," << A_n << std::endl; std::complex f = A_p * std::exp(std::complex(0., -w_p * t)) + A_n * std::exp(std::complex(0., -w_n * t)); @@ -211,7 +212,7 @@ void potential_resonance_narrow_sweep() double amplitudes[]{.1, .4, .7}; - double freq_start = 1.; + double freq_start = 1.1; double freq_end = 1.7; double freq_increment = .002; size_t freq_iterations = @@ -268,7 +269,7 @@ void potential_resonance_narrow_sweep_interaction() double amplitudes[]{.1, .4, .7}; - double freq_start = 1.; + double freq_start = 1.1; double freq_end = 1.7; double freq_increment = .002; size_t freq_iterations = diff --git a/src/scripts/plot_particles_left.py b/src/scripts/plot_particles_left.py index 125b73d..5116be7 100644 --- a/src/scripts/plot_particles_left.py +++ b/src/scripts/plot_particles_left.py @@ -23,12 +23,16 @@ def main(): files = [ "output/time_dependent_potential/wide_sweep.txt", "output/time_dependent_potential/narrow_sweep.txt", + "output/time_dependent_potential/narrow_sweep_fine.txt", "output/time_dependent_potential/narrow_sweep_interactions.txt", + "output/time_dependent_potential/narrow_sweep_interactions_fine.txt", ] outputs = [ "../latex/images/particles_left_wide_sweep.pdf", "../latex/images/particles_left_narrow_sweep.pdf", + "../latex/images/particles_left_narrow_sweep_fine.pdf", "../latex/images/particles_left_narrow_sweep_interactions.pdf", + "../latex/images/particles_left_narrow_sweep_interactions_fine.pdf", ] for file, output in zip(files, outputs): with open(file) as f: