OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

add offset to Domain for swc

Browse Source
This commit is contained in:
James McClure 2022-05-25 17:59:08 -04:00
parent df27167212
commit dc78491a9c
2 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
common/Domain.cpp
View File

@ -101,14 +101,13 @@ void Domain::read_swc(const std::string &Filename) {
if (value_y < min_cy) min_cy = value_y;
if (value_z < min_cz) min_cz = value_z;
}
/* shift the swc data
/* shift the swc data */
printf(" shift swc data by %f, %f, %f \n",min_cx,min_cy, min_cz);
for (count=0; count<number_of_lines; count++){
List_cx[count] -= min_cx;
List_cy[count] -= min_cy;
List_cz[count] -= min_cz;
List_cx[count] -= offset_x*voxel_length;
List_cy[count] -= offset_y*voxel_length;
List_cz[count] -= offset_z*voxel_length;
}
*/
}
/* everybody gets the swc file */
Comm.bcast(List_cx,number_of_lines,0);
@ -309,6 +308,7 @@ void Domain::initialize(std::shared_ptr<Database> db) {
int nx = n[0];
int ny = n[1];
int nz = n[2];
offset_x = offset_y = offset_z = 0;
if (d_db->keyExists("InletLayers")) {
auto InletCount = d_db->getVector<int>("InletLayers");
@ -510,6 +510,9 @@ void Domain::Decomp(const std::string &Filename) {
xStart = offset[0];
yStart = offset[1];
zStart = offset[2];
offset_x = xStart;
offset_y = yStart;
offset_z = zStart;
}
if (database->keyExists("InletLayers")) {
auto InletCount = database->getVector<int>("InletLayers");

1
common/Domain.h
View File

@ -134,6 +134,7 @@ public: // Public variables (need to create accessors instead)
int Nx, Ny, Nz, N;
int inlet_layers_x, inlet_layers_y, inlet_layers_z;
int outlet_layers_x, outlet_layers_y, outlet_layers_z;
int offset_x, offset_y, offset_z;
int inlet_layers_phase; //as usual: 1->n, 2->w
int outlet_layers_phase;
double porosity;