Block-smoothers with variable block size

[Nicolas Tardieu]

Dear PETSc experts,

I am currently designing a multigrid PC that needs to operate on variable
block size matrices. I know PETSc does not support this feature.
So, I will do it "by hand".  Neverthless I would like to be able to use some
smoothers like Jacobi or SOR in their block version.

What is the best solution that would work both in sequential and parallel?

- define some MATSHELL that handle variable block size
    -> will I be able to use in some way a block Jacobi or SOR?
    -> would it be possible to use CSR inside the MATSHELL and manage the
variable blocks.

- use a AIJ matrix and manage the blocks using some indirection mechanism
and reprogram Jacobi or SOR



Thank you,

Nicolas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080320/85514329/attachment.htm>