コード例 #1
0
ファイル: postprocess.py プロジェクト: nthakkar/emclaw
 def get_num_cells(self):
     solution = IO()
     solution.path = self.path
     solution.file_prefix = self.file_prefix
     solution.frame = 0
     solution.read_aux = False
     solution.read_petsc()
     self.__setattr__('num_cells',solution.x.centers.size)
     return self
コード例 #2
0
ファイル: postprocess.py プロジェクト: nthakkar/emclaw
 def _PlotFrame(self,frame):
     solution = IO()
     solution.path = self.path
     solution.file_prefix = self.file_prefix
     solution.frame = frame
     solution.read_aux = True
     solution.read_petsc()
     if solution.num_dim==1:
         self.PlotQ1D(solution)
     if solution.num_dim==2:
         self.PlotQ2D(solution)
コード例 #3
0
ファイル: postprocess.py プロジェクト: nthakkar/emclaw
 def PlotFrames(self):
     solution = IO()
     solution.path = self.path
     solution.file_prefix = self.file_prefix
     for frame in self.plot_frames:
         solution.frame = frame
         solution.read_aux = True
         solution.read_petsc()
         if solution.num_dim==1:
             self.PlotQ1D(solution)
         if solution.num_dim==2:
             self.PlotQ2D(solution)
コード例 #4
0
ファイル: postprocess.py プロジェクト: nthakkar/emclaw
 def post_frame_calc(self,frame):
     solution = IO()
     solution.path = self.path
     solution.file_prefix = self.file_prefix
     solution.frame = frame
     solution.read_aux = True
     solution.read_petsc()
     if solution.num_dim==1:
         self.QuadField1D(solution,self.quad)
         self.EMEnergy(solution)
     elif solution.num_dim==2:
         self.QuadField2D(solution,self.quad)
コード例 #5
0
ファイル: postprocess.py プロジェクト: nthakkar/emclaw
 def PeakSampleFrame(self,frame,n):
     solution = IO()
     solution.path = self.path
     solution.file_prefix = self.file_prefix
     solution.frame = frame
     solution.read_petsc()
     x = solution.x.centers.copy()
     # print x.shape
     for k,mode in enumerate(self.sample_mode):
         if self.num_dim==1:
             self.FieldIntensity1D(solution)
             self.Poyinting1D(solution)
             self.q_sampled[k,0,n] = solution.t
             self.S_sampled[k,0,n] = solution.t
             self.I_sampled[k,0,n] = solution.t
             maxif = 0.0
             if mode=='peak':
                 self.q_sampled[k,1,n],self.q_sampled[k,3,n] = self.PeakSample1D(solution.q[0])
                 if self.q_sampled[k,1,n]==-99:
                     self.q_sampled[k,1,n] = self.q_sampled[k,1,n-1]
                     self.q_sampled[k,3,n] = self.q_sampled[k,3,n-1]
                 self.q_sampled[k,2,n],self.q_sampled[k,4,n] = self.PeakSample1D(solution.q[1])
                 if self.q_sampled[k,2,n]==-99:
                     self.q_sampled[k,2,n] = self.q_sampled[k,2,n-1]
                     self.q_sampled[k,4,n] = self.q_sampled[k,4,n-1]
                 self.I_sampled[k,1,n],self.I_sampled[k,2,n] = self.PeakSample1D(self.I)
                 if self.q_sampled[k,1,n]==-99:
                     self.I_sampled[k,1,n] = self.I_sampled[k,1,n-1]
                     self.I_sampled[k,2,n] = self.I_sampled[k,2,n-1]
                 self.S_sampled[k,1,n],self.S_sampled[k,2,n] = self.PeakSample1D(self.S)
                 if self.q_sampled[k,1,n]==-99:
                     self.S_sampled[k,1,n] = self.S_sampled[k,1,n-1]
                     self.S_sampled[k,2,n] = self.S_sampled[k,2,n-1]
             if mode=='width':
                 self.q_sampled[k,1,n],self.q_sampled[k,3,n] = self.PeakWidth1D(solution.q[0],x)
                 if self.q_sampled[k,1,n]==-99:
                     self.q_sampled[k,1,n] = self.q_sampled[k,1,n-1]
                     self.q_sampled[k,3,n] = self.q_sampled[k,3,n-1]
                 self.q_sampled[k,2,n],self.q_sampled[k,4,n] = self.PeakWidth1D(solution.q[1],x)
                 if self.q_sampled[k,2,n]==-99:
                     self.q_sampled[k,2,n] = self.q_sampled[k,2,n-1]
                     self.q_sampled[k,4,n] = self.q_sampled[k,4,n-1]
                 self.I_sampled[k,1,n],self.I_sampled[k,2,n] = self.PeakWidth1D(self.I,x)
                 if self.q_sampled[k,1,n]==-99:
                     self.I_sampled[k,1,n] = self.I_sampled[k,1,n-1]
                     self.I_sampled[k,2,n] = self.I_sampled[k,2,n-1]
                 self.S_sampled[k,1,n],self.S_sampled[k,2,n] = self.PeakWidth1D(self.S,x)
                 if self.q_sampled[k,1,n]==-99:
                     self.S_sampled[k,1,n] = self.S_sampled[k,1,n-1]
                     self.S_sampled[k,2,n] = self.S_sampled[k,2,n-1]
コード例 #6
0
ファイル: postprocess.py プロジェクト: nthakkar/emclaw
    def get_num_dim(self):
        solution = IO()
        solution.path = self.path
        solution.file_prefix = self.file_prefix
        solution.frame = 50
        solution.read_aux = True
        solution.read_petsc()
        self.__setattr__('num_dim',solution.num_dim)
        if self.num_dim==1:
            self.RefIndex1D(solution)
        if self.num_dim==2:
            self.RefIndex2D(solution)

        self.__setattr__('cmax',self.co/self.n.min())
        self.__setattr__('cmin',self.co/self.n.max())

        return self