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>