SPARTA (23 Dec 2017)
# advect particles via free molecular flow on a uniform grid
# particles reflect off global box boundaries

variable            x index 10
variable            y index 10
variable            z index 10

variable            lx equal $x*1.0e-5
variable            lx equal 64*1.0e-5
variable            ly equal $y*1.0e-5
variable            ly equal 40*1.0e-5
variable            lz equal $z*1.0e-5
variable            lz equal 40*1.0e-5

variable            n equal 10*$x*$y*$z
variable            n equal 10*64*$y*$z
variable            n equal 10*64*40*$z
variable            n equal 10*64*40*40

seed	    	    12345
dimension   	    3
global              gridcut 1.0e-5

boundary	    rr rr rr

create_box  	    0 ${lx} 0 ${ly} 0 ${lz}
create_box  	    0 0.00064 0 ${ly} 0 ${lz}
create_box  	    0 0.00064 0 0.0004 0 ${lz}
create_box  	    0 0.00064 0 0.0004 0 0.0004
Created orthogonal box = (0 0 0) to (0.00064 0.0004 0.0004)
create_grid 	    $x $y $z
create_grid 	    64 $y $z
create_grid 	    64 40 $z
create_grid 	    64 40 40
WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (../grid.cpp:381)
Created 102400 child grid cells
  parent cells = 1
  CPU time = 0.0328851 secs
  create/ghost percent = 69.789 30.211

balance_grid        rcb part
Balance grid migrated 99200 cells
  CPU time = 0.150901 secs
  reassign/sort/migrate/ghost percent = 47.6352 0.0875301 43.0553 9.22195

species		    ar.species Ar
mixture		    air Ar vstream 0.0 0.0 0.0 temp 273.15

global              nrho 7.07043E22
global              fnum 7.07043E6

create_particles    air n $n
create_particles    air n 1024000
Created 1024000 particles
  CPU time = 0.0147841 secs

stats		    100
compute             temp temp
stats_style	    step cpu np nattempt ncoll c_temp

# equilibrate with large timestep to unsort particles
# then benchmark with normal timestep

timestep 	    7.00E-8
run                 30
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 1.6875 1.6875 1.6875
  grid      (ave,min,max) = 1.51388 1.51388 1.51388
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 3.20138 3.20138 3.20138
Step CPU Np Natt Ncoll c_temp 
       0            0  1024000        0        0    273.67084 
      30   0.21286297  1024000        0        0    273.67084 
Loop time of 0.212881 on 64 procs for 30 steps with 1024000 particles

Particle moves    = 30720000 (30.7M)
Cells touched     = 151340692 (151M)
Particle comms    = 7559306 (7.56M)
Boundary collides = 2688228 (2.69M)
Boundary exits    = 0 (0K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Surf reactions    = 0 (0K)
Collide attempts  = 0 (0K)
Collide occurs    = 0 (0K)
Reactions         = 0 (0K)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 2.25478e+06
Particle-moves/step: 1.024e+06
Cell-touches/particle/step: 4.92645
Particle comm iterations/step: 3.03333
Particle fraction communicated: 0.246071
Particle fraction colliding with boundary: 0.0875074
Particle fraction exiting boundary: 0
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Surf-reactions/particle/step: 0
Collision-attempts/particle/step: 0
Collisions/particle/step: 0
Reactions/particle/step: 0

Move  time (%) = 0.124245 (58.3636)
Coll  time (%) = 0 (0)
Sort  time (%) = 0 (0)
Comm  time (%) = 0.0873736 (41.0434)
Modfy time (%) = 0 (0)
Outpt time (%) = 0.00125042 (0.587379)
Other time (%) = 1.19433e-05 (0.00561031)

Particles: 16000 ave 16212 max 15779 min
Histogram: 3 6 6 4 11 11 8 5 7 3
Cells:      1600 ave 1600 max 1600 min
Histogram: 64 0 0 0 0 0 0 0 0 0
GhostCell: 1488.69 ave 2320 max 992 min
Histogram: 8 8 8 18 0 14 0 6 0 2
EmptyCell: 774.312 ave 1328 max 535 min
Histogram: 18 0 16 14 0 8 6 0 0 2
timestep 	    7.00E-9
run 		    $t
run 		    100
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 1.71387 1.6875 3.375
  grid      (ave,min,max) = 1.51388 1.51388 1.51388
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 3.22775 3.20138 4.88888
Step CPU Np Natt Ncoll c_temp 
      30            0  1024000        0        0    273.67084 
     100   0.15473104  1024000        0        0    273.67084 
     130   0.22153616  1024000        0        0    273.67084 
Loop time of 0.22155 on 64 procs for 100 steps with 1024000 particles

Particle moves    = 102400000 (102M)
Cells touched     = 142466937 (142M)
Particle comms    = 2664745 (2.66M)
Boundary collides = 895520 (0.896M)
Boundary exits    = 0 (0K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Surf reactions    = 0 (0K)
Collide attempts  = 0 (0K)
Collide occurs    = 0 (0K)
Reactions         = 0 (0K)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 7.22185e+06
Particle-moves/step: 1.024e+06
Cell-touches/particle/step: 1.39128
Particle comm iterations/step: 1
Particle fraction communicated: 0.0260229
Particle fraction colliding with boundary: 0.00874531
Particle fraction exiting boundary: 0
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Surf-reactions/particle/step: 0
Collision-attempts/particle/step: 0
Collisions/particle/step: 0
Reactions/particle/step: 0

Move  time (%) = 0.192399 (86.8424)
Coll  time (%) = 0 (0)
Sort  time (%) = 0 (0)
Comm  time (%) = 0.0263642 (11.8999)
Modfy time (%) = 0 (0)
Outpt time (%) = 0.00274794 (1.24032)
Other time (%) = 3.84934e-05 (0.0173746)

Particles: 16000 ave 16310 max 15621 min
Histogram: 1 1 3 3 17 15 14 4 4 2
Cells:      1600 ave 1600 max 1600 min
Histogram: 64 0 0 0 0 0 0 0 0 0
GhostCell: 1488.69 ave 2320 max 992 min
Histogram: 8 8 8 18 0 14 0 6 0 2
EmptyCell: 774.312 ave 1328 max 535 min
Histogram: 18 0 16 14 0 8 6 0 0 2