source: bacteria.py @ 22:69a4724e7cac

importˆàrandom
fromˆàmathˆàimportˆà*
importˆàvector

classˆàBacteriaˆà(object):
ˆà ˆà defˆà__init__(self,ˆàposition,ˆàvelocity,ˆàlifetime,ˆàmutated):
ˆà ˆà ˆà ˆà self.position =ˆàposition
ˆà ˆà ˆà ˆà self.velocity =ˆàvelocity
ˆà ˆà ˆà ˆà self.lifetime =ˆàlifetime
ˆà ˆà ˆà ˆà self.mutated =ˆàmutated

ˆà ˆà defˆà__eq__(self,other):
ˆà ˆà ˆà ˆà ifˆàself.position==other.position and\
ˆà ˆà ˆà ˆà ˆà ˆàself.velocity==other.velocity and\
ˆà ˆà ˆà ˆà ˆà ˆàself.lifetime==other.lifetime and\
ˆà ˆà ˆà ˆà ˆà ˆàself.mutated==other.mutated:
ˆà ˆà ˆà ˆà ˆà ˆà returnˆàTrue
ˆà ˆà ˆà ˆà returnˆàFalse

ˆà ˆà defˆà__ne__(self,other):
ˆà ˆà ˆà ˆà returnˆànot(self==other)

ˆà ˆà defˆàcheck_collision(self,ˆàbacteria,ˆàdelta,ˆàbact_r,delta_t):
ˆà ˆà ˆà ˆà forˆàbact inˆàbacteria:
ˆà ˆà ˆà ˆà ˆà ˆà ifˆà0ˆà<ˆàabs(self.position -ˆàbact.position)ˆà<ˆàdelta +ˆà2*bact_r:
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà returnˆàbact
ˆà ˆà ˆà ˆà returnˆàFalse

ˆà ˆà defˆàcheck_dish_wall(self,position,radius,delta,bact_r):
ˆà ˆà ˆà ˆà ifˆàabs(position)ˆà>ˆàradius -ˆà(bact_r +ˆàdelta):
ˆà ˆà ˆà ˆà ˆà ˆà returnˆàTrue
ˆà ˆà ˆà ˆà returnˆàFalse

ˆà ˆà defˆàmove(self,ˆàbacteria,ˆàdelta_t,ˆàradius,ˆàdelta,ˆàbact_r):
ˆà ˆà ˆà ˆà ifˆàself.check_dish_wall(self.position,radius,delta,bact_r):
ˆà ˆà ˆà ˆà ˆà ˆà self.velocity =ˆàself.velocity.angleToCoord(-(pi/2ˆà-
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.velocity.angle()))
ˆà ˆà ˆà ˆà self.position +=ˆàself.velocity*delta_t

ˆà ˆà defˆàcollision(self,ˆàother,ˆàbacteria,ˆàdelta_t,ˆàradius,ˆàdelta,ˆàbact_r):
ˆà ˆà ˆà ˆà self.velocity,ˆàother.velocity =ˆàother.velocity,ˆàself.velocity 
ˆà ˆà ˆà ˆà self.move(bacteria,ˆàdelta_t,ˆàradius,ˆàdelta,ˆàbact_r)
ˆà ˆà ˆà ˆà other.move(bacteria,ˆàdelta_t,ˆàradius,ˆàdelta,ˆàbact_r)
ˆà ˆà 
ˆà ˆà defˆàreprod(self,ˆàtype,ˆàenv,ˆàp_max):
ˆà ˆà ˆà ˆà ifˆàtypeˆà==ˆà's':
ˆà ˆà ˆà ˆà ˆà ˆà prob =ˆà1ˆà-ˆàp_max*(env/100)
ˆà ˆà ˆà ˆà else:
ˆà ˆà ˆà ˆà ˆà ˆà prob =ˆà0.001
ˆà ˆà ˆà ˆà returnˆàrandom.random()ˆà<ˆàprobˆà ˆà 

ˆà ˆà defˆàasexual(self,ˆàbacteria,ˆàdelta_t,ˆàradius,ˆàdelta,ˆàbact_r,ˆàfull_lifetime):
ˆà ˆà ˆà ˆà ifˆàself.find_place_a(bacteria,ˆàdelta,ˆàbact_r)\
ˆà ˆà ˆà ˆà ˆà ˆàandˆànot(self.check_dish_wall(self.position,radius,delta,3*bact_r)):
ˆà ˆà ˆà ˆà ˆà ˆà pos1=self.position +ˆàvector.Vector(-(bact_r +ˆà0.75*delta),0)
ˆà ˆà ˆà ˆà ˆà ˆà vel1=vector.Vector(1,0).angleToCoord(self.velocity.angle())
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.append(Bacteria(pos1,ˆàself.rnd_velocity(pi),
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàfull_lifetime,ˆàself.if_mutated_a()))
ˆà ˆà ˆà ˆà ˆà ˆà pos2=self.position +ˆàvector.Vector((bact_r +ˆà0.75*delta),0)
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.append(Bacteria(pos2,ˆàself.rnd_velocity(0),
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàfull_lifetime,ˆàself.if_mutated_a()))
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.remove(self)
ˆà ˆà ˆà ˆà ˆà ˆà returnˆà1
ˆà ˆà ˆà ˆà returnˆà0


