using lattice boltzmann -it's very slow


I have just begun to work with lb/fluid in order to add hydrodynamic forces according to the lattice boltzmann method but the problem is that it is vert slow.It took about 1 hour to reach step 2!

I have a 606060 box and just 9 particles

neighbor 3 bin
neigh_modify every 10 one 10000

dielectric 3.0
pair_style coul/debye 1.4 3.1
pair_coeff * *

fix 21 all lb/fluid 1 1 3 0.864 setGamma 20 dx 1.0
fix vflow all lb/viscous

and just this!I had to reduce my data file to 9 particles and also reduce my box size but it wasn’t quite helpful.
Do you have any idea what is going on here?

I suggest you first try to run the several

examples in examples/USER/lb and examine

those scripts to see what is different.

If you still have Qs, then email the authors

of the package, listed in src/USER-LB/README.