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

build passed

Browse Source
This commit is contained in:
Zhe Li 2022-04-28 16:55:58 +10:00
parent a00a3606f7
commit bdd0efd36e
2 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
common/ScaLBL.h
View File

@ -373,6 +373,7 @@ extern "C" void ScaLBL_D3Q7_AAeven_Poisson_ElectricPotential(int *Map, double *d
* @param finish - lattice node to finish loop
* @param Np - size of local sub-domain (derived from Domain structure)
*/
extern "C" void ScaLBL_D3Q7_Poisson_Init(int *Map, double *dist, double *Psi, int start, int finish, int Np);
extern "C" void ScaLBL_D3Q19_Poisson_Init(int *Map, double *dist, double *Psi, int start, int finish, int Np);

5
models/PoissonSolver.cpp
View File

@ -686,11 +686,9 @@ void ScaLBL_Poisson::Initialize(double time_conv_from_Study){
void ScaLBL_Poisson::Run(double *ChargeDensity, bool UseSlippingVelBC, int timestep_from_Study){
double error = 1.0;
if (lattice_scheme.compare("D3Q7")==0){
if (rank==0) printf("LB-Poisson Solver: Use D3Q7 lattice structure.\n");
timestep=0;
double error = 1.0;
while (timestep < timestepMax && error > tolerance) {
//************************************************************************/
// *************ODD TIMESTEP*************//
@ -762,7 +760,6 @@ void ScaLBL_Poisson::Run(double *ChargeDensity, bool UseSlippingVelBC, int times
host_Error = new double [Np];
timestep=0;
double error = 1.0;
while (timestep < timestepMax && error > tolerance) {
//************************************************************************/
// *************ODD TIMESTEP*************//