[petsc-users] question about the BuildGradientReconstruction

[Rongliang Chen]

Dear all,

I have a question about the gradient reconstruction for the dmplexfvm.

Why the ghost cells are ignored during the gradient reconstruction in 
the function BuildGradientReconstruction_Internal?

For the tetrahedron mesh, if a cell (on the corner) whose three faces 
are on the boundary of the computational domain, then only one cell can 
be used to reconstruct the gradient (it is deficient for the least 
square). I found that, for this situation, the accuracy of the gradient 
reconstructed by the least square is a problem. Do you have any 
suggestions for deal with this situation?

Best regards,
Rongliang