Esempio n. 1
0
	def showPhase(self, ppart, kpic, plottype=None):  #data is the particle data, ppart.  kpic is array of num particles per tile
		if not self.isGraphing():
			return
		#shape is the bounds of the histogram, [[xmin,xmax], [ymin,ymax]]
		numPart = np.sum(kpic) #number of particles
		numTiles = np.size(kpic) #number of tiles

		xvInTile = np.zeros( (2,numPart), ppart.dtype )
		isum = 0
		for k, kk in enumerate(kpic):  #De-tile particles in to one big fat array
			xvInTile[0,isum:kk+isum] = ppart[0, 0:kk, k]
			xvInTile[1,isum:kk+isum] = ppart[1, 0:kk, k]
			isum += kk

		dv1 = Graphs.DrawPhase( xvInTile ) #copy the data
		if plottype != None:
			dv1.plottype = plottype
		self._sendplot(dv1)
Esempio n. 2
0
    def showPhase(
        self,
        ppart,
        kpic,
        plottype=None,
        title=None,
        early=None,
        twophase=None,
        plotlabels=["Position", "Velocity"]
    ):  # data is the particle data, ppart.  kpic is array of num particles per tile
        if self.norun:
            return
        pt = plottype
        if pt == None:
            pt = "DRAWPHASE"
        if not self.isGraphing(pt):
            return
        if early is not None:
            self.graphBeforeEndOfFF(pt, early)
        # shape is the bounds of the histogram, [[xmin,xmax], [ymin,ymax]]
        numPart = np.sum(kpic)  # number of particles
        numTiles = np.size(kpic)  # number of tiles

        xvInTile = np.zeros((2 + 1, numPart), ppart.dtype)
        isum = 0
        for k, kk in enumerate(
                kpic):  # De-tile particles in to one big fat array
            xvInTile[0, isum:kk + isum] = ppart[0, 0:kk, k]
            xvInTile[1, isum:kk + isum] = ppart[1, 0:kk, k]
            if twophase is not None:
                xvInTile[twophase, isum:kk + isum] = ppart[twophase, 0:kk, k]
            isum += kk

        dv1 = Graphs.DrawPhase(xvInTile,
                               title=title,
                               twophase=twophase,
                               plotlabels=plotlabels)  # copy the data
        if plottype != None:
            dv1.plottype = plottype
        self._sendplot(dv1)