پردازش موازی (Parallel Computing)، در متلب (MATLAB)
نرم افزار متلب (MATLAB)، این امکان را به وجود آورده است که در کامپیوترهای دارای پردازنده (CPU) چند هسته ای، از پردازش موازی (Parallel Computing) استفاده نماییم. در این مبحث آموزشی، قصد داریم که نحوه نوشتن کدهای متلب، برای پردازش موازی کدها را شرح بدهیم.
اگر بخواهیم کدهای متلب مورد نظرمان، به صورت موازی پردازش شوند (Parallel Computing)، باید شکل کلی کدها به صورت زیر باشد :
که در آن، به جای عبارت % your code ، باید کدهای متلب مورد نظر خود را بنویسید. به عدد 2 که درون پرانتز اولین دستور matlabpool نوشته شده است، دقت کنید. این عدد، تعداد هسته هایی (cores) را مشخص کرده است که قصد دارید توسط آنها، کدها به صورت موازی پردازش شوند. بنابراین مثلا اگر کامپیوتر شما دارای پردازنده (CPU) چهار هسته ای است، شما باید، حداکثر، عدد 4 را به جای آن بنویسید.
همچنین توجه داشته باشید که شما می توانید تعیین کنید که همه کدها یا یک یا چند بخش دلخواه از کدهای برنامه متلب خود را، به صورت موازی، پردازش نمایید، بنابراین تنها کافی است که آن بخش های مورد نظرتان از کدها را در میان دو عبارتی که در بالا به کار بردیم، قرار بدهید.
به مثال زیر توجه کنید :
در ابتدا کدی را به صورت معمولی (بدون پردازش موازی) می نویسیم و سپس همان را با پردازش موازی اجرا می کنیم و همچنین زمان اجرای دستورات را هم با دو دستور tic و toc محاسبه می کنیم تا ببینیم که پردازش موازی، تا چه حد در سریعتر اجرا شدن کدهای برنامه، تاثیر دارد.
برنامه ای شامل دو حلقه for متوالی می نویسیم، البته به صورت معمولی و بدون استفاده از پردازش موازی :
نتیجه :
مشاهده می کنید که زمان پردازش دو حلقه for ، برابر 0.005328 ثانیه بوده است.
حالا این بار، از پردازش موازی استفاده می کنیم تا ببینیم تا چه حد، زمان پردازش کدها کاهش می یابد (البته شروع و پایان خود پردازش موازی، چند ثانیه ای طول می کشد، ولی این میزان، برای برنامه های بزرگ که مد نظر اکثر کاربران است، بسیار ناچیز می باشد) :
نتیجه :
این بار، زمان پردازش دو حلقه for ، برابر با 0.002837 ثانیه شده است که از عدد قبلی (پردازش معمولی، پردازش سری)، کمتر است (تقریبا نصف).
دقت کنید که در مثال بالا، تنها دو حلقه for ، به طور موازی با هم اجرا می شوند، اما کدهای درون هر حلقه for ، به طور موازی اجرا نخواهند شد (منظور این است که مثلا اگر حلقه for دارای پارامتری از 1 تا 1000 است، این 1000 بار اجرا شدن حلقه، به صورت موازی نخواهد بود، بلکه به صورت پشت سرهم و سری است).
اما شاید بخواهیم که دستورات درون حلقه for نیز به صورت موازی اجرا گردند، برای این منظور، دستور parfor در متلب در نظر گرفته شده است. دقت کنید که برای استفاده از این دستور، باید هر گام (مرحله) از اجرای حلقه، از دیگر گام ها (مراحل)، مستقل باشد. مثلا اگر در همان مثال قبل، دستورات for را به parfor تبدیل کنیم، جوابی غیر منتظره خواهیم گرفت، زیرا در هر گام (مرحله)، مقدار جدیدی برای متغیر x یا y در نظر گرفته می شود، بنابراین یک متغیر، در گام های مختلف اجرای حلقه، مقادیر متفاوت خواهد داشت، بنابراین گام ها از هم مستقل نیستند. پس نباید برای آن، پردازش موازی به کار ببریم.
برای درک بهتر این موضوع، به مثال زیر توجه کنید :
در همان مثال قبل، به جای دستور for ، از دستور parfor استفاده می کنیم :
نتیجه :
مطمئنا انتظار چنین پاسخی را نداشتید. این پاسخ به این دلیل است که گام های اجرای حلقه، از یکدیگر مستقل نمی باشند.
اکنون، مثالی را برای parfor می نویسیم که گام های حلقه، از هم مستقل باشند (تعداد گام ها را کمتر کردیم تا نتیجه راحتتر قابل نمایش باشد) :
نتیجه :
همین نتیجه را انتظار داشتیم. در این مثال، x و y ، هر کدام دارای 10 عنصر هستند، بنابراین، مثلا اگر دستورات حلقه x ، به صورت موازی اجرا گردند، 10 گام آن، در هم تداخل ندارند، زیرا هر بار، یک عنصر متفاوت از x ، مقداردهی می شود (هیچ دو گامی نیست که در آنها، یک عنصر مشترک، مقداردهی شود).
یک سوال مهم، این است که آیا اگر یک تابع را به صورت m-file تعریف کنیم و سپس در یک برنامه متلب، دو بار به آن ارجاع دهیم و برنامه به صورت موازی اجرا شود، آیا عملکرد برنامه صحیح خواهد بود ؟
حال پاسخ این سوال را، با یک مثال، برای شما روشن خواهم کرد :
ابتدا یک تابع با نام my_function می سازیم (تابع به صورت m-file)، که کدهای آن به صورت زیر می باشد :
این کدها را درون یک m-file نوشته ایم و سپس آن را با نام my_function ذخیره نموده ایم. همچنین فایل این تابع و فایل برنامه ای که خواهیم نوشت، باید هر دو درون یک فولدر قرار داشته باشند (یا مسیر هر دو فایل را به مسیر های نرم افزار متلب اضافه کنیم).
سپس یک برنامه متلب می نویسم که در آن، دو بار به تابع my_function ، ارجاع داده شده است :
نتیجه :
مشاهده می کنید که نتایج صحیح است (البته این پاسخ کلی برای پرسش مورد نظر ما است، اما امکان دارد که کاربر به گونه ای کدهای تابع را بنویسد که تداخل به وجود بیاید، مفهوم تداخل را قبلا برایتان شرح دادم، بنابراین همیشه چک کنید که آیا مراحل عملیات های مختلف، با هم تداخل دارند یا خیر).
بنابراین، پاسخ کلی، هم ((بله)) است و هم ((نه)). یعنی اینکه، اگر کدهای نوشته شده در تابع، تداخل در فرآیند موازی بودن، ایجاد نکنند، مشکلی وجود ندارد (مثل مثال ساده بالا که از کدهای تابع، معلوم است که تداخل ایجاد نمی کند) و اگر کدها، تداخل به وجود بیاورند، نتیجه ای غیرمنتظره، به دست می آید. پس همیشه، مستقل بودن پردازش ها را بررسی کنید.
برای آنکه خیالتان راحت شود، ابتدا کدها را به صورت سری اجرا کنید و بعد یک بار به طور موازی و نتیجه ها را مقایسه کنید. نگویید که خوب اینکه شد دو بار اجرا ! ، خود من برای شبیه سازی انتشار امواج الکترومغناطیسی، مجبور بودم که یک برنامه سنگین بنویسم و بعد آن را صدها بار اجرا نمایم و حالت های مختلف ناشی از ورودی های مختلف را مقایسه کنم، بنابراین تنها کافی بود که دو بار برنامه را تست کنم و بعد صدها بار، از برنامه بهینه شده استفاده نمایم، پروژه تعداد زیادی از دانشجویان نیز اینگونه می باشد و در این موارد، این دو بار اجرا کردن، ارزش صرف وقت را دارد.
akiko thanks
پارسا سلام. ممنون از اطلاعات مفیدتون
سوالی که دارم این است که من کدی دارم که در هر مرحله عملیاتی را انجام میدهد و نتایج در مرحله بعدی بکار میرود. زمان انجام آنالیز در حدود 2 هفته می شود. آیا میتوان با پردازش موازی زمان را کاهش داد؟ ممنون
admin سلام.
با توجه به اینکه باید نتایج هر مرحله برای استفاده در مرحله بعد به کار رود، پس باید هر مرحله را با پردازش موازی اجرا کنید و پس از پایان، به سراغ مرحله بعد بروید (نمی توان مراحل را به صورت موازی اجرا کرد).
نکته مهم این است که کدهای هر مرحله هم مهم است و باید ببینید که می توان آنها را به صورت موازی پردازش کرد یا خیر. بنابراین نمیشه نظر کلی داد و باید خود برنامه نویس، با بررسی کدها، این مورد را چک کند.
sina سلام.خسته نباشید.بابت سایت خوب و مفیدتون بسیار متشکرم.
سوالی داشتم
برنامه نوشتم که باید حلقه های آن از یک تا 2000 و بیشتر بچرخد. مشکل این است که از دور 499 به بعد یک خط از حلقه که عمل به توان رساندن را انجام می دهد inf می شود
برای حل این مشکل (برطرف کردن inf و دیدن جواب) چه راه کاری را پیشنهاد می کنید؟
سیاوش سلام دمتون گرم
مه سلام
ممنون از اطلاعات خوبتون
من کدی نوشتم که نیاز به شبکه بندی بسیار ریز دارد. وقتی تعداد شبکه رو بالا تعیین میکنم متلب پیغام out of memory میده.
آیا با موازی سازی ممکنه این مشکل حل بشه یا تنها راهش استفاده از کامپیوتر قویتره؟
ممنون
وحید مرسی از مطلبتون
سوال من اینه که ما از متلب برای بهینه سازی استفاده میکنیم و برای تابع فیتنسمون از یه نرم افزار استفاده میشه
یعنی این نرم افزار باید ران بشه به دفعات. این نرم افزار از متلب توسط بتچ فایل فراخوانی میشه.
ایا راهی هست که همزمان چند بتچ فایل فعال بشن و سپس بدون تداخل باهم چند نرم افزارو به صورت موازی ران کنن؟
hitman این دستور دیگه کار نمیکنه ممنون میشم اگه دستور parpool رو توضیح بدین
فاطمه سلام
ممنون از سایت خوبتون
من میخواستم نمودار هم براش رسم کنم باید چیکار کنم
اگه میشه لطف کنید راهنماییم کنید
meisamfeizollahzadeh سلام و خسته نباشید لطفا مرا در اجرای این کد در متلب راهنمایی فرمایید میگن کد تست شده و مشکلی نداره : کد زیر کد روش درونیابی نیوتن در متلب است که در زمان اجرا داده ها رو به شکل زیر وارد کردم اما این ارور رو میده که خیلی لازمه باید چهارشنبه ارائه بدم.از زحمات و راهنماییهایتان نهایت تشکر رادارم.
x=[300,400,500,600,700,800]
lengh(x)
y=[.024,.035,.046,.058,.067,0.083]
a=350
newtondd(x,y,a)
n=lengh(x); Error in ==> newtondd at 11
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------
کد:
<?php
function newtondd(x,y,a)
% Newton divided difference
disp(' Newton divided difference')
disp('______________________________________________________________________')
disp(' x y f[,] f[,,] f[,,,] ... ')
disp('______________________________________________________________________')
n=lengh(x);
for k=1:n-1
% Compute first divided difference
d(k,1)=(y(k+1)-y(k))/(x(k+1)-x(k));
end;
for i=2:n-1
for k=1:n-i
% Compute ith divided difference
d(k,i)=(d(k+1,i-1)-d(k,i-1))/(x(k+i)-x(k));
end;
end;
% print results
if (rem(n,2)==0)
p=n/2;
m=n/2;
else
p=fix(n/2);
m=fix(n/2)+1;
end;
for i=1:p
fprintf('%8.2f %8.2f',x(i),y(i));
for k=1:i-1
fprintf(' %8.5f',d(i-k,2*k));
end;
fprintf('\n ')
for k=1:i
fprintf(' %8.5f',d(i-k+1,2*k-1));
end;
fprintf('\n')
end;
j=p;
for i=m:-1:1
j=j+1;
fprintf('%8.2f %8.2f',x(j),y(j));
for k=1:i-1
fprintf(' %8.5f',d(j-k,2*k));
end;
fprintf('\n ')
for k=1:i-1
fprintf(' %8.5f',d(j-k+1,2*k-1));
end;
fprintf('\n')
end;
%%
%%
p=0;
for ii=1:length(x)
if ii==1
p=p+y(1);
else
temp=1;
for jj=2:ii
temp=temp*(a-x(jj-1));
end
p=p+d(1,ii-1)*temp;
end
end
fprintf('Pn(a=%.2f)=%.2f\n',a,p);
ati 
mary 
مهدی 
مهمان 
فیروزه 
حسین 
سیاوش 
ابوعلی زاده