[petsc-users] MatZeroRows question

[Smith, Barry F.]

  They call the function and provide zero rows


> On Dec 6, 2017, at 10:52 AM, Felipe Giacomelli <fe.wallner at gmail.com> wrote:
> 
> 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