## sample code for stokes_res_ewald_3f()
## Copyright (C) 2006-2007 Kengo Ichiki <kichiki@users.sourceforge.net>
## $Id: manual.tex,v 1.5 2008/10/12 20:16:53 kichiki Exp $
# lattice vector
l = [10.0, 10.0, 10.0]
# position of particles (np = 8)
pos = [0.0, 0.0, 0.0,\
0.0, 0.0, 5.0,\
0.0, 5.0, 0.0,\
0.0, 5.0, 5.0,\
5.0, 0.0, 0.0,\
5.0, 0.0, 5.0,\
5.0, 5.0, 0.0,\
5.0, 5.0, 5.0]
# velocity of particles
u = [1.0, 1.0, 1.0,\
1.0, 1.0, 1.0,\
1.0, 1.0, 1.0,\
1.0, 1.0, 1.0,\
1.0, 1.0, 1.0,\
1.0, 1.0, 1.0,\
1.0, 1.0, 1.0,\
1.0, 1.0, 1.0]
# solve resistance problem and obtain force
f = stokes_res_3f(pos, u, l, 60.25)