apis/examples_python/transport2_8py_source.html

가져오기os

가져오기sys

from게임가져오기GamsWorkspace


GAMS_DATA ="""

설정

i '통조림 공장' / 시애틀, 샌디에이고 /

j 'markets' / 뉴욕, 시카고, 토피카 /;


매개변수

a(i) '경우에 따른 식물 i의 용량'

/시애틀 350

샌디에고 600 /


b(j) '경우에 따라 시장 j의 수요'

/뉴욕 325

시카고 300

토피카 275 /;


테이블 d(i,j) '거리(천 마일)'

뉴욕 시카고 토피카

시애틀 2.5 1.7 1.8

산디에고 2.5 1.8 1.4;


스칼라 f '1,000마일당 케이스당 운임(달러)' / 90 /;

"""


GAMS_MODEL ="""

설정

i '통조림 식물'

j '시장';


매개변수

a(i) '경우에 따른 식물 i의 용량'

b(j) '경우에 따라 시장 j의 수요'

d(i,j) '거리(천 마일)';


$if not set incname $abort '제공된 데이터 파일에 대한 포함 파일 이름이 없습니다'

$include %incname%


매개변수 c(i,j) '케이스당 수천 달러의 운송 비용';

c(i,j) = f*d(i,j)/1000;


변수

x(i,j) '케이스 내 배송 수량'

z '총 운송 비용(수천 달러)';


양수 변수 x;


수식

비용 '목적 함수 정의'

supply(i) '공장 i의 공급 제한을 준수하세요'

demand(j) '시장 j의 수요 충족';


비용.. z =e= sum((i,j), c(i,j)*x(i,j));


공급(i).. 합계(j, x(i,j)) =l= a(i);


수요(j).. 합계(i, x(i,j)) =g= b(j);


모델 전송 /all/;


z를 최소화하는 lp를 사용하여 전송을 해결합니다.


x.l, x.m 표시;

"""


if__name__ =="__main__":

sys_dir = sys.argv[1]iflen(sys.argv) > 1그밖에 없음

work_dir = sys.argv[2]iflen(sys.argv) > 2그밖에 없음

ws = GamsWorkspace(system_directory=sys_dir,working_directory=work_dir)


    함께오픈(os.path.join(ws.working_directory,"tdata.gms"), "w") as파일:

파일.쓰기(GAMS_DATA)


작업 = ws.add_job_from_string(GAMS_MODEL)

opt = ws.add_options()

opt.defines["inc이름"] ="tdata"

job.run(옵션)

    forrec에job.out_db["x"]:

인쇄(

f"x(rec.key(0),rec.key(1)): 레벨=rec.level 한계=rec.marginal"

        )