apis/examples_python/transport3_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) '거리(천 마일)';


스칼라 f '1,000마일당 케이스당 달러 단위의 화물';


$if 설정되지 않음 gdxincname $abort '제공된 데이터 파일에 대한 포함 파일 이름 없음'

$gdxIn %gdxincname%

$load i j a b d f

$gdxIn


매개변수 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)


작업 = ws.add_job_from_string(GAMS_DATA)

job.run()

job.out_db.export(os.path.join(ws.working_directory,"tdata.gdx"))

작업 = ws.add_job_from_string(GAMS_MODEL)


opt = ws.add_options()

opt.defines["gdxincname"] ="tdata"

opt.all_model_types ="익스프레스"

job.run(옵션)

    forrec에job.out_db["x"]:

인쇄(

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

        )


job_data = ws.add_job_from_string(GAMS_DATA)

job_model = ws.add_job_from_string(GAMS_MODEL)

job_data.run()

opt.defines["gdxincname"] = job_data.out_db.name

job_model.run(opt, 데이터베이스=job_data.out_db)

    forrec에job_model.out_db["x"]:

인쇄(

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

        )