{"id":164499,"date":"2012-12-31T20:27:04","date_gmt":"2012-12-31T16:27:04","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=164499"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=164499","title":{"rendered":"<span class=\"post_title\">\u041c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0435 \u044d\u043b\u0435\u043a\u0442\u0440\u0438\u0447\u0435\u0441\u043a\u043e\u0433\u043e \u043f\u043e\u043b\u044f \u0441\u0440\u0435\u0434\u0441\u0442\u0432\u0430\u043c\u0438 CUDA<\/span>"},"content":{"rendered":"<div class=\"content html_format\"> \t\t\t\u0414\u0430\u043d\u043d\u0430\u044f \u0441\u0442\u0430\u0442\u044c\u044f \u043d\u0430\u043f\u0438\u0441\u0430\u043d\u0430 \u0441 \u0446\u0435\u043b\u044c\u044e \u043f\u0440\u043e\u0434\u0435\u043c\u043e\u043d\u0441\u0442\u0440\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u043a\u0430\u043a \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0442\u0435\u0445\u043d\u043e\u043b\u043e\u0433\u0438\u0438 CUDA \u043c\u043e\u0436\u043d\u043e \u0441\u043c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u043f\u0440\u043e\u0441\u0442\u043e\u0435 \u0432\u0437\u0430\u0438\u043c\u043e\u0434\u0435\u0439\u0441\u0442\u0432\u0438\u0435 \u0437\u0430\u0440\u044f\u0436\u0435\u043d\u044b\u0445 \u0447\u0430\u0441\u0442\u0438\u0446 (\u0441\u043c. <a href=\"http:\/\/en.wikipedia.org\/wiki\/Coulomb%27s_law\">\u0417\u0430\u043a\u043e\u043d \u041a\u0443\u043b\u043e\u043d\u0430<\/a>). \u0414\u043b\u044f \u0432\u044b\u0432\u043e\u0434\u0430 \u0441\u0442\u0430\u0442\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438 \u044f \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043b \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0443 <b>freeglut<\/b>.<br \/>  \u041a\u0430\u043a \u043f\u0438\u0448\u0443\u0442 \u0447\u0430\u0441\u0442\u0435\u043d\u044c\u043a\u043e \u043d\u0430 \u0425\u0430\u0431\u0440\u0435: <a name=\"habracut\"><\/a><\/p>\n<p>  \u0415\u0449\u0451 \u0432\u043e \u0432\u0440\u0435\u043c\u0435\u043d\u0430 Turbo Pascal \u0438 Borland C 5.02 \u043c\u0435\u043d\u044f \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043e\u0432\u0430\u043b \u043f\u0440\u043e\u0446\u0435\u0441\u0441 \u0441\u0438\u043c\u0443\u043b\u044f\u0446\u0438\u0438 \u043f\u0440\u043e\u0441\u0442\u044b\u0445 \u0444\u0438\u0437\u0438\u0447\u0435\u0441\u043a\u0438\u0445 \u044f\u0432\u043b\u0435\u043d\u0438\u0439 \u043f\u0440\u0438\u0440\u043e\u0434\u044b \u2014 \u0431\u0443\u0434\u044c \u0442\u043e \u0441\u0438\u043b\u044b \u0433\u0440\u0430\u0432\u0438\u0442\u0430\u0446\u0438\u043e\u043d\u043d\u043e\u0433\u043e \u043f\u0440\u0438\u0442\u044f\u0436\u0435\u043d\u0438\u044f \u0438\u043b\u0438 \u043a\u0443\u043b\u043e\u043d\u043e\u0432\u0441\u043a\u043e\u0433\u043e \u0432\u0437\u0430\u0438\u043c\u043e\u0434\u0435\u0439\u0441\u0442\u0432\u0438\u044f. \u041a \u0441\u043e\u0436\u0430\u043b\u0435\u043d\u0438\u044e, \u043c\u043e\u0449\u043d\u043e\u0441\u0442\u0438 \u043e\u0434\u043d\u043e\u0433\u043e \u043f\u0440\u043e\u0446\u0435\u0441\u0441\u043e\u0440\u0430 \u0441 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u0438\u043c\u0438 \u0441\u043e\u0442\u043d\u044f\u043c\u0438 \u041c\u0413\u0446 \u043d\u0435 \u0445\u0432\u0430\u0442\u0430\u043b\u043e \u0434\u043b\u044f \u043f\u0440\u043e\u0441\u0442\u044b\u0445 \u043c\u0430\u0442\u0435\u043c\u0430\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0445 \u043e\u043f\u0435\u0440\u0430\u0446\u0438\u0439, \u043a\u043e\u0438\u0445 \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u043e \u0431\u044b\u043b\u043e \u0432\u0435\u043b\u0438\u043a\u043e.<\/p>\n<p>  \u0425\u043e\u0442\u044c \u044f \u0438 \u043d\u0435 \u043f\u043e\u0433\u0440\u043e\u043c\u043c\u0438\u0441\u0442, \u043c\u043d\u0435 \u0437\u0430\u0445\u043e\u0442\u0435\u043b\u043e\u0441\u044c \u043e\u0434\u043d\u0430\u0436\u0434\u044b \u0432\u0441\u043f\u043e\u043c\u043d\u0438\u0442\u044c \u0434\u0430\u043b\u0451\u043a\u043e\u0435 \u0434\u0435\u0442\u0441\u0442\u0432\u043e, \u0434\u0430 \u0438 \u043f\u043e\u0434\u0443\u0447\u0438\u0442\u044c \u043e\u0441\u043d\u043e\u0432\u044b JavaScript, \u0438 \u044f \u043a\u0430\u043a-\u0442\u043e \u0441\u043c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u043b \u0433\u0440\u0430\u0432\u0438\u0442\u0430\u0446\u0438\u043e\u043d\u043d\u043e\u0435 \u0432\u0437\u0430\u0438\u043c\u043e\u0434\u0435\u0439\u0441\u0442\u0432\u0438\u0435 \u043b\u044e\u0431\u043e\u0433\u043e \u0440\u0430\u0437\u0443\u043c\u043d\u043e\u0433\u043e \u0447\u0438\u0441\u043b\u0430 \u043c\u0430\u0442\u0435\u0440\u0438\u0430\u043b\u044c\u043d\u044b\u0445 \u0442\u043e\u0447\u0435\u043a. \u0422\u0430\u043a, \u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u0447\u0430\u0441\u0442\u0438\u0447\u043a\u0430 \u043d\u0438\u0436\u0435 \u043a\u0440\u0443\u0442\u0438\u0442\u0441\u044f \u0432\u043e\u043a\u0440\u0443\u0433 \u00ab\u0441\u043e\u043b\u043d\u0446\u0430\u00bb \u043f\u043e\u043b\u044c\u0437\u0443\u044f\u0441\u044c \u041d\u044c\u044e\u0442\u043e\u043d\u043e\u0432\u0441\u043a\u0438\u043c \u0437\u0430\u043a\u043e\u043d\u043e\u043c:<br \/>  <img decoding=\"async\" src=\"http:\/\/habrastorage.org\/storage2\/a99\/1d7\/d8c\/a991d7d8cad963868cfa2b535f0870a6.png\"\/><br \/>  \u0415\u0441\u043b\u0438 \u0432\u0430\u043c \u0431\u0443\u0434\u0435\u0442 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u043e, \u0438 \u043c\u0435\u043d\u044f \u043d\u0435 \u043e\u0441\u0442\u0430\u0432\u0438\u0442 \u044d\u043d\u0442\u0443\u0437\u0438\u0430\u0437\u043c, \u0442\u043e \u044f \u043e\u043f\u0438\u0448\u0443 \u0432\u0430\u043c \u0441\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u0442\u0430\u043a\u043e\u0439 \u0432\u043e\u0442 \u0434\u0435\u043c\u043a\u0438, \u0447\u0442\u043e\u0431\u044b \u0432\u044b \u0441\u043c\u043e\u0433\u043b\u0438 \u043d\u0430\u043f\u0438\u0441\u0430\u0442\u044c \u0441\u0432\u043e\u0438\u0445 \u0437\u043b\u044b\u0445 \u043f\u0442\u0438\u0446 \u0432 \u0432\u0430\u043a\u0443\u0443\u043c\u0435.<\/p>\n<p>  \u0418\u0442\u0430\u043a, \u0445\u0432\u0430\u0442\u0438\u0442 \u043d\u043e\u0441\u0442\u0430\u043b\u044c\u0433\u0438\u0440\u043e\u0432\u0430\u0442\u044c, \u0434\u0430\u0432\u0430\u0439\u0442\u0435 <s>\u043f\u043e\u0433\u0440\u043e\u043c<\/s>\u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u0438\u0440\u043e\u0432\u0430\u0442\u044c. \u041c\u043e\u044f \u043f\u043b\u0430\u0442\u0444\u043e\u0440\u043c\u0430 \u2014 \u044d\u0442\u043e \u043d\u043e\u0443\u0442\u0431\u0443\u043a \u0441 Ubuntu \u0438 \u0432\u0438\u0434\u0435\u043e\u043a\u0430\u0440\u0442\u043e\u0439 GeForce GT 540m \u0438 \u0443\u0441\u0442\u0430\u043d\u043e\u0432\u043b\u0435\u043d\u043d\u044b\u043c CUDA Toolkit 5.0. \u0415\u0441\u043b\u0438 \u0432\u044b \u043d\u0430\u0447\u0438\u043d\u0430\u0435\u0442\u0435 \u0437\u043d\u0430\u043a\u043e\u043c\u0441\u0442\u0432\u043e \u0441 CUDA, \u0442\u043e \u043c\u043e\u0433\u0443 \u0432\u0430\u043c \u043f\u043e\u0441\u043e\u0432\u0435\u0442\u043e\u0432\u0430\u0442\u044c \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u0443\u044e \u043a\u043d\u0438\u0433\u0443 <a href=\"https:\/\/developer.nvidia.com\/content\/cuda-example-introduction-general-purpose-gpu-programming-0\">CUDA by Example<\/a>.<\/p>\n<p>  \u0414\u043b\u044f \u043d\u0430\u0447\u0430\u043b\u0430 \u0441\u043e\u0437\u0434\u0430\u0439\u0442\u0435 \u043f\u0440\u043e\u0435\u043a\u0442 \u0438\u0437 \u0448\u0430\u0431\u043b\u043e\u043d\u0430 Nsight \u0438 \u043f\u043e\u0434\u043a\u043b\u044e\u0447\u0438\u0442\u0435 \u043a \u043d\u0435\u043c\u0443 freeglut \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044f \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0438 <i>GLU, GL, glut<\/i>. \u042d\u0442\u043e \u043f\u043e\u0437\u0432\u043e\u043b\u0438\u0442 \u043d\u0430\u043c \u043e\u0442\u043e\u0431\u0440\u0430\u0437\u0438\u0442\u044c \u043f\u043e\u043b\u0443\u0447\u0435\u043d\u043d\u043e\u0435 \u0432\u0435\u043a\u0442\u043e\u0440\u043d\u043e\u0435 \u043f\u043e\u043b\u0435, \u043a\u043e\u0442\u043e\u0440\u043e\u0435 \u043c\u044b \u0440\u0430\u0441\u0441\u0447\u0438\u0442\u0430\u0435\u043c. \u041f\u0440\u0438\u0431\u043b\u0438\u0437\u0438\u0442\u0435\u043b\u044c\u043d\u0430\u044f \u0441\u0442\u0440\u0443\u043a\u0442\u0443\u0440\u0430 \u043d\u0430\u0448\u0435\u0433\u043e \u0444\u0430\u0439\u043b\u0430 \u0441 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439 <code>main<\/code> \u0431\u0443\u0434\u0435\u0442 \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0430\u044f:  <\/p>\n<pre><code class=\"cpp\">#include &lt;stdio.h&gt; #include &lt;stdlib.h&gt;  int main(int argc, char** argv) { \t\/\/ Initialize freeglut \t\/\/ ... \t \t\/\/ Initialize Host and Device variables \t\/\/ ... \t \t\/\/ Launch Kernel to render the image \t \t\/\/ Display Image \t \t\/\/ Free resources \t \treturn 0; } <\/code><\/pre>\n<p>  \u0414\u043b\u044f \u043d\u0430\u0447\u0430\u043b\u0430, \u0432\u0441\u0451 \u0432\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u043f\u0440\u043e\u0441\u0442\u043e. \u0414\u0430\u0432\u0430\u0439\u0442\u0435, \u0441\u043f\u0435\u0440\u0432\u0430, \u043e\u0442\u043e\u0431\u0440\u0430\u0437\u0438\u043c \u043e\u043a\u043d\u043e \u043d\u0443\u0436\u043d\u043e\u0433\u043e \u043d\u0430\u043c \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u043d\u0443 \u0438 \u0441\u0440\u0435\u043d\u0434\u0435\u0440\u0438\u043c \u0447\u0451\u0440\u043d\u044b\u0439 \u044d\u043a\u0440\u0430\u043d. \u041a\u0441\u0442\u0430\u0442\u0438, \u0441\u043a\u0430\u0436\u0443 \u0432\u0430\u043c \u043e\u0447\u0435\u0432\u0438\u0434\u043d\u0443\u044e \u0432\u0435\u0449\u044c: <i>\u0435\u0441\u043b\u0438 \u043a\u0430\u043a\u0430\u044f-\u0442\u043e \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0432\u0430\u043c \u043d\u0435\u0437\u043d\u0430\u043a\u043e\u043c\u0430, \u0442\u043e \u043f\u0440\u043e\u0441\u0442\u043e \u043f\u043e\u0433\u0443\u0433\u043b\u0438\u0442\u0435<\/i>.  <\/p>\n<h4>\u0418\u043d\u0438\u0446\u0438\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u044f \u043e\u043a\u043d\u0430<\/h4>\n<p>  \u0422\u0430\u043a \u043a\u0430\u043a \u044f, \u043a\u043e\u0433\u0434\u0430 \u0437\u0430\u0434\u0443\u043c\u044b\u0432\u0430\u043b \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u0435, \u043d\u0435 \u043f\u043b\u0430\u043d\u0438\u0440\u043e\u0432\u0430\u043b \u0434\u0438\u043d\u0430\u043c\u0438\u0447\u0435\u0441\u043a\u043e\u0435 \u0438\u0437\u043c\u0435\u043d\u0435\u043d\u0438\u0435 \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u043e\u043a\u043d\u0430, \u0442\u043e \u044f \u0440\u0435\u0448\u0438\u043b \u0443\u043a\u0430\u0437\u0430\u0442\u044c \u0435\u0433\u043e \u0434\u0435\u043a\u043b\u0430\u0440\u0430\u0442\u0438\u0432\u043d\u043e, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044f \u0433\u043b\u043e\u0431\u0430\u043b\u044c\u043d\u044b\u0435 \u043a\u043e\u043d\u0441\u0442\u0430\u043d\u0442\u044b:  <\/p>\n<pre><code class=\"cpp\">const int width = 1280; const int height = 720; <\/code><\/pre>\n<p>  \u0410\u0445 \u0434\u0430! \u041d\u0435 \u0441\u0443\u0434\u0438\u0442\u0435 \u0441\u0442\u0440\u043e\u0433\u043e \u0437\u0430 \u0441\u043c\u0435\u0448\u0438\u0432\u0430\u043d\u0438\u0435 \u0441\u0442\u0438\u043b\u0435\u0439 C++ \u0438 C99. \u041d\u0435\u043a\u043e\u0442\u043e\u0440\u044b\u0435 \u0432\u0435\u0449\u0438, \u0442\u0438\u043f\u0430 \u043f\u0435\u0440\u0435\u0433\u0440\u0443\u0437\u043a\u0438, \u044f \u043f\u043e\u0440\u043e\u0439 \u043d\u0430\u0445\u043e\u0436\u0443 \u043f\u043e\u043b\u0435\u0437\u043d\u044b\u043c, \u043d\u043e \u0438 \u043e\u0442 malloc \u044f \u043d\u0435 \u043e\u0442\u043a\u0430\u0437\u0430\u043b\u0441\u044f, \u0447\u0442\u043e\u0431\u044b \u0445\u043e\u0442\u044c \u043a\u0430\u043a-\u043d\u0438\u0431\u0443\u0434\u044c \u0441\u043b\u0435\u0434\u043e\u0432\u0430\u0442\u044c \u0441\u0442\u0438\u043b\u044e \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u044f\u043c \u0441 CUDA.<\/p>\n<p>  \u0421\u0430\u043c\u043e \u043e\u0442\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435 \u043e\u043a\u043d\u0430, \u043e\u0431\u0440\u0430\u0431\u043e\u0442\u043a\u0430 \u043a\u043b\u0430\u0432\u0438\u0448\u0438 Escape \u0438 \u043f\u0440\u043e\u0441\u0442\u043e\u0439 \u0440\u0435\u043d\u0434\u0435\u0440 \u043c\u043e\u0436\u043d\u043e \u043e\u0441\u0443\u0449\u0435\u0441\u0442\u0432\u0438\u0442\u044c \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u043c \u043a\u043e\u0434\u043e\u043c:  <\/p>\n<pre><code class=\"cpp\">#include &lt;GL\/freeglut.h&gt; #include &lt;GL\/gl.h&gt; #include &lt;GL\/glext.h&gt;  \/\/ ...  void key(unsigned char key, int x, int y) { \tswitch (key) { \tcase 27: \t\tprintf(&quot;Exit application\\n&quot;); \t\tglutLeaveMainLoop(); \t\tbreak; \t} }  void draw(void) { \tglClearColor(0.0, 0.0, 0.0, 1.0); \tglClear(GL_COLOR_BUFFER_BIT); \t\/\/ glDrawPixels(width, height, GL_RGBA, GL_UNSIGNED_BYTE, screen); \tglFlush(); }  \/\/ ... \t\/\/ Initialize freeglut \tglutInit(&argc, argv); \tglutInitDisplayMode(GLUT_SINGLE | GLUT_RGBA); \tglutInitWindowSize(width, height); \tglutCreateWindow(&quot;Electric field&quot;); \tglutDisplayFunc(draw); \tglutKeyboardFunc(key); \tglutSetOption(GLUT_ACTION_ON_WINDOW_CLOSE, GLUT_ACTION_CONTINUE_EXECUTION);  \t\/\/ ...  \t\/\/ Display Image \tglutMainLoop(); <\/code><\/pre>\n<p>  \u041f\u0440\u0438 \u0437\u0430\u043f\u0443\u0441\u043a\u0435 \u044d\u0442\u043e\u0433\u043e \u043a\u043e\u0434\u0430, \u0432\u044b \u043d\u0438\u0447\u0435\u0433\u043e \u043a\u0440\u043e\u043c\u0435 \u0447\u0451\u0440\u043d\u043e\u0433\u043e \u044d\u043a\u0440\u0430\u043d\u0430 \u043d\u0435 \u0443\u0432\u0438\u0434\u0438\u0442\u0435. \u0422\u0430\u043a \u0447\u0442\u043e \u044f \u043e\u043f\u0443\u0449\u0443 \u0437\u0434\u0435\u0441\u044c \u0441\u043a\u0440\u0438\u043d\u0448\u043e\u0442. \u041e\u0431\u0440\u0430\u0442\u0438\u0442\u0435 \u0432\u043d\u0438\u043c\u0430\u043d\u0438\u0435 \u043d\u0430 \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0435 \u043c\u043e\u043c\u0435\u043d\u0442\u044b:  <\/p>\n<ul>\n<li><code>glutLeaveMainLoop<\/code> \u2014 \u043a\u043e\u0433\u0434\u0430 \u043d\u0430\u0436\u0430\u0442 Escape, \u0442\u043e freeglut \u043f\u043e\u0437\u0432\u043e\u043b\u044f\u0435\u0442 \u043f\u0440\u0435\u0440\u0432\u0430\u0442\u044c \u0433\u043b\u0430\u0432\u043d\u044b\u0439 \u0446\u0438\u043a\u043b, \u0447\u0442\u043e\u0431\u044b \u043c\u044b \u0441\u043c\u043e\u0433\u043b\u0438 \u043e\u0441\u0432\u043e\u0431\u043e\u0434\u0438\u0442\u044c \u0440\u0435\u0441\u0443\u0440\u0441\u044b \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u044f. \u0415\u0441\u043b\u0438 \u0445\u043e\u0447\u0435\u0442\u0441\u044f, \u0442\u043e \u043c\u043e\u0436\u043d\u043e \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u0442\u044c \u0438 <code>atexit<\/code>;<\/li>\n<li><code>glDrawPixels<\/code> \u2014 \u043f\u043e\u0437\u0436\u0435 \u044f \u0441\u043e\u0437\u0434\u0430\u043c \u0431\u0443\u0444\u0435\u0440 \u0438\u0437 <code>unsigned byte<\/code>, \u0433\u0434\u0435 \u044f \u0431\u0443\u0434\u0443 \u0445\u0440\u0430\u043d\u0438\u0442\u044c \u0432\u0441\u0435 \u043f\u0438\u043a\u0441\u0435\u043b\u0438 (\u0438\u0445 \u0446\u0432\u0435\u0442\u0430) \u0438 \u043f\u043e\u0442\u043e\u043c \u043e\u0442\u043e\u0431\u0440\u0430\u0436\u0430\u0442\u044c \u0432 \u043e\u043a\u043d\u0435;<\/li>\n<li><code>glutSetOption<\/code> \u2014 \u044d\u0442\u043e \u043c\u043d\u0435 \u043f\u043e\u0437\u0432\u043e\u043b\u0438\u0442\u044c \u043f\u0440\u043e\u0434\u043e\u043b\u0436\u0438\u0442\u044c \u0432\u044b\u043f\u043e\u043b\u043d\u0435\u043d\u0438\u0435 \u043a\u043e\u0434\u0430 \u0432 <code>main<\/code>, \u0447\u0442\u043e\u0431\u044b \u043e\u0447\u0438\u0441\u0442\u0438\u0442\u044c \u0440\u0435\u0441\u0443\u0440\u0441\u044b \u0432 \u0432\u0438\u0434\u0435 \u0434\u0438\u043d\u0430\u043c\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u043f\u0430\u043c\u044f\u0442\u0438.<\/li>\n<\/ul>\n<p>  <\/p>\n<h4>\u0420\u0430\u0437\u0431\u0438\u0435\u043d\u0438\u0435 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u0435 \u043d\u0430 \u0431\u043b\u043e\u043a\u0438 \u0438 \u043f\u043e\u0442\u043e\u043a\u0438<\/h4>\n<p>  \u041f\u043e\u0431\u0430\u043b\u043e\u0432\u0430\u0432\u0448\u0438\u0441\u044c \u0441 \u043f\u0440\u0438\u043c\u0435\u0440\u0430\u043c\u0438, \u044f \u0441\u043c\u043e\u0433 \u0434\u043b\u044f \u0441\u0435\u0431\u044f \u0443\u044f\u0441\u043d\u0438\u0442\u044c, \u0447\u0442\u043e \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u044f\u0435\u0442 \u0441\u043e\u0431\u043e\u0439 \u0442\u0430\u043a\u043e\u0435 \u043f\u043e\u043d\u044f\u0442\u0438\u0435 \u043a\u0430\u043a \u0431\u043b\u043e\u043a \u0438 \u043f\u043e\u0442\u043e\u043a, \u0438 \u0447\u0442\u043e \u0434\u0435\u043b\u0430\u0442\u044c \u0441 \u0438\u0445 \u0440\u0430\u0437\u043c\u0435\u0440\u043d\u043e\u0441\u0442\u044c\u044e. \u0412 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435, \u0435\u0441\u043b\u0438 \u043e\u043f\u0443\u0441\u0442\u0438\u0442\u044c \u0442\u0435\u0445\u043d\u0438\u0447\u0435\u0441\u043a\u0438\u0435 \u0434\u0435\u0442\u0430\u043b\u0438, \u0442\u043e \u0431\u043b\u043e\u043a \u2014 \u044d\u0442\u043e \u043a\u043b\u0435\u0442\u043a\u0430, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0435\u0441\u0442\u044c \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u043f\u043e\u0442\u043e\u043a\u043e\u0432. \u041a\u043e\u0433\u0434\u0430 \u043f\u0440\u043e\u0438\u0441\u0445\u043e\u0434\u0438\u0442 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u0435 \u0446\u0435\u043b\u043e\u0439 \u0441\u0435\u0442\u043a\u0438 \u0431\u043b\u043e\u043a\u043e\u0432, \u0442\u043e \u0432 \u043a\u0430\u043a\u043e\u0439-\u0442\u043e \u043c\u043e\u043c\u0435\u043d\u0442 \u0432\u0440\u0435\u043c\u0435\u043d\u0438 \u043a\u0430\u0436\u0434\u044b\u0439 \u0431\u043b\u043e\u043a \u0432\u044b\u0441\u0442\u0440\u0435\u043b\u0438\u0432\u0430\u0435\u0442 \u0441\u0432\u043e\u0438 \u043f\u043e\u0442\u043e\u043a\u0438, \u0434\u043e\u0436\u0438\u0434\u0430\u0435\u0442\u0441\u044f \u0438\u0445 \u0437\u0430\u0432\u0435\u0440\u0448\u0435\u043d\u0438\u044f \u0438 \u043f\u0435\u0440\u0435\u0434\u0430\u0435\u0442 \u0443\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u0438\u0435 \u0434\u0440\u0443\u0433\u043e\u043c\u0443 \u0431\u043b\u043e\u043a\u0443. \u0422\u0430\u043a, \u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u0432\u044b \u043c\u043e\u0436\u0435\u0442\u0435 \u043d\u0430\u0433\u0435\u043d\u0435\u0440\u0438\u0442\u044c \u0431\u043b\u043e\u043a\u0438 \u0432 \u043e\u0434\u043d\u043e\u043c\u0435\u0440\u043d\u0443\u044e \u043f\u043e\u043b\u043e\u0441\u043a\u0443, \u0438\u043b\u0438 \u0432 \u0434\u0432\u0443\u0445\u043c\u0435\u0440\u043d\u0443\u044e \u0441\u0435\u0442\u043a\u0443.<\/p>\n<p>  \u041f\u043e\u0442\u043e\u043a\u0438, \u0432 \u0441\u0432\u043e\u044e \u043e\u0447\u0435\u0440\u0435\u0434\u044c, \u0442\u043e\u0436\u0435 \u043c\u043e\u0433\u0443\u0442 \u0440\u0430\u0441\u043f\u043e\u043b\u0430\u0433\u0430\u0442\u044c\u0441\u044f \u0432 \u0431\u043b\u043e\u043a\u0435 \u0432 \u0432\u0438\u0434\u0435 \u043e\u0434\u043d\u043e\u043c\u0435\u0440\u043d\u043e\u0439 \u043f\u043e\u043b\u043e\u0441\u044b, \u0438\u043b\u0438 \u0434\u0432\u0443\u0445\u043c\u0435\u0440\u043d\u043e\u0439 \u0442\u0430\u0431\u043b\u0438\u0446\u044b. \u041a\u0441\u0442\u0430\u0442\u0438, \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u044b \u0432\u0430\u0440\u0438\u0430\u043d\u0442\u044b \u0441 \u0442\u0440\u0435\u043c\u044f \u0438\u0437\u043c\u0435\u0440\u0435\u043d\u0438\u044f\u043c\u0438, \u043d\u043e \u044f \u0438\u0445 \u0443\u043f\u0443\u0449\u0443, \u0442\u0430\u043a \u043a\u0430\u043a \u043e\u043d\u0438 \u043f\u043e\u043a\u0430 \u043d\u0435 \u0432\u043f\u0438\u0441\u044b\u0432\u0430\u044e\u0442\u0441\u044f \u0432 \u043c\u043e\u044e \u0437\u0430\u0434\u0430\u0447\u0443.<\/p>\n<p>  \u042f \u043c\u043e\u0433\u0443 \u0440\u0430\u0437\u0431\u0438\u0442\u044c \u044d\u043a\u0440\u0430\u043d \u043d\u0430 \u0431\u043e\u043b\u044c\u0448\u0438\u0435 \u043a\u0432\u0430\u0434\u0440\u0430\u0442\u044b (\u0431\u043b\u043e\u043a\u0438), \u0432\u043d\u0443\u0442\u0440\u0438 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u043f\u0438\u043a\u0441\u0435\u043b\u0438 \u0431\u0443\u0434\u0443\u0442 \u0440\u0435\u043d\u0434\u0435\u0440\u0438\u0442\u044c\u0441\u044f \u043e\u0442\u0434\u0435\u043b\u044c\u043d\u044b\u043c \u043f\u043e\u0442\u043e\u043a\u043e\u043c. \u0422\u0430\u043a\u0430\u044f \u0437\u0430\u0434\u0430\u0447\u0430 \u2014 \u043e\u043f\u0442\u0438\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043d\u0438\u0435 \u0433\u0440\u0430\u0444\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u043a\u0430\u0440\u0442\u044b, \u043d\u043e \u0447\u0442\u043e\u0431\u044b \u0440\u0430\u0437\u0431\u0438\u0442\u044c \u044d\u043a\u0440\u0430\u043d \u043d\u0430 \u0431\u043b\u043e\u043a\u0438, \u043f\u0440\u0438\u0434\u0435\u0442\u0441\u044f \u043f\u0440\u0438\u043c\u0435\u043d\u0438\u0442\u044c \u043f\u0440\u043e\u0441\u0442\u0443\u044e \u0430\u0440\u0438\u0444\u043c\u0435\u0442\u0438\u043a\u0443 \u2014 \u044f \u0443\u0437\u043d\u0430\u044e \u0447\u0438\u0441\u043b\u043e \u043f\u043e\u0442\u043e\u043a\u043e\u0432 \u0434\u043e\u0437\u0432\u043e\u043b\u0435\u043d\u043d\u043e\u0433\u043e \u043c\u043d\u0435 \u0432 \u043f\u0440\u0435\u0434\u0435\u043b\u0430\u0445 \u0431\u043b\u043e\u043a\u0430, \u0438 \u0437\u0430\u0442\u0435\u043c \u043f\u043e\u0434\u0435\u043b\u044e \u0440\u0430\u0437\u043c\u0435\u0440 \u043e\u043a\u043d\u0430 \u043d\u0430 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u043f\u043e\u0442\u043e\u043a\u043e\u0432 \u043f\u043e \u0433\u043e\u0440\u0438\u0437\u043e\u043d\u0442\u0430\u043b\u0438 \u0438 \u0432\u0435\u0440\u0442\u0438\u043a\u0430\u043b\u0438.<\/p>\n<p>  \u0427\u0442\u043e\u0431\u044b \u0443\u0437\u043d\u0430\u0442\u044c \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u043f\u043e\u0442\u043e\u043a\u043e\u0432 \u0432\u043d\u0443\u0442\u0440\u0438 \u0431\u043b\u043e\u043a\u0430, \u0432\u043e\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u043c\u0441\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439 <code>cudaGetDeviceProperties<\/code>:  <\/p>\n<pre><code class=\"cpp\">int threadsCount(void) { \tint deviceCount = 0; \tcudaGetDeviceCount(&deviceCount);  \tif (deviceCount &lt;= 0) { \t\tprintf(&quot;No CUDA devices\\n&quot;); \t\texit(-1); \t}  \tcudaDeviceProp properties; \tcudaGetDeviceProperties(&properties, 0);  \treturn properties.maxThreadsPerBlock; } <\/code><\/pre>\n<p>  \u041c\u043e\u044f \u043a\u0430\u0440\u0442\u0430 \u0432\u043e\u0437\u0432\u0440\u0430\u0449\u0430\u0435\u0442 \u0447\u0438\u0441\u043b\u043e 1024. \u041a \u0441\u043e\u0436\u0430\u043b\u0435\u043d\u0438\u044e, \u044f \u043d\u0435 \u043c\u043e\u0433\u0443 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u0442\u044c \u0432\u0441\u0435 \u043f\u043e\u0442\u043e\u043a\u0438, \u043f\u043e\u0442\u043e\u043c\u0443 \u0447\u0442\u043e \u0432 \u0440\u0435\u0436\u0438\u043c\u0435 \u043e\u0442\u043b\u0430\u0434\u043a\u0438, \u043f\u0440\u0438 \u0434\u0435\u043b\u0435\u043d\u0438\u0438 <code>float<\/code> \u0447\u0438\u0441\u0435\u043b, \u044f \u043f\u043e\u043b\u0443\u0447\u0430\u044e \u043e\u0448\u0438\u0431\u043a\u0443 <code>cudaErrorLaunchOutOfResources<\/code>. \u042f \u0442\u0430\u043a \u0438 \u043d\u0435 \u0441\u043c\u043e\u0433 \u043e\u0442\u044b\u0441\u043a\u0430\u0442\u044c, \u043f\u043e\u0447\u0435\u043c\u0443 \u043c\u043d\u0435 \u043d\u0435 \u0445\u0432\u0430\u0442\u0438\u043b\u043e \u0440\u0435\u0441\u0443\u0440\u0441\u043e\u0432 (\u0438\u043b\u0438 \u0440\u0435\u0433\u0438\u0441\u0442\u0440\u043e\u0432, \u043a\u0430\u043a \u043f\u0438\u0448\u0435\u0442 \u0438\u043d\u0442\u0435\u0440\u043d\u0435\u0442), \u043f\u043e\u044d\u0442\u043e\u043c\u0443 \u0435\u0434\u0438\u043d\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0435 \u0442\u0440\u0438 \u0440\u0435\u0448\u0435\u043d\u0438\u044f, \u0447\u0442\u043e \u044f \u0441\u043c\u043e\u0433 \u043d\u0430\u0439\u0442\u0438 \u2014 \u044d\u0442\u043e \u0441\u043d\u0438\u0436\u0435\u043d\u0438\u0435 \u0447\u0438\u0441\u043b\u0430 \u043f\u043e\u0442\u043e\u043a\u043e\u0432 \u043d\u0430 \u0431\u043b\u043e\u043a, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043d\u0438\u0435 release \u0440\u0435\u0436\u0438\u043c\u0430 \u043f\u0440\u0438 \u043a\u043e\u043c\u043f\u0438\u043b\u044f\u0446\u0438\u0438, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043d\u0438\u0435 \u043e\u043f\u0446\u0438\u0438 <code>maxrregcount<\/code> \u0441 \u043b\u044e\u0431\u044b\u043c \u0447\u0438\u0441\u043b\u043e\u043c \u0431\u043e\u043b\u044c\u0448\u0435 \u043d\u0443\u043b\u044f.<br \/>  \u0415\u0441\u043b\u0438 \u0432\u044b \u0434\u043e\u0447\u0438\u0442\u0430\u043b\u0438 \u043c\u043e\u044e \u0441\u0442\u0430\u0442\u044c\u044e \u0434\u043e \u044d\u0442\u043e\u0433\u043e \u043c\u043e\u043c\u0435\u043d\u0442\u0430, \u0442\u043e \u044f \u0440\u0430\u0434, \u0447\u0442\u043e \u043d\u0435 \u0437\u0440\u044f \u0441\u0438\u0436\u0443 \u043d\u043e\u0447\u044c\u044e \u043f\u0435\u0440\u0435\u0434 \u043a\u043e\u043c\u043f\u043e\u043c. \u0410 \u0435\u0441\u043b\u0438 \u0432\u044b \u0437\u043d\u0430\u0435\u0442\u0435 \u043e\u0442\u0432\u0435\u0442, \u043a\u0430\u043a \u0443\u0437\u043d\u0430\u0442\u044c \u0432 \u0447\u0435\u043c \u0438\u043c\u0435\u043d\u043d\u043e \u0431\u044b\u043b\u0430 \u043c\u043e\u044f \u043e\u0448\u0438\u0431\u043a\u0430, \u0438 \u0441\u043e\u043e\u0431\u0449\u0438\u0442\u0435 \u043c\u043d\u0435 \u0435\u0433\u043e, \u0442\u043e \u044f \u0431\u0443\u0434\u0443 \u0440\u0430\u0434 \u0432\u0434\u0432\u043e\u0439\u043d\u0435.<\/p>\n<p>  \u0422\u0435\u043f\u0435\u0440\u044c \u0440\u0430\u0437\u043e\u0431\u044c\u0451\u043c \u043d\u0430\u0448\u0435 \u043f\u043e\u043b\u0435 \u043d\u0430 \u0431\u043b\u043e\u043a\u0438 \u0438 \u043f\u043e\u0442\u043e\u043a\u0438:  <\/p>\n<pre><code class=\"cpp\">void setupGrid(dim3* blocks, dim3* threads, int maxThreads) { \tthreads-&gt;x = 32; \tthreads-&gt;y = maxThreads \/ threads-&gt;x - 2; \/\/ to avoid cudaErrorLaunchOutOfResources error  \tblocks-&gt;x = (width + threads-&gt;x - 1) \/ threads-&gt;x; \tblocks-&gt;y = (height + threads-&gt;y - 1) \/ threads-&gt;y; }  \/\/ ... \t\/\/ Initialize Host \tdim3 blocks, threads; \tsetupGrid(&blocks, &threads, threadsCount()); \tprintf(&quot;Debug: blocks(%d, %d), threads(%d, %d)\\nCalculated Resolution: %d x %d\\n&quot;, \t\tblocks.x, blocks.y, threads.x, threads.y, blocks.x * threads.x, \t\tblocks.y * threads.y); <\/code><\/pre>\n<p>  \u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u0431\u0443\u0434\u0435\u0442 \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u043c (\u044f \u0437\u0430\u0432\u0435\u0434\u043e\u043c\u043e \u0443\u0432\u0435\u043b\u0438\u0447\u0438\u043b \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u0431\u043b\u043e\u043a\u043e\u0432 \u043d\u0430 \u0435\u0434\u0438\u043d\u0438\u0446\u0443 \u0432 \u0442\u0435\u0445 \u0441\u043b\u0443\u0447\u0430\u044f\u0445, \u043a\u043e\u0433\u0434\u0430 \u0440\u0430\u0437\u043c\u0435\u0440 \u044d\u043a\u0440\u0430\u043d\u0430 \u043d\u0435 \u043a\u0440\u0430\u0442\u0435\u043d \u0447\u0438\u0441\u043b\u0443 \u043f\u043e\u0442\u043e\u043a\u043e\u0432):  <\/p>\n<pre><code>Debug: blocks(40, 24), threads(32, 30) Calculated Resolution: 1280 x 720 <\/code><\/pre>\n<p>  <\/p>\n<h4>\u0422\u0435\u043e\u0440\u0438\u044f \u0438 \u043f\u0440\u0430\u043a\u0442\u0438\u043a\u0430<\/h4>\n<p>  \u0417\u0430\u0440\u044f\u0434\u044b, \u0440\u0430\u0441\u0441\u0442\u043e\u044f\u043d\u0438\u044f, \u043a\u0443\u043b\u043e\u043d\u044b \u0438 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u044b \u2014 \u044d\u0442\u043e \u0437\u0432\u0443\u0447\u0438\u0442 \u0445\u043e\u0440\u043e\u0448\u043e. \u041d\u043e \u043a\u043e\u0433\u0434\u0430 \u0434\u0435\u043b\u043e \u0434\u043e\u0445\u043e\u0434\u0438\u0442 \u0434\u043e \u043c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u044f, \u0442\u043e \u0442\u0443\u0442 \u043c\u044b \u043d\u0435 \u043c\u043e\u0436\u0435\u043c \u0443\u0436\u0435 \u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u0442\u044c\u0441\u044f \u0441\u0443\u043f\u0435\u0440-\u043c\u0430\u043b\u0435\u043d\u044c\u043a\u0438\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438, \u043f\u043e\u0442\u043e\u043c\u0443 \u0447\u0442\u043e \u043e\u043d\u0438 \u043c\u043e\u0433\u0443\u0442 \u043f\u0440\u043e\u043f\u0430\u0441\u0442\u044c \u0432 \u043f\u0440\u043e\u0446\u0435\u0441\u0441\u0435 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f. \u041f\u043e\u044d\u0442\u043e\u043c\u0443 \u043c\u044b \u0443\u043f\u0440\u043e\u0441\u0442\u0438\u043c \u0437\u0430\u0434\u0430\u0447\u0443:<\/p>\n<p>  \u0417\u0430\u0440\u044f\u0434 \u2014 \u044d\u0442\u043e \u0441\u0442\u0440\u0443\u043a\u0442\u0443\u0440\u0430 \u0441 \u043a\u043e\u043e\u0440\u0434\u0438\u043d\u0430\u0442\u0430\u043c\u0438 (\u043a\u0430\u043a \u043d\u0430 \u044d\u043a\u0440\u0430\u043d\u0435).  <\/p>\n<pre><code class=\"cpp\">struct charge { \tint x, y, q; }; <\/code><\/pre>\n<p>  \u041d\u0430\u043f\u0440\u044f\u0436\u0451\u043d\u043d\u043e\u0441\u0442\u044c \u044d\u043b\u0435\u043a\u0442\u0440\u0438\u0447\u0435\u0441\u043a\u043e\u0433\u043e \u043f\u043e\u043b\u044f \u0432 \u0442\u043e\u0447\u043a\u0435 \u2014 \u044d\u0442\u043e \u0432\u0435\u043a\u0442\u043e\u0440, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043c\u044b \u0432\u044b\u0440\u0430\u0437\u0438\u043c \u043f\u0440\u043e\u0435\u043a\u0446\u0438\u044f\u043c\u0438 \u0438 \u0444\u0443\u043d\u043a\u0446\u0438\u0438, \u0432\u044b\u0447\u0438\u0441\u043b\u044f\u044e\u0449\u0435\u0439 \u0438\u0445 (\u0432\u043e\u0442 \u0438\u043c\u0435\u043d\u043d\u043e \u0437\u0434\u0435\u0441\u044c \u044f \u0438 \u043f\u043e\u043b\u0443\u0447\u0438\u043b \u0441\u0432\u043e\u044e \u043e\u0448\u0438\u0431\u043a\u0443):  <\/p>\n<pre><code class=\"cpp\">const float k = 50.0f; const float minDistance = 0.9f; \/\/ not to divide by zero const float maxForce = 1e7f;  struct force { \tfloat fx, fy;  \t__device__ force() : \t\t\tfx(0.0f), fy(0.0f) { \t}  \t__device__ force(int fx, int fy) : \t\t\tfx(fx), fy(fy) { \t}  \t__device__ float length2() const { \t\treturn (fx * fx + fy * fy); \t}  \t__device__ float length() const { \t\treturn sqrtf(length2()); \t}  \t__device__ void calculate(const charge& q, int probe_x, int probe_y) { \t\t\/\/ F(1-&gt;2) = k * q1 * q2 \/ r(1-&gt;2)^2 * vec_r(1-&gt;2) \/ abs(vec_r(1-&gt;2)) \t\t\/\/ e = vec_F \/ q2 \t\tfx = probe_x - q.x; \t\tfy = probe_y - q.y;  \t\tfloat l = length(); \t\tif (l &lt;= minDistance) { \t\t\treturn; \t\t}  \t\tfloat e = k * q.q \/ (l * l * l); \t\tif (e &gt; maxForce) { \t\t\tfx = fy = maxForce; \t\t} else { \t\t\tfx *= e; \t\t\tfy *= e; \t\t} \t}  \t__device__ force operator +(const force& f) const { \t\treturn force(fx + f.fx, fy + f.fy); \t}  \t__device__ force operator -(const force& f) const { \t\treturn force(fx - f.fx, fy - f.fy); \t}  \t__device__ force& operator +=(const force& f) { \t\tfx += f.fx; \t\tfy += f.fy; \t\treturn *this; \t}  \t__device__ force& operator -=(const force& f) { \t\tfx -= f.fx; \t\tfy -= f.fy; \t\treturn *this; \t} }; <\/code><\/pre>\n<p>  \u0417\u0430\u0434\u0443\u043c\u043a\u0430 \u043c\u043e\u044f \u043f\u0440\u043e\u0441\u0442\u0430\u044f \u2014 \u044f \u0441\u043e\u0437\u0434\u0430\u043c \u0434\u0432\u0443\u0445\u043c\u0435\u0440\u043d\u044b\u0439 \u043c\u0430\u0441\u0441\u0438\u0432 \u0438\u0437 \u0432\u0435\u043a\u0442\u043e\u0440\u043e\u0432 \u043d\u0430\u043f\u0440\u044f\u0436\u0451\u043d\u043d\u043e\u0441\u0442\u0435\u0439. \u0412 \u043a\u0430\u0436\u0434\u043e\u0439 \u0442\u043e\u0447\u043a\u0435 \u044d\u043a\u0440\u0430\u043d\u0430, \u043d\u0430\u043f\u0440\u044f\u0436\u0451\u043d\u043d\u043e\u0441\u0442\u044c \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0432\u043d\u0430 \u0441\u0443\u043c\u043c\u0435 \u0432\u0435\u043a\u0442\u043e\u0440\u043e\u0432 \u043a\u0443\u043b\u043e\u043d\u043e\u0432\u0441\u043a\u0438\u0445 \u0441\u0438\u043b. \u0417\u0432\u0443\u0447\u0438\u0442 \u0437\u0430\u0443\u043c\u043d\u043e, \u043d\u043e \u044f \u0438\u043d\u0430\u0447\u0435 \u043d\u0435 \u0441\u043c\u043e\u0433 \u0441\u0444\u043e\u0440\u043c\u0443\u043b\u0438\u0440\u043e\u0432\u0430\u0442\u044c. \u041f\u0440\u043e\u0449\u0435 \u0441\u043a\u0430\u0437\u0430\u0442\u044c, \u043a\u043e\u0434 \u0431\u0443\u0434\u0435\u0442 \u0432\u043e\u0442 \u0442\u0430\u043a\u0438\u043c (\u0432\u044b \u0435\u0433\u043e \u043f\u0440\u043e\u0447\u0442\u0438\u0442\u0435, \u043d\u043e \u043d\u0435 \u043f\u0438\u0448\u0438\u0442\u0435 \u0435\u0449\u0451 \u0432 \u0441\u0432\u043e\u0439 \u043f\u0440\u043e\u0435\u043a\u0442):  <\/p>\n<pre><code class=\"cpp\">\tforce temp_f; \tfor (int i = 0; i &lt; chargeCount; i++) { \t\ttemp_f.calculate(charges[i], x, y); \t\t*f += temp_f; \t} <\/code><\/pre>\n<p>  \u0415\u0441\u043b\u0438 \u043c\u044b \u043f\u043e\u043f\u0430\u0434\u0451\u043c \u0432 \u0437\u0430\u0440\u044f\u0434, \u0442\u043e \u044f \u043f\u0440\u043e\u043f\u0443\u0441\u043a\u0430\u044e \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u0435 \u043d\u0430\u043f\u0440\u044f\u0436\u0435\u043d\u043d\u043e\u0441\u0442\u0438 \u0432 \u044d\u0442\u043e\u0439 \u0442\u043e\u0447\u043a\u0435 (\u0438\u0431\u043e \u0434\u0435\u043b\u0438\u0442\u044c \u043d\u0430 \u043d\u043e\u043b\u044c \u043d\u0435\u043b\u044c\u0437\u044f). \u0422\u0430\u043a \u0447\u0442\u043e \u043d\u0430 \u044d\u043a\u0440\u0430\u043d\u0435 \u0432\u044b \u0434\u043e\u043b\u0436\u043d\u044b \u0443\u0432\u0438\u0434\u0435\u0442\u044c \u0447\u0451\u0440\u043d\u0443\u044e \u0442\u043e\u0447\u043a\u0443 \u2014 \u0437\u0430\u0440\u044f\u0434.<\/p>\n<p>  \u0422\u0435\u043f\u0435\u0440\u044c, \u043a\u043e\u0433\u0434\u0430 \u043c\u044b \u0441\u043e\u0437\u0434\u0430\u043b\u0438 \u0438\u043d\u0444\u0440\u0430\u0441\u0442\u0440\u0443\u043a\u0442\u0443\u0440\u0443 \u2014 \u0434\u0430\u0432\u0430\u0439\u0442\u0435 \u043a\u043e\u0434\u0438\u0442\u044c \u043f\u043e\u0434 CUDA. \u0414\u043b\u044f \u043d\u0430\u0447\u0430\u043b\u0430 \u0441\u043e\u0437\u0434\u0430\u0434\u0438\u043c \u043d\u0435\u043e\u0431\u0445\u043e\u0434\u0438\u043c\u044b\u0435 \u043f\u0435\u0440\u0435\u043c\u0435\u043d\u043d\u044b\u0435 \u0445\u043e\u0441\u0442\u0430 \u0438 \u0443\u0441\u0442\u0440\u043e\u0439\u0441\u0442\u0432\u0430 \u0438 \u043e\u0441\u0432\u043e\u0431\u043e\u0434\u0438\u043c \u0438\u0445!  <\/p>\n<h4>\u0421\u043e\u0437\u0434\u0430\u0435\u043c \u0438 \u043e\u0441\u0432\u043e\u0431\u043e\u0436\u0434\u0430\u0435\u043c \u0440\u0435\u0441\u0443\u0440\u0441\u044b<\/h4>\n<p>  \u041d\u0430\u043c \u043f\u043e\u043d\u0430\u0434\u043e\u0431\u044f\u0442\u0441\u044f \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0435 \u043f\u0435\u0440\u0435\u043c\u0435\u043d\u043d\u044b\u0435:  <\/p>\n<ul>\n<li>\u044d\u043a\u0440\u0430\u043d \u2014 \u0445\u043e\u0441\u0442 \u2014 \u0446\u0432\u0435\u0442\u043d\u044b\u0435 \u043f\u0438\u043a\u0441\u0435\u043b\u0438;<\/li>\n<li>\u044d\u043a\u0440\u0430\u043d \u2014 \u0443\u0441\u0442\u0440\u043e\u0439\u0441\u0442\u0432\u043e \u2014 \u043e\u043d\u0438 \u0436\u0435;<\/li>\n<li>\u0437\u0430\u0440\u044f\u0434\u044b \u2014 \u0443\u0441\u0442\u0440\u043e\u0439\u0441\u0442\u0432\u043e \u2014 \u043a\u043e\u043e\u0440\u0434\u0438\u043d\u0430\u0442\u044b \u0438 \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u044b \u0437\u0430\u0440\u044f\u0434\u043e\u0432 \u0432 \u043f\u043e\u0441\u0442\u043e\u044f\u043d\u043d\u043e\u0439 \u043f\u0430\u043c\u044f\u0442\u0438 \u0443\u0441\u0442\u0440\u043e\u0439\u0441\u0442\u0432\u0430.<\/li>\n<\/ul>\n<p>  \u0421\u043e\u0437\u0434\u0430\u0451\u043c:  <\/p>\n<pre><code class=\"cpp\">uchar4* screen = NULL;  \/\/ ...  \tscreen = (uchar4*) malloc(width * height * sizeof(uchar4)); \tmemset(screen, 0, width * height * sizeof(uchar4));  \tuchar4 *dev_screen = NULL; \tcudaMalloc((void**) &dev_screen, width * height * sizeof(uchar4)); \tcudaMemset(dev_screen, 0, width * height * sizeof(uchar4));  \/\/ ...  \t\/\/ Free resources \tfree(screen); \tscreen = NULL;  \tcudaFree(dev_screen); \tdev_screen = NULL; <\/code><\/pre>\n<p>  \u0417\u0430\u0440\u044f\u0434\u044b \u043f\u0440\u043e\u0438\u043d\u0438\u0446\u0438\u0430\u043b\u0438\u0437\u0438\u0440\u0443\u0435\u043c \u043f\u0440\u043e\u0441\u0442\u043e. \u042f \u043d\u0435 \u0441\u0431\u0440\u0430\u0441\u044b\u0432\u0430\u044e \u0433\u0435\u043d\u0435\u0440\u0430\u0442\u043e\u0440 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0445 \u0447\u0438\u0441\u0435\u043b, \u0447\u0442\u043e\u0431\u044b \u043f\u043e\u043b\u0443\u0447\u0438\u0442\u044c \u043e\u0434\u043d\u0443 \u0438 \u0442\u0443 \u0436\u0435 \u043a\u0430\u0440\u0442\u0438\u043d\u0443. \u0415\u0441\u043b\u0438 \u0432\u044b \u0445\u043e\u0442\u0438\u0442\u0435 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u043e\u0435 \u043f\u043e\u043b\u0435 \u043a\u0430\u0436\u0434\u044b\u0439 \u0440\u0430\u0437, \u0442\u043e \u043c\u043e\u0436\u0435\u0442\u0435 \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u0432\u044b\u0437\u043e\u0432 <code>srand<\/code>:  <\/p>\n<pre><code class=\"cpp\">const int chargeCount = 10; __constant__ charge dev_charges[chargeCount]; const int maxCharge = 1000; const int minCharge = -1000;  \/\/ ...  void prepareCharges(void) { \tcharge* charges = (charge*) malloc(chargeCount * sizeof(charge)); \tfor (int i = 0; i &lt; chargeCount; i++) { \t\tcharges[i].x = rand() % width; \t\tcharges[i].y = rand() % height; \t\tcharges[i].q = rand() % (maxCharge - minCharge) + minCharge;  \t\tprintf(&quot;Debug: Charge #%d (%d, %d, %d)\\n&quot;, i, charges[i].x, \t\t\t\tcharges[i].y, charges[i].q); \t} \tcudaMemcpyToSymbol(dev_charges, charges, chargeCount * sizeof(charge)); }  \/\/ ... \tprepareCharges(); <\/code><\/pre>\n<p>  \u041f\u043e\u043c\u043d\u0438\u0442\u0435 \u0437\u0430\u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0438\u0440\u043e\u0432\u0430\u043d\u043d\u0443\u044e \u0441\u0442\u0440\u043e\u043a\u0443, \u0440\u0435\u043d\u0434\u0435\u0440\u044f\u0449\u0443\u044e \u044d\u043a\u0440\u0430\u043d? \u0422\u0435\u043f\u0435\u0440\u044c \u0432\u044b \u043c\u043e\u0436\u0435\u0442\u0435 \u0435\u0451 \u0440\u0430\u0441\u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0438\u0440\u043e\u0432\u0430\u0442\u044c:  <\/p>\n<pre><code class=\"cpp\">\tglDrawPixels(width, height, GL_RGBA, GL_UNSIGNED_BYTE, screen); <\/code><\/pre>\n<p>  <\/p>\n<h4>\u0421\u043e\u0437\u0434\u0430\u0434\u0438\u043c \u044f\u0434\u0440\u043e<\/h4>\n<p>  \u0422\u0435\u043f\u0435\u0440\u044c \u0441\u043e\u0437\u0434\u0430\u0434\u0438\u043c \u0444\u0443\u043d\u043a\u0446\u0438\u044e, \u043a\u043e\u0442\u043e\u0440\u0443\u044e \u043c\u044b \u0431\u0443\u0434\u0435\u043c \u0432\u044b\u0437\u044b\u0432\u0430\u0442http:\/\/habrahabr.ru\/sandbox\/add\/#\u044c \u0432 \u043a\u0430\u0436\u0434\u043e\u043c \u043f\u043e\u0442\u043e\u043a\u0435, \u0447\u0442\u043e\u0431\u044b \u0441\u0440\u0435\u043d\u0434\u0435\u0440\u0438\u0442\u044c \u043a\u0430\u0436\u0434\u044b\u0439 \u043f\u0438\u043a\u0441\u0435\u043b\u044c \u043d\u0430 \u044d\u043a\u0440\u0430\u043d\u0435. \u042d\u0442\u0430 \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0431\u0443\u0434\u0435\u0442 \u0432\u044b\u043f\u043e\u043b\u043d\u044f\u0442\u044c \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0435\u0435: \u043e\u043d\u0430 \u0432\u044b\u0447\u0438\u0441\u043b\u0438\u0442 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043d\u0430\u043f\u0440\u044f\u0436\u0451\u043d\u043d\u043e\u0441\u0442\u0438 \u043f\u043e\u043b\u044f, \u0438 \u0432\u044b\u0437\u043e\u0432\u0435\u0442 \u0434\u0440\u0443\u0433\u0443\u044e \u0444\u0443\u043d\u043a\u0446\u0438\u044e, \u0447\u0442\u043e\u0431\u044b \u0441\u043e\u043f\u043e\u0441\u0442\u0430\u0432\u0438\u0442\u044c \u043a\u0430\u043a\u043e\u0439-\u043d\u0438\u0431\u0443\u0434\u044c \u0446\u0432\u0435\u0442 \u0442\u043e\u0447\u043a\u0435.<\/p>\n<pre><code class=\"cpp\">__device__ uchar4 getColor(const force& f) { \tuchar4 color; \tcolor.x = color.y = color.z = color.w = 0;  \tfloat l = f.length(); \tcolor.x = (l &gt; maxLengthForColor ? 255 : l * 256 \/ maxLengthForColor);  \treturn color; }  __global__ void renderFrame(uchar4* screen) { \t\/\/ build the field and render the frame } <\/code><\/pre>\n<p>  \u0427\u0442\u043e\u0431\u044b \u0432\u044b\u0437\u0432\u0430\u0442\u044c \u044f\u0434\u0440\u043e\u2026 \u0430 \u0432\u043f\u0440\u043e\u0447\u0435\u043c \u0438 \u0442\u0430\u043a \u0432\u0441\u0451 \u044f\u0441\u043d\u043e:  <\/p>\n<pre><code class=\"cpp\">\t\/\/ Launch Kernel to render the image \trenderFrame&lt;&lt;&lt;blocks, threads&gt;&gt;&gt;(dev_screen);  \tcudaMemcpy(screen, dev_screen, width * height * sizeof(uchar4), \t\t\tcudaMemcpyDeviceToHost); <\/code><\/pre>\n<p>  \u0422\u0435\u043f\u0435\u0440\u044c \u0434\u043e\u0431\u0430\u0432\u0438\u043c \u043d\u0435\u043c\u043d\u043e\u0433\u043e \u043a\u043e\u0434\u0430 \u0432 <code>renderFrame<\/code> \u0438 \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u0442\u043e, \u0447\u0442\u043e \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u043e:  <\/p>\n<pre><code class=\"cpp\">__global__ void renderFrame(uchar4* screen) { \tint x = blockIdx.x * blockDim.x + threadIdx.x; \tint y = blockIdx.y * blockDim.y + threadIdx.y;  \tforce f, temp_f; \tfor (int i = 0; i &lt; chargeCount; i++) { \t\ttemp_f.calculate(dev_charges[i], x, y); \t\tf += temp_f; \t}  \tscreen[x + y * width] = getColor(f); } <\/code><\/pre>\n<p>  <\/p>\n<h4>\u0412\u044b\u0432\u043e\u0434<\/h4>\n<p>  \u0427\u0430\u0441\u0442\u043e \u043d\u0430\u0440\u043e\u0434 \u043f\u0440\u043e\u043f\u0443\u0441\u043a\u0430\u0435\u0442 \u0441\u0442\u0430\u0442\u044c\u044e, \u043f\u043e\u0442\u043e\u043c\u0443 \u0447\u0442\u043e \u0432\u0430\u0436\u0435\u043d \u0432\u044b\u0432\u043e\u0434. \u041c\u043d\u0435 \u0437\u0434\u0435\u0441\u044c \u043d\u0435\u0447\u0435\u0433\u043e \u0441\u043a\u0430\u0437\u0430\u0442\u044c, \u043a\u0440\u043e\u043c\u0435 \u043a\u0430\u043a \u043e\u043f\u0443\u0431\u043b\u0438\u043a\u043e\u0432\u0430\u0442\u044c \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0443, \u0447\u0442\u043e \u0443 \u043c\u0435\u043d\u044f \u043f\u043e\u043b\u0443\u0447\u0438\u043b\u0430\u0441\u044c. \u0421\u0441\u044b\u043b\u043a\u0443 \u043d\u0430 \u043f\u043e\u0447\u0442\u0438 \u043f\u0443\u0441\u0442\u043e\u0439 \u0440\u0435\u043f\u043e\u0437\u0438\u0442\u043e\u0440\u0438\u0439 \u043d\u0430\u0439\u0434\u0451\u0442\u0435 <a href=\"https:\/\/github.com\/Vest\/ElectroField\">\u0437\u0434\u0435\u0441\u044c<\/a>.<br \/>  <img decoding=\"async\" src=\"http:\/\/habrastorage.org\/storage2\/3cd\/e28\/d58\/3cde28d5846677bd7ee66a0bb20d0730.png\"\/><br \/>  \u0412 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435, \u0435\u0441\u0442\u044c \u043a\u0443\u0434\u0430 \u0440\u0430\u0437\u0432\u0438\u0432\u0430\u0442\u044c\u0441\u044f \u2014 \u043c\u043e\u0436\u043d\u043e \u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440 \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u0430\u043d\u0438\u043c\u0430\u0446\u0438\u044e, \u043f\u043e\u043a\u0430\u0437\u044b\u0432\u0430\u044e\u0449\u0443\u044e \u043a\u0430\u043a \u043c\u0435\u043d\u044f\u0435\u0442\u0441\u044f \u043f\u043e\u043b\u0435 \u043f\u0440\u0438 \u043f\u0435\u0440\u0435\u043c\u0435\u0449\u0435\u043d\u0438\u0438 \u0437\u0430\u0440\u044f\u0434\u043e\u0432. \u041d\u043e \u043e\u043f\u044f\u0442\u044c \u0436\u0435, \u0435\u0441\u043b\u0438 \u044d\u0442\u043e \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u043e. \t\t\t \t\t\t<\/p>\n<div class=\"clear\"><\/div>\n<\/p><\/div>\n<p> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"http:\/\/habrahabr.ru\/post\/164499\/\"> http:\/\/habrahabr.ru\/post\/164499\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<div class=\"content html_format\"> \t\t\t\u0414\u0430\u043d\u043d\u0430\u044f \u0441\u0442\u0430\u0442\u044c\u044f \u043d\u0430\u043f\u0438\u0441\u0430\u043d\u0430 \u0441 \u0446\u0435\u043b\u044c\u044e \u043f\u0440\u043e\u0434\u0435\u043c\u043e\u043d\u0441\u0442\u0440\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u043a\u0430\u043a \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0442\u0435\u0445\u043d\u043e\u043b\u043e\u0433\u0438\u0438 CUDA \u043c\u043e\u0436\u043d\u043e \u0441\u043c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u043f\u0440\u043e\u0441\u0442\u043e\u0435 \u0432\u0437\u0430\u0438\u043c\u043e\u0434\u0435\u0439\u0441\u0442\u0432\u0438\u0435 \u0437\u0430\u0440\u044f\u0436\u0435\u043d\u044b\u0445 \u0447\u0430\u0441\u0442\u0438\u0446 (\u0441\u043c. <a href=\"http:\/\/en.wikipedia.org\/wiki\/Coulomb%27s_law\">\u0417\u0430\u043a\u043e\u043d \u041a\u0443\u043b\u043e\u043d\u0430<\/a>). \u0414\u043b\u044f \u0432\u044b\u0432\u043e\u0434\u0430 \u0441\u0442\u0430\u0442\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438 \u044f \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043b \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0443 <b>freeglut<\/b>.<br \/>  \u041a\u0430\u043a \u043f\u0438\u0448\u0443\u0442 \u0447\u0430\u0441\u0442\u0435\u043d\u044c\u043a\u043e \u043d\u0430 \u0425\u0430\u0431\u0440\u0435: <\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[],"tags":[],"class_list":["post-164499","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/164499","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=164499"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/164499\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=164499"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=164499"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=164499"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}