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

fix bad ref to D3Q7

Browse Source
This commit is contained in:
James McClure 2022-04-22 14:30:03 -04:00
parent 481a258fbd
commit d3dc35c018
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
models/IonModel.cpp
View File

@ -809,9 +809,9 @@ void ScaLBL_IonModel::AssignIonConcentrationMembrane( double *Ci, int ic) {
int idx = Map(i, j, k); int idx = Map(i, j, k);
if (!(idx < 0)) { if (!(idx < 0)) {
if (MembraneDistance(i,j,k) < 0.0) { if (MembraneDistance(i,j,k) < 0.0) {
VALUE = MembraneIonConcentration[ic]* (h * h * h * 1.0e-18); VALUE = MembraneIonConcentration[ic];//* (h * h * h * 1.0e-18);
} else { } else {
VALUE = IonConcentration[ic]* (h * h * h * 1.0e-18); VALUE = IonConcentration[ic];//* (h * h * h * 1.0e-18);
} }
Ci[idx] = VALUE; Ci[idx] = VALUE;

4
models/PoissonSolver.cpp
View File

@ -721,9 +721,9 @@ void ScaLBL_Poisson::SolvePoissonAAodd(double *ChargeDensity, bool UseSlippingVe
void ScaLBL_Poisson::SolvePoissonAAeven(double *ChargeDensity, bool UseSlippingVelBC){ void ScaLBL_Poisson::SolvePoissonAAeven(double *ChargeDensity, bool UseSlippingVelBC){
ScaLBL_Comm->SendD3Q19AA(fq); //READ FROM NORMAL ScaLBL_Comm->SendD3Q19AA(fq); //READ FROM NORMAL
ScaLBL_D3Q7_AAeven_Poisson(dvcMap, fq, ChargeDensity, Psi, ElectricField, tau, epsilon_LB, UseSlippingVelBC, ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np); ScaLBL_D3Q19_AAeven_Poisson(dvcMap, fq, ChargeDensity, Psi, ElectricField, tau, epsilon_LB, UseSlippingVelBC, ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np);
ScaLBL_Comm->RecvD3Q19AA(fq); //WRITE INTO OPPOSITE ScaLBL_Comm->RecvD3Q19AA(fq); //WRITE INTO OPPOSITE
ScaLBL_D3Q7_AAeven_Poisson(dvcMap, fq, ChargeDensity, Psi, ElectricField, tau, epsilon_LB, UseSlippingVelBC, 0, ScaLBL_Comm->LastExterior(), Np); ScaLBL_D3Q19_AAeven_Poisson(dvcMap, fq, ChargeDensity, Psi, ElectricField, tau, epsilon_LB, UseSlippingVelBC, 0, ScaLBL_Comm->LastExterior(), Np);
ScaLBL_Comm->Barrier(); ScaLBL_Comm->Barrier();
//ScaLBL_Comm->SolidDirichletAndNeumannD3Q7(fq, Psi, Psi_BCLabel); //ScaLBL_Comm->SolidDirichletAndNeumannD3Q7(fq, Psi, Psi_BCLabel);
} }