these modifications compute the virtual-casing self-consistent vacuum…

[SRHudson]

… field,

and this means that the Picard iterations for computing the self-consistent external field is not required. This work is incomplete: at present, the derivative of the modified matrix w.r.t. interface geometry is not yet implemented, which means that Lfindzero = 1 must be used.