def _calculate_ptdf(self): ptdf_r, ldf, ldf_c = tx_calc.calculate_ptdf_ldf(self._branches,self._buses,self.branches_keys,self.buses_keys,self._reference_bus,self._base_point,\ mapping_bus_to_idx=self._busname_to_index_map) self.PTDFM = ptdf_r self.LDF = ldf self.LDF_C = ldf_c ## protect the arrays using numpy self.PTDFM.flags.writeable = False self.LDF.flags.writeable = False self.LDF_C.flags.writeable = False
def _calculate_ptdf(self): ptdf_r, ldf, ldf_c = tx_calc.calculate_ptdf_ldf( self._branches, self._buses, self.branches_keys, self.buses_keys, self._reference_bus, self._base_point) ## protect the arrays using numpy ptdf_r.flags.writeable = False ldf.flags.writeable = False ldf_c.flags.writeable = False self.PTDFM = ptdf_r self.LDF = ldf self.LDF_C = ldf_c
def create_dicts_of_ptdf_losses(md,base_point=BasePointType.SOLUTION): branches = dict(md.elements(element_type='branch')) buses = dict(md.elements(element_type='bus')) branch_attrs = md.attributes(element_type='branch') bus_attrs = md.attributes(element_type='bus') reference_bus = md.data['system']['reference_bus'] ptdf_r, ldf, ldf_c = tx_calc.calculate_ptdf_ldf(branches,buses,branch_attrs['names'],bus_attrs['names'],reference_bus,base_point) phi_from, phi_to = tx_calc.calculate_phi_constant(branches,branch_attrs['names'],bus_attrs['names'],ApproximationType.PTDF_LOSSES) phi_loss_from, phi_loss_to = tx_calc.calculate_phi_loss_constant(branches,branch_attrs['names'],bus_attrs['names'],ApproximationType.PTDF_LOSSES) _len_branch = len(branch_attrs['names']) _mapping_branch = {i: branch_attrs['names'][i] for i in list(range(0,_len_branch))} _len_bus = len(bus_attrs['names']) _mapping_bus = {i: bus_attrs['names'][i] for i in list(range(0,_len_bus))} for idx,branch_name in _mapping_branch.items(): branch = md.data['elements']['branch'][branch_name] _row_ptdf_r = {bus_attrs['names'][i]: ptdf_r[idx,i] for i in list(range(0,_len_bus))} branch['ptdf_r'] = _row_ptdf_r _row_ldf = {bus_attrs['names'][i]: ldf[idx,i] for i in list(range(0,_len_bus))} branch['ldf'] = _row_ldf branch['ldf_c'] = ldf_c[idx] for idx, bus_name in _mapping_bus.items(): bus = md.data['elements']['bus'][bus_name] _row_phi_from = {branch_attrs['names'][i]: phi_from[idx,i] for i in list(range(0,_len_branch)) if phi_from[idx, i] != 0.} bus['phi_from'] = _row_phi_from _row_phi_to = {branch_attrs['names'][i]: phi_to[idx,i] for i in list(range(0,_len_branch)) if phi_to[idx, i] != 0.} bus['phi_to'] = _row_phi_to _row_phi_loss_from = {branch_attrs['names'][i]: phi_loss_from[idx,i] for i in list(range(0,_len_branch)) if phi_loss_from[idx, i] != 0.} bus['phi_loss_from'] = _row_phi_loss_from _row_phi_loss_to = {branch_attrs['names'][i]: phi_loss_to[idx,i] for i in list(range(0,_len_branch)) if phi_loss_to[idx, i] != 0.} bus['phi_loss_to'] = _row_phi_loss_to