"Daralagodu Dattatreya Jois, Sathwik Bharadw" <[email protected]> writes:
> I am using AIJ matrix to solve Laplace problem in finite element > framework. To apply Neumann boundary conditions I need to obtain > values of first and last few columns and subtract it with the > corresponding right hand side vector. I understand that > MatGetColumnVector and MatGetValues are not collective. Is there any > other alternative petsc calls to achieve this in parallel? Inhomogeneous Neumann conditions are normally applied by integration over the boundary. But if you need a linear combination of some columns of the matrix, that's a MatMult.
signature.asc
Description: PGP signature
