Dear PETSc experts,<br><br>I am currently designing a multigrid PC that needs to operate on variable block size matrices. I know PETSc does not support this feature.<br>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.<br>
<br>What is the best solution that would work both in sequential and parallel?<br><br>- define some MATSHELL that handle variable block size<br> -> will I be able to use in some way a block Jacobi or SOR?<br> -> would it be possible to use CSR inside the MATSHELL and manage the variable blocks.<br>
<br>- use a AIJ matrix and manage the blocks using some indirection mechanism and reprogram Jacobi or SOR<br><br><br><br>Thank you,<br><br>Nicolas<br><br>