Hello, According to PETSc documentation, all processes that share the matrix MUST call this routine for the parallel case. However, it is possible that, after the domain decomposition, some subdomains wouldn’t have Dirichlet boundary conditions. Hence, MatZeroRows wouldn’t be called by all processes. Is there a standard (or proposed) approach for this scenario?
Thank you, -- Felipe M Wallner Giacomelli
