2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
4 // + This file is part of enGrid. +
6 // + Copyright 2008,2009 Oliver Gloth +
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. +
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. +
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/>. +
21 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
24 #include "foamobject.h"
27 #include "guimainwindow.h"
29 FoamObject::FoamObject()
34 int FoamObject::deleteBetween(int i
, QString str1
, QString str2
)
36 int i1
= m_Buffer
.indexOf(str1
, i
);
38 int i2
= m_Buffer
.indexOf(str2
, i1
);
40 m_Buffer
= m_Buffer
.remove(i1
, i2
- i1
+ str2
.size());
47 void FoamObject::stripBuffer()
49 while (deleteBetween(0, "/*", "*/")) {};
50 while (deleteBetween(0, "//", "\n")) {};
51 int i
= m_Buffer
.indexOf("FoamFile", 0);
53 deleteBetween(i
, "{", "}");
55 m_Buffer
.remove("FoamFile");
56 m_Buffer
= m_Buffer
.replace("{", " ");
57 m_Buffer
= m_Buffer
.replace("}", " ");
58 m_Buffer
= m_Buffer
.replace("(", " ");
59 m_Buffer
= m_Buffer
.replace(")", " ");
60 m_Buffer
= m_Buffer
.replace(";", " ");
61 m_Buffer
= m_Buffer
.simplified();
64 void FoamObject::readFile(QString file_name
)
66 file_name
= m_CaseDir
+ "/" + file_name
;
67 QFile
file(file_name
);
68 if (!file
.open(QIODevice::ReadOnly
)) {
69 EG_ERR_RETURN(QString("error loading file:\n") + file_name
);
72 m_Buffer
= f
.readAll();
76 void FoamObject::buildMaps()
79 readFile("constant/polyMesh/points");
81 QTextStream
f(getBuffer());
85 readFile("constant/polyMesh/neighbour");
86 QTextStream
f(getBuffer());
90 m_FirstBoundaryFace
= 0;
91 while (m_FirstBoundaryFace
< num_neigh
) {
96 ++m_FirstBoundaryFace
;
100 m_VolToSurfMap
.fill(-1, num_foam_nodes
);
101 readFile("constant/polyMesh/faces");
102 int num_surf_nodes
= 0;
106 QTextStream
f(getBuffer());
108 for (int i
= 0; i
< num_foam_faces
; ++i
) {
111 for (int j
= 0; j
< num_nodes
; ++j
) {
114 max_node
= max(node
, max_node
);
115 if (i
>= m_FirstBoundaryFace
) {
116 if (m_VolToSurfMap
[node
] == -1) {
117 m_VolToSurfMap
[node
] = num_surf_nodes
;
124 m_SurfToVolMap
.fill(-1, num_surf_nodes
);
125 for (int i
= 0; i
< m_VolToSurfMap
.size(); ++i
) {
126 if (m_VolToSurfMap
[i
] != -1) {
127 if (m_VolToSurfMap
[i
] > m_SurfToVolMap
.size()) {
130 m_SurfToVolMap
[m_VolToSurfMap
[i
]] = i
;
133 for (int i
= 0; i
< m_SurfToVolMap
.size(); ++i
) {
134 if (m_SurfToVolMap
[i
] == -1) {
140 void FoamObject::setCaseDir(QString case_dir
)
142 m_CaseDir
= case_dir
;
143 GuiMainWindow::pointer()->setXmlSection("openfoam/CaseDir",m_CaseDir
);