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
