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

update color

Browse Source
This commit is contained in:
James E McClure 2022-02-03 09:14:33 -05:00
parent a181dfc85d
commit 3adbd04d38
1 changed files with 13 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
models/ColorModel.cpp
View File

@ -115,8 +115,7 @@ void ScaLBL_ColorModel::ReadParams(string filename) {
inletB = 0.f; inletB = 0.f;
outletA = 0.f; outletA = 0.f;
outletB = 1.f; outletB = 1.f;
BoundaryCondition = 0; BoundaryCondition = 0;
if (color_db->keyExists("BC")) { if (color_db->keyExists("BC")) {
BoundaryCondition = color_db->getScalar<int>("BC"); BoundaryCondition = color_db->getScalar<int>("BC");
@ -483,12 +482,22 @@ void ScaLBL_ColorModel::Create() {
// copy the neighbor list // copy the neighbor list
ScaLBL_CopyToDevice(NeighborList, neighborList, neighborSize); ScaLBL_CopyToDevice(NeighborList, neighborList, neighborSize);
ScaLBL_Comm->Barrier();
delete[] neighborList; delete[] neighborList;
// initialize phi based on PhaseLabel (include solid component labels) // initialize phi based on PhaseLabel (include solid component labels)
double *PhaseLabel; double *PhaseLabel;
PhaseLabel = new double[N]; PhaseLabel = new double[Nx*Ny*Nz];
ScaLBL_Comm->Barrier();
AssignComponentLabels(PhaseLabel); AssignComponentLabels(PhaseLabel);
ScaLBL_CopyToDevice(Phi, PhaseLabel, N * sizeof(double)); ScaLBL_Comm->Barrier();
ScaLBL_CopyToDevice(Phi, PhaseLabel, Nx*Ny*Nz * sizeof(double));
ScaLBL_Comm->Barrier();
if (rank == 0)
printf("Model created \n");
delete[] PhaseLabel; delete[] PhaseLabel;
} }