convert_geom_file.cc File Reference

#include <iostream>
#include <fstream>
#include <vector>
#include <math.h>

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

 {
 
 // Convert argument to strings that specify the input file name
 string mesh_file_name(argv[1]);
  
// Read the output mesh file to find informations about the nodes
// and elements of the mesh 
 
ifstream infile(mesh_file_name.c_str(), ios_base::in);
unsigned n_node;
infile>>n_node;
vector<double> x(n_node);
vector<double> y(n_node);
vector<int> vertinfo(n_node);
 for(unsigned i=0;i<n_node;i++)
  { 
   infile>>x[i];
   infile>>y[i];
   infile>>vertinfo[i];
  }
unsigned n_vx;
infile>>n_vx;
vector<int> nodecode(n_vx);
vector<int> icurv(n_vx);
vector<double> ucurv(n_vx);
 for(unsigned i=0;i<n_vx;i++)
  {
   infile>>nodecode[i];
   infile>>icurv[i];
   infile>>ucurv[i];
  }
unsigned n_local_node;
infile>>n_local_node;
unsigned n_element;
infile>>n_element;
unsigned b=n_local_node*n_element;
vector<int> v(b);
vector<int> edgeinfo(b);
unsigned k=0;
for(unsigned i=0;i<n_element;i++)
 {
  for(unsigned j=0;j<n_local_node;j++)
   {
    infile>>v[k];
    k++;
   }
 }
unsigned l=0;
for(unsigned i=0;i<n_element;i++)
 {
  for(unsigned j=0;j<n_local_node;j++)
   {
    infile>>edgeinfo[l];
    l++;
   }
 }
 
infile.close();
 
// Create a file of type ".quad" to visualize the mesh with geomview
 
 
 unsigned nn=0;
 char result[20];
 sprintf(result,"%s","mesh.quad");
 ofstream outfile(result,ios_base::out);
 outfile<<"QUAD"<<'\n';
 for(unsigned i=0;i<n_element;i++)
   {
    for(unsigned j=0;j<n_local_node;j++) 
     {
      outfile<<x[v[nn]-1]<<" "<<y[v[nn]-1]<<" 0 ";
      nn++; 
     }
    outfile<<'\n';
   } 
 outfile.close(); 
 
 
 
}  //end of main

References b, i, j, k, out(), v, plotDoE::x, and y.