Skip to content
Snippets Groups Projects
Commit ef0369fe authored by Henry's avatar Henry
Browse files

Test-fvc: minor enhancement

parent 63e8be78
Branches
Tags
No related merge requests found
...@@ -35,17 +35,18 @@ Description ...@@ -35,17 +35,18 @@ Description
int main(int argc, char *argv[]) int main(int argc, char *argv[])
{ {
#include "setRootCase.H"
# include "setRootCase.H" #include "createTime.H"
#include "createMesh.H"
# include "createTime.H" volScalarField fx(pow(mesh.C().component(vector::X), 1));
# include "createMesh.H"
volScalarField fx(pow(mesh.C().component(vector::X), 2));
fx.write(); fx.write();
volScalarField gradx4(fvc::grad(fx)().component(vector::X)); volScalarField gradx4(fvc::grad(fx)().component(vector::X));
gradx4.write(); gradx4.write();
Info<< fvc::snGrad(fx) - (mesh.Sf()/mesh.magSf())().component(vector::X) << endl;
//volVectorField curlC(fvc::curl(1.0*mesh.C())); //volVectorField curlC(fvc::curl(1.0*mesh.C()));
//curlC.write(); //curlC.write();
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment