The weighted_jacobi function seems to do the job of being an iterative solver.
I know that if omega=1 is effectively removes the weighting but I think it would be nice to have a method called jacobi that performs the unweighted iteration.
Make sure all new functions include a docstring explaining how to use them.
The
weighted_jacobifunction seems to do the job of being an iterative solver.I know that if
omega=1is effectively removes the weighting but I think it would be nice to have a method calledjacobithat performs the unweighted iteration.Make sure all new functions include a docstring explaining how to use them.