On Dec 12, 2007, at 12:40 PM, [EMAIL PROTECTED] wrote:
. We have run into a problem in solving the
Poisson equation for a potential field with a varying charge
density. I thought
that someone might be able to explain why this script is
disregarding the charge
density when it solves for the field.
What are you seeing that leads you to that conclusion? Given your
boundary conditions, unless the solution you are seeing is POT=0
everywhere, then the charge density *is* having an effect. When I run
your code, POT is definitely not zero everywhere and has undulations
that seem consistent with your imposed sinusoidal charged density.
What version of FiPy are you running?
What platform?
What are you seeing?
What did you expect to see?
Have you plotted C_li or Rho to see if the charge distribution is
what you think it is?
