Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  1. MATLAB のホームタブで、並列メニューを選択します。[クラスタの作成と管理]を選択します。

    Image RemovedImage Added
  2. クラスタプロファイルの追加 > MATLAB Job Scheduler をクリックします。

    Image RemovedImage Added
  3. プロファイル名をダブルクリックして、MATLAB Job Scheduler プロファイルの名前を変更します。
    今回はMJSProfile1

  4. プロファイルを選択し、ツールバーの[編集]をクリックして編集します。

    Image RemovedImage Added

  5. [完了]をクリックします。二枚目の画像は、設定後のMATLAB Job Schedulerクラスタ・プロファイルを示します。

    Image AddedImage RemovedImage Removed
  6. Image Added

    [規定に設定] を選択して、このプロファイルをデフォルトにします。

    Image Removed
  7. クラスタープロファイルの検証をクリックします。クラスタの検証が成功すると、以下のようになります。
    ※プロパティの隣の検証のタブに移動する。

    Image Removed

検証のエラー

...

  1. Image Added

  2. クラスタープロファイルの検証をクリックします。その後、ユーザー名とパスワードが聞かれます。
    クラスタの検証が成功すると、以下のようになります。

    Image AddedImage Added

演習

  • 準備:Coding files>Block Process On Large Imageファイルのinput_img.ngをC:\Users\ユーザー名\Documents\MATLABに移す。

  • localとサーバーの切り替えは、クラスタープロファイルマネージャーでの既定の設定を切り替えることで可能。

    Image Added
  • 並列処理ありとなしの切り替えは、コマンドにおける'UseParallel',trueを'UseParallel',falseに書き換える。

演習1:大きな画像に対するブロック処理

  1. 並列処理ありの場合

    Code Block
    tic % Start stopwatch timer
    % Input image
    input_img = "input_img.jpg"; % Add image path
    % Initialize Edge detection function
    fun = @(block_struct) edge(block_struct.data,"canny");
    % Covert source image from RGB to GRAY
    input_image= rgb2gray(imread(input_img));
    % Perform Parallel Block Process
    result = blockproc(input_image,[25 25],fun, ...
       'UseParallel',true);
    toc % Terminate stopwatch timer 
    % Show ouput image
    imshow(result)
    Image Added

  2. 並列処理なしの場合

    Code Block
    tic % Start stopwatch timer
    % Input image
    input_img = "input_img.jpg"; % Add image path
    % Initialize Edge detection function
    fun = @(block_struct) edge(block_struct.data,"canny");
    % Covert source image from RGB to GRAY
    input_image= rgb2gray(imread(input_img));
    % Perform Parallel Block Process
    result = blockproc(input_image,[25 25],fun, ...
       'UseParallel',false);
    toc % Terminate stopwatch timer 
    % Show ouput image
    imshow(result
    Image AddedImage Added

計算時間の比較(単位:秒)

並列処理なし

並列処理あり

windows

(Local)

windows

(Local)

4 workers

CentOS

(a9 Server)

12 workers

11.97

2.73

2.79

演習2:グローバルミニマムの探索

  1. 並列処理ありの場合

    Code Block
    tic % Start stopwatch timer
    % Consider a function with several local minima.
    fun = @(x) x.^2 + 4*sin(5*x);
    fplot(fun,[-10,10])
    rng default % For reproducibility
    opts = optimoptions(@fmincon,'Algorithm','sqp');
    problem = createOptimProblem('fmincon','objective',...
        fun,'x0',3,'lb',-5,'ub',5,'options',opts);
    ms = MultiStart('UseParallel', true);
    %To search for the global minimum, run MultiStart on 2000 instances of the problem using the fmincon 'sqp' algorithm.
    [x,f] = run(ms,problem,2000)
    toc % Terminate stopwatch timer
    Image Added

  2. 並列処理なしの場合

    Code Block
    tic % Start stopwatch timer
    % Consider a function with several local minima.
    fun = @(x) x.^2 + 4*sin(5*x);
    fplot(fun,[-10,10])
    rng default % For reproducibility
    opts = optimoptions(@fmincon,'Algorithm','sqp');
    problem = createOptimProblem('fmincon','objective',...
        fun,'x0',3,'lb',-5,'ub',5,'options',opts);
    ms = MultiStart('UseParallel', false);
    %To search for the global minimum, run MultiStart on 2000 instances of the problem using the fmincon 'sqp' algorithm.
    [x,f] = run(ms,problem,2000)
    toc % Terminate stopwatch timer

並列処理なし

並列処理あり

windows

(Local)

windows

(Local)

8 workers

CentOS

(a9 Server)

12 workers

6.20

1.77

1.52

演習3:SVM分類器の最適化

  1. 並列処理ありの場合

    Code Block
    tic % Start stopwatch timer 
    load ionosphere % Load the ionosphere data set.
    rng default
    % Find hyperparameters that minimize five-fold cross-validation loss by using automatic hyperparameter optimization. For reproducibility, set the random seed and use the 'expected-improvement-plus' acquisition function.
    Mdl = fitcsvm(X,Y,'OptimizeHyperparameters','auto', ...
    'HyperparameterOptimizationOptions',struct('UseParallel',true))
    toc % Terminate stopwatch timer
    Image AddedImage Added

  2. 並列処理なしの場合

    Code Block
    tic % Start stopwatch timer 
    load ionosphere % Load the ionosphere data set.
    rng default
    % Find hyperparameters that minimize five-fold cross-validation loss by using automatic hyperparameter optimization. For reproducibility, set the random seed and use the 'expected-improvement-plus' acquisition function.
    Mdl = fitcsvm(X,Y,'OptimizeHyperparameters','auto', ...
    'HyperparameterOptimizationOptions',struct('UseParallel',false))
    toc % Terminate stopwatch timer

並列処理なし

並列処理あり

windows

(Local)

windows

(Local)

8 workers

CentOS

(a9 Server)

12 workers

32.38

7.78

7.86

演習4:並列処理によるデータのクラスタリング

  1. 並列処理ありの場合

    Code Block
    Mu = bsxfun(@times,ones(20,300),(1:20)'); % Gaussian mixture mean
    rn300 = randn(300,300);
    Sigma = rn300'*rn300; % Symmetric and positive-definite covariance
    Mdl = gmdistribution(Mu,Sigma); % Define the Gaussian mixture distribution
    
    rng(1); % For reproducibility
    X = random(Mdl,10000);
    % Specify the options for parallel computing.
    stream = RandStream('mlfg6331_64');  % Random number stream
    options = statset('UseParallel',1,'UseSubstreams',1,...
        'Streams',stream);
    % Cluster the data using k-means clustering. Specify that there are k = 200 clusters in the data and increase the number of iterations. Typically, the objective function contains local minima. Specify 10 replicates to help find a lower, local minimum.
    
    tic; % Start stopwatch timer
    [idx,C,sumd,D] = kmeans(X,200,'Options',options,'MaxIter',10000,...
        'Display','final','Replicates',10);
    toc % Terminate stopwatch timer
    Image Added

  2. 並列処理なしの場合

    Code Block
    Mu = bsxfun(@times,ones(20,300),(1:20)'); % Gaussian mixture mean
    rn300 = randn(300,300);
    Sigma = rn300'*rn300; % Symmetric and positive-definite covariance
    Mdl = gmdistribution(Mu,Sigma); % Define the Gaussian mixture distribution
    
    rng(1); % For reproducibility
    X = random(Mdl,10000);
    % Specify the options for parallel computing.
    stream = RandStream('mlfg6331_64');  % Random number stream
    options = statset('UseParallel',false,'UseSubstreams',1,...
        'Streams',stream);
    % Cluster the data using k-means clustering. Specify that there are k = 200 clusters in the data and increase the number of iterations. Typically, the objective function contains local minima. Specify 10 replicates to help find a lower, local minimum.
    
    tic; % Start stopwatch timer
    [idx,C,sumd,D] = kmeans(X,200,'Options',options,'MaxIter',10000,...
        'Display','final','Replicates',10);
    toc % Terminate stopwatch timer

並列処理なし

並列処理あり

windows

(Local)

windows

(Local)

8 workers

CentOS

(a9 Server)

12 workers

8.45

9.90

11.00

並列処理ありと並列処理なしの違い

  • 並列処理あり

    • コマンドにおいて'UseParallel',true

    • 処理中の際に左下のアイコンが緑色になる。

      Image Added
  • 並列処理なし

    • コマンドにおいて'UseParallel',false

    • 処理中の際に左下のアイコンが青色になる。

      Image Added