GopherProxy

	tadd explicit call for periodic boundaries - sphere - GPU-based 3D discrete ele…
	git clone git://src.adamsgaard.dk/sphere
	Log
	Files
	Refs
	LICENSE
	---
	commit 58668c68b1fce324d35d1c86991240dabd48591e
	parent ad617bec05a87782ca2fc5a55b50c8765bb216ff
	Author: Anders Damsgaard <[email protected]>
	Date: Fri, 20 Jun 2014 21:40:33 +0200

	add explicit call for periodic boundaries

	Diffstat:
	M python/cube-init.py \| 2 ++
	M python/fluidshear.py \| 1 +
	M python/segregation.py \| 2 ++
	M python/shear-test.py \| 9 +++++++++

	4 files changed, 14 insertions(+), 0 deletions(-)
	---
	diff --git a/python/cube-init.py b/python/cube-init.py
	t@@ -5,6 +5,8 @@ init = sphere.sim('cube-init', np=1e2)

	init.generateRadii(psd='uni', radius_mean=0.01, radius_variance=0.002)

	+init.periodicBoundariesXY()
	+
	# Initialize positions in random grid (also sets world size)
	init.initRandomGridPos(gridnum=(6, 6, 1e12))

	diff --git a/python/fluidshear.py b/python/fluidshear.py
	t@@ -9,6 +9,7 @@ sid = 'fluidshear'
	sim = sphere.sim(sid + '-init', np = 24000, fluid = False)
	#sim.cleanup()
	sim.radius[:] = 0.05
	+sim.periodicBoundariesXY()
	sim.initRandomGridPos(gridnum = [20, 20, 9000])
	sim.initTemporal(total = 10.0, file_dt = 0.05)
	sim.g[2] = -9.81
	diff --git a/python/segregation.py b/python/segregation.py
	t@@ -32,6 +32,8 @@ init.generateRadii(radius_mean = 0.08)
	# Use default params
	init.defaultParams(gamma_n = 100.0, mu_s = 0.4, mu_d = 0.4)

	+init.periodicBoundariesXY()
	+
	# Initialize positions in random grid (also sets world size)
	hcells = np*(1.0/3.0)0.6
	init.initRandomGridPos(gridnum = numpy.array([hcells, hcells, 1e9]))
	diff --git a/python/shear-test.py b/python/shear-test.py
	t@@ -34,6 +34,9 @@ init.defaultParams(gamma_n = 100.0, mu_s = 0.6, mu_d = 0.6)
	# Add gravity
	init.g[2] = -9.81

	+# Periodic x and y boundaries
	+init.periodicBoundariesXY()
	+
	# Initialize positions in random grid (also sets world size)
	hcells = np**(1.0/3.0)
	init.initRandomGridPos(gridnum = numpy.array([hcells, hcells, 1e9]), periodic …
	t@@ -71,6 +74,9 @@ for devs in devslist:
	lastf = status(sim_id + "-init")
	cons.readbin("../output/" + sim_id + "-init.output{:0=5}.bin".format(lastf…

	+ # Periodic x and y boundaries
	+ cons.periodicBoundariesXY()
	+
	# Setup consolidation experiment
	cons.consolidate(deviatoric_stress = devs, periodic = init.periodic)

	t@@ -117,6 +123,9 @@ for devs in devslist:
	lastf = status(sim_id + "-cons-devs{}".format(devs))
	shear.readbin("../output/" + sim_id + "-cons-devs{}.output{:0=5}.bin".form…

	+ # Periodic x and y boundaries
	+ shear.periodicBoundariesXY()
	+
	# Setup shear experiment
	shear.shear(shear_strain_rate = 0.05, periodic = init.periodic)