<div dir="ltr">Thanks both of you.<div>Great suggestions.<br><div><br></div><div>Ling</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 3, 2017 at 12:55 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  The TS/SNES/KSP/PC allocate various work vectors and matrices which they keep around for additional solvers. Thus if you want to change the size of your problem and/or number of degrees of freedom per process you need to do one of 2 things<br>
<br>
   1) destroy the old TS/SNES/ or KSP and create a new one each time you change your vectors to a different size or<br>
<br>
   2) call TSReset/SNESReset/or KSPReset each time you change your vectors, this causes all the internal work vectors and matrices to be freed so they can get replaced but keeps all the solver options etc.<br>
<br>
   2) is preferred to 1.<br>
<br>
   Barry<br>
<div><div class="h5"><br>
> On Sep 3, 2017, at 1:44 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
><br>
> On Sun, Sep 3, 2017 at 1:02 PM, Zou, Ling <<a href="mailto:ling.zou@inl.gov">ling.zou@inl.gov</a>> wrote:<br>
> Hi All,<br>
><br>
> I wonder if there is a good example on how PETSc should work with mesh refinement. I have done some online search/PETSc-user email search, but did not find a good one.<br>
><br>
> Here is my case. Currently, assuming NO mesh refinement (so I know exactly the problem size, mesh connectivity, non-zero pattern of matrices, etc.)<br>
> In a transient simulation, my PETSc code looks like:<br>
><br>
> ==============================<wbr>==============================<wbr>============<br>
> SetupVectors(); // unknown vec, old unknown, etc.<br>
> SetupResFunc(); // residual function, basically pointing to the discretized PDE<br>
> SetupMatrices(); // allocate memory for preconditioning matrix based on non-zero pattern (so it is mesh dependent)<br>
>                            // I also use finite differencing to get preconditioning matrix, so there is another MatFDColoring object associated with the matrix<br>
> for (int step = 0; step < max_step; step++)<br>
> {<br>
>   SNESSolve();<br>
> }<br>
> ==============================<wbr>==============================<wbr>============<br>
><br>
> This setup works perfectly for my problems.<br>
> (I understand it is low efficient to use FD to get the preconditioning matrix, but it is not a priority to me at this moment).<br>
><br>
> However, now I would like to do some mesh refinement, so the previous setup would not work, because the problem size, non-zero pattern of the matrix and MatFDColoring change with mesh size (and connectivity).<br>
> I suppose my code should be changed like this:<br>
><br>
> ==============================<wbr>==============================<wbr>============<br>
> SetupVectors();<br>
> SetupResFunc();<br>
> SetupMatrices();<br>
><br>
> for (int step = 0; step < max_step; step++)<br>
> {<br>
>   MeshRefinement(); // I will manage my mesh<br>
><br>
>   AdjustVectors();   // resize vector<br>
>   AdjustMatrices(); // resize matrix, re-allocate memory, adjust non-zero pattern, etc.<br>
><br>
>   SNESSolve();<br>
> }<br>
> ==============================<wbr>==============================<wbr>============<br>
><br>
> The issue is that I understand one should not re-allocate memory to both PETSc vector and matrix dynamically during simulation, which has been discussed here:<br>
</div></div>> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__scicomp.stackexchange.com_questions_373_is-2Dit-2Dpossible-2Dto-2Ddynamically-2Dresize-2Da-2Dsparse-2Dmatrix-2Din-2Dthe-2Dpetsc-2Dlibrary&d=DwIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=kuHHom1yjd94zUrBWecnYg&m=lEB-2RrZ7_xUBhwaeQnezW0fuTVkBbMZVcRtx8Fntz0&s=6JH1-QnmGFstCWaXbdRsDRn-lSHYB-gqGOlsCFyhDVQ&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=https-3A__<wbr>scicomp.stackexchange.com_<wbr>questions_373_is-2Dit-<wbr>2Dpossible-2Dto-2Ddynamically-<wbr>2Dresize-2Da-2Dsparse-<wbr>2Dmatrix-2Din-2Dthe-2Dpetsc-<wbr>2Dlibrary&d=DwIFAg&c=<wbr>54IZrppPQZKX9mLzcGdPfFD1hxrcB_<wbr>_aEkJFOKJFd00&r=<wbr>kuHHom1yjd94zUrBWecnYg&m=lEB-<wbr>2RrZ7_<wbr>xUBhwaeQnezW0fuTVkBbMZVcRtx8Fn<wbr>tz0&s=6JH1-QnmGFstCWaXbdRsDRn-<wbr>lSHYB-gqGOlsCFyhDVQ&e=</a><br>
<span class="">><br>
> So, the final questions are:<br>
> 1) Any good example on PETSc working with mesh refinement?<br>
> 2) Any suggestion on my proposed code structure?<br>
><br>
> 1) No<br>
><br>
> 2) Yes. We advocate using the DM to describe your data layout. This allows all Vec and Mats to be created automatically, and<br>
> also tells the solver (SNESSetDM, etc.) the sizes and potential decomposition.<br>
><br>
> We are working on a mesh adaptivity interface which will allow adaptivity to be used with all DM problems without any extra<br>
> code. We should have some examples with this soon.<br>
><br>
>   Thanks,<br>
><br>
>      Matt<br>
><br>
> Best,<br>
><br>
> Ling<br>
><br>
><br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
> -- Norbert Wiener<br>
><br>
</span>> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.caam.rice.edu_-7Emk51_&d=DwIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=kuHHom1yjd94zUrBWecnYg&m=lEB-2RrZ7_xUBhwaeQnezW0fuTVkBbMZVcRtx8Fntz0&s=URW-NglPj9JSNGOxqvL8qv9WilvxNzmYIw8A6IEn5tg&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.<wbr>caam.rice.edu_-7Emk51_&d=<wbr>DwIFAg&c=<wbr>54IZrppPQZKX9mLzcGdPfFD1hxrcB_<wbr>_aEkJFOKJFd00&r=<wbr>kuHHom1yjd94zUrBWecnYg&m=lEB-<wbr>2RrZ7_<wbr>xUBhwaeQnezW0fuTVkBbMZVcRtx8Fn<wbr>tz0&s=URW-<wbr>NglPj9JSNGOxqvL8qv9WilvxNzmYIw<wbr>8A6IEn5tg&e=</a><br>
<br>
</blockquote></div><br></div>