% /**************************************************************** % Copyright (C) 1997 Lucent Technologies % All Rights Reserved % % Permission to use, copy, modify, and distribute this software and % its documentation for any purpose and without fee is hereby % granted, provided that the above copyright notice appear in all % copies and that both that the copyright notice and this % permission notice and warranty disclaimer appear in supporting % documentation, and that the name of Lucent or any of its entities % not be used in advertising or publicity pertaining to % distribution of the software without specific, written prior % permission. % % LUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, % INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS. % IN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY % SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES % WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER % IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, % ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF % THIS SOFTWARE. % ****************************************************************/ function [g,A,B] = evalg(x) global clb cub clu ceq mp0 n xlc xuc [g,cjac] = amplfunc(x, 1); A = [ cjac(ceq,:) ]; B = [ cjac(clb,:); -cjac(cub,:); kron(cjac(clu,:),[1;-1]);... zeros(size(xlc,2)+size(xuc,2), n) ]; j = 0; for i = xlc, j = j + 1; B(j,i) = 1; end for i = xuc, j = j + 1; B(j,i) = -1; end .