On Wed, Nov 7, 2012 at 9:31 AM, Arun Kamath <span dir="ltr"><<a href="mailto:arun.kamath.m@gmail.com" target="_blank">arun.kamath.m@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi!<div>I am working with a CFD code  where I need to solve the Poisson Pressure equation. We solve a two phase flow problem which results in a jump in the coefficients and a non symmetric system.</div><div><br></div><div>

The boundary conditions are fulfilled using immersed boundary method ghost cells at the solid boundaries of the domain.</div><div>I read about the DAs to deal with the boundary conditions, but they are applicable to only "IJK" matrices as per my understanding.</div>

<div>So when I am using the "Aij" (CSR) matrix, how would you suggest I deal with the ghost cells of the physical domain while forming the matrix and the vector? </div></blockquote><div><br></div><div>Are you saying that you have an unstructured grid with ghost cells, and you would like PETSc to produce</div>
<div>a Vec of the correct size and Mat with the correct preallocation for that grid and discretization, just as DA</div><div>does for structured grids?</div><div><br></div><div>  Thanks,</div><div><br></div><div>      Matt</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>Thanks a lot.</div><span class="HOEnZb"><font color="#888888">
-- <br><div dir="ltr">Arun Kamath <div>PhD Stipendiat</div><div>NTNU, Trondheim</div></div><br>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <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>
</div>