[petsc-users] Some papers for additive schwarz and blocked jacobi?

[Dürrwang, Jürgen]

@ Jed Brown

Yes, I tried some PETSC examples and I modified one for my stuff. It works very well on my Xeon quadcore, but my intention is to mix CPU and GPU code. I want a paralell domain decomposition using jacobi block method for runing ILU(0) on each block(number of blocks = number of CPU cores). Then I want to take the results of each blocksolution as a preconditioner for a cg solver on GPU.

At the moment I can decompose my matrix in  four jacobi block matrices. I compared my results with petsc and they are the same. But now I don’t know if I have to run my cg solver on each block or could I put the results of each blocked-ILU together and the use this as preconditioner for the non blocked matrix(my large input matrix).

When I run my PETSC program with –ksp_monitor/-ksp_view, I can only see that the matrix will be blocked and then ILU(0) will applied and after that there comes a cg solver. But there is no information about how to bring the results of each solved block together for using in the cg computing….thats my issue



Jed Brown wrote:

I don't understand what you're asking. Have you looked at some PETSc examples? You should be able to do this stuff with run-time options. If you are trying something that doesn't work, send the options you ran with and the output. Running with -ksp_view is helpful to see what method is being used.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110606/8f84aae1/attachment-0001.htm>