apis/examples_matlab_control/transport3_8m_source.html

function transport3(varargin)


    % create a directory

    workingDirectory = fullfile(pwd, 'transport3');

    mkdir(workingDirectory);


    % create a workspace

    wsInfo = gams.control.WorkspaceInfo();

    wsInfo.workingDirectory = workingDirectory;

    if nargin > 0

        wsInfo.systemDirectory = varargin{1};

    end

    ws = gams.control.Workspace(wsInfo);


    data = {

    'Sets                                                             '

    '  i   canning plants   / seattle, san-diego /                    '

    '  j   markets          / new-york, chicago, topeka / ;           '

    'Parameters                                                       '

    '                                                                 '

    '  a(i)  capacity of plant i in cases                             '

    '                     /    seattle     350                        '

    '                          san-diego   600  /                     '

    '                                                                 '

    '  b(j)  demand at market j in cases                              '

    '                     /    new-york    325                        '

    '                          chicago     300                        '

    '                          topeka      275  / ;                   '

    '                                                                 '

    'Table d(i,j)  distance in thousands of miles                     '

    '               new-york       chicago      topeka                '

    '  seattle          2.5           1.7          1.8                '

    '  san-diego        2.5           1.8          1.4  ;             '

    '                                                                 '

    'Scalar f  freight in dollars per case per thousand miles  /90/ ;'};

    data = sprintf('%s\n', data{:});


    model = {

    'Sets                                                                       '

    '      i   canning plants                                                   '

    '      j   markets                                                          '

    '                                                                           '

    ' Parameters                                                                '

    '      a(i)   capacity of plant i in cases                                  '

    '      b(j)   demand at market j in cases                                   '

    '      d(i,j) distance in thousands of miles                                '

    ' Scalar f  freight in dollars per case per thousand miles;                 '

    '                                                                           '

    '$if not set gdxincname $abort ''no include file name for data file provided'''

    '$gdxin %gdxincname%                                                        '

    '$load i j a b d f                                                          '

    '$gdxin                                                                     '

    '                                                                           '

    ' Parameter c(i,j)  transport cost in thousands of dollars per case ;       '

    '                                                                           '

    '            c(i,j) = f * d(i,j) / 1000 ;                                   '

    '                                                                           '

    ' Variables                                                                 '

    '       x(i,j)  shipment quantities in cases                                '

    '       z       total transportation costs in thousands of dollars ;        '

    '                                                                           '

    ' Positive Variable x ;                                                     '

    '                                                                           '

    ' Equations                                                                 '

    '                                                                           '

    '      cost        define objective function                                '

    '      supply(i)   observe supply limit at plant i                          '

    '      demand(j)   satisfy demand at market j ;                             '

    '                                                                           '

    '  cost ..        z  =e=  sum((i,j), c(i,j)*x(i,j)) ;                       '

    '                                                                           '

    '  supply(i) ..   sum(j, x(i,j))  =l=  a(i) ;                               '

    '                                                                           '

    '  demand(j) ..   sum(i, x(i,j))  =g=  b(j) ;                               '

    '                                                                           '

    ' Model transport /all/ ;                                                   '

    '                                                                           '

    ' Solve transport using lp minimizing z ;                                   '

    '                                                                           '

    ' Display x.l, x.m ;                                                        '};

    model = sprintf('%s\n', model{:});


    % create and run a job from a data file, then explicitly export to a GDX file

    t3 = ws.addJobFromString(data);

    t3.run();

    t3.outDB.export(fullfile(ws.workingDirectory, 'tdata.gdx'));


    % clear database

    t3.outDB.dispose();


    % run a job using an instance of Options that defines the data include file

    t3 = ws.addJobFromString(model);

    opt = ws.addOptions();

    opt.defines('gdxincname', 'tdata');

    t3.run(opt);


    % retrieve Variable 'x' from Job's output databases

    for x = t3.outDB.getVariable('x').records

        fprintf('x(%s,%s): level=%g marginal=%g\n', x{1}.keys{:}, x{1}.level, x{1}.marginal);

    end


    % clear database

    t3.outDB.dispose();


    % similar to the previous run but without exporting database into a file

    t3a = ws.addJobFromString(data);

    t3b = ws.addJobFromString(model);

    opt = ws.addOptions();

    t3a.run();

    opt.defines('gdxincname', t3a.outDB.name);

    t3b.run(opt, t3a.outDB);


    % retrieve Variable 'x' from Job's output databases

    for x = t3b.outDB.getVariable('x').records

        fprintf('x(%s,%s): level=%g marginal=%g\n', x{1}.keys{:}, x{1}.level, x{1}.marginal);

    end


    % clear option and database

    opt.dispose();

    t3a.outDB.dispose();

    t3b.outDB.dispose();


    % remove working directory

    rmdir(ws.workingDirectory, 's');


end