|
[Sponsors] |
May 31, 2008, 06:11 |
Help to solve 2 eqns simultaneously in MPI
|
#1 |
Guest
Posts: n/a
|
Hi,
I obtain 2 linear eqns (u and v velocity) from the momentum eqn in my CFD code. Instead of solving eqn 1 in parallel, and then subsequently eqn 2 in parallel, I am thinking of solving the 2 eqns at the same time, using half the number of processors on each eqn. In other words, when using 4 processors, I use 2 processors for eqn 1 and 2 processors for eqn 2. Sonmeone told to use the MPI_Group features but after looking on the net, I still can't really understand how to do it. Any experts in MPI around here? Or if someone can point me to some examples, tt'll be great. Btw, I'm using the parallel solver PETSc. Thanks alot! |
|
May 31, 2008, 10:01 |
Re: Help to solve 2 eqns simultaneously in MPI
|
#2 |
Guest
Posts: n/a
|
You can create separate PETSc objects on separate communicators, but I doubt this is something you really want to do. In particular, it only possibly makes sense if you are absolutely sure that there will never be coupling between the components, even indirectly. My recommendation would be to solve them simultaneously as a coupled system (i.e. in the same Krylov iteration), but use a split preconditioner (PCFIELDSPLIT will make this easy, but you can have more manual control with PCSHELL). A normal preconditioner may work very well too. This is far more flexible and should scale well. Also, if you move to nonlinear problems, you will generally get much faster convergence with Newton-Krylov iterations on the whole coupled system. In general, anything `special' can and should go in the preconditioner.
|
|
June 1, 2008, 21:21 |
Re: Help to solve 2 eqns simultaneously in MPI
|
#3 |
Guest
Posts: n/a
|
Thanks Jed.
It's obtained from a semi-implicit scheme so that u/v are definitely not coupled together. Anyway, I'm still not able to find a suitable mpi sample code which suits my problem. Anyone can help? Thanks |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Sgimpi | pere | OpenFOAM | 27 | September 24, 2011 08:57 |
Solve Flow or VOF simultaneously ? | Ramsey | FLUENT | 1 | February 16, 2011 14:16 |
Need help to solve N-S eqns with F.V.M | vasilis | Main CFD Forum | 4 | October 7, 2008 03:56 |
Error using LaunderGibsonRSTM on SGI ALTIX 4700 | jaswi | OpenFOAM | 2 | April 29, 2008 11:54 |
Is Testsuite on the way or not | lakeat | OpenFOAM Installation | 6 | April 28, 2008 12:12 |