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]
