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

fix dumb bug of asssigning inverse of tau

Browse Source
This commit is contained in:
Rex Zhe Li 2020-09-29 15:48:39 -04:00
parent 471e78703a
commit a02288631a
1 changed files with 3 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
models/IonModel.cpp
View File

@ -693,15 +693,9 @@ void ScaLBL_IonModel::Run(double *Velocity, double *ElectricField){
//2. ElectricField is from Poisson model //2. ElectricField is from Poisson model
//LB-related parameter //LB-related parameter
vector<double> rlx(tau.begin(),tau.end()); vector<double> rlx;
for (double item : rlx){ for (unsigned int ic=0;ic<tau.size();ic++){
item = 1.0/item; rlx.push_back(1.0/tau[ic]);
}
//**debug
if (rank==0){
for (unsigned int ic=0;ic<rlx.size();ic++){
printf("**debug: Ion %i rlx = %.5g; LB diffusivity = %.5g \n",ic+1,rlx[ic],IonDiffusivity[ic]);
}
} }
//.......create and start timer............ //.......create and start timer............