fixed density choice: instead of taking first match, the one with smallest length...
[engrid.git] / src / polydatareader.cpp
blobdf3ea822ce4caaefd7ec333667bddf2910b79baa
1 //
2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
3 // + +
4 // + This file is part of enGrid. +
5 // + +
6 // + Copyright 2008,2009 Oliver Gloth +
7 // + +
8 // + enGrid is free software: you can redistribute it and/or modify +
9 // + it under the terms of the GNU General Public License as published by +
10 // + the Free Software Foundation, either version 3 of the License, or +
11 // + (at your option) any later version. +
12 // + +
13 // + enGrid is distributed in the hope that it will be useful, +
14 // + but WITHOUT ANY WARRANTY; without even the implied warranty of +
15 // + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +
16 // + GNU General Public License for more details. +
17 // + +
18 // + You should have received a copy of the GNU General Public License +
19 // + along with enGrid. If not, see <http://www.gnu.org/licenses/>. +
20 // + +
21 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
23 #include "polydatareader.h"
24 #include "vtkEgPolyDataToUnstructuredGridFilter.h"
26 #include <vtkXMLPolyDataReader.h>
28 #include <QFileInfo>
29 #include "guimainwindow.h"
31 PolyDataReader::PolyDataReader()
33 setFormat("VTK poly data files(*.vtp)");
34 setExtension(".vtp");
37 void PolyDataReader::operate()
39 try {
40 QFileInfo file_info(GuiMainWindow::pointer()->getFilename());
41 readInputFileName(file_info.completeBaseName() + ".vtp");
42 if (isValid()) {
43 EG_VTKSP(vtkXMLPolyDataReader,vtp);
44 EG_VTKSP(vtkEgPolyDataToUnstructuredGridFilter,pd2ug);
45 vtp->SetFileName(qPrintable(getFileName()));
46 pd2ug->SetInput(vtp->GetOutput());
47 pd2ug->Update();
48 grid->DeepCopy(pd2ug->GetOutput());
49 createBasicFields(grid, grid->GetNumberOfCells(), grid->GetNumberOfPoints());
50 UpdateNodeIndex(grid);
51 UpdateCellIndex(grid);
52 EG_VTKDCC(vtkIntArray, bc, grid, "cell_code");
53 for (vtkIdType id_cell = 0; id_cell < grid->GetNumberOfPoints(); ++id_cell) {
54 bc->SetValue(id_cell,99);
57 } catch (Error err) {
58 err.display();