1. Loading a Mesh
The next step is to load a mesh.
The loadMesh
function has a _name
option set by default as the
default value of the --gmsh.filename
option that point either to a
.geo
, either to a .msh
, or a .h5
file. Meshes in general are
more detail into this section.
auto mesh=loadMesh( _mesh=new Mesh<Simplex<2>> );
1.1. Exporting the Mesh for visualisation
See this section for more details about exporting and visualizing meshes.
1.2. Implementation
An implementation reads as follows:
#include <feel/feelfilters/loadmesh.hpp>
#include <feel/feelfilters/exporter.hpp>
#include <feel/feelvf/vf.hpp>
int main( int argc, char** argv )
{
using namespace Feel;
using namespace Feel::vf;
// initialize Feel++ Environment
Environment env( _argc=argc, _argv=argv,
_about=about( _name="mymesh" ,
_author="Feel++ Consortium",
_email="feelpp-devel@feelpp.org" ) );
// tag::mesh[]
// create a mesh with GMSH using Feel++ geometry tool
auto mesh = loadMesh(_mesh=new Mesh<Simplex<2>>);
// end::mesh[]
LOG(INFO) << "mesh " << soption(_name="gmsh.filename") << " loaded";
LOG(INFO) << "volume =" << integrate( _range=elements( mesh ), _expr=cst( 1. ) ).evaluate();
LOG(INFO) << "surface = " << integrate( _range=boundaryfaces( mesh ), _expr=cst( 1. ) ).evaluate();
}
and the associated config file
** xref:{examplesdir}/03-mymesh.cfg[]