[sharonyue]

Hi All,

In mathematics, I think laplacian(A,p) and div(A*fvc::grad(p)) are the same,

so can I turn this code into:

fvScalarMatrix pEqn
            (
                fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
            );

fvScalarMatrix pEqn
            (
                fvm::div(A*fvc::grad(p)) == fvc::div(phiHbyA)
            );

Thanks for your attention.