Ellipsoid Meshing Issue in C++

#include <iostream> #include <fstream> namespace nglib { #include <nglib.h> } int main(int argc, char** argv) { using namespace nglib; char file_name[256]; strcpy(file_name, argv[1]); std::cout << file_name << std::endl; // Initialise the Netgen Core library nglib::Ng_Init(); long nvertices = 0; long ntriangles = 0; long ntetrahedrons = 0; double* vertexlist = NULL; int* trianglelist = NULL; int* tetralist = NULL; nglib::Ng_CSG_GenerateMesh(file_name, &nvertices, &ntriangles, &ntetrahedrons, vertexlist, trianglelist, tetralist); // deconstruct Netgen library. nglib::Ng_Exit(); return 0; } DLL_HEADER Ng_Result Ng_CSG_GenerateMesh (const char * filename, long *nvertices, long *ntriangles, long *ntetrahedrons, double * vertexlist, int * trianglelist, int * tetralist);

# ## An elliptic cylinder # algebraic3d solid test = ellipsoid (0, 0, 0; 2, 0, 0; 0, 1, 0; 0, 0, 1) -maxh=0.5; tlo test;

define solid test Calc Triangle Approximation Object 0 has 882 triangles Start Findpoints main-solids: 1 Found points 0 Analyze spec points done find points done Find edges Add line segment to smooth surface 0 edges found Check intersecting edges CalcLocalH: 2 Points 0 Elements 0 Surface Elements Start Findpoints main-solids: 1 Found points 0 Analyze spec points done Find edges Add line segment to smooth surface 0 edges found Check intersecting edges Start Findpoints main-solids: 1 Found points 0 Analyze spec points done Find edges Add line segment to smooth surface 0 edges found Check intersecting edges Surface 1 / 1 load internal triangle rules [b]faces = 66 trials = 1000 elements = 65 els/sec = 650000 faces = 113 trials = 2000 elements = 111 els/sec = 1.11e+06 faces = 141 trials = 3000 elements = 139 els/sec = 1.39e+06 faces = 174 trials = 4000 elements = 172 els/sec = 1.72e+06[/b] ........