trimesh.cpp

Demonstration of Creation an Artificial Triangular Mesh.

Demonstration of Creation an Artificial Triangular Mesh

#include <iostream>
#include <bundle/Bundle.hpp>
#include <grid/CartesianChart.hpp>
 
using namespace Fiber;
using namespace std;
using namespace MemCore;
using namespace Eagle;
 
int main()
{
BundlePtr B; 
 
double  time = 0.0; 
 
Slice   &T0 = (*B)[time]; 
 
Grid    &G = T0[ "MySpecialTriangularMesh" ]; 
 
Skeleton&MeshVertices     = G[ SkeletonID(3,0) ]; 
Skeleton&MeshConnectivity = G[ SkeletonID(3,1) ]; 
 
RefPtr<Chart> myCartesianChart = G.makeChart( typeid(Fiber::CartesianChart3D) ); 
 
Representation&CartesianVertices = MeshVertices[ myCartesianChart ]; 
Representation&CellsAsVertices   = MeshConnectivity[ MeshVertices ]; 
 
RefPtr<MemArray<1, Vector<double,3> > >    Vertices = new MemArray<1, Vector<double,3> >(5); 
RefPtr<MemArray<1, Vector<int,3> > >   Connectivity = new MemArray<1, Vector<int,3> >(3); 
 
// Define the vertices of the mesh
        (*Vertices)[0] = 1.0, 3.4, 5.1;
        (*Vertices)[1] = 5.0, 1.3, 2.2;
        (*Vertices)[2] = 3.0, 2.8, 4.3;
        (*Vertices)[3] = 1.4, 8.9, 7.6;
        (*Vertices)[4] = 5.3, 2.1, 6.12; 
 
// Define the triangles of the mesh
        (*Connectivity)[0] = 0,1,2; 
        (*Connectivity)[1] = 3,2,4;
        (*Connectivity)[2] = 0,3,2; 
 
        CartesianVertices [ FIBER_POSITIONS ] = new Field(Vertices, Cache::MemCache() ); 
        CellsAsVertices.setPositions( new Field(Connectivity, Cache::MemCache() ) );
 
        B.save("MyMesh.f5");
 
        return 0;
}