Implement total_force_external

2023-09-29 14:41:04 +02:00 · 2023-09-29 14:41:04 +02:00 · ad10bf1948
commit ad10bf1948
parent 8e9f3db316
1 changed files with 4 additions and 0 deletions
--- 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)