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

example for how to structure ion comm loop

Browse Source
This commit is contained in:
James McClure 2020-08-10 12:54:55 -04:00
parent 28988ef6ba
commit b4d50ee821
1 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
models/IonModel.cpp
View File

@ -275,17 +275,13 @@ void ScaLBL_IonModel::Run(double *Velocity, double *ElectricField){
//Update ion concentration and charge density
for (int ic=0; ic<number_ion_species; ic++){
ScaLBL_Comm->SendD3Q7AA(fq, ic); //READ FROM NORMAL
}
for (int ic=0; ic<number_ion_species; ic++){
//TODO should this loop be merged with Send & Recv ?
//Sum up distribution to get ion concentration
ScaLBL_D3Q7_AAodd_IonConcentration(NeighborList, &fq[ic*Np*7],&Ci[ic*Np],ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np);
ScaLBL_Comm->RecvD3Q7AA(fq, ic); //WRITE INTO OPPOSITE
/* Update exterior ion concentration */
}
ScaLBL_D3Q7_IonChargeDensity(Ci, ChargeDensity, IonValence, number_ion_species, ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np);
for (int ic=0; ic<number_ion_species; ic++){
ScaLBL_Comm->RecvD3Q7AA(fq, ic); //WRITE INTO OPPOSITE
}
for (int ic=0; ic<number_ion_species; ic++){
//TODO should this loop be merged with Send & Recv ?
//Sum up distribution to get ion concentration