mlp_lpvect.html

国外MPI教材

<h2> C Version of Laplace Program (vector) </h2>   <pre><br /><br />#include &lt;stdio.h&gt;<br />#include &lt;stdlib.h&gt;<br />#include &lt;math.h&gt;<br /></pre><pre>#define PI 3.14159<br /><br />main(int argc, char* argv[]) {<br />        const int imax=2001,jmax=2001;<br />        const int itmax=100;<br />        int it;<br />        double u[imax][jmax], du[imax][jmax],dumax;<br />        const double umax=10.0;<br /><br />        for(int j=0;j&lt;jmax;++j) {<br />          for(int i=0;i&lt;(imax-1);++i) {<br />            u[i][j]=0.0;<br />           du[i][j]=0.0;<br />          }<br />          u[imax-1][j]=umax*sin(PI*(double)j/(double)jmax-1);<br />        }<br /><br />        // Main Computation Loop<br />        for(it=0;it&lt;itmax;++it) {<br />          dumax=0.0; <br /><br />          for (int j=1;j&lt;(jmax-1);++j) {<br />            for (int i=1;i&lt;(imax-1);++i) {<br />              du[i][j]=0.25*(u[i-1][j]+u[i+1][j]+u[i][j-1]+u[i][j+1])<br />                             -u[i][j];<br />            }<br />          }<br />           <br />          for (int j=1;j&lt;(jmax-1);++j) {<br />            for (int i=1;i&lt;(imax-1);++i) {<br />               dumax=dumax &lt; fabs(du[i][j]) ? fabs(du[i][j]) : dumax;<br />            }<br />          }<br /><br />          for (int j=1;j&lt;(jmax-1);++j) {<br />            for (int i=1;i&lt;(imax-1);++i) {<br />               u[i][j]=u[i][j]+du[i][j];<br />            }<br />          }<br />     <br />          printf(&quot;%d %g\n&quot;,it,dumax);<br />        }<br />}<br /></pre>