Implement total_force

2023-09-29 14:48:32 +02:00
parent c611fa70a2
commit eb989e0293
1 changed files with 1 additions and 1 deletions
--- a/src/PenningTrap.cpp
+++ b/src/PenningTrap.cpp
@@ -105,7 +105,7 @@ arma::vec PenningTrap::total_force_particles(int i)

 arma::vec PenningTrap::total_force(int i)
 {
-
+    return this->total_force_external(i) - this->total_force_particles(i);
 }

 void PenningTrap::evolve_RK4(double dt)