an>
0049   <span class="comment">%  pi += rNu*rTmp;</span>
0050   <span class="comment">%  printf(&quot;pi: %e, rTmp: %e\n&quot;, pi, rTmp);</span>
0051   <span class="comment">%}</span>
0052   
0053   <span class="comment">%if(Uequal) {</span>
0054   <span class="comment">%  v_copy(Uequal, u_it);</span>
0055   <span class="comment">%  for(i=0; i&lt;Uequal-&gt;dim; i++)</span>
0056   <span class="comment">%    u_d-&gt;ve[i] = d-&gt;ve[X-&gt;dim+i];</span>
0057   <span class="comment">%    // smer pro multip