From ad10bf19487deb5778d028e37048727b2c92e0c9 Mon Sep 17 00:00:00 2001 From: Cory Date: Fri, 29 Sep 2023 14:41:04 +0200 Subject: [PATCH] Implement total_force_external --- src/PenningTrap.cpp | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/PenningTrap.cpp b/src/PenningTrap.cpp index 538bb4c..f83576c 100644 --- a/src/PenningTrap.cpp +++ b/src/PenningTrap.cpp @@ -76,7 +76,11 @@ arma::vec PenningTrap::force_on_particle(int i, int j) arma::vec PenningTrap::total_force_external(int i) { + Particle p = this->particles.at(i); + arma::vec force = p.q + *(this->external_E_field(p.r_vec) + this->external_B_field(p.r_vec)); + return force; } arma::vec PenningTrap::total_force_particles(int i)