[petsc-users] Does Petsc have any feature to handle the mirror points?

[Li, Zhisong (lizs)]

 Hi, Petsc Team,

In finite difference computations, we often utilize schemes of five or more points stencil in one direction. To calculate the interior point next to the boundary, we need to assume one point out of the domain by extrapolation (such as a mirror point). As we are required to include all interior points and boundary points in the system matrix, does Petsc have any feature to handle these mirror points ? In assembling the matrix for a 3D system of multiple variables and equations, besides setting up the boundary conditions, I think this will cost lots of extra effort to specially manage the points next to the boundary, if Petsc doesn't have any easy way to do this.

Thank you.


Zhisong Li
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100810/7921ee8f/attachment.htm>