<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" id="owaParaStyle">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
P {margin-top:0;margin-bottom:0;}</style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi Matthew<br>
<br>
It seems like the problem is not fully fixed. I have changed the code to now run on with both 2,3 and 4 cells. When I run the code using NP = 1..3 I get different result both for NP=1 to NP=2/3 when cell count is larger than 2.<br>
<br>
Kind regards,<br>
Morten<br>
____<br>
mpiexec -np 1 ./ex18k<br>
cells 2<br>
Loc size: 36<br>
Trace of matrix: 132.000000<br>
cells 3<br>
Loc size: 48<br>
Trace of matrix: 192.000000<br>
cells 4<br>
Loc size: 60<br>
Trace of matrix: 258.000000<br>
mpiexec -np 2 ./ex18k<br>
cells 2<br>
Loc size: 24<br>
Loc size: 24<br>
Trace of matrix: 132.000000<br>
cells 3<br>
Loc size: 36<br>
Loc size: 24<br>
Trace of matrix: 198.000000<br>
cells 4<br>
Loc size: 36<br>
Loc size: 36<br>
Trace of matrix: 264.000000<br>
mpiexec -np 3 ./ex18k<br>
cells 2<br>
Loc size: 24<br>
Loc size: 24<br>
Loc size: 0<br>
Trace of matrix: 132.000000<br>
cells 3<br>
Loc size: 24<br>
Loc size: 24<br>
Loc size: 24<br>
Trace of matrix: 198.000000<br>
cells 4<br>
Loc size: 36<br>
Loc size: 24<br>
Loc size: 24<br>
Trace of matrix: 264.000000<br>
<br>
<br>
<br>
<br>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF16481" style="direction:ltr"><font size="2" face="Tahoma" color="#000000"><b>From:</b> petsc-users-bounces@mcs.anl.gov [petsc-users-bounces@mcs.anl.gov] on behalf of Morten Nobel-Jørgensen [mono@dtu.dk]<br>
<b>Sent:</b> Sunday, September 25, 2016 11:15 AM<br>
<b>To:</b> Matthew Knepley<br>
<b>Cc:</b> PETSc [petsc-users@mcs.anl.gov]<br>
<b>Subject:</b> Re: [petsc-users] DMPlex problem<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi Matthew
<div><br>
</div>
<div>Thank you for the bug-fix :) I can confirm that it works :)</div>
<div><br>
</div>
<div>And thanks for your hard work on PETSc - your work is very much appreciated!</div>
<div><br>
</div>
<div>Kind regards,</div>
<div>Morten<br>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF211326" style="direction:ltr"><font size="2" face="Tahoma" color="#000000"><b>From:</b> Matthew Knepley [knepley@gmail.com]<br>
<b>Sent:</b> Friday, September 23, 2016 2:46 PM<br>
<b>To:</b> Morten Nobel-Jørgensen<br>
<b>Cc:</b> PETSc [petsc-users@mcs.anl.gov]<br>
<b>Subject:</b> Re: [petsc-users] DMPlex problem<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Fri, Sep 23, 2016 at 7:45 AM, Matthew Knepley <span dir="ltr">
<<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote"><span class="">On Fri, Sep 23, 2016 at 3:48 AM, Morten Nobel-Jørgensen
<span dir="ltr"><<a href="mailto:mono@dtu.dk" target="_blank">mono@dtu.dk</a>></span> wrote:<br>
</span>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:tahoma; color:rgb(0,0,0); font-size:10pt">
<div><span style="font-size:10pt">Dear PETSc developers</span></div>
<div><br>
</div>
<span class="">
<div>Any update on this issue regarding DMPlex? Or is there any obvious workaround that we are unaware of?</div>
</span></div>
</div>
</blockquote>
<div><br>
</div>
<div>I have fixed this bug. It did not come up in nightly tests because we are not using MatSetValuesLocal(). Instead we</div>
<div>use MatSetValuesClosure() which translates differently.</div>
<div><br>
</div>
<div>Here is the branch</div>
<div><br>
</div>
<div> <a href="https://bitbucket.org/petsc/petsc/branch/knepley/fix-dm-ltog-bs" target="_blank">https://bitbucket.org/petsc/<wbr>petsc/branch/knepley/fix-dm-<wbr>ltog-bs</a></div>
<div><br>
</div>
<div>and I have merged it to next. It will go to master in a day or two.</div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Also, here is the cleaned up source with no memory leaks.</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote"><span class="">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:tahoma; color:rgb(0,0,0); font-size:10pt">
<div>Also should we additionally register the issue on Bitbucket or is reporting the issue on the mailing list enough?</div>
</div>
</div>
</blockquote>
<div><br>
</div>
</span>
<div>Normally we are faster, but the start of the semester was hard this year.</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div>
<div class="h5">
<div> </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:tahoma; color:rgb(0,0,0); font-size:10pt">
<div>Kind regards,</div>
<div>Morten</div>
<br>
<div style="font-family:"times new roman"; color:rgb(0,0,0); font-size:16px"><span>
<hr>
<div style="direction:ltr"><font size="2" face="Tahoma" color="#000000"><b>From:</b> Matthew Knepley [<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>]<br>
<b>Sent:</b> Friday, September 09, 2016 12:21 PM<br>
<b>To:</b> Morten Nobel-Jørgensen<br>
<b>Cc:</b> PETSc [<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>]<br>
<b>Subject:</b> Re: [petsc-users] DMPlex problem<br>
</font><br>
</div>
<div></div>
</span>
<div>
<div>
<div>
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Fri, Sep 9, 2016 at 4:04 AM, Morten Nobel-Jørgensen <span dir="ltr">
<<a href="mailto:mono@dtu.dk" target="_blank">mono@dtu.dk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:tahoma; color:rgb(0,0,0); font-size:10pt">
Dear PETSc developers and users,<br>
<br>
Last week we posted a question regarding an error with DMPlex and multiple dofs and have not gotten any feedback yet. This is uncharted waters for us, since we have gotten used to an extremely fast feedback from the PETSc crew. So - with the chance of sounding
impatient and ungrateful - we would like to hear if anybody has any ideas that could point us in the right direction?<br>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>This is my fault. You have not gotten a response because everyone else was waiting for me, and I have been</div>
<div>slow because I just moved houses at the same time as term started here. Sorry about that.</div>
<div><br>
</div>
<div>The example ran for me and I saw your problem. The local-tp-global map is missing for some reason.</div>
<div>I am tracking it down now. It should be made by DMCreateMatrix(), so this is mysterious. I hope to have</div>
<div>this fixed by early next week.</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:tahoma; color:rgb(0,0,0); font-size:10pt">
We have created a small example problem that demonstrates the error in the matrix assembly.<br>
<div><font size="1"><span style="font-size:13px"><font size="2"><br>
Thanks,<br>
Morten</font></span></font></div>
<div><font size="1"><span style="font-size:13px"><font size="2"><br>
</font></span></font></div>
<div><font size="1"><span style="font-size:13px"><font size="2"><br>
</font></span></font></div>
</div>
</div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div>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</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
<div>
<div class="h5"><br>
<br clear="all">
<div><br>
</div>
-- <br>
<div>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</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="gmail_signature">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</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>