## 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)