[petsc-users] Petsc has generated inconsistent data - Divide by zero!

[TAY wee-beng]

Hi Jed,

I think I understand what you mean. Are you saying that I have 2 
separate nxn matrix to solve.

I combine them to form a 2n x 2n matrix and solve. I think it will also 
improve the scalability with regards to MPI too.

Is that so?

Yours sincerely,

TAY wee-beng


On 10/7/2011 1:02 AM, Jed Brown wrote:
> On Sat, Jul 9, 2011 at 17:57, TAY wee-beng <zonexo at gmail.com 
> <mailto:zonexo at gmail.com>> wrote:
>
>     Ok thanks. Btw, is it better to solve my x/y eqns together? How
>     can I do that?
>
>
> Sure, interlace the x and y degrees of freedom and assemble the matrix 
> with x equations in the even row/column numbers and y equations in the 
> odd row/column numbers.
>
> Note that the equations are not decoupled for some boundary conditions 
> (e.g. slip on a curved surface or free surface).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110710/b3312c9a/attachment.htm>