An add-on question. Member function *mg::SmootherRelaxation::smooth* needs
to call *RelaxationType::step* to actually do iterations, which in my case
is *PreconditionJacobi*. And *PreconditionJacobi::step* is an interface
calling *Jacobi_step* from its MatrixType inside.
However this is not
Hi all,
As described in the title, I am looking for an appropriate preconditioner
for an asymmetric matrix implemented in matrix-free style.
I've gone though all the matrix-free tutorials and none of them talks about
the asymmetric case. (Almost all of them use the Chebyshev smoother with