FYS3150/Project-3
3
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Coryab/code #11

Merged
coryab merged 5 commits from coryab/code into develop 2023-10-24 18:35:52 +00:00
Conversation 0 Commits 5 Files Changed 14 +7 -2
7 changed files with 7 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files
Showing only changes of commit 9bc6038af0 - Show all commits

BIN
latex/images/particles_left_narrow_sweep.pdf
View File

Binary file not shown.

BIN
latex/images/particles_left_narrow_sweep_fine.pdf Normal file
View File

Binary file not shown.

BIN
latex/images/particles_left_narrow_sweep_interactions.pdf
View File

Binary file not shown.

BIN
latex/images/particles_left_narrow_sweep_interactions_fine.pdf Normal file
View File

Binary file not shown.

BIN
latex/images/particles_left_wide_sweep.pdf
View File

Binary file not shown.

5
src/main.cpp
View File

@ -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 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_p = (25. + w_n * 20.) / (w_n - w_p);
double A_n = -(25. + w_p * 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<double> f = std::complex<double> f =
A_p * std::exp(std::complex<double>(0., -w_p * t)) A_p * std::exp(std::complex<double>(0., -w_p * t))
+ A_n * std::exp(std::complex<double>(0., -w_n * t)); + A_n * std::exp(std::complex<double>(0., -w_n * t));
@ -211,7 +212,7 @@ void potential_resonance_narrow_sweep()
double amplitudes[]{.1, .4, .7}; double amplitudes[]{.1, .4, .7};
double freq_start = 1.; double freq_start = 1.1;
double freq_end = 1.7; double freq_end = 1.7;
double freq_increment = .002; double freq_increment = .002;
size_t freq_iterations = size_t freq_iterations =
@ -268,7 +269,7 @@ void potential_resonance_narrow_sweep_interaction()
double amplitudes[]{.1, .4, .7}; double amplitudes[]{.1, .4, .7};
double freq_start = 1.; double freq_start = 1.1;
double freq_end = 1.7; double freq_end = 1.7;
double freq_increment = .002; double freq_increment = .002;
size_t freq_iterations = size_t freq_iterations =

4
src/scripts/plot_particles_left.py
View File

@ -23,12 +23,16 @@ def main():
files = [ files = [
"output/time_dependent_potential/wide_sweep.txt", "output/time_dependent_potential/wide_sweep.txt",
"output/time_dependent_potential/narrow_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.txt",
"output/time_dependent_potential/narrow_sweep_interactions_fine.txt",
] ]
outputs = [ outputs = [
"../latex/images/particles_left_wide_sweep.pdf", "../latex/images/particles_left_wide_sweep.pdf",
"../latex/images/particles_left_narrow_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.pdf",
"../latex/images/particles_left_narrow_sweep_interactions_fine.pdf",
] ]
for file, output in zip(files, outputs): for file, output in zip(files, outputs):
with open(file) as f: with open(file) as f: