Hi Matt, 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.
Yan For the Dirichlet B.C.s, I did not understand the coefficients **<http://www.physicsforums.com/showthread.php?t=51412>given >> below. Isn't correct to set the v[0]=Hx*Hy here? >> >> if (i==0 || j==0 || i==mx-1 || j==my-1) { >> >> if (user->bcType == DIRICHLET) { >> v[0] = 2.0*rho*(HxdHy + HydHx); >> } >> } >> > > This is the proper scaling. > > the Dirichlet B.C.s: Which scaling do you think is proper, "v[0]=Hx*Hy" or 2.0*rho*(HxdHy + HydHx)? 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? Matt > > > > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20091211/0d316c6d/attachment.htm>