Merge pull request #11 from gkoolstra/dev-gerwin

Fixes another error in the equations of motion
gkoolstra · Jul 17, 2024 · 1e34bc8 · 1e34bc8
2 parents 7efa299 + 964899b
commit 1e34bc8
Show file tree

Hide file tree

Showing 3 changed files with 7,751 additions and 235 deletions.
diff --git a/examples/97_transmission_line.ipynb b/examples/97_transmission_line.ipynb
diff --git a/examples/98_electron_capacitance.ipynb b/examples/98_electron_capacitance.ipynb
diff --git a/quantum_electron/eom_solver.py b/quantum_electron/eom_solver.py
@@ -107,14 +107,14 @@ def setup_eom_coupled_lc(self, ri: ArrayLike,
                               [Cdot / D, (C1 + Cdot) / D]])
 
         K[2:num_electrons + 2, 0] = K[0, 2:num_electrons + 2] = q_e / D * \
-            ((C2 + Cdot) * self.Ex_up(xe, ye) - Cdot * self.Ex_down(xe, ye))
+            ((C2 + Cdot) * self.Ex_up(xe, ye) + Cdot * self.Ex_down(xe, ye))
         K[2:num_electrons + 2, 1] = K[1, 2:num_electrons + 2] = q_e / D * \
-            ((C1 + Cdot) * self.Ex_down(xe, ye) - Cdot * self.Ex_up(xe, ye))
+            ((C1 + Cdot) * self.Ex_down(xe, ye) + Cdot * self.Ex_up(xe, ye))
 
         K[num_electrons + 2:2 * num_electrons + 2, 0] = K[0, num_electrons + 2:2 * num_electrons +
-                                                          2] = q_e / D * ((C2 + Cdot) * self.Ey_up(xe, ye) - Cdot * self.Ey_down(xe, ye))
+                                                          2] = q_e / D * ((C2 + Cdot) * self.Ey_up(xe, ye) + Cdot * self.Ey_down(xe, ye))
         K[num_electrons + 2:2 * num_electrons + 2, 1] = K[1, num_electrons + 2:2 * num_electrons +
-                                                          2] = q_e / D * ((C1 + Cdot) * self.Ey_down(xe, ye) - Cdot * self.Ey_up(xe, ye))
+                                                          2] = q_e / D * ((C1 + Cdot) * self.Ey_down(xe, ye) + Cdot * self.Ey_up(xe, ye))
 
         kij_plus = np.zeros((num_electrons, num_electrons))
         kij_minus = np.zeros((num_electrons, num_electrons))