ksp ex29.c B.C.s and Forcing terms.

[Ryan Yan]

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>