예제 #1
0
파일: PM.py 프로젝트: fiskpralin/Tota
	def __init__(self, G=None, vis=True, anim=False, mtype='2a', ttype='random'):
		SimExtend.__init__(self, G, vis, anim, animDelay=0.3) #animDelay==standard value for animations
		self.testvar=range(30000000)
		self.initialize()
		self.stats={'plant attempts':0, 'mound attempts':0, 'remound attempts':0, 'stumps in WA':None, 'stumps in WA sum diameter':0, 'immobile boulder struck':0, 'immobile vol sum':0, 'number of dibble disr stones in mound':0, 'dibble distr stone vol cum':0, 'queue percent':0,'work percent':0, 'work time':0,'rest time':0  }
		if not self.G.terrain:
			self.G.terrain=PlantMTerrain(G=self.G, ttype=ttype)
		tic=time.clock()
		if self.G.automatic=='undefined':
			if mtype[0:2]=='1a':
				self.G.automatic={'mound': False, 'plant': True, 'automove': False, 'micrositeSelection': False, 'moveToMicro': False,'haltMound':False, 'clearForOtherHead': False}
			else:
				self.G.automatic={'mound': True, 'plant': True, 'automove': True, 'micrositeSelection': False, 'moveToMicro': False, 'haltMound':True, 'clearForOtherHead': True}
		#self.test=testProcess1(name='sdf', sim=self)
		#self.activate(self.test, self.test.run())
		#self.simulate(until=5000)
		self.m=PlantMachine(name="planter", sim=self, G=self.G, mtype=mtype)
		self.activate(self.m,self.m.run())
		self.simulate(until=self.G.maxtime)
		return None
		#postprocessing
		self.productivity=len(self.m.treesPlanted)/(self.now()/3600.) #trees/hour
		print "result: %d trees in %d time, productivity: %f trees/hour"%(len(self.m.treesPlanted), self.now(), self.productivity)		
		self.setQueuePerc()
		print "stumps in WA sum diameter:", self.stats['stumps in WA sum diamater']
		print "number of remounds", self.stats['remound attempts']
		if self.p: #if we have a plotter instance
			self.p.terminate()
			plt.show() #nothing is processed in paralell with this, thus last:
예제 #2
0
	def __init__(self, G=None, vis=True, anim=False):
		SimExtend.__init__(self,G, vis=vis, anim=anim, animDelay=5) #does some common stuff, i.e seed
		if not self.G.terrain:
			self.G.terrain=setUpTerrain(self.G)
		h=Harvester(self.G) #does all the harvester stuffi n the constructor.
		self.trees=0 #total number of picked up trees
		self.m=Forwarder(name="forwy", sim=self, G=self.G)
		self.activate(self.m,self.m.run())
		self.simulate(until=1000000)
		if self.p:
			self.p.terminate()
			#postprocessing
			plt.show()
		print "simulation is done."
예제 #3
0
	def __init__(self, G=None, vis=True, anim=False, head='BC',nCranes=2,series=None, bundler=False, twigCrack=False, observer=False):
		SimExtend.__init__(self,G, vis, anim, animDelay=1.2,series=series) #does some common stuff, e.g. seed

		#if no file to read the simulationsparameters from:
		setDefaultThinningParams(self.G.simParam) #sets the parameters to default values
			
		if not self.G.terrain:
			self.G.areaPoly=[(0,0), (25,0), (25,40), (0,40)] #default for thinning files.
			self.G.terrain=Terrain(G=self.G)
			self.G.terrain.readTrees(thinning=True)
		craneMax=self.G.simParam['maxCraneLength']
		startPos=[random.uniform(craneMax, 25-craneMax), -4]
		self.m=ThinningMachine(name="thinny", sim=self, G=self.G, head=head, nCranes=nCranes,startPos=startPos, bundler=bundler,twigCrack=twigCrack)
		self.treeStats() #do some statistics on the trees
		self.activate(self.m,self.m.run())
		if observer:
			self.o=Observer(name="obsy", sim=self, G=self.G)
			self.activate(self.o, self.o.run())
		self.simulate(until=10000)
		#some simulation information:
		if observer:
			self.timeStats()
			"""
			print self.o.tstep*sum([c[1] for c in self.o.lAMoni]), 'left head active time via monitor'
			print self.o.tstep*sum([c[1] for c in self.o.lWBMoni]), 'left head waiting for Bundler'
			print self.stats['oneCraneWorkTime'], 'one crane work time'
			print self.stats['oneCraneWaitDriverTime'], 'one crane wait for driver time'
			if self.m.hasBundler:
				print self.stats['oneCraneWaitBundlerTime'], 'one crane wait for bundler time'
			if len(self.m.heads)==2:
				print self.o.tstep*sum([c[1] for c in self.o.rAMoni]), 'right head active time via monitor'
				print self.o.tstep*sum([c[1] for c in self.o.rWBMoni]), 'right head waiting for Bundler'
				print self.stats['twoCranesWorkTime'], 'two cranes work time'
				print self.stats['twoCranesWaitDriverTime'], 'two cranes wait for driver time'
				if self.m.hasBundler:
					print self.stats['twoCranesWaitBundlerTime'], 'two cranes wait for bundler time'
			"""
		self.stats['productivity']=len(self.m.trees)/self.now()*3600
		self.stats['outTake']=100*len(self.m.trees)/float(len(self.G.terrain.trees))
		self.stats['groundAreaRatio']=sum([t.dbh**2*pi*0.25 for t in self.m.trees])/max(1,sum([t.dbh**2*pi*0.25 for t in self.G.terrain.trees]))
		
		print "Simulation ended. Productivity: %.1f trees/hour. %.1f %% of the trees were harvested"%(len(self.m.trees)/self.now()*3600, self.stats['outTake'])
		print "average trees per corridor: %f average log weight in corridor: %f trees in main road: %f"%(self.avTreesPerCorr, self.avCorrTreeWeight, self.treesPerRoad)
		self.setQueuePerc()
		print self.getGeneralStatsString()
		if self.p: #plot
			self.p.terminate()
			if vis: plt.show()
예제 #4
0
파일: PM.py 프로젝트: fiskpralin/pibsgraph
	def __init__(self, G=None, vis=True, anim=False, mtype='2a', ttype='random'):
		self.initialize()
		SimExtend.__init__(self, G, vis, anim, animDelay=0.3) #animDelay==standard value for animations
		if not self.G.simParam:
			self.G.simParam=paramsForSensAn()
			assert self.G.simParam
		self.stats={'plant attempts':0, 'mound attempts':0, 'remound attempts':0, 'stumps in WA':None, 'noStumps in WA':0,'stumps in WA sum diameter':0,'meanSurfBoulderDiamWA':0,'noSurfBouldersWA':0,'immobile boulder struck':0, 'immobile vol sum':0, 'number of dibble disr stones in mound':0, 'dibble distr stone vol cum':0, 'queue percent':0,'work percent':0, 'work time':0,'rest time':0, 'humus depths':[], 'inverting attempts':0, 'reinverting attempts':0}
		if not self.G.terrain:
			self.G.terrain=PlantMTerrain(ttype=ttype)
		else:
			if ttype!='random':
				raise Exception('Terrain is already specified. The input ttype is only used if G.terrain=None')
		self.stats['noSurfBoulders']=self.G.terrain.noSBoulders
	   	if self.stats['noSurfBoulders']==0: 
			self.stats['meanSurfBoulderDiam']=0
	   	else:
			self.stats['meanSurfBoulderDiam']=sum([2*sstone.radius for sstone in self.G.terrain.surfaceBoulders])/float(self.stats['noSurfBoulders']) #gives mean diameter
			print 'average surface boulder diameter in total area:',self.stats['meanSurfBoulderDiam']
		if self.G.automatic=='undefined':
			if mtype[0:2]=='1a':
				self.G.automatic={'mound': False, 'plant': True, 'automove': False, 'micrositeSelection': False, 'moveToMicro': False,'haltMound':False, 'clearForOtherHead': False}
			else:
				self.G.automatic={'mound': True, 'plant': True, 'automove': True, 'micrositeSelection': False, 'moveToMicro': False, 'haltMound':True, 'clearForOtherHead': True}
		self.m=PlantMachine(name="planter", sim=self, G=self.G, mtype=mtype)
		self.activate(self.m,self.m.run())
		self.simulate(until=self.G.maxtime)
		#postprocessing
		self.productivity=len(self.m.treesPlanted)/(self.now()/3600.) #trees/hour
		print "result: %d trees in %d time, productivity: %f trees/hour"%(len(self.m.treesPlanted), self.now(), self.productivity)		
		self.setQueuePerc()
		print "stumps in WA sum diameter:", self.stats['stumps in WA sum diameter']
		print 'Mean diameter of surfaceboulders in WA',self.stats['meanSurfBoulderDiamWA']
		print 'Number of SurfaceBoulders in WA',self.stats['noSurfBouldersWA']
		print "time to search for microsite", self.m.timeConsumption['searchTime']
		print "number of remounds", self.stats['remound attempts']
		if self.p: #if we have a plotter instance
			self.p.terminate()
			plt.show() #nothing is processed in parallell with this, thus last: