parfeval
을 사용하여 파라미터 스윕 중에 플로팅하기
이 예제에서는 parfeval
을 사용하여 병렬 파라미터 스윕을 수행하고 DataQueue
객체를 사용하여 계산 중에 결과를 다시 보내는 방법을 보여줍니다.
parfeval
은 MATLAB을 차단하지 않으므로 계산이 수행되는 동안 계속해서 작업할 수 있습니다.
이 예제는 로렌츠 연립상미분방정식에서 파라미터 및 에 대한 파라미터 스윕을 수행하고 이 연립방정식의 혼돈적인 성질을 보여줍니다.
파라미터 그리드 만들기
파라미터 스윕에서 탐색할 파라미터의 범위를 정의합니다.
gridSize = 40; sigma = linspace(5, 45, gridSize); rho = linspace(50, 100, gridSize); beta = 8/3;
meshgrid
함수를 사용하여 파라미터의 2차원 그리드를 만듭니다.
[rho,sigma] = meshgrid(rho,sigma);
figure 객체를 만들고 'Visible'
을 true
로 설정하면 라이브 스크립트 밖의 새 창에서 열립니다. 파라미터 스윕의 결과를 시각화하려면 곡면 플롯을 만드십시오. NaN
을 사용하여 곡면의 Z
구성요소를 초기화하면 빈 플롯이 만들어집니다.
figure('Visible',true); surface = surf(rho,sigma,NaN(size(sigma))); xlabel('\rho','Interpreter','Tex') ylabel('\sigma','Interpreter','Tex')
병렬 환경 설정하기
parpool
함수를 사용하여 병렬 워커의 풀을 만듭니다.
parpool;
Starting parallel pool (parpool) using the 'Processes' profile ... Connected to the parallel pool (number of workers: 6).
워커에서 데이터를 보내려면 DataQueue
객체를 만드십시오. 워커가 afterEach
함수를 사용하여 데이터를 보낼 때마다 곡면 플롯을 업데이트하는 함수를 설정합니다. updatePlot
함수는 예제의 마지막 부분에 정의된 지원 함수입니다.
Q = parallel.pool.DataQueue; afterEach(Q,@(data) updatePlot(surface,data));
병렬 파라미터 스윕 수행하기
파라미터를 정의한 후 병렬 파라미터 스윕을 수행할 수 있습니다.
작업량을 분산하면 parfeval
은 더 효율적으로 작동합니다. 작업량을 분산하려면 탐색하고자 하는 파라미터를 파티션 단위로 그룹화하십시오. 이 예제에서는 콜론 연산자(:
)를 사용하여 균일한 크기(step
)의 파티션으로 분할합니다. 결과로 생성되는 배열 partitions
에는 파티션 경계가 담겨 있습니다. 이때 마지막 파티션의 끝점을 추가해야 함을 명심하십시오.
step = 100; partitions = [1:step:numel(sigma), numel(sigma)+1]
partitions = 1×17
1 101 201 301 401 501 601 701 801 901 1001 1101 1201 1301 1401 1501 1601
최상의 성능을 구현하려면 다음과 같은 파티션 크기로 분할하십시오.
파티션을 나누는데 드는 오버헤드에 비해 계산 시간이 더 길 정도로 충분히 큰 크기.
모든 워커에서 계산을 수행할 수 있을 정도로 충분히 작은 크기.
병렬 워커에서 함수 실행을 표현하고 결과를 유지하려면 Future 객체를 사용하십시오.
f(1:numel(partitions)-1) = parallel.FevalFuture;
parfeval
함수를 사용하여 병렬 워커에 계산을 분담할 수 있습니다. parameterSweep
은 탐색할 파라미터의 파티션에 대해 로렌츠 연립방정식을 푸는 아래 스크립트의 마지막에 정의된 헬퍼 함수입니다. 출력 인수가 하나이므로 parfeval
에서 출력값 개수로 1
을 지정해야 합니다.
for ii = 1:numel(partitions)-1 f(ii) = parfeval(@parameterSweep,1,partitions(ii),partitions(ii+1),sigma,rho,beta,Q); end
parfeval
은 MATLAB을 차단하지 않으므로 계산이 수행되는 동안 계속해서 작업할 수 있습니다. 워커는 병렬로 계산하며 중간 결과가 나오는 대로 DataQueue
를 통해 전송합니다.
parfeval
이 완료될 때까지 MATLAB을 차단하려면 Future 객체에 wait
함수를 사용하십시오. wait
함수는 후속 코드가 parfeval
계산의 완료 여부에 종속되는 경우에 유용하게 사용할 수 있습니다.
wait(f);
parfeval
계산을 마치면 wait
가 완료된 후에 나머지 코드를 실행할 수 있습니다. 예를 들어, 결과 곡면의 등고선을 플로팅해 보겠습니다. fetchOutputs
함수를 사용하여 Future 객체에 저장된 결과를 가져오십시오.
results = reshape(fetchOutputs(f),gridSize,[]); contourf(rho,sigma,results) xlabel('\rho','Interpreter','Tex') ylabel('\sigma','Interpreter','Tex')
파라미터 스윕에 계산 리소스가 더 필요하고 클러스터에 액세스할 수 있는 경우 parfeval
계산을 확장할 수 있습니다. 자세한 내용은 Scale Up from Desktop to Cluster 항목을 참조하십시오.
헬퍼 함수 정의하기
탐색할 파라미터의 파티션에 대해 로렌츠 연립방정식을 푸는 헬퍼 함수를 정의합니다. DataQueue
객체의 send
함수를 사용하여 중간 결과를 MATLAB 클라이언트로 전송합니다.
function results = parameterSweep(first,last,sigma,rho,beta,Q) results = zeros(last-first,1); for ii = first:last-1 lorenzSystem = @(t,a) [sigma(ii)*(a(2) - a(1)); a(1)*(rho(ii) - a(3)) - a(2); a(1)*a(2) - beta*a(3)]; [t,a] = ode45(lorenzSystem,[0 100],[1 1 1]); result = a(end,3); send(Q,[ii,result]); results(ii-first+1) = result; end end
새 데이터가 도착하면 곡면 플롯을 업데이트하는 또다른 헬퍼 함수를 정의합니다.
function updatePlot(surface,data) surface.ZData(data(1)) = data(2); drawnow('limitrate'); end
참고 항목
parpool
| parallel.pool.DataQueue
| afterEach
| parfeval