def __init__(self, **kwargs): """Constructor""" args = (Set(),) if 'expr' in kwargs: raise ValueError( "ConstraintList does not accept the 'expr' keyword") Constraint.__init__(self, *args, **kwargs)
def __init__(self, **kwargs): """Constructor""" args = (Set(), ) self._nobjectives = 0 if 'expr' in kwargs: raise ValueError( "ObjectiveList does not accept the 'expr' keyword") Objective.__init__(self, *args, **kwargs)
def __init__(self, **kwds): #kwds['dense'] = False args = (Set(), ) IndexedVar.__init__(self, *args, **kwds)
def __init__(self, **kwds): kwds['dense'] = False args = (Set(),) self._nvars = 0 IndexedVar.__init__(self, *args, **kwds)
def declare_parameters(model): model.Alpha = Param(within=NonNegativeReals, doc="factor for area cost" ) model.Beta = Param(within=PercentFraction , doc="exponent for the area cost" ) model.Cost_hx = Param(within=NonNegativeReals, doc="fixed charge for exchangers" ) model.Cost_cu = Param(within=NonNegativeReals, doc="utility cost coeffiecient for cooling utility j" ) model.Cost_hu = Param(within=NonNegativeReals, doc="utility cost coeffiecient for heating utility k" ) model.Delta_t_min = Param(within=NonNegativeReals, doc="minimum temperature approach" ) model.Number_stages = Param(within=PositiveIntegers, doc="number of stages" ) model.Number_hot_stream = Param(within=PositiveIntegers, doc="number of hot streams" ) model.Number_cold_stream = Param(within=PositiveIntegers, doc="number of cold streams" ) model.First_stage = Param(within=PositiveIntegers, default=1, doc='Index of the first stage') model.Last_stage = Param(within=PositiveIntegers, default=model.Number_stages+1, doc='Index of the last stage' ) model.HP = RangeSet(1, model.Number_hot_stream, doc="set of hot process streams i" ) model.CP = RangeSet(1, model.Number_cold_stream, doc="set of cold process streams j" ) model.ST = RangeSet(model.First_stage, model.Number_stages, doc="set of stages in the superstructure" ) model.K = RangeSet(model.First_stage, model.Last_stage, doc="set of temperature locations" ) model.First_Stage_Set = RangeSet(model.First_stage, model.First_stage) model.K_Take_First_Stage = model.K - model.First_Stage_Set model.Fh = Param(model.HP, doc="flow capicity of hot stream i" ) model.Fc = Param(model.CP, doc="flow capacity of cold stream j" ) model.Hh = Param(model.HP, doc="heat transfer coefficient for hot stream i" ) model.Hc = Param(model.CP, doc="heat transfer coefficient for cold stream j" ) model.H_cu = Param( doc="heat transfer coefficient for cold utility" ) model.H_hu = Param( doc="heat transfer coefficient for hot utility" ) model.Th_in = Param(model.HP, doc="inlet temperature of hot stream i" ) model.Tc_in = Param(model.CP, doc="inlet temperature of cold stream j" ) model.Th_out = Param(model.HP, doc="outlet temperature of hot stream i" ) model.Tc_out = Param(model.CP, doc="outlet temperature of cold stream j" ) model.T_cu_in = Param(within=PositiveReals, doc="inlet temperature of cold utility" ) model.T_cu_out = Param(within=PositiveReals, doc="outlet temperature of cold utility" ) model.T_hu_in = Param(within=PositiveReals, doc="inlet temperature of hot utility" ) model.T_hu_out = Param(within=PositiveReals, doc="outlet temperature of hot utility" ) ############################### # Initialised Parameters # ############################### model.Reclmtd_gradient_points = Set(model.HP, model.CP, model.ST, dimen=2, initialize=[]) model.Reclmtd_cu_gradient_points = Set(model.HP, dimen=1, initialize=[]) model.Reclmtd_hu_gradient_points = Set(model.CP, dimen=1, initialize=[]) model.U = Param(model.HP, model.CP, within=PositiveReals, initialize=u_init) model.U_cu = Param(model.HP, within=PositiveReals, initialize=u_cu_init) model.U_hu = Param(model.CP, within=PositiveReals, initialize=u_hu_init) model.Ech = Param(model.HP, within=PositiveReals, initialize=ech_init) model.Ecc = Param(model.CP, within=PositiveReals, initialize=ecc_init) model.Omega_ij = Param(model.HP, model.CP, within=PositiveReals, initialize=Omega_ij_init) model.Omega_i = Param(model.HP, within=PositiveReals, initialize=Omega_i_init) model.Omega_j = Param(model.CP, within=PositiveReals, initialize=Omega_j_init) model.Gamma = Param(model.HP, model.CP, within=NonNegativeReals, initialize=Gamma_init) model.Th_breakpoints = Set(model.HP, model.ST, dimen=1, ordered=True, initialize=lambda model, i, k: th_breakpoints_init(model, i)) model.Thx_breakpoints = Set(model.HP, model.CP, model.ST, dimen=1, ordered=True, initialize=lambda model, i, j, k: thx_breakpoints_init(model, i, j, k)) model.Tc_breakpoints = Set(model.CP, model.K_Take_First_Stage, dimen=1, ordered=True, initialize=lambda model, j, k: tc_breakpoints_init(model, j)) model.Tcx_breakpoints = Set(model.HP, model.CP, model.ST, dimen=1, ordered=True, initialize=lambda model, i, j, k: tcx_breakpoints_init(model, i, j, k)) model.Q_breakpoints = Set(model.HP, model.CP, model.ST, dimen=1, ordered=True, initialize=lambda model, i, j, k: q_breakpoints_init(model, i, j)) model.Q_cu_breakpoints = Set(model.HP, dimen=1, ordered=True, initialize=lambda model, i: q_cu_breakpoints_init(model, i)) model.Q_hu_breakpoints = Set(model.CP, dimen=1, ordered=True, initialize=lambda model, j: q_hu_breakpoints_init(model, j)) model.Area_beta_breakpoints = Set(model.HP, model.CP, model.ST, dimen=1, ordered=True, initialize=lambda model,i,j,k: area_beta_breakpoints_init(model,i,j)) model.Area_beta_exp = Set(model.HP, model.CP, model.ST, dimen=1, ordered=True, initialize=lambda model,i,j,k: map(lambda A: pow(A, model.Beta), model.Area_beta_breakpoints[i,j,k])) model.Area_beta_gradients = Set(model.HP, model.CP, model.ST, dimen=1, ordered=True, initialize=area_beta_gradients_init) model.Area_cu_beta_breakpoints = Set(model.HP, dimen=1, ordered=True, initialize=lambda model,i: area_cu_beta_breakpoints_init(model,i)) model.Area_cu_beta_exp = Set(model.HP, dimen=1, ordered=True, initialize=lambda model,i: map(lambda A: pow(A, model.Beta), model.Area_cu_beta_breakpoints[i])) model.Area_cu_beta_gradients = Set(model.HP, dimen=1, ordered=True, initialize=area_cu_beta_gradients_init) model.Area_hu_beta_breakpoints = Set(model.CP, dimen=1, ordered=True, initialize=lambda model,j: area_hu_beta_breakpoints_init(model,j)) model.Area_hu_beta_exp = Set(model.CP, dimen=1, ordered=True, initialize=lambda model,j: map(lambda A: pow(A, model.Beta), model.Area_hu_beta_breakpoints[j])) model.Area_hu_beta_gradients = Set(model.CP, dimen=1, ordered=True, initialize=area_hu_beta_gradients_init)