xfygogo/phasicFlow
1
0
Fork 0
mirror of https://github.com/PhasicFlow/phasicFlow.git synced 2025-07-08 03:07:03 +00:00
Code Issues Packages Projects Releases Wiki Activity

some modification for write rectMesh results in vtk format

Browse Source
This commit is contained in:
wanqing0421
2025-06-26 13:21:28 +08:00
parent a18936c8ec
commit be086ffb67
2 changed files with 41 additions and 25 deletions
Download Patch File Download Diff File Expand all files Collapse all files

63
src/PostprocessData/region/regionPoints/rectMeshRegionPoints/rectMeshRegionPoints.cpp
View File

@ -2,10 +2,8 @@
#include "fieldsDataBase.hpp" #include "fieldsDataBase.hpp"
#include "numericConstants.hpp" #include "numericConstants.hpp"
namespace pFlow::postprocessData
{
rectMeshRegionPoints::rectMeshRegionPoints pFlow::postprocessData::rectMeshRegionPoints::rectMeshRegionPoints
( (
const dictionary &dict, const dictionary &dict,
fieldsDataBase &fieldsDataBase fieldsDataBase &fieldsDataBase
@ -13,16 +11,13 @@ rectMeshRegionPoints::rectMeshRegionPoints
: :
regionPoints(dict, fieldsDataBase) regionPoints(dict, fieldsDataBase)
{ {
const auto& rectMeshInfo = dict.subDict("rectMeshInfo"); const auto& rectMeshInfo = dict.subDict("rectMeshInfo");
boxRegion_ = box(rectMeshInfo.subDict("boxInfo"));
auto minP = rectMeshInfo.getVal<realx3>("min"); nx = rectMeshInfo.getValMax<uint32>("nx", 1);
auto maxP = rectMeshInfo.getVal<realx3>("max"); ny = rectMeshInfo.getValMax<uint32>("ny", 1);
auto nx = rectMeshInfo.getValMax<uint32>("nx", 1); nz = rectMeshInfo.getValMax<uint32>("nz", 1);
auto ny = rectMeshInfo.getValMax<uint32>("ny", 1);
auto nz = rectMeshInfo.getValMax<uint32>("nz", 1);
boxRegion_ = box(minP, maxP); // boxRegion_ = box(minP, maxP);
cells_ = uint32x3(nx, ny, nz); cells_ = uint32x3(nx, ny, nz);
uint32 nCells = nx * ny * nz; uint32 nCells = nx * ny * nz;
@ -52,7 +47,7 @@ rectMeshRegionPoints::rectMeshRegionPoints
} }
} }
bool rectMeshRegionPoints::update() bool pFlow::postprocessData::rectMeshRegionPoints::update()
{ {
const auto points = database().updatePoints(); const auto points = database().updatePoints();
for (auto& elem : selectedPoints_) for (auto& elem : selectedPoints_)
@ -78,21 +73,39 @@ bool rectMeshRegionPoints::update()
return true; return true;
} }
bool rectMeshRegionPoints::write(iOstream &os) const bool pFlow::postprocessData::rectMeshRegionPoints::write(iOstream &os) const
{ {
os << "# rect mesh region points\n"; os << "# vtk DataFile Version 3.0" << endl;
os <<"# min point: "<< boxRegion_.minPoint() << endl; os << "postProcessData" << endl;
os <<"# max point: "<< boxRegion_.maxPoint() << endl; os << "ASCII" << endl;
os <<"# nx: "<< cells_.x() <<endl; os << "DATASET RECTILINEAR_GRID" << endl;
os <<"# ny: "<< cells_.y() <<endl; os << "DIMENSIONS " << nx + 1 << " " << ny + 1 << " " << nz + 1 << endl;
os <<"# nz: "<< cells_.z() <<endl;
os << "time/No. "; auto [x, y , z] = boxRegion_.minPoint();
for (uint32 i = 0; i < cells_.x() * cells_.y() * cells_.z(); ++i) auto [dx, dy, dz] = (boxRegion_.maxPoint() - boxRegion_.minPoint()) / realx3(nx, ny, nz);
os << "X_COORDINATES " << nx + 1 << " float\n";
for(int32 i = 0; i < nx + 1; i++)
{ {
os << i << " "; os << x << "\n";
x += dx;
} }
os << endl;
os << "Y_COORDINATES " << ny + 1 << " float\n";
for(int32 j = 0; j < ny + 1; j++)
{
os << y << "\n";
y += dy;
}
os << "Z_COORDINATES " << nz + 1 << " float\n";
for(int32 j = 0; j < nz + 1; j++)
{
os << z << "\n";
z += dz;
}
os << "CELL_DATA " << nx * ny * nz << endl;
return true; return true;
} }
} // End namespace pFlow::postprocessData

3
src/PostprocessData/region/regionPoints/rectMeshRegionPoints/rectMeshRegionPoints.hpp
View File

@ -51,6 +51,9 @@ private:
/// box object defining the region for point selection /// box object defining the region for point selection
box boxRegion_; box boxRegion_;
/// Number of cells in each direction
uint32 nx, ny, nz;
/// store the cells that are inside the box region /// store the cells that are inside the box region
uint32x3 cells_; uint32x3 cells_;