def __init__(self, particle):
Collider.__init__(self, particle)
posVariance = 40
radiusVariance = 15
c1 = Vec(random.uniform(-1, 1), random.uniform(-1, 1)) * posVariance
dist1 = c1.magnitude()
r1 = dist1 + 2 + random.uniform(0, radiusVariance)
c2 = -c1.copy() / c1.magnitude() * random.uniform(r1 - dist1, posVariance)
dist2 = c2.magnitude()
r2 = dist2 + (r1 - dist1)
maximum = (r1 + r2 - (dist1 + dist2)) * 0.5
r3 = random.uniform(2, maximum)
r4 = random.uniform(2, maximum)
c1c2 = c2 - c1
c1c2Dir = math.atan2(c1c2.y, c1c2.x)
c1c2Dist = c1c2.magnitude()
c1c3Angle = calcAngle(r1, r3, r2, c1c2Dist)
c1c4Angle = calcAngle(r1, r4, r2, c1c2Dist)
c2c3Angle = calcAngle(r2, r3, r1, c1c2Dist)
c2c4Angle = calcAngle(r2, r4, r1, c1c2Dist)
a1 = c1c2Dir - c1c3Angle
a2 = c1c2Dir + c1c4Angle
a3 = c1c2Dir + math.pi - c2c4Angle
a4 = c1c2Dir + math.pi + c2c3Angle
n1 = makeVec(a1)
n2 = makeVec(a2)
n3 = makeVec(a3)
n4 = makeVec(a4)
c3 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r3) + c1
c4 = Vec(math.cos(a2), math.sin(a2)) * (r1 - r4) + c1
self.arcs = [Arc(c1, r1, n1, n2), Arc(c4, r4, n2, n3), Arc(c2, r2, n3, n4), Arc(c3, r3, n4, n1)]
def __init__(self, particle, n):
Collider.__init__(self, particle)
posVariance = 100
radiusVariance = 40
self.n = n
'''angle = 2 * math.pi / (n / 2)
c = []
dist = []
r = []
# starting point
start = Vec(random.uniform(-1, 1), random.uniform(-1, 1))
# rotation angle
angle = 2 * math.pi / (n / 2)
# create circles
for i in range(n/2):
c += [start.rotated(-angle * i) * 40]
r += [n ** 2]
dist += [c[i].magnitude()]
# small circle radii
for i in range(n/2):
r += [n]
cc = {}
ccDir = {}
ccDist = {}
ccAngle = {}
# circle-circle vectors and vars
for i in range(n/2):
j = (i+1) % (n/2)
cc[(i, j)] = c[j] - c[i]
ccDir[(i, j)] = math.atan2(cc[(i, j)].x, cc[(i, j)].y)
ccDist[(i, j)] = cc[(i, j)].magnitude()
ccAngle[(i, i + n/2)] = calcAngle(r[i], r[i+n/2], r[j], ccDist[(i, j)])
ccAngle[(j, i + n/2)] = calcAngle(r[j], r[i+n/2], r[i], ccDist[(i, j)])
# WORKS TO HEAR!!!!
a = []
a += [ccDir[(0, 1)] - ccAngle[(0, n/2)]]
for i in range(n/2 - 1):
j = (n/2 - 1 - i)
k = (n/2 - i) % (n/2)
a += [ccDir[(j, k)] + math.pi + ccAngle[(k, n-i-1)]]
a += [ccDir[(j, k)] - ccAngle[(j, n-i-1)]]
a += [ccDir[(0, 1)] + math.pi + ccAngle[1, n/2]]
# WORKS TO HEAR
nv = []
for i in range(n):
nv += [makeVec(a[i])]
print "n:" + str(n)
c += [Vec(math.cos(a[0]), math.sin(a[0])) * (r[0] - r[n/2]) + c[0]]
print "a[0], a[0], r[0], r["+str(n/2)+"], c[0]"
for i in range(1, n/2):
m = n - 2 * i
c += [Vec(math.cos(a[m]), math.sin(a[m])) * (r[i] - r[i+n/2]) + c[i]]
print "a[" + str(m) + "], a[" + str(m) + "], r[" + str(i) + "], r[" + str(i+n/2) + "], c[" + str(i) + "]"
c5 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r5) + c1
c6 = Vec(math.cos(a7), math.sin(a7)) * (r2 - r6) + c2
c7 = Vec(math.cos(a5), math.sin(a5)) * (r3 - r7) + c3
c8 = Vec(math.cos(a3), math.sin(a3)) * (r4 - r8) + c4
self.arcs = []
#self.arcs += [Arc(c[0], r[0], n[0], n[1])]
for i in range(n/2):
v = n - 1 - i
w = n/2 - 1 - i
x = (1 + i * 2) % 8
y = (2 + i * 2) % 8
z = (3 + i * 2) % 8
self.arcs += [Arc(c[v], r[v], nv[x], nv[y])]
self.arcs += [Arc(c[w], r[w], nv[y], nv[z])]'''
if n == 4:
self.c1 = Vec(random.uniform(-1, 1), random.uniform(-1, 1)) * posVariance
c1 = self.c1
dist1 = c1.magnitude()
self.r1 = dist1 + 2 + random.uniform(0, radiusVariance)
r1 = self.r1
self.c2 = -c1.copy() / c1.magnitude() * random.uniform(r1 - dist1, posVariance)
c2 = self.c2
dist2 = c2.magnitude()
self.r2 = dist2 + 2 + random.uniform(0, radiusVariance)
r2 = self.r2
maximum = (r1 + r2 - (dist1 + dist2)) * 0.5
r3 = random.uniform(2, maximum)
r4 = random.uniform(2, maximum)
c1c2 = c2 - c1
c1c2Dir = math.atan2(c1c2.y, c1c2.x)
c1c2Dist = c1c2.magnitude()
c1c3Angle = calcAngle(r1, r3, r2, c1c2Dist)
c1c4Angle = calcAngle(r1, r4, r2, c1c2Dist)
c2c3Angle = calcAngle(r2, r3, r1, c1c2Dist)
c2c4Angle = calcAngle(r2, r4, r1, c1c2Dist)
a1 = c1c2Dir - c1c3Angle
a2 = c1c2Dir + c1c4Angle
a3 = c1c2Dir + math.pi - c2c4Angle
a4 = c1c2Dir + math.pi + c2c3Angle
n1 = makeVec(a1)
n2 = makeVec(a2)
n3 = makeVec(a3)
n4 = makeVec(a4)
c3 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r3) + c1
c4 = Vec(math.cos(a2), math.sin(a2)) * (r1 - r4) + c1
self.arcs = [Arc(c1, r1, n1, n2),
Arc(c4, r4, n2, n3),
Arc(c2, r2, n3, n4),
Arc(c3, r3, n4, n1)]
elif n == 6:
posVariance = 80
radiusVariance = 50
#c1 = Vec(random.uniform(-1, 1), random.uniform(-1, 1)) * posVariance
c1 = Vec(0, 1) * 40
dist1 = c1.magnitude()
#r1 = dist1 + 2 + random.uniform(0, radiusVariance)
r1 = 60
angle = 2 * math.pi / (n / 2)
#c2 = c1.rotated(angle) / c1.magnitude() * random.uniform(0, posVariance / 2)
c2 = c1.rotated(-angle) / c1.magnitude() * 40
dist2 = c2.magnitude()
#r2 = dist2 + 2 + random.uniform(r1, radiusVariance)
r2 = 60
c3 = c2.rotated(-angle) / c2.magnitude() * 40
dist3 = c3.magnitude()
#r3 = dist3 + 2 + random.uniform(r2, radiusVariance)
r3 = 60
#maximum = (r1 + r2 + r3 - (dist1 + dist2 + dist3)) / 6
c1c2 = c2 - c1
c2c3 = c3 - c2
c3c1 = c1 - c3
c1c2Dir = math.atan2(c1c2.y, c1c2.x)
c1c2Dist = c1c2.magnitude()
c2c3Dir = math.atan2(c2c3.y, c2c3.x)
c2c3Dist = c2c3.magnitude()
c3c1Dir = math.atan2(c3c1.y, c3c1.x)
c3c1Dist = c3c1.magnitude()
#r4 = random.uniform(0, (r1 + r2 - c1c2Dist) / 2)
#r5 = random.uniform(0, (r2 + r3 - c2c3Dist) / 2)
#r6 = random.uniform(0, (r3 + r1 - c3c1Dist) / 2)
r4 = 5
r5 = 5
r6 = 5
# Calculate Angles
c1c4Angle = calcAngle(r1, r4, r2, c1c2Dist)
c2c4Angle = calcAngle(r2, r4, r1, c1c2Dist)
c2c5Angle = calcAngle(r2, r5, r3, c2c3Dist)
c3c5Angle = calcAngle(r3, r5, r2, c2c3Dist)
c3c6Angle = calcAngle(r3, r6, r1, c3c1Dist)
c1c6Angle = calcAngle(r1, r6, r3, c3c1Dist)
# Get local angles
a1 = c1c2Dir - c1c4Angle
a2 = c3c1Dir + math.pi + c1c6Angle
a3 = c3c1Dir - c3c6Angle
a4 = c2c3Dir + math.pi + c3c5Angle
a5 = c2c3Dir - c2c5Angle
a6 = c1c2Dir + math.pi + c2c4Angle
n1 = makeVec(a1)
n2 = makeVec(a2)
n3 = makeVec(a3)
n4 = makeVec(a4)
n5 = makeVec(a5)
n6 = makeVec(a6)
c4 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r4) + c1
c5 = Vec(math.cos(a5), math.sin(a5)) * (r2 - r5) + c2
c6 = Vec(math.cos(a3), math.sin(a3)) * (r3 - r6) + c3
self.arcs = [Arc(c1, r1, n1, n2),
Arc(c6, r6, n2, n3),
Arc(c3, r3, n3, n4),
Arc(c5, r5, n4, n5),
Arc(c2, r2, n5, n6),
Arc(c4, r4, n6, n1)]
elif n == 8:
posVariance = 80
radiusVariance = 50
#c1 = Vec(random.uniform(-1, 1), random.uniform(-1, 1)) * posVariance
c1 = Vec(0, 1) * 20
dist1 = c1.magnitude()
#r1 = dist1 + 2 + random.uniform(0, radiusVariance)
r1 = 60
angle = 2 * math.pi / (n / 2)
#c2 = c1.rotated(angle) / c1.magnitude() * random.uniform(0, posVariance / 2)
c2 = c1.rotated(-angle) / c1.magnitude() * 40
dist2 = c2.magnitude()
#r2 = dist2 + 2 + random.uniform(r1, radiusVariance)
r2 = 60
c3 = c2.rotated(-angle) / c2.magnitude() * 40
dist3 = c3.magnitude()
#r3 = dist3 + 2 + random.uniform(r2, radiusVariance)
r3 = 60
c4 = c3.rotated(-angle) / c3.magnitude() * 40
dist4 = c4.magnitude()
#r3 = dist3 + 2 + random.uniform(r2, radiusVariance)
r4 = 60
#maximum = (r1 + r2 + r3 - (dist1 + dist2 + dist3)) / 6
c1c2 = c2 - c1
c2c3 = c3 - c2
c3c4 = c4 - c3
c4c1 = c1 - c4
c1c2Dir = math.atan2(c1c2.y, c1c2.x)
c1c2Dist = c1c2.magnitude()
c2c3Dir = math.atan2(c2c3.y, c2c3.x)
c2c3Dist = c2c3.magnitude()
c3c4Dir = math.atan2(c3c4.y, c3c4.x)
c3c4Dist = c3c4.magnitude()
c4c1Dir = math.atan2(c4c1.y, c4c1.x)
c4c1Dist = c4c1.magnitude()
#r4 = random.uniform(0, (r1 + r2 - c1c2Dist) / 2)
#r5 = random.uniform(0, (r2 + r3 - c2c3Dist) / 2)
#r6 = random.uniform(0, (r3 + r1 - c3c1Dist) / 2)
r5 = 4
r6 = 4
r7 = 4
r8 = 4
# Calculate Angles
c1c5Angle = calcAngle(r1, r5, r2, c1c2Dist)
c2c5Angle = calcAngle(r2, r5, r1, c1c2Dist)
c2c6Angle = calcAngle(r2, r6, r3, c2c3Dist)
c3c6Angle = calcAngle(r3, r6, r2, c2c3Dist)
c3c7Angle = calcAngle(r3, r7, r4, c3c4Dist)
c4c7Angle = calcAngle(r4, r7, r3, c3c4Dist)
c4c8Angle = calcAngle(r4, r8, r1, c4c1Dist)
c1c8Angle = calcAngle(r1, r8, r4, c4c1Dist)
# Get local angles
a1 = c1c2Dir - c1c5Angle
a2 = c4c1Dir + math.pi + c1c8Angle
a3 = c4c1Dir - c4c8Angle
a4 = c3c4Dir + math.pi + c4c7Angle
a5 = c3c4Dir - c3c7Angle
a6 = c2c3Dir + math.pi + c3c6Angle
a7 = c2c3Dir - c2c6Angle
a8 = c1c2Dir + math.pi + c2c5Angle
n1 = makeVec(a1)
n2 = makeVec(a2)
n3 = makeVec(a3)
n4 = makeVec(a4)
n5 = makeVec(a5)
n6 = makeVec(a6)
n7 = makeVec(a7)
n8 = makeVec(a8)
c5 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r5) + c1
c6 = Vec(math.cos(a7), math.sin(a7)) * (r2 - r6) + c2
c7 = Vec(math.cos(a5), math.sin(a5)) * (r3 - r7) + c3
c8 = Vec(math.cos(a3), math.sin(a3)) * (r4 - r8) + c4
self.arcs = [Arc(c1, r1, n1, n2),
Arc(c8, r8, n2, n3),
Arc(c4, r4, n3, n4),
Arc(c7, r7, n4, n5),
Arc(c3, r3, n5, n6),
Arc(c6, r6, n6, n7),
Arc(c2, r2, n7, n8),
Arc(c5, r5, n8, n1)]
self.boundRadius = 0
for arc in self.arcs:
radius = abs(self.particle.p - arc.pos) + arc.radius
if radius > self.boundRadius:
self.boundRadius = radius
self.boundRadius = 150
def __init__(self, particle, n):
Collider.__init__(self, particle)
posVariance = 40
radiusVariance = 15
self.n = n
if n == 4:
self.c1 = Vec(random.uniform(-1, 1), random.uniform(-1, 1)) * posVariance
dist1 = c1.magnitude()
self.r1 = dist1 + 2 + random.uniform(0, radiusVariance)
self.c2 = -c1.copy() / c1.magnitude() * random.uniform(r1 - dist1, posVariance)
dist2 = c2.magnitude()
self.r2 = dist2 + 2 + random.uniform(0, radiusVariance)
maximum = (r1 + r2 - (dist1 + dist2)) * 0.5
self.r3 = random.uniform(2, maximum)
self.r4 = random.uniform(2, maximum)
c1c2 = c2 - c1
c1c2Dir = math.atan2(c1c2.y, c1c2.x)
c1c2Dist = c1c2.magnitude()
c1c3Angle = calcAngle(r1, r3, r2, c1c2Dist)
c1c4Angle = calcAngle(r1, r4, r2, c1c2Dist)
c2c3Angle = calcAngle(r2, r3, r1, c1c2Dist)
c2c4Angle = calcAngle(r2, r4, r1, c1c2Dist)
a1 = c1c2Dir - c1c3Angle
a2 = c1c2Dir + c1c4Angle
a3 = c1c2Dir + math.pi - c2c4Angle
a4 = c1c2Dir + math.pi + c2c3Angle
n1 = makeVec(a1)
n2 = makeVec(a2)
n3 = makeVec(a3)
n4 = makeVec(a4)
c3 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r3) + c1
c4 = Vec(math.cos(a2), math.sin(a2)) * (r1 - r4) + c1
self.arcs = [Arc(c1, r1, n1, n2),
Arc(c4, r4, n2, n3),
Arc(c2, r2, n3, n4),
Arc(c3, r3, n4, n1)]
elif n == 6:
posVariance = 80
radiusVariance = 50
#c1 = Vec(random.uniform(-1, 1), random.uniform(-1, 1)) * posVariance
c1 = Vec(0, 1) * 40
dist1 = c1.magnitude()
#r1 = dist1 + 2 + random.uniform(0, radiusVariance)
r1 = 60
angle = 2 * math.pi / (n / 2)
print str(angle)
#c2 = c1.rotated(angle) / c1.magnitude() * random.uniform(0, posVariance / 2)
c2 = c1.rotated(-angle) / c1.magnitude() * 40
dist2 = c2.magnitude()
#r2 = dist2 + 2 + random.uniform(r1, radiusVariance)
r2 = 60
c3 = c2.rotated(-angle) / c2.magnitude() * 40
dist3 = c3.magnitude()
#r3 = dist3 + 2 + random.uniform(r2, radiusVariance)
r3 = 60
#maximum = (r1 + r2 + r3 - (dist1 + dist2 + dist3)) / 6
c1c2 = c2 - c1
c2c3 = c3 - c2
c3c1 = c1 - c3
c1c2Dir = math.atan2(c1c2.y, c1c2.x)
c1c2Dist = c1c2.magnitude()
c2c3Dir = math.atan2(c2c3.y, c2c3.x)
c2c3Dist = c2c3.magnitude()
c3c1Dir = math.atan2(c3c1.y, c3c1.x)
c3c1Dist = c3c1.magnitude()
#r4 = random.uniform(0, (r1 + r2 - c1c2Dist) / 2)
#r5 = random.uniform(0, (r2 + r3 - c2c3Dist) / 2)
#r6 = random.uniform(0, (r3 + r1 - c3c1Dist) / 2)
r4 = 5
r5 = 5
r6 = 5
# Calculate Angles
c1c4Angle = calcAngle(r1, r4, r2, c1c2Dist)
c2c4Angle = calcAngle(r2, r4, r1, c1c2Dist)
c2c5Angle = calcAngle(r2, r5, r3, c2c3Dist)
c3c5Angle = calcAngle(r3, r5, r2, c2c3Dist)
c3c6Angle = calcAngle(r3, r6, r1, c3c1Dist)
c1c6Angle = calcAngle(r1, r6, r3, c3c1Dist)
print "c1c2Dir:" + str(c1c2Dir)
print "c1c4Angle:" + str(c1c4Angle)
# Get local angles
a1 = c1c2Dir - c1c4Angle
a2 = c3c1Dir + math.pi + c1c6Angle
a3 = c3c1Dir - c3c6Angle
a4 = c2c3Dir + math.pi + c3c5Angle
a5 = c2c3Dir - c2c5Angle
a6 = c1c2Dir + math.pi + c2c4Angle
print "a1:" + str(a1)
n1 = makeVec(a1)
n2 = makeVec(a2)
n3 = makeVec(a3)
n4 = makeVec(a4)
n5 = makeVec(a5)
n6 = makeVec(a6)
c4 = Vec(math.cos(a1), math.sin(a1)) * (r1 - r4) + c1
c5 = Vec(math.cos(a5), math.sin(a5)) * (r2 - r5) + c2
c6 = Vec(math.cos(a3), math.sin(a3)) * (r3 - r6) + c3
self.arcs = [Arc(c1, r1, n1, n2),
Arc(c6, r6, n2, n3),
Arc(c3, r3, n3, n4),
Arc(c5, r5, n4, n5),
Arc(c2, r2, n5, n6),
Arc(c4, r4, n6, n1)]
self.boundRadius = 0
for arc in self.arcs:
radius = abs(self.particle.p - arc.pos) + arc.radius
if radius > self.boundRadius:
self.boundRadius = radius
self.boundRadius = 150