From f2e9004f96c9e76e7146f648dcb17cc9af21d870 Mon Sep 17 00:00:00 2001 From: Cory Date: Fri, 29 Sep 2023 14:07:27 +0200 Subject: [PATCH] Implement external_B_field --- src/PenningTrap.cpp | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/src/PenningTrap.cpp b/src/PenningTrap.cpp index f1c62df..b445f21 100644 --- a/src/PenningTrap.cpp +++ b/src/PenningTrap.cpp @@ -48,7 +48,13 @@ arma::vec PenningTrap::external_E_field(arma::vec r) arma::vec PenningTrap::external_B_field(arma::vec r) { + double x = r(0), y = r(1); + arma::vec::fixed<3> res; + res(0) = y*this->B_0; + res(1) = -x*this->B_0; + res(2) = 0.; + return res; } arma::vec PenningTrap::force_on_particle(int i, int j)