int main (int argc, char *argv[])
{
int i, n, ib, nb, nz, nv, celldim, phydim;
int nn, type, *elems = 0, idata[5];
cgsize_t ne;
char *p, basename[33], title[65];
float value, *var;
SOLUTION *sol;
FILE *fp;
if (argc < 2)
print_usage (usgmsg, NULL);
ib = 0;
basename[0] = 0;
while ((n = getargs (argc, argv, options)) > 0) {
switch (n) {
case 'a':
ascii = 1;
break;
case 'b':
ib = atoi (argarg);
break;
case 'B':
strncpy (basename, argarg, 32);
basename[32] = 0;
break;
case 'w':
weighting = 1;
break;
case 'S':
usesol = atoi (argarg);
break;
}
}
if (argind > argc - 2)
print_usage (usgmsg, "CGNSfile and/or Tecplotfile not given");
if (!file_exists (argv[argind]))
FATAL (NULL, "CGNSfile does not exist or is not a file");
/* open CGNS file */
printf ("reading CGNS file from %s\n", argv[argind]);
nb = open_cgns (argv[argind], 1);
if (!nb)
FATAL (NULL, "no bases found in CGNS file");
if (*basename && 0 == (ib = find_base (basename)))
FATAL (NULL, "specified base not found");
if (ib > nb) FATAL (NULL, "base index out of range");
cgnsbase = ib ? ib : 1;
if (cg_base_read (cgnsfn, cgnsbase, basename, &celldim, &phydim))
FATAL (NULL, NULL);
if (celldim != 3 || phydim != 3)
FATAL (NULL, "cell and physical dimension must be 3");
printf (" using base %d - %s\n", cgnsbase, basename);
if (NULL == (p = strrchr (argv[argind], '/')) &&
NULL == (p = strrchr (argv[argind], '\\')))
strncpy (title, argv[argind], sizeof(title));
else
strncpy (title, ++p, sizeof(title));
title[sizeof(title)-1] = 0;
if ((p = strrchr (title, '.')) != NULL)
*p = 0;
read_zones ();
if (!nZones)
FATAL (NULL, "no zones in the CGNS file");
/* verify dimensions fit in an integer */
for (nz = 0; nz < nZones; nz++) {
if (Zones[nz].nverts > CG_MAX_INT32)
FATAL(NULL, "zone size too large to write with integers");
if (Zones[nz].type == CGNS_ENUMV(Unstructured)) {
count_elements (nz, &ne, &type);
if (ne > CG_MAX_INT32)
FATAL(NULL, "too many elements to write with integers");
}
}
nv = 3 + check_solution ();
/* open Tecplot file */
printf ("writing %s Tecplot data to <%s>\n",
ascii ? "ASCII" : "binary", argv[++argind]);
if (NULL == (fp = fopen (argv[argind], ascii ? "w+" : "w+b")))
FATAL (NULL, "couldn't open Tecplot output file");
/* write file header */
if (ascii)
fprintf (fp, "TITLE = \"%s\"\n", title);
else {
fwrite ("#!TDV75 ", 1, 8, fp);
i = 1;
write_ints (fp, 1, &i);
write_string (fp, title);
}
/* write variables */
if (ascii) {
fprintf (fp, "VARIABLES = \"X\", \"Y\", \"Z\"");
if (usesol) {
sol = Zones->sols;
for (n = 0; n < sol->nflds; n++)
fprintf (fp, ",\n\"%s\"", sol->flds[n].name);
}
}
else {
write_ints (fp, 1, &nv);
write_string (fp, "X");
write_string (fp, "Y");
write_string (fp, "Z");
if (usesol) {
sol = Zones->sols;
for (n = 0; n < sol->nflds; n++)
write_string (fp, sol->flds[n].name);
}
}
/* write zones */
if (!ascii) {
for (nz = 0; nz < nZones; nz++) {
if (Zones[nz].type == CGNS_ENUMV(Structured)) {
idata[0] = 0; /* BLOCK */
idata[1] = -1; /* color not specified */
idata[2] = (int)Zones[nz].dim[0];
idata[3] = (int)Zones[nz].dim[1];
idata[4] = (int)Zones[nz].dim[2];
}
else {
count_elements (nz, &ne, &type);
idata[0] = 2; /* FEBLOCK */
idata[1] = -1; /* color not specified */
idata[2] = (int)Zones[nz].dim[0];
idata[3] = (int)ne;
idata[4] = type;
}
value = 299.0;
write_floats (fp, 1, &value);
write_string (fp, Zones[nz].name);
write_ints (fp, 5, idata);
}
value = 357.0;
write_floats (fp, 1, &value);
}
for (nz = 0; nz < nZones; nz++) {
printf (" zone %d...", nz+1);
fflush (stdout);
read_zone_grid (nz+1);
ne = 0;
type = 2;
nn = (int)Zones[nz].nverts;
var = (float *) malloc (nn * sizeof(float));
if (NULL == var)
FATAL (NULL, "malloc failed for temp float array");
if (Zones[nz].type == CGNS_ENUMV(Unstructured))
elems = volume_elements (nz, &ne, &type);
if (ascii) {
if (Zones[nz].type == CGNS_ENUMV(Structured))
fprintf (fp, "\nZONE T=\"%s\", I=%d, J=%d, K=%d, F=BLOCK\n",
Zones[nz].name, (int)Zones[nz].dim[0],
(int)Zones[nz].dim[1], (int)Zones[nz].dim[2]);
else
fprintf (fp, "\nZONE T=\"%s\", N=%d, E=%d, F=FEBLOCK, ET=%s\n",
Zones[nz].name, nn, (int)ne, type == 2 ? "TETRAHEDRON" : "BRICK");
}
else {
value = 299.0;
write_floats (fp, 1, &value);
i = 0;
write_ints (fp, 1, &i);
i = 1;
for (n = 0; n < nv; n++)
write_ints (fp, 1, &i);
}
for (n = 0; n < nn; n++)
var[n] = (float)Zones[nz].verts[n].x;
write_floats (fp, nn, var);
for (n = 0; n < nn; n++)
var[n] = (float)Zones[nz].verts[n].y;
write_floats (fp, nn, var);
for (n = 0; n < nn; n++)
var[n] = (float)Zones[nz].verts[n].z;
write_floats (fp, nn, var);
if (usesol) {
read_solution_field (nz+1, usesol, 0);
sol = &Zones[nz].sols[usesol-1];
if (sol->location != CGNS_ENUMV(Vertex))
cell_vertex_solution (nz+1, usesol, weighting);
for (nv = 0; nv < sol->nflds; nv++) {
for (n = 0; n < nn; n++)
var[n] = (float)sol->flds[nv].data[n];
write_floats (fp, nn, var);
}
}
free (var);
if (Zones[nz].type == CGNS_ENUMV(Unstructured)) {
if (!ascii) {
i = 0;
write_ints (fp, 1, &i);
}
nn = 1 << type;
for (i = 0, n = 0; n < ne; n++, i += nn)
write_ints (fp, nn, &elems[i]);
free (elems);
}
puts ("done");
}
fclose (fp);
cg_close (cgnsfn);
return 0;
}
void cCellMesh::get_mesh(std::string file_name){
// local variables
std::ifstream cell_file(file_name.c_str()); // open the mesh file
std::string line; // file line buffer
std::vector <std::string> tokens; // tokenized line
// get the mesh nodes
std::cout << "<MESH> getting the mesh nodes..." << std::endl;
while(getline(cell_file, line)){
if(line != "$Nodes") continue;
getline(cell_file, line);
nodes_count = std::atof(line.c_str());
coordinates.resize(nodes_count, Eigen::NoChange);
for(tElement n = 0; n < nodes_count; n++){
getline(cell_file, line);
boost::split(tokens, line, boost::is_any_of(", "), boost::token_compress_on);
for(int m = 0; m < 3; m++) coordinates(n, m) = atof(tokens[m + 1].c_str());
}
break;
}
surface_node.resize(nodes_count, Eigen::NoChange);
surface_node.setZero();
// get the mesh elements
std::cout << "<MESH> getting the mesh elements..." << std::endl;
while(getline(cell_file, line)){
if(line != "$Elements") continue;
getline(cell_file, line);
total_elements_count = std::atof(line.c_str());
surface_elements.resize(total_elements_count, Eigen::NoChange); // overkill for now
volume_elements.resize(total_elements_count, Eigen::NoChange); //
for(tElement n = 0; n < total_elements_count; n++){
getline(cell_file, line);
boost::split(tokens, line, boost::is_any_of(", "), boost::token_compress_on);
int element_type = atoi(tokens[1].c_str());
if(element_type == 2){ // surface triangles
for(int m = 0; m < 3; m++){ // index from zero
tElement i = atol(tokens[m + 5].c_str()) - 1;
surface_elements(surface_elements_count, m) = i;
surface_node(i) = true;
}
surface_elements_count++;
}
if(element_type == 4){ // volume tetrahedrons
for(int m = 0; m < 4; m++) // index from zero
volume_elements(volume_elements_count, m) = atol(tokens[m + 5].c_str()) - 1;
volume_elements_count++;
}
}
surface_elements.conservativeResize(surface_elements_count, Eigen::NoChange); // correct the size
volume_elements.conservativeResize(volume_elements_count, Eigen::NoChange); //
break;
}
// get the node data
std::cout << "<MESH> getting the mesh node data..." << std::endl;
while(getline(cell_file, line)){
if(line != "\"distance to nearest lumen\"") continue;
else for(int i = 0; i < 6; i++) getline(cell_file, line); // skip six lines
node_data.resize(nodes_count, Eigen::NoChange);
for(tElement n = 0; n < nodes_count; n++){
getline(cell_file, line);
boost::split(tokens, line, boost::is_any_of(", "), boost::token_compress_on);
node_data(n, dist_lumen) = atof(tokens[1].c_str());
}
break;
}
cell_file.close();
}
