apis/examples_python/cutstock_8py_source.html

가져오기sys

from감스가져오기GamsModifier, GamsWorkspace


GAMS_MASTER_MODEL ="""

'너비'를 설정합니다;


매개변수

w(i) '너비'

d(i) '수요';


스칼라

r '원시 너비';


$gdxIn csdata

$load i w d r

$gdxIn


$pmax를 설정하지 않은 경우 $set pmax 1000

설정

p '가능한 패턴' / 1*%pmax% /

pp(p) 'p의 동적 하위 집합';


매개변수

aip(i,p) 'p에서 성장하는 패턴의 너비 i 수';


변수

xp(p) '사용된 패턴'

z '객관적 변수';


정수 변수 xp;

xp.up(p) = sum(i, d(i));


수식

numpat

수요(i);


numpat.. z =e= sum(pp, xp(pp));


수요(i).. sum(pp, aip(i,pp)*xp(pp)) =g= d(i);


모델 마스터 / numpat, 수요 /;

"""


GAMS_SUB_MODEL ="""

'너비'를 설정합니다;


매개변수 w(i) '너비';

스칼라 r '원시 너비';


$gdxIn csdata

$load i w r

$gdxIn


매개변수

demdual(i) '마스터 수요 제약의 이중' / #i eps /;


변수

z

y(i) '새 패턴';

정수 변수 y;

y.up(i) = ceil(r/w(i));


수식

defobj

배낭;


defobj.. z =e= 1 - sum(i, demdual(i)*y(i));


배낭.. sum(i, w(i)*y(i)) =l= r;


모델 가격 / defobj, 배낭 /;

"""


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)


opt = ws.add_options()

cutstock_data = ws.add_database("csdata")

opt.all_model_types ="복합단지"

opt.optcr = 0.0# 최적성을 해결

max_pattern = 35


opt.defines["pmax"] = str(max_pattern)

opt.defines["solveMasterAs"] ="RMIP"


r = 100# 원시 너비

raw_width = cutstock_data.add_parameter("r", 0,"원시 너비")

raw_width.add_record().value = 100


d ="i1": 97,"i2": 610,"i3": 395,"i4": 211}

수요 = cutstock_data.add_parameter("d", 1,"수요")

너비 = cutstock_data.add_set("나", 1,"너비")

    fork, v에d.items():

widths.add_record(k)

demand.add_record(k).value = v


w ="i1": 47,"i2": 36,"i3": 31,"i4": 14}

폭 = cutstock_data.add_parameter("w", 1,"너비")

    fork, v에w.items():

width.add_record(k).value = v


cp_master = ws.add_checkpoint()

job_master_init = ws.add_job_from_string(GAMS_MASTER_MODEL)

job_master_init.run(opt, cp_master, 데이터베이스=cutstock_data)

job_master = ws.add_job_from_string(

        "execute_load 'csdata', aip, pp; %solveMasterAs%를 사용하여 마스터 최소 z를 해결합니다.",

cp_master,

    )


패턴 = cutstock_data.add_set("pp", 1,"패턴 색인")

pattern_data = cutstock_data.add_parameter("aip", 2,"패턴 데이터")


    # 초기 패턴: 패턴은 너비 i를 유지합니다.

pattern_count = 0

    fork, vinw.items():

pattern_count += 1

pattern_data.add_record(

(k, 패턴.add_record(str(pattern_count)).key(0))

).value = (int)(r / v)


cp_sub = ws.add_checkpoint()

job_sub = ws.add_job_from_string(GAMS_SUB_MODEL)

job_sub.run(opt, cp_sub, 데이터베이스=cutstock_data)

mi_sub = cp_sub.add_modelinstance()


    # 수정자 demdual 정의

demand_dual = mi_sub.sync_db.add_parameter(

        "뎀듀얼", 1,"마스터의 이중 요구"

    )

mi_sub.instantiate("mip를 사용하여 min z 가격 책정", GamsModifier(demand_dual), opt)


    # 새로운 패턴 찾기

    그동안 참:

job_master.run(opt, cp_master, 데이터베이스=cutstock_data)

        # 이중을 mi_sub.sync_db DB에 복사

demand_dual.clear()

        fordem에job_master.out_db["수요"]:

demand_dual.add_record(dem.key(0)).value = dem.marginal

mi_sub.solve()

        ifmi_sub.sync_db["z"].first_record().level < -0.00001:

            ifpattern_count == max_pattern:

인쇄(

f"패턴이 벗어났습니다. max_pattern(현재 max_pattern)을 늘리세요."

                )

                휴식

            그밖에:

인쇄(f"새 패턴! 값: mi_sub.sync_db['z'].first_record().level")

pattern_count += 1

s = 패턴.add_record(str(pattern_count))

                foryinmi_sub.sync_db["y"]:

                    ify.레벨 > 0.5:

pattern_data.add_record((y.key(0), s.key(0))).value = round(

y.level

                        )

        그밖에:

            휴식


    # 최종 MIP 해결

opt.defines["solveMasterAs"] ="MIP"

job_master.run(opt, 데이터베이스=cutstock_data)

인쇄(f"최적의 해결책: job_master.out_db['z'].first_record().level")

    forxp에job_master.out_db["xp"]:

        ifxp.레벨 > 0.5:

인쇄(f" 패턴 xp.key(0) xp.level 번:")

aip = job_master.out_db["아이피"].first_record((" ", xp.key(0)))

            그동안 참:

인쇄(f" aip.key(0): aip.value")

                if 아님aip.move_next():

                    휴식