def run(self): if not self._kill: self.dataobj = casio() self.dataobj.readinp(self.parent._filename) self.dataobj.readcax() self.emit(SIGNAL('progressbar_update(int)'),90) if not self._kill: self.dataobj.calcbtf() if not self._kill: self.parent.dataobj = self.dataobj
def load_pickle(self,filename): self.statusBar().showMessage('Importing data from %s' % filename, 2000) self.dataobj = casio() self.dataobj.loadpic(filename) self.init_pinobjects() #fuetype = 'SVEA-96' #self.dataobj.btf = btf(self.dataobj,fuetype) #self.setpincoords() #self.draw_fuelmap() #self.set_pinvalues() # Update case number list box ncases = len(self.dataobj.cases) for i in range(1,ncases+1): self.case_cbox.addItem(str(i)) self.connect(self.case_cbox, SIGNAL('currentIndexChanged(int)'), self.fig_update) self.fig_update()
SK = min([SK1,SK2]) # Sum weighting factors SWJ = (WP[i-1,j] + WP[i+1,j] + WP[i,j-1] + WP[i,j+1])*WJ # Side rods SWK = (WP[i-1,j-1] + WP[i+1,j-1] + WP[i-1,j+1] + WP[i+1,j+1])*WK # Diagonal rods DOW[i-1,j-1] = (RP[i,j] + SJ + SK)/(1.0 + SWJ + SWK)*np.sqrt(Ntotrods/float(Nhotrods)) + AC[i-1,j-1] return DOW def calc_btf(fuetype,POW3): naxial_nodes = 25 DOW = np.zeros((naxial_nodes,POW3.shape[1],POW3.shape[2])) # Calculate lattice K-factors for all nodes for z in range(naxial_nodes): DOW[z,:,:] = rfact_axial(fuetype,POW3[z,:,:]) # Integrate along z-direction to get pinwise R-factors DOX = sum(DOW,0)/naxial_nodes return DOX if __name__ == '__main__': casobj = casio() casobj.loadpic('caxfiles_a10xm.p') POW3 = casobj.pow3(casobj)