Hi Matt,<br>Thank you very much for the reply. Now, I got the Neumann part. But I am still a bit confused about the Dirichlet part. Please see the following quote.<br><br>Yan<br><br><div class="gmail_quote"><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="gmail_quote"><div class="im"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
For the Dirichlet B.C.s, I did not understand the coefficients <a href="http://www.physicsforums.com/showthread.php?t=51412" target="_blank"><i></i></a>given below. Isn't correct to set the v[0]=Hx*Hy here?<br>
<pre width="80"><font size="4">if (i==0 || j==0 || i==mx-1 || j==my-1) {</font><br></pre>if (user->bcType == DIRICHLET) {<br> v[0] = 2.0*rho*(HxdHy + HydHx);<br> }<br>}<br></blockquote></div><div><br>This is the proper scaling.<br>
</div></div></blockquote><div>the Dirichlet B.C.s:<br>Which scaling do you think is proper, "v[0]=Hx*Hy" or 2.0*rho*(HxdHy + HydHx)? <br>If it is 2.0*rho*(HxdHy + HydHx), can you say a little bit more about why is this one? I only see a factor of Hx*Hy when we set up the RHS. Did I miss something?<br>
<br><br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="gmail_quote"><div> Matt<br> </div><br></div><font color="#888888"><br>
<br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>
</font></blockquote></div><br>