Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'Poisson Solver' topicStudent Support Forum > General > "Poisson Solver"

Next Comment >Help | Reply To Topic
Author Comment/Response
Joel D
10/08/12 8:06pm

Hello All,

I am looking for some help with a Poisson Solver I am writing in Mathematica. The code is quite long with Arrays plugged in, so I posted the full details can be found at http://pastebin.com/u/JoelD

I am calculating voltages given charge densities using the central difference method derived from Poisson's Eqn. After calculating the voltage, I test the data set for convergence. I am setting convergence thresholds on the order of 10^-1000+. I have the loop set up to kick out after 10000 iterations incase something goes awry, as a fail safe. I have a loop counter in place for sanity. The program seems to run fine as long as the convergence threshold is set to 10^-100.

My question is this: No matter what I update the threshold too, ex, 10^-1000, 10^-10000, the computation stops after 633 iterations and kicks out of the loop. I would appreciate any help with this, I am completely stuck. Again, I know this description is limited, so please see the attached url http://pastebin.com/u/JoelD
for the full program.

Thank you in advance

Joel

URL: http://pastebin.com/u/JoelD,

Subject (listing for 'Poisson Solver')
Author Date Posted
Poisson Solver Joel D 10/08/12 8:06pm
Re: Poisson Solver Bill Simpson 10/09/12 12:22pm
Re: Re: Poisson Solver Joel D 10/09/12 6:11pm
Re: Re: Re: Poisson Solver Bill Simpson 10/10/12 01:41am
Re: Re: Re: Re: Poisson Solver Joel Dale 10/10/12 6:49pm
Re: Re: Re: Re: Poisson Solver Joel D 10/10/12 8:10pm
Next Comment >Help | Reply To Topic