function [bus, gen, branch] = opfsoln(baseMVA, bus0, gen0, branch0, ...
								Ybus, Yf, Yt, V, Sg, lambda, ref, pv, pq, mpopt)
%OTOPFSOL   Updates bus, gen, branch data structures to match opf soln.
%   [bus, gen, branch] = opfsoln(baseMVA, bus0, gen0, branch0, ...
%                              Ybus, Yf, Yt, V, Sg, lambda, ref, pv, pq, mpopt)

%   MATPOWER Version 2.0
%   by Ray Zimmerman, PSERC Cornell    12/24/97
%   Copyright (c) 1996, 1997 by Power System Engineering Research Center (PSERC)
%   See http://www.pserc.cornell.edu/ for more info.

%% constants
j = sqrt(-1);
nl = size(branch0, 1);		%% number of lines

%% define named indices into data matrices
[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
	VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
	GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;
[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, ...
	RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST] = idx_brch;

%% initialize return values
bus		= bus0;
gen		= gen0;
branch	= branch0;

%%----- update bus voltages -----
bus(:, VM) = abs(V);
bus(:, VA) = angle(V) * 180 / pi;

%%----- update Pg and Qg for all gens -----
%% generator info
on = find(gen(:, GEN_STATUS));				%% which generators are on?
gbus = gen(on, GEN_BUS);					%% what buses are they at?
refgen = find(gen(:, GEN_BUS) == ref);		%% which is the reference gen?

%% copy back Pg & Vg
gen(:, PG) = zeros(size(gen, 1), 1);
gen(on, PG) = real(Sg) * baseMVA;
gen(on, VG) = bus(gbus, VM);

%% compute Qg if not passed in with Sg
if ~any(imag(Sg))
	%% This is slow in Matlab 5 ...
	% Sg = V(gbus) .* conj(Ybus(gbus, :) * V);
	%% ... so we do this instead ...
	temp = Ybus.';
	Sg = V(gbus) .* conj(temp(:, gbus).' * V);
	Qg = imag(Sg) + bus(gbus, QD) / baseMVA;	%% inj Q + local Qd
else
	Qg = imag(Sg);
end

%% update Qg for all generators
gen(:, QG) = zeros(size(gen, 1), 1);		%% zero out all Qg
gen(on, QG) = Qg * baseMVA;					%% except for on-line generators

%%----- update/compute branch power flows -----
Sf = V(branch(:, F_BUS)) .* conj(Yf * V) * baseMVA;	%% complex power at "from" bus
St = V(branch(:, T_BUS)) .* conj(Yt * V) * baseMVA;	%% complex power injected at "to" bus
branch(:, [PF, QF, PT, QT]) = [real(Sf) imag(Sf) real(St) imag(St)];

%%----- update lambda's and mu's -----
%% sizes of things
nb = size(bus, 1);
nl = size(branch, 1);
npv	= length(pv);
npq	= length(pq);
ng = npv + 1;			%% number of generators that are turned on

%% initialize with all zeros
bus(:, [LAM_P, LAM_Q, MU_VMIN, MU_VMAX]) = zeros(nb, 4);
gen(:, [MU_PMIN, MU_PMAX, MU_QMIN, MU_QMAX]) = zeros(size(gen, 1), 4);
branch(:, [MU_SF, MU_ST]) = zeros(nl, 2);

%% set up indexing for lambda
i1 = 1;			i2 = nb;			%% i1:i2 - P mismatch, all buses
i3 = i2 + 1;	i4 = i2 + nb;		%% i3:i4 - Q mismatch, all buses
i5 = i4 + 1;	i6 = i4 + nb;		%% i5:i6 - Vmin, all buses
i7 = i6 + 1;	i8 = i6 + nb;		%% i7:i8 - Vmax, all buses
i9 = i8 + 1;	i10 = i8 + ng;		%% i9:i10 - Pmin, gen buses
i11 = i10 + 1;	i12 = i10 + ng;		%% i11:i12 - Pmax, gen buses
i13 = i12 + 1;	i14 = i12 + ng;		%% i13:i14 - Qmin, gen buses
i15 = i14 + 1;	i16 = i14 + ng;		%% i15:i16 - Qmax, gen buses
i17 = i16 + 1;	i18 = i16 + nl;		%% i17:i18 - |Sf| line limit
i19 = i18 + 1;	i20 = i18 + nl;		%% i19:i20 - |St| line limit

%% copy multipiers to bus, gen, branch
bus(:, LAM_P)		= lambda(i1:i2) / baseMVA;
bus(:, LAM_Q)		= lambda(i3:i4) / baseMVA;
bus(:, MU_VMIN)		= lambda(i5:i6);
bus(:, MU_VMAX)		= lambda(i7:i8);

gen(on, MU_PMIN)	= lambda(i9:i10) / baseMVA;
gen(on, MU_PMAX)	= lambda(i11:i12) / baseMVA;
gen(on, MU_QMIN)	= lambda(i13:i14) / baseMVA;
gen(on, MU_QMAX)	= lambda(i15:i16) / baseMVA;

branch(:, MU_SF)	= lambda(i17:i18) / baseMVA;
branch(:, MU_ST)	= lambda(i19:i20) / baseMVA;

return;