Skip to content
Snippets Groups Projects

Fix sign of dipolar vector potential for phase

Merged Franco Bonafé requested to merge fix_phase_vecpot_sign_mxll_coupling into release-14.0
All threads resolved!
2 files
+ 9
11
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -323,7 +323,8 @@ contains
case (VELOCITY_GAUGE_DIPOLE)
call hm_base%allocate_field(mesh, FIELD_UNIFORM_VECTOR_POTENTIAL, .false.)
hm_base%uniform_vector_potential(1:space%dim) = hm_base%uniform_vector_potential(1:space%dim) + &
! The minus sign is here is for the wrong convention of Octopus for the phase
+2
hm_base%uniform_vector_potential(1:space%dim) = hm_base%uniform_vector_potential(1:space%dim) - &
this%vec_pot_dip(1:space%dim)
case (FULL_MINIMAL_COUPLING)
Loading