Python API to evaluate ewald Coulomb force

[JZ Most]

Thanks for reply. I also wrote a PDE to simulate the setting, it gives completely different results. I expect the middle region without fixed charged should not be empty. Also, the free ions should not stack around the fixed charges, it should form a relatively smooth decreasing curve away from the fixed charges, for the reason that the diffusion also contributes to the distribution.