<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class=""><br class=""></div> You need to do a mapping of your global numbering to the standard PETSc numbering and use the PETSc numbering for all access to vectors and matrices. <div class=""><br class=""></div><div class=""> <a href="https://petsc.org/release/docs/manualpages/AO/AOCreate.html" class="">https://petsc.org/release/docs/manualpages/AO/AOCreate.html</a> provides one approach to managing the renumbering.</div><div class=""><br class=""></div><div class=""> Barry</div><div class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Mar 9, 2022, at 3:42 PM, Nicolás Barnafi <<a href="mailto:nabw91@gmail.com" class="">nabw91@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hi community, <div class=""><br class=""></div><div class="">I have an application with polytopal meshes (elements of arbitrary shape) where the distribution of dofs is not PETSc friendly, meaning that it is not true that cpu0 owns dofs [0,a), then cpu1 owns [a,b) and so on, but instead the distribution is in fact random. Another important detail is that boundary dofs are shared, meaning that if dof 150 is on the boundary, each subdomain vector has dof 150. </div><div class=""><br class=""></div><div class="">Under this considerations: </div><div class=""><br class=""></div><div class="">i) Is it possible to give an arbitrary mapping to the matrix structure or is the blocked distribution hard coded? </div><div class="">ii) Are the repeated boundary dofs an issue when computing a Fieldsplit preconditioner in parallel?</div><div class=""><br class=""></div><div class="">Best regards,<br clear="all" class=""><div class="">Nicolas</div></div></div>
</div></blockquote></div><br class=""></div></body></html>