apis/examples_python/transport8a_8py_source.html

가져오기os

가져오기sys

from게임가져오기GamsWorkspace, GamsModifier

from다중 처리가져오기잠금, 프로세스, 대기열


GAMS_MODEL ="""

세트

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

j 마켓 / 뉴욕, 시카고, 토피카 / ;


매개변수

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 ;


수식

비용 정의 목적 함수

공급(i) 공장 i의 공급 제한을 준수합니다.

수요(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/ ;

"""


defscen_solve(cp_file, bmult_queue, queue_lock, io_lock):

ws = GamsWorkspace()

cp = ws.add_checkpoint(cp_file)

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 값을 동적으로 가져옵니다.

queue_lock.acquire()

        ifbmult_queue.empty():

queue_lock.release()

            반환

b = bmult_queue.get()

queue_lock.release()

bmult.first_record().value = b

mi.solve()


        # 보고서 정보가 엉망이 되는 것을 방지하기 위해 출력을 중요한 섹션으로 만들어야 합니다.

io_lock.acquire()

인쇄("시나리오 bmult="+ 문자열(b) +":")

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

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

인쇄(" 개체: "+ str(mi.sync_db.get_variable("z").find_record().level))

io_lock.release()


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)

cp_file = os.path.join(ws.working_directory, cp.name) +".g00"


bmult_queue = 대기열()

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

    forb에b다중:

bmult_queue.put(b)


    # 여러 모델 인스턴스를 병렬로 해결

io_lock = 잠금()

queue_lock = 잠금()


nr_workers = 4

프로세스 =

    fori에범위(nr_workers):

프로세스[i] = 프로세스(

target=scen_solve, args=(cp_file, bmult_queue, queue_lock, io_lock)

        )

프로세스[i].start()

    foriin범위(nr_workers):

프로세스[i].join()