I find both sor and mglin depend on zero boundary values to succeed. This make them useless in practice. Text in chapter 19 is also vague on boundary conditions and should have offer more discussions on how to apply Dirichlet/Neumann/mixed boundary conditions to the routine. Is there anyone there managed to add boundary conditions to sor/mglin?