I am not sure if this is the right place for the post.
Would it be useful to have a relaxation parameter in compute_fp? That is, v(t+1) = (1-eta)v(t) + eta Tv(t), where T is the operator whose fp we are after? Not so much in Bellman iteration, but say, to find fixed points of Aiyagari models, a relaxation parameter is useful.
Seems easy to do.