apis/examples_python/transport12_8py_source.html

가져오기sys

from게임가져오기 (

GamsWorkspace,

GamsModifier,

GamsEquation,

슬롯 사이트변수,

GamsException,

업데이트액션,

VarType,

SymbolUpdateType,

)


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/;

"""


#은 인스턴스화되지 않은 GamsModelInstance로 호출되어야 합니다.

defguss_call(scenario_dict, mi,solv_statement, opt=없음, mi_opt=없음, 출력=없음):

    ifscenario_dict.dimension != 3:

        올림GamsException("시나리오 사전은 3차원이어야 합니다")


scen_name = 시나리오_dict.first_record((" ", "시나리오", " ")).key(0)

scen_symbol = 시나리오_dict.database[scen_name]

modifier_list = []

    forrec에scenario_dict:

        ifrec.key(1).lower() =="시나리오":

            계속

        ifrec.key(1).lower() =="매개변수":

modifier_dim = (

scenario_dict.database[rec.key(2)].dimension - scen_symbol.dimension

            )

            ifmodifier_dim < 0:

                올림GamsException("크기"+ Rec.key(2) +"너무 작습니다")

modifier_list.append(

                (

GamsModifier(

mi.sync_db.add_parameter(rec.key(0), modifier_dim,"")

                    ),

scenario_dict.database[rec.key(2)],

                )

            )

        엘리프rec.key(1).lower()in["낮은", "상단", "고정됨"]:

modifier_dim = (

scenario_dict.database[rec.key(2)].dimension - scen_symbol.dimension

            )

            ifmodifier_dim < 0:

                올림GamsException("크기"+ Rec.key(2) +"너무 작습니다")

modifier_var =없음

            시도해 보세요:

modifier_var = 시나리오_dict.database[rec.key(0)]

            제외:

modifier_var = mi.sync_db.add_variable(

rec.key(0), modifier_dim, VarType.Free,""

                )

            ifrec.key(1).lower() =="낮은":

update_action = UpdateAction.Lower

            엘리프rec.key(1).lower() =="상단":

update_action = UpdateAction.Upper

            그밖에:# 수정됨

update_action = UpdateAction.Fixed

modifier_list.append(

                (

GamsModifier(

modifier_var,

update_action,

mi.sync_db.add_parameter(rec.key(2), modifier_dim,""),

                    ),

scenario_dict.database[rec.key(2)],

                )

            )

        엘리프rec.key(1).lower()in["레벨", "한계"]:

            # GamsDatabase에 매개변수가 있는지 확인하고 없으면 예외가 발생합니다.

scenario_dict.database[rec.key(2)]

        그밖에:

            올림GamsException("UpdateAction을 처리할 수 없습니다"+ Rec.key(1))


ml = [m[0]formin수정자_목록]

mi.instantiate(solve_statement, ml, opt)

out_list = []

    fors에scen_symbol:

        form, p_scen에수정자_목록:

p = m.gams_symbolifm.data_symbolis 없음 그밖에m.data_symbol

p.clear()# SymbolUpdateType=BaseCase를 구현함


rec =없음

필터 = [""] * p_scen.dimension

            foriin범위(scen_symbol.dimension):

필터[i] = s.key(i)

            foriin범위(scen_symbol.dimension, p_scen.dimension):

필터[i] =" "

            시도해 보세요:

rec = p_scen.first_record(필터)

            제외:

                계속


            그동안 참:

my_keys = [

rec.key(scen_symbol.dimension + i)foriin범위(p.dimension)

]

p.add_record(my_keys).value = Rec.value

                if 아님rec.move_next():

                    휴식


mi.solve(SymbolUpdateType.BaseCase, 출력, mi_opt)

        if 아님out_list:

            forrec에scenario_dict:

                ifrec.key(1).lower()in["레벨", "한계"]:

out_list.append(

                        (

mi.sync_db[rec.key(0)],

scenario_dict.database[rec.key(2)],

rec.key(1).lower(),

                        )

                    )


        fortupinout_list:

my_keys = [""] * (scen_symbol.dimension + len(tup[0].first_record().keys))

            fori에범위(scen_symbol.dimension):

my_keys[i] = s.key(i)


            ifisinstance(tup[0], (GamsVariable, GamsEquation)):

                if트업[2]in["레벨", "한계"]:

                    forrecintup[0]:

                        fori에범위(len(rec.keys)):

my_keys[scen_symbol.dimension + i] = s.key(i)

tup[1].add_record(my_keys).value = (

rec.leveliftup[2] =="레벨" 그밖에rec.marginal

                        )


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_list = [0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3]

db = ws.add_database()

슬롯 사이트 = db.add_set("장면", 1,"")

bmult = db.add_parameter_dc("b다중 목록", [장면])

zscen = db.add_parameter_dc("zscen", [장면])


    fori, bin열거(bmult_list):

bmult.add_record(f"si").값 = b

슬롯 사이트add_record(f"si")


scenario_dict = db.add_set("scenario_dict", 3,"")

scenario_dict.add_record((슬롯 사이트name,"시나리오", ""))

scenario_dict.add_record(("bmult", "매개변수", bmult.name))

scenario_dict.add_record(("z", "레벨", zscen.name))


guss_call(scenario_dict, mi,"운송은 lp min z를 사용합니다")


    forrec에db[zscen.name]:

인쇄(f"rec.key(0) obj: rec.value")


mi2 = cp.add_modelinstance()

db2 = ws.add_database()


scen2 = db2.add_set("장면", 1,"")

zscen2 = db2.add_parameter_dc("zscen", [scen2])

xup = db2.add_parameter("xup", 3,"")


    forjin범위(4):

        forrec_iinjob.out_db["나"]:

            forrec_j에job.out_db["j"]:

xup.add_record((f"sj", rec_i.key(0), rec_j.key(0))).value = j + 1

scen2.add_record(f"sj")


scenario_dict2 = db2.add_set("scenario_dict", 3,"")

scenario_dict2.add_record((scen2.name,"시나리오", ""))

scenario_dict2.add_record(("x", "낮은", xup.name))

scenario_dict2.add_record(("z", "레벨", zscen2.name))


guss_call(scenario_dict2, mi2,"운송은 lp min z를 사용합니다", 출력=sys.stdout)


    forrecindb2[zscen2.name]:

인쇄(f"rec.key(0) obj: rec.value")