Sang pham van <[email protected]> writes: > Yes, It does converge toward 0-flow when I refine the mesh ! > > In the code, coordinates are set by > ierr = > DMDASetUniformCoordinates(da_prop,0.0+0.5*dx,1.0-0.5*dx,0.0+0.5*dy,1.0-0.5*dy,0.,0);CHKERRQ(ierr); > Does it set coordinate for elements center here? > When I want to try non-uniform grid, should I just modify coordinates > attached with the DM, or I need more implementation?
You can modify coordinates. It looks like the code does the correct coordinate transformations. (Dave's real code runs on non-affine meshes for science reasons.)
pgpP3u59FtF6F.pgp
Description: PGP signature
