Re: [petsc-users] MatZeroRowsColumns eliminates incorrectly in parallel

This is incorrect and will not work. for (j = 0; j < daInfo.my; ++j) { dirichlet_rows[n_dirichlet_rows++] = j * daInfo.mx * ncom

[petsc-users] MatZeroRowsColumns eliminates incorrectly in parallel

Dear all, I tried to use MatZeroRowsColumns to eliminate Dirichlet boundary nodes. However, it cannot eliminate correctly in parallel. Please see the attached code which uses DMDA to create the matrix. When I used one process, it works as expected. For two processes, the domain is split in the x