hi all, i just wonder the code , mglin, written in this chapter, if it works in this way: input the rhs of the PDE and it will return the solution in one step. i.e. if mglin runs only one time? or do i need to add a convergent for-loop until iu converges and output iu?