apis/examples_python/transport7_8py_source.html

가져오기sys

에서게임가져오기GamsWorkspace, GamsModifier, UpdateAction, VarType


GAMS_MODEL ="""

설정

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 /

bmult '수요 승수' / 1 /;


매개변수 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).. sum(i, x(i,j)) =g= bmult*b(j);


모델 전송 /all/;

"""


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)


cp = ws.add_checkpoint()


    # GamsJob을 실행하여 GamsCheckpoint를 초기화합니다.

작업 = ws.add_job_from_string(GAMS_MODEL)

job.run(체크포인트=cp)


    # GamsModelInstance를 생성하고 다른 스칼라 bmult를 사용하여 여러 번 해결합니다.

mi = cp.add_modelinstance()

bmult = mi.sync_db.add_parameter("bmult", 0,"수요 승수")

opt = ws.add_options()

opt.all_model_types ="복합체"


    # GAMSModelInstance를 인스턴스화하고 모델 정의 및 GAMSModifier를 전달하여 bmult 변경 가능 선언

mi.instantiate("운송은 lp min z를 사용합니다", GamsModifier(bmult), opt)


bmult.add_record().value = 1.0

bmult_list = [0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3]


    forb에bmult_list:

bmult.first_record().value = b

mi.solve()

인쇄(f"시나리오 bmult=b:")

인쇄(f" 모델 상태: mi.model_status")

인쇄(f" 해결 상태: mi.solver_status")

인쇄(f" Obj: mi.sync_db['z'].first_record().level")


    # GamsModelInstance를 생성하고 차단된 네트워크의 단일 링크로 해결

mi = cp.add_modelinstance()

x = mi.sync_db.add_variable("x", 2, VarType.Positive)

xup = mi.sync_db.add_parameter("xup", 2,"x의 상한")


    # GamsModelInstance를 인스턴스화하고 모델 정의와 GamsModifier를 전달하여 x의 상한을 변경 가능으로 선언

mi.instantiate("운송은 lp min z를 사용합니다", GamsModifier(x, UpdateAction.Upper, xup))

mi.solve()


    fori에job.out_db["나"]:

        forjinjob.out_db["j"]:

xup.clear()

xup.add_record((i.key(0), j.key(0))).value = 0

mi.solve()

인쇄(f"시나리오 링크가 차단되었습니다: i.key(0) - j.key(0)")

인쇄(f" 모델 상태: mi.model_status")

인쇄(f" 솔베상태: mi.solver_status")

인쇄(f" 개체: mi.sync_db['z'].find_record().level")