On Thu, Apr 21, 2011 at 7:07 AM, TAY wee-beng <zonexo at gmail.com> wrote:
> Hi, > > Since there is a similar topic on improving performance of CFD using PETSc > earlier, I hope to improve my current CFD solver performance too. > > I wrote it in Fortran90 with MPI. It is a Immersed Boundary Method (IBM) > Navier-Stokes Cartesian grid solver in 2D, although I hope to extend it to > 3D in the future. I solve the NS equations using fractional step which > results in 2 equations - the momentum and Poisson equations. They are then > linearized into systems of equations. > > I currently solve the momentum solver using PETSc with KSPBCGS. For the > Poisson equation, I was using hypre's BoomerAMG. But I have changed to using > the geometric multigrid solver from hypre since it's slightly faster. > > Currently, I am dividing my grid along the y direction for MPI into equal > size for each processor. I guess this is not very efficient since beyond 4 > processors, the scaling factor drops. > > I think implementing the distributed array should increase performance, is > that so? I wonder how difficult it is because most examples are in C and I > am not so used to that. I am also using staggered grid but I will most > likely changed to a collocated grid arrangement. > It will definitely improve scalability. I don't think conversion should be that hard. Matt What other suggestions do you have to improve the solver's performance using > PETSc? > > Thank you very much. > > Yours sincerely, > > TAY wee-beng > > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110421/b73a3235/attachment.htm>
