Hello PETSc users I was wondering if you've got any suggestion for the issue mentioned in the topic.
I'm solving a linear system on my unstructured grid that has its own partitioning. When I set up the linear system in PETSc I use the model partitioning whose grid points might not be evenly balanced (generally is not) over the processors. So far I use VecScatter and IS to re-order the solution vector with the application numbering. What I aim at, now, is to let PETSc to decide the linear system partitioning with PETSC_DECIDE local dimension and see if I get some gain in performances, but in this case I have to scatter between vectors that, at any rate, have different parallel layouts (one of the requirements of VecScatter is to have 2 Vecs with same layouts) Does anyone of you have ever faced a similar problem? any suggestion about how to circumvent this? Feel free to ask for further information if my explanation is not clear enough Thanks in advance Ivano
