Thank you Jed, I guess the thing that makes it complicated application of boundary conditions. Since XYZGHOSTED is not available in fortran, I am using XYZPERIODIC in order to implement bc's, as it is done in ex31.c in SNES Combining XYZPERIODIC with larger block size is still not much clear for me.
By the way, what is the current status of DMDACreate3D and other DMDA routines and their fortran support ? According to ex11f90.F and ex22f.F in the manual page of DMDACreate3d , these routine(s) provide ghost cell support for fortran in 3D ? If ghost cell support is available, implementation would be relatively easy Cheers, Ilyas 2011/4/28 Jed Brown <jed at 59a2.org> > On Thu, Apr 28, 2011 at 15:41, ilyas ilyas <ilyascfd at gmail.com> wrote: > >> Would you please explain more the first option? > > > You create a system with larger block size. If you have two node-centered > values and one cell-centered, you would use a block size of 3. Each center > would be associated with the node to its lower left, for example. There will > be a fringe of "cell-centers" that extend out of domain on the right and > top, you use x = 0 for these equations and the Jacobian for those equations > will just have a 1 on the diagonal. > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110430/42626df4/attachment.htm>
