I2E_DATA Converts data from internal to external bus numbering.
VAL = I2E_DATA(MPC, VAL, OLDVAL, ORDERING)
VAL = I2E_DATA(MPC, VAL, OLDVAL, ORDERING, DIM)
For a case struct using internal indexing, this function can be
used to convert other data structures as well by passing in 3 or 4
extra parameters in addition to the case struct. If the value passed
in the 2nd argument (VAL) is a column vector, it will be converted
according to the ordering specified by the 4th argument (ORDERING,
described below). If VAL is an n-dimensional matrix, then the
optional 5th argument (DIM, default = 1) can be used to specify
which dimension to reorder. The 3rd argument (OLDVAL) is used to
initialize the return value before converting VAL to external
indexing. In particular, any data corresponding to off-line gens
or branches or isolated buses or any connected gens or branches
will be taken from OLDVAL, with VAL supplying the rest of the
returned data.
The ORDERING argument is used to indicate whether the data
corresponds to bus-, gen- or branch-ordered data. It can be one
of the following three strings: 'bus', 'gen' or 'branch'. For
data structures with multiple blocks of data, ordered by bus,
gen or branch, they can be converted with a single call by
specifying ORDERING as a cell array of strings.
Any extra elements, rows, columns, etc. beyond those indicated
in ORDERING, are not disturbed.
Examples:
A_ext = i2e_data(mpc, A_int, A_orig, {'bus','bus','gen','gen'}, 2);
Converts an A matrix for user-supplied OPF constraints from
internal to external ordering, where the columns of the A
matrix correspond to bus voltage angles, then voltage
magnitudes, then generator real power injections and finally
generator reactive power injections.
gencost_ext = i2e_data(mpc, gencost_int, gencost_orig, {'gen','gen'}, 1);
Converts a GENCOST matrix that has both real and reactive power
costs (in rows 1--ng and ng+1--2*ng, respectively).
See also E2I_DATA, I2E_FIELD, INT2EXT.0001 function newval = i2e_data(mpc, val, oldval, ordering, dim) 0002 %I2E_DATA Converts data from internal to external bus numbering. 0003 % 0004 % VAL = I2E_DATA(MPC, VAL, OLDVAL, ORDERING) 0005 % VAL = I2E_DATA(MPC, VAL, OLDVAL, ORDERING, DIM) 0006 % 0007 % For a case struct using internal indexing, this function can be 0008 % used to convert other data structures as well by passing in 3 or 4 0009 % extra parameters in addition to the case struct. If the value passed 0010 % in the 2nd argument (VAL) is a column vector, it will be converted 0011 % according to the ordering specified by the 4th argument (ORDERING, 0012 % described below). If VAL is an n-dimensional matrix, then the 0013 % optional 5th argument (DIM, default = 1) can be used to specify 0014 % which dimension to reorder. The 3rd argument (OLDVAL) is used to 0015 % initialize the return value before converting VAL to external 0016 % indexing. In particular, any data corresponding to off-line gens 0017 % or branches or isolated buses or any connected gens or branches 0018 % will be taken from OLDVAL, with VAL supplying the rest of the 0019 % returned data. 0020 % 0021 % The ORDERING argument is used to indicate whether the data 0022 % corresponds to bus-, gen- or branch-ordered data. It can be one 0023 % of the following three strings: 'bus', 'gen' or 'branch'. For 0024 % data structures with multiple blocks of data, ordered by bus, 0025 % gen or branch, they can be converted with a single call by 0026 % specifying ORDERING as a cell array of strings. 0027 % 0028 % Any extra elements, rows, columns, etc. beyond those indicated 0029 % in ORDERING, are not disturbed. 0030 % 0031 % Examples: 0032 % A_ext = i2e_data(mpc, A_int, A_orig, {'bus','bus','gen','gen'}, 2); 0033 % 0034 % Converts an A matrix for user-supplied OPF constraints from 0035 % internal to external ordering, where the columns of the A 0036 % matrix correspond to bus voltage angles, then voltage 0037 % magnitudes, then generator real power injections and finally 0038 % generator reactive power injections. 0039 % 0040 % gencost_ext = i2e_data(mpc, gencost_int, gencost_orig, {'gen','gen'}, 1); 0041 % 0042 % Converts a GENCOST matrix that has both real and reactive power 0043 % costs (in rows 1--ng and ng+1--2*ng, respectively). 0044 % 0045 % See also E2I_DATA, I2E_FIELD, INT2EXT. 0046 0047 % MATPOWER 0048 % $Id: i2e_data.m,v 1.1 2011/11/09 21:32:13 cvs Exp $ 0049 % by Ray Zimmerman, PSERC Cornell 0050 % Copyright (c) 2009-2011 by Power System Engineering Research Center (PSERC) 0051 % 0052 % This file is part of MATPOWER. 0053 % See http://www.pserc.cornell.edu/matpower/ for more info. 0054 % 0055 % MATPOWER is free software: you can redistribute it and/or modify 0056 % it under the terms of the GNU General Public License as published 0057 % by the Free Software Foundation, either version 3 of the License, 0058 % or (at your option) any later version. 0059 % 0060 % MATPOWER is distributed in the hope that it will be useful, 0061 % but WITHOUT ANY WARRANTY; without even the implied warranty of 0062 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 0063 % GNU General Public License for more details. 0064 % 0065 % You should have received a copy of the GNU General Public License 0066 % along with MATPOWER. If not, see <http://www.gnu.org/licenses/>. 0067 % 0068 % Additional permission under GNU GPL version 3 section 7 0069 % 0070 % If you modify MATPOWER, or any covered work, to interface with 0071 % other modules (such as MATLAB code and MEX-files) available in a 0072 % MATLAB(R) or comparable environment containing parts covered 0073 % under other licensing terms, the licensors of MATPOWER grant 0074 % you additional permission to convey the resulting work. 0075 0076 if ~isfield(mpc, 'order') 0077 error('i2e_data: mpc does not have the ''order'' field required for conversion back to external numbering.'); 0078 end 0079 o = mpc.order; 0080 if o.state ~= 'i' 0081 error('i2e_data: mpc does not appear to be in internal order'); 0082 end 0083 if nargin < 5 0084 dim = 1; 0085 end 0086 if ischar(ordering) %% single set 0087 if strcmp(ordering, 'gen') 0088 v = get_reorder(val, o.(ordering).i2e, dim); 0089 else 0090 v = val; 0091 end 0092 newval = set_reorder(oldval, v, o.(ordering).status.on, dim); 0093 else %% multiple sets 0094 be = 0; %% base, external indexing 0095 bi = 0; %% base, internal indexing 0096 for k = 1:length(ordering) 0097 ne = size(o.ext.(ordering{k}), 1); 0098 ni = size(mpc.(ordering{k}), 1); 0099 v = get_reorder(val, bi+(1:ni), dim); 0100 oldv = get_reorder(oldval, be+(1:ne), dim); 0101 new_v{k} = i2e_data(mpc, v, oldv, ordering{k}, dim); 0102 be = be + ne; 0103 bi = bi + ni; 0104 end 0105 ni = size(val, dim); 0106 if ni > bi %% the rest 0107 v = get_reorder(val, bi+1:ni, dim); 0108 new_v{length(new_v)+1} = v; 0109 end 0110 newval = cat(dim, new_v{:}); 0111 end