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 128*1.0e-5
variable            ly equal $y*1.0e-5
variable            ly equal 80*1.0e-5
variable            lz equal $z*1.0e-5
variable            lz equal 80*1.0e-5

variable            n equal 10*$x*$y*$z
variable            n equal 10*128*$y*$z
variable            n equal 10*128*80*$z
variable            n equal 10*128*80*80

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.00128 0 ${ly} 0 ${lz}
create_box  	    0 0.00128 0 0.0008 0 ${lz}
create_box  	    0 0.00128 0 0.0008 0 0.0008
Created orthogonal box = (0 0 0) to (0.00128 0.0008 0.0008)
create_grid 	    $x $y $z
create_grid 	    128 $y $z
create_grid 	    128 80 $z
create_grid 	    128 80 80
WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (../grid.cpp:381)
Created 819200 child grid cells
  parent cells = 1
  CPU time = 0.0295031 secs
  create/ghost percent = 42.27 57.73

balance_grid        rcb part
Balance grid migrated 818280 cells
  CPU time = 0.368648 secs
  reassign/sort/migrate/ghost percent = 80.5769 0.0377048 8.81329 10.5721

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 8192000
Created 8192000 particles
  CPU time = 0.016037 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  8192000        0        0    273.20841 
      30   0.40395379  8192000        0        0    273.20841 
Loop time of 0.404128 on 1024 procs for 30 steps with 8192000 particles

Particle moves    = 245760000 (246M)
Cells touched     = 1221147106 (1.22B)
Particle comms    = 86483587 (86.5M)
Boundary collides = 10747126 (10.7M)
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: 593871
Particle-moves/step: 8.192e+06
Cell-touches/particle/step: 4.96886
Particle comm iterations/step: 3.83333
Particle fraction communicated: 0.351903
Particle fraction colliding with boundary: 0.0437302
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.259532 (64.2202)
Coll  time (%) = 0 (0)
Sort  time (%) = 0 (0)
Comm  time (%) = 0.143727 (35.5648)
Modfy time (%) = 0 (0)
Outpt time (%) = 0.000742621 (0.183759)
Other time (%) = 0.000126316 (0.0312563)

Particles: 8000 ave 8286 max 7696 min
Histogram: 3 15 70 140 264 226 174 90 33 9
Cells:      800 ave 800 max 800 min
Histogram: 1024 0 0 0 0 0 0 0 0 0
GhostCell: 1290.78 ave 1552 max 640 min
Histogram: 8 12 42 56 93 50 283 25 130 325
EmptyCell: 632.784 ave 912 max 0 min
Histogram: 1 11 10 49 124 93 246 138 192 160
timestep 	    7.00E-9
run 		    $t
run 		    100
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 
      30            0  8192000        0        0    273.20841 
     100   0.18394709  8192000        0        0    273.20841 
     130   0.26232004  8192000        0        0    273.20841 
Loop time of 0.262513 on 1024 procs for 100 steps with 8192000 particles

Particle moves    = 819200000 (819M)
Cells touched     = 1143018746 (1.14B)
Particle comms    = 31480212 (31.5M)
Boundary collides = 3578620 (3.58M)
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: 3.04747e+06
Particle-moves/step: 8.192e+06
Cell-touches/particle/step: 1.39529
Particle comm iterations/step: 1
Particle fraction communicated: 0.038428
Particle fraction colliding with boundary: 0.00436843
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.177874 (67.7582)
Coll  time (%) = 0 (0)
Sort  time (%) = 0 (0)
Comm  time (%) = 0.0827971 (31.5402)
Modfy time (%) = 0 (0)
Outpt time (%) = 0.00156433 (0.595907)
Other time (%) = 0.000277504 (0.105711)

Particles: 8000 ave 8232 max 7717 min
Histogram: 2 22 56 117 201 220 220 97 68 21
Cells:      800 ave 800 max 800 min
Histogram: 1024 0 0 0 0 0 0 0 0 0
GhostCell: 1290.78 ave 1552 max 640 min
Histogram: 8 12 42 56 93 50 283 25 130 325
EmptyCell: 632.784 ave 912 max 0 min
Histogram: 1 11 10 49 124 93 246 138 192 160
