int main(int argc, char* argv[]) {
if (argc != 2) {
std::cerr << "Usage: " << argv[0] << " [run_config.json]";
return -1;
}
Config c;
c.LoadFromJson(path(argv[1]));
ProcessMergesDump(c);
ProcessNewWordsDump(c);
ProcessVertices(c);
ProcessEdges(c);
ProcessQueue(c);
ProcessPairsClasses(c);
}
RawModel OBJLoader::LoadObjModel(const std::string& fileName, Loader& loader)
{
clock_t startTime = clock();
// Open the file as read only
FILE* file;
if (fopen_s(&file, ("../res/models/" + fileName + ".obj").c_str(), "r") != 0)
{
printf("Failed to open: %s\n", fileName);
}
// Storage variables
std::vector<glm::vec2> textures, tempTextures;
std::vector<glm::vec3> vertices, normals, tempNormals;
std::vector<int> indices;
char *type, *token, *stop = 0;
double x, y, z;
char line[256];
while (fgets(line, 256, file) != NULL)
{
token = NULL;
type = strtok_s(line, " ", &token);
// V is vertex points
if (type[0] == 'v' && type[1] == NULL)
{
x = strtod(token, &stop);
token = stop + 1; // Move to the next value
y = strtod(token, &stop);
token = stop + 1; // Move to the next value
z = strtod(token, &stop);
// Store a new vertex
vertices.push_back(glm::vec3(x, y, z));
}
// VT is vertex texture coordinates
else if (type[0] == 'v' && type[1] == 't')
{
x = strtod(token, &stop);
token = stop + 1; // Move to the next value
y = 1 - strtod(token, &stop);
// Store a new texture
tempTextures.push_back(glm::vec2(x, y));
}
else if (type[0] == 'v' && type[1] == 'n')
{
x = strtod(token, &stop);
token = stop + 1; // Move to the next value
y = strtod(token, &stop);
token = stop + 1; // Move to the next value
z = strtod(token, &stop);
// Store a new normal
tempNormals.push_back(glm::vec3(x, y, z));
}
// F is the index list for faces
else if (type[0] == 'f')
{
if (indices.size() == 0)
{
// Set the size of the array
textures.resize(vertices.size());
normals.resize(vertices.size());
}
// Process set of vertex data
ProcessVertices(token, indices, tempTextures, textures, tempNormals, normals);
}
}
fclose(file);
printf("Load time: %dms\n", clock() - startTime);
return loader.LoadToVAO(vertices, textures, normals, indices);
}
/**
* Gmsh file contains a list of nodes and their coordinates, along with
* a list of elements and those nodes which define them. We read in and
* store the list of nodes in #m_node and store the list of elements in
* #m_element. Each new element is supplied with a list of entries from
* #m_node which defines the element. Finally some mesh statistics are
* printed.
*
* @param pFilename Filename of Gmsh file to read.
*/
void InputVtk::Process()
{
if (m_mesh->m_verbose)
{
cout << "InputVtk: Start reading file..." << endl;
}
vtkPolyDataReader *vtkMeshReader = vtkPolyDataReader::New();
vtkMeshReader->SetFileName(m_config["infile"].as<string>().c_str());
vtkMeshReader->Update();
vtkPolyData *vtkMesh = vtkMeshReader->GetOutput();
vtkPoints *vtkPoints = vtkMesh->GetPoints();
const int numCellTypes = 3;
vtkCellArray* vtkCells[numCellTypes];
LibUtilities::ShapeType vtkCellTypes[numCellTypes];
int vtkNumPoints[numCellTypes];
vtkCells[0] = vtkMesh->GetPolys();
vtkCells[1] = vtkMesh->GetStrips();
vtkCells[2] = vtkMesh->GetLines();
vtkCellTypes[0] = LibUtilities::eTriangle;
vtkCellTypes[1] = LibUtilities::eTriangle;
vtkCellTypes[2] = LibUtilities::eSegment;
vtkNumPoints[0] = 3;
vtkNumPoints[1] = 3;
vtkNumPoints[2] = 2;
vtkIdType npts;
vtkIdType *pts = 0;
double p[3];
for (int i = 0; i < vtkPoints->GetNumberOfPoints(); ++i)
{
vtkPoints->GetPoint(i, p);
if ((p[0] * p[0]) > 0.000001 && m_mesh->m_spaceDim < 1)
{
m_mesh->m_spaceDim = 1;
}
if ((p[1] * p[1]) > 0.000001 && m_mesh->m_spaceDim < 2)
{
m_mesh->m_spaceDim = 2;
}
if ((p[2] * p[2]) > 0.000001 && m_mesh->m_spaceDim < 3)
{
m_mesh->m_spaceDim = 3;
}
m_mesh->m_node.push_back(boost::shared_ptr<Node>(new Node(i, p[0], p[1], p[2])));
}
for (int c = 0; c < numCellTypes; ++c)
{
vtkCells[c]->InitTraversal();
for (int i = 0; vtkCells[c]->GetNextCell(npts, pts); ++i)
{
for (int j = 0; j < npts - vtkNumPoints[c] + 1; ++j)
{
// Create element tags
vector<int> tags;
tags.push_back(0); // composite
tags.push_back(vtkCellTypes[c]); // element type
// Read element node list
vector<NodeSharedPtr> nodeList;
for (int k = j; k < j + vtkNumPoints[c]; ++k)
{
nodeList.push_back(m_mesh->m_node[pts[k]]);
}
// Create element
ElmtConfig conf(vtkCellTypes[c],1,false,false);
ElementSharedPtr E = GetElementFactory().
CreateInstance(vtkCellTypes[c],
conf,nodeList,tags);
// Determine mesh expansion dimension
if (E->GetDim() > m_mesh->m_expDim) {
m_mesh->m_expDim = E->GetDim();
}
m_mesh->m_element[E->GetDim()].push_back(E);
}
}
}
ProcessVertices();
ProcessEdges();
ProcessFaces();
ProcessElements();
ProcessComposites();
}
