At present the library supports reading only gmsh grids consisting of linear elements of triangles and quadrilaterals. The mesh can consist of both types of cells. It uses the older version of Gmsh grid format, so set

Mesh.MshFileVersion = 2.2;

in your geo file. Or specify the format on the command line

gmsh -2 -format msh22 foo.geo -o foo.msh

• src/main.cc : test grid preprocessor
• src/ccfv1.cc: first order FV code
• src/ccfv2.cc: second order FV code

cd src
make
cd ../tests
gmsh -2 cylinder.geo
../src/main

See the output files cylinder.vtk and test_msh.msh.

visit -o cylinder.vtk
gmsh test_msh.msh

cd tests
gmsh -2 square_tri.geo -o ccfv1.msh
../src/ccfv1
visit -o sol*.vtk

This code reads some parameters from an input file, see tests/input.txt file.

Run astyle on source code before committing it

astyle --options=./astyle.rc grid.cc

To understand the methods used in this code, you can refer to the following two books.

• Rainald Lohner: Applied Computational Fluid Dynamics Techniques
• Jiri Blazek: Computational Fluid Dynamics

• Shashwat Tiwari
• Praveen Chandrashekar