# -*- encoding=utf-8 -*-
O.engines = [
        ForceResetter(),
        InsertionSortCollider([
                Bo1_Sphere_Aabb(),
                Bo1_Box_Aabb(),
        ]),
        InteractionLoop([Ig2_Sphere_Sphere_ScGeom(), Ig2_Box_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()]),
        NewtonIntegrator(damping=.2, exactAsphericalRot=True, gravity=[0, 0, -9.81])
]
# support
O.bodies.append(box([0, 0, -1.5], [3, 3, .2], fixed=True))
# stand-alone sphere
O.bodies.append(sphere([0, 0, 0], .5))
# clumps
relPos = [(0, -.5, -.5), (0, .5, 0), (.5, 0, 0), (0, 0, .5)]
coords = [(-2, 0, 0), (2, 0, 0), (0, 2, 0), (0, -2, 0)]
for i, cc in enumerate(coords):
	# This shorthand command does something like this:
	# O.bodies.appendClumped([sphere(...),sphere(...),sphere(...)])
	# and returns tuple of clumpId,[bodyId1,bodyId2,bodyId3]
	clump, spheres = O.bodies.appendClumped(
	        [sphere([relPos[j][0] + coords[i][0], relPos[j][1] + coords[i][1], relPos[j][2] + coords[i][2]], .5) for j in range(0, i + 1)]
	)
	print(clump, spheres)
O.dt = PWaveTimeStep()
O.saveTmp('init')