ˆà ˆà defˆàsexualˆà(self,ˆàother,ˆàbacteria,ˆàdelta_t,ˆàradius,ˆàdelta,ˆàbact_r,
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà full_lifetime):
ˆà ˆà ˆà ˆà num_of_mut=1
ˆà ˆà ˆà ˆà center =ˆà(self.position+other.position)*0.5
ˆà ˆà ˆà ˆà ifˆàself.find_place_s(other,ˆàbacteria,ˆàdelta,ˆàbact_r)ˆàand\
ˆà ˆà ˆà ˆà ˆà ˆànot(self.check_dish_wall(center,radius,delta,bact_r*3)):
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.append(Bacteria(center +
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàvector.Vector(-(bact_r*sqrt(2)ˆà+ˆà0.75*delta),0),
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.rnd_velocity(pi),ˆàfull_lifetime,
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.if_mutated_s(other,num_of_mut)))
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.append(Bacteria(center +
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàvector.Vector((bact_r*sqrt(2)ˆà+ˆà0.75*delta),0),
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.rnd_velocity(0),ˆàfull_lifetime,
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.if_mutated_s(other,num_of_mut)))
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.append(Bacteria(center +
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàvector.Vector(0,ˆà-(bact_r*sqrt(2)ˆà+ˆà0.75*delta)),
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.rnd_velocity(3*pi/2),ˆàfull_lifetime,
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.if_mutated_s(other,num_of_mut)))
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.append(Bacteria(center +
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàvector.Vector(0,ˆà(bact_r*sqrt(2)ˆà+ˆà0.75*delta)),
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.rnd_velocity(pi/2),ˆàfull_lifetime,
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆàself.if_mutated_s(other,num_of_mut)))
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.remove(other)
ˆà ˆà ˆà ˆà ˆà ˆà bacteria.remove(self)
ˆà ˆà ˆà ˆà ˆà ˆà returnˆànum_of_mut
ˆà ˆà ˆà ˆà else:
ˆà ˆà ˆà ˆà ˆà ˆà self.collision(other,bacteria,delta_t,radius,delta,bact_r)
ˆà ˆà ˆà ˆà ˆà ˆà returnˆà0

ˆà ˆà defˆàfind_place_a(self,ˆàbacteria,ˆàdelta,ˆàbact_r):ˆà#searches for place for children
ˆà ˆà ˆà ˆà forˆàbact inˆàbacteria:
ˆà ˆà ˆà ˆà ˆà ˆà ifˆà0ˆà<ˆàabs(bact.position -ˆàself.position)ˆà<ˆà(bact_r*3ˆà+ˆàdelta*2):
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà returnˆàFalse
ˆà ˆà ˆà ˆà returnˆàTrue

ˆà ˆà ˆà ˆà ˆà ˆà 
ˆà ˆà defˆàfind_place_s(self,ˆàother,ˆàbacteria,ˆàdelta,ˆàbact_r):
ˆà ˆà ˆà ˆà forˆàbact inˆàbacteria:
ˆà ˆà ˆà ˆà ˆà ˆà ifˆàbact!=selfˆàandˆàbact!=other:
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ifˆà0ˆà<ˆàabs(bact.position -ˆà(self.position +ˆàother.position)*0.5)\
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà<ˆà((sqrt(2)ˆà+ˆà1)*bact_r*2ˆà+ˆàdelta*2):
ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà ˆà returnˆàFalse
ˆà ˆà ˆà ˆà returnˆàTrue


ˆà ˆà defˆàrnd_velocity(self,angle):
ˆà ˆà ˆà ˆà vel=vector.Vector(1,0)
ˆà ˆà ˆà ˆà returnˆàvel.angleToCoord(angle +ˆàpi*random.randrange(-4,5)/16)

ˆà ˆà defˆàif_mutated_a(self):
ˆà ˆà ˆà ˆà returnˆàself.mutated

ˆà ˆà defˆàif_mutated_s(self,other,n):
ˆà ˆà ˆà ˆà ifˆàself.mutated==Falseˆàandˆàother.mutated ==ˆàFalse:
ˆà ˆà ˆà ˆà ˆà ˆà returnˆàFalse
ˆà ˆà ˆà ˆà r =ˆà(random.randint(0,ˆà3)ˆà==ˆà0)
ˆà ˆà ˆà ˆà ifˆàr:
ˆà ˆà ˆà ˆà ˆà ˆà n+=1
ˆà ˆà ˆà ˆà returnˆàr