intpt - diag(s) zeros(n,m) diag(x) ]; if min(x)...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function [dx,dy,ds]=intpt(A,b,c,x,y,s) % Finds a affine scaling search direction for LP in SEF % Input: % A,b,c data for LP problem max c'x s.t. Ax=b, x>=0 % x,y,s current (int) feasible primal-dual pair, x>0, s>0 % [m,n]=size(A); M=[ zeros(n) A' -eye(n) A zeros(m) zeros(m,n)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: diag(s) zeros(n,m) diag(x) ]; if min(x) < 0 | min(s) < 0, disp('WARNING: x,y does not yield an interior point') end %RHS=[-(A'*y-s-c) ; -(A*x-b);-x.*s]; RHS=[zeros(n,1) ; zeros(m,1);-x.*s]; sol=M\RHS; dx=sol(1:n);dy=sol(n+1:n+m);ds=sol(n+m+1:end);...
View Full Document

Ask a homework question - tutors are online