
<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8" /></head><body style='font-size: 10pt; font-family: Verdana,Geneva,sans-serif'>

<p>Dear all, <br /><br />I hope this email finds you well. <br />We are currently working on a solver which is employing DMDA with SNES. The jacobian is computed via FDColoring, ie: </p>

<p>call DMDACreate1D(PETSC_COMM_WORLD, DM_BOUNDARY_GHOSTED, NY, NC, NGc, PETSC_NULL_INTEGER, dmF, ierr)</p>

<p>! ----- Some steps ... </p>

<p>call DMCreateColoring(dmF, IS_COLORING_GLOBAL, iscoloring, ierr)<br />call MatFDColoringCreate(Jac,iscoloring, matfdcoloring, ierr)<br />call MatFDColoringSetFunction(matfdcoloring, FormFunction, CTX, ierr)<br />call MatFDColoringSetUp(Jac ,iscoloring,matfdcoloring, ierr)<br />call SNESSetJacobian(snes, Jac, Jac, SNESComputeJacobianDefaultColor, matfdcoloring, ierr)</p>

<p>Everything is running smoothly. <br /><span style="font-size: 10pt;">Recently, we modified the boundary conditions such as to use periodic BC: </span></p>

<p>call DMDACreate1D(PETSC_COMM_WORLD, DM_BOUNDARY_PERIODIC, NY, NC, NGc, PETSC_NULL_INTEGER, dmF, ierr)</p>

<p>We then encountered frequent crashes when calling MatFDColoringSetUp, depending on the number of cells NY. After looking for an solution, I found this old thread: <a href="https://lists.mcs.anl.gov/pipermail/petsc-users/2013-May/017449.html">https://lists.mcs.anl.gov/pipermail/petsc-users/2013-May/017449.html</a> <br />It appears that when using periodic BC, FDColoring can only be used if the number of cells is divisible by 2*NGc+1. Even though this is only a slight annoyance, I was wondering if you were working on this matter / if you had a quick fix at hand? At any rate, I think it would be nice if a warning was displayed in the FDColoring documentation?  <br /><br />Respectfully, <br />Pierre Bernigaud </p>


</body></html>