I have been looking at the example problem vp1_np.py. I wanted to tweak the initial condition a little, so I inserted this line of code just before the "problem = ..." line:
u_ = interpolate(Constant(0.0),V)
Then I called the program with the command
python vp1_np.py m l 1 10 5
The consequence was that it claimed convergence in 1 iteration (not 5 as the unchanged file yields), and the answer was substantially different from the exact solution. I observed "convergence" in 1 iteration also if I used
u_ = interpolate(Constant(10.0),V)
again to the wrong solution.
I am using version 1.3 from the repositories on Ubuntu 13.10.
What am I doing wrong?