This is output from the run_it_simple_PT.jl file (see 7295f64, note that it is in a little extra branch).
After 10^4 iterations, the residuals are low everywhere except at the upper x-boundary, so we were wondering whether it has something to do with the boundary conditions.
After 10^5 iterations, it already looks a bit different:
One idea was to set Dirichlet B.C. also at the upper boundaries. Then the residuals after 10^5 iterations look like this:
In any case the model is far away from bringing the residuals to machine precision.
This is output from the
run_it_simple_PT.jlfile (see 7295f64, note that it is in a little extra branch).After 10^4 iterations, the residuals are low everywhere except at the upper x-boundary, so we were wondering whether it has something to do with the boundary conditions.
After 10^5 iterations, it already looks a bit different:
One idea was to set Dirichlet B.C. also at the upper boundaries. Then the residuals after 10^5 iterations look like this:
In any case the model is far away from bringing the residuals to machine precision.