Vacuum Boundary Conditions in Dielectric Planar should be part of 3d boundary condition code.
From my understanding, the vacuum boundary condition for dielectric planar should be accessible in the is_3d
branch of the conclude
function not in the non-three-dimensional (= 2d) branch.
if not self.is_3d:
self.results.eps_perp = -pref * cov_perp
self.results.eps_perp_self = -pref * cov_perp_self
self.results.eps_perp_coll = -pref * cov_perp_coll
self.results.deps_perp = pref * dcov_perp
if self.vac:
self.results.eps_perp *= 2.0 / 3.0
self.results.eps_perp_self *= 2.0 / 3.0
self.results.eps_perp_coll *= 2.0 / 3.0
self.results.deps_perp *= 2.0 / 3.0
else:
self.results.eps_perp = -cov_perp / (pref**-1 + var_perp / self.results.V)
self.results.deps_perp = pref * dcov_perp
self.results.eps_perp_self = (-pref * cov_perp_self) / (
1 + pref / self.results.V * var_perp
)
self.results.eps_perp_coll = (-pref * cov_perp_coll) / (
1 + pref / self.results.V * var_perp
)