Dear all
In stokes_iterative demo:
A, bb = assemble_system(a, L, bcs)
P, btmp = assemble_system(b, L, bcs)
solver.set_operators(A, P)
solver.solve(U.vector(), bb)
solver.set_operators(A, P) is a little bit unclear.
Does this mean the system P^{-1} A X=P^{-1}bb will be solved?
Thanks!
