Hi Yan Zhang,
Thanks for the toolbox. i have a question. If i want to apply your code to flux maximizing flow, where should i change? i think what i have to do is to change u_t = div(\nabla u) to u_t = div(\nabla I) or u_t = div(V), where I is some input image and V is vector field.