Hi! 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.
The boundary conditions are fulfilled using immersed boundary method ghost cells at the solid boundaries of the domain. I read about the DAs to deal with the boundary conditions, but they are applicable to only "IJK" matrices as per my understanding. 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? Thanks a lot. -- Arun Kamath PhD Stipendiat NTNU, Trondheim -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121107/2290453f/attachment.html>
