def init_matPara(self): if not self.mat_f: if PARAS.ACTION_PARAS[MAT_MODE]: raise ('can not file mat file with mat mode on') return print 'self.mat_f=', self.mat_f matIn = f_load_mat(self.mat_f) if not self.weight_raw_in: self.weight_raw_in = [matIn['Weight']] # use mat file weight if no weight input self.Threshold = matIn['Threshold'].flatten()[0] self.Iteration = matIn['Iteration'] self.NormMode = [] self.Outloop = 0 self.GameMode = matIn['Mode'][0] self.constraint_mode = 'lag' #o_file = 'fmin_iteration_out_sdn.txt' if 0: self.debug_inputPara()
def init_matPara(self): if not self.mat_f: if PARAS.ACTION_PARAS[MAT_MODE]: raise ('can not file mat file with mat mode on') return print 'self.mat_f=', self.mat_f matIn = f_load_mat(self.mat_f) if not self.weight_raw_in: self.weight_raw_in = [matIn['Weight'] ] # use mat file weight if no weight input self.Threshold = matIn['Threshold'].flatten()[0] self.Iteration = matIn['Iteration'] self.NormMode = [] self.Outloop = 0 self.GameMode = matIn['Mode'][0] self.constraint_mode = 'lag' #o_file = 'fmin_iteration_out_sdn.txt' if 0: self.debug_inputPara()
for c_idx in range(colNum): qmat.append([key, list(qmat_in[:, c_idx])]) for c_idx in range(colNum, var_num): qmat.append([[c_idx], [0.0]]) return qmat def f_recalculate(self): # called by external entity self.f_check() self.f_init_miqp_m() self.f_init() self.f_simulate() if __name__ == "__main__": #matIn = f_load_mat('USNET_2.mat') matIn = f_load_mat('case2a.mat') #matIn = f_load_mat('case2.mat') ''' D,C,PMAX,Beta,I,H,L,EM,RM, Weight,Threshold,Iteration,NormMode,Outloop, GameMode,RES_ratio,constraint_mode, i_sdiag = None, normScale=None, i_fileName=TEM_OUT, plot_en=False,srecord_en=False ''' #print matIn D = matIn['D'].flatten() C = matIn['S'].flatten() PMAX = matIn['PMAX'].flatten() Beta = matIn['Beta'].flatten() I = matIn['I'].flatten()[0]
self.m_miqp_rhs = np.concatenate((bineq,beq),axis=1) #print np.shape(Amat), np.shape(self.m_miqp_rhs),self.m_miqp_rhs self.m_miqp_cols = self.f_get_cols(Amat) #print self.m_miqp_cols sense_L = 'L'*(2*self.m_ResNum+self.m_RtNum) # dim: (2N+R) sense_E = 'E'*self.m_EnbNum self.m_miqp_sense = sense_L + sense_E #print self.m_miqp_sense # super(miqpSolverDelay, self).f_miqp_constraints() if __name__ == "__main__": #matIn = f_load_mat('USNET_2.mat') #matIn = f_load_mat('case2a.mat') matIn = f_load_mat('USNET_11.mat') # matIn = f_load_mat('USNET_12.mat') #matIn = f_load_mat('case2.mat') ''' D,C,PMAX,Beta,I,H,L,EM,RM, Weight,Threshold,Iteration,NormMode,Outloop, GameMode,RES_ratio,constraint_mode, i_sdiag = None, normScale=None, i_fileName=TEM_OUT, plot_en=False,srecord_en=False ''' #print matIn D = matIn['D'].flatten();C = matIn['S'].flatten();PMAX = matIn['PMAX'].flatten();Beta = matIn['Beta'].flatten(); I = matIn['I'].flatten()[0];H = matIn['H'].flatten()[0];L = matIn['L'].flatten()[0];EM = matIn['EM'];RM = matIn['RM']; Weight = matIn['Weight'] # 0.5 # Threshold = matIn['Threshold'].flatten()[0];Iteration = matIn['Iteration'];
qmat.append([key,list(qmat_in[:,c_idx])]) for c_idx in range(colNum,var_num): qmat.append([[c_idx],[0.0]]) return qmat def f_recalculate(self): # called by external entity self.f_check(); self.f_init_miqp_m() self.f_init(); self.f_simulate(); if __name__ == "__main__": #matIn = f_load_mat('USNET_2.mat') matIn = f_load_mat('case2a.mat') #matIn = f_load_mat('case2.mat') ''' D,C,PMAX,Beta,I,H,L,EM,RM, Weight,Threshold,Iteration,NormMode,Outloop, GameMode,RES_ratio,constraint_mode, i_sdiag = None, normScale=None, i_fileName=TEM_OUT, plot_en=False,srecord_en=False ''' #print matIn D = matIn['D'].flatten();C = matIn['S'].flatten();PMAX = matIn['PMAX'].flatten();Beta = matIn['Beta'].flatten(); I = matIn['I'].flatten()[0];H = matIn['H'].flatten()[0];L = matIn['L'].flatten()[0];EM = matIn['EM'];RM = matIn['RM']; Weight = 0.75 #matIn['Weight'] # 0.5 # Threshold = matIn['Threshold'].flatten()[0];Iteration = matIn['Iteration']; NormMode = [];Outloop = 0;
self.m_rescap_pen = self.m_rescap_pen + pen_dif rescap_pen = self.m_rescap_pen # print max(self.m_rescap_pen) # print res_cap_indic return rescap_pen def f_recalculate(self): # called by external entity self.f_check() self.f_init() self.f_simulate() if __name__ == "__main__": # matIn = f_load_mat('USNET_2.mat') matIn = f_load_mat("case2a.mat") """ D,C,PMAX,Beta,I,H,L,EM,RM, Weight,Threshold,Iteration,NormMode,Outloop, GameMode,RES_ratio,constraint_mode, i_sdiag = None, normScale=None, i_fileName=TEM_OUT, plot_en=False,srecord_en=False """ # print matIn D = matIn["D"].flatten() C = matIn["S"].flatten() PMAX = matIn["PMAX"].flatten() Beta = matIn["Beta"].flatten() I = matIn["I"].flatten()[0] H = matIn["H"].flatten()[0]