<pre style="white-space:pre-wrap;color:rgb(0,0,0)"><pre style="white-space:pre-wrap">My previous post gave a completely wrong meaning I guess, with a missing punctuation.</pre><pre style="white-space:pre-wrap">Matt, you got it right. I wanted to know how do I deal with the ghost cells arising from the physical domain, when I am using an unstructured grid. Since DA works with structured grids only, as far as I understood.</pre>
<pre style="white-space:pre-wrap">Thanks.</pre><pre style="white-space:pre-wrap"><br></pre><pre style="white-space:pre-wrap">Are you saying that you have an unstructured grid with ghost cells, and you
would like PETSc to produce
a Vec of the correct size and Mat with the correct preallocation for that
grid and discretization, just as DA
does for structured grids?

  Thanks,

      Matt</pre></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><br></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)">On Wed, Nov 7, 2012 at 9:31 AM, Arun Kamath <<a href="https://lists.mcs.anl.gov/mailman/listinfo/petsc-users">arun.kamath.m at gmail.com</a>> wrote:

><i> Hi!
</i>><i> I am working with a CFD code  where I need to solve the Poisson Pressure
</i>><i> equation. We solve a two phase flow problem which results in a jump in the
</i>><i> coefficients and a non symmetric system.
</i>><i>
</i>><i> The boundary conditions are fulfilled using immersed boundary method ghost
</i>><i> cells at the solid boundaries of the domain.
</i>><i> I read about the DAs to deal with the boundary conditions, but they are
</i>><i> applicable to only "IJK" matrices as per my understanding.
</i>><i> So when I am using the "Aij" (CSR) matrix, how would you suggest I deal
</i>><i> with the ghost cells of the physical domain while forming the matrix and
</i>><i> the vector?
</i>><i>
</i>



><i> Thanks a lot.
</i>><i> --
</i>><i> Arun Kamath
</i>><i> PhD Stipendiat
</i>><i> NTNU, Trondheim
</i>><i>
</i>></pre><div><br></div>-- <br><div dir="ltr">Arun Kamath <div>PhD Stipendiat</div><div>NTNU, Trondheim</div></div><br>