> What you imply that this is a correct usage for block matrix if I use
> gmres+any preconditioner (in this case, block jacobi) :
>
> PETScWrappers::PreconditionBlockJacobi
> preconditioner(precon_matrix.block(0,0));
> gmres.solve (system_matrix.block(0,0), solution.block(0),
> system_rhs.block(0) , preconditioner);
>
>
> NOT :
>
> PETScWrappers::PreconditionBlockJacobi preconditioner(precon_matrix);
> gmres.solve (system_matrix, solution, system_rhs , preconditioner);
Yes.
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [EMAIL PROTECTED]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
