Hi, when solving my coupled nonlinear PDE systems I need to extract the boundary of the domain, compute its curvature and the force there, and move mesh. I wanted to know whether I can write a subroutine for these treatments and call this subroutine. I did not find a good example in the tutorial py files. Thanks!
BTW, Maybe I shall say "define/call functions" but function has its specific meaning in Fenics so I use subroutines instead.
