You may take a look at ~petsc/src/snes/examples/tutorials/ex19.c which solves a Nonlinear driven cavity over structured 2d grid.
Hong On Wed, 21 Jan 2009, Etienne Ahusborde wrote: > Hello, > > I am a new user of PETSc. > > My goal is to solve linear systems coming from the Navier-Stokes equations. > Components of the velocity field > are coupled each other, so that, in 2D, the matrix is composed of four blocks > (two diagonal blocks, one for each component, > and two extra-diagonal blocks for the coupling between u and v). > Discretization is made on a staggered and structured grid > (MAC grid, one grid for each component). > > After having read the user manual, I an wondering if there is a way to take > benefit of the structured grid as the Semi-Structured-Grid System interface > of Hypre, or if I have to use the parallel AIJ Sparse matrices. We have > understood that the distributed arrays are > rather adapted to scalar equations (in comparison with coupled vectorial > equations). > > In your opinion, what is the best way to treat our problem? > > Kind regards > > Etienne Ahusborde > > > >
