[petsc-users] DMPlex problem

[Matthew Knepley]

On Fri, Sep 23, 2016 at 7:45 AM, Matthew Knepley <knepley at gmail.com> wrote:

> On Fri, Sep 23, 2016 at 3:48 AM, Morten Nobel-Jørgensen <mono at dtu.dk>
> wrote:
>
>> Dear PETSc developers
>>
>> Any update on this issue regarding DMPlex? Or is there any obvious
>> workaround that we are unaware of?
>>
>
> I have fixed this bug. It did not come up in nightly tests because we are
> not using MatSetValuesLocal(). Instead we
> use MatSetValuesClosure() which translates differently.
>
> Here is the branch
>
>   https://bitbucket.org/petsc/petsc/branch/knepley/fix-dm-ltog-bs
>
> and I have merged it to next. It will go to master in a day or two.
>

Also, here is the cleaned up source with no memory leaks.

  Matt


> Also should we additionally register the issue on Bitbucket or is
>> reporting the issue on the mailing list enough?
>>
>
> Normally we are faster, but the start of the semester was hard this year.
>
>   Thanks,
>
>      Matt
>
>
>> Kind regards,
>> Morten
>>
>> ------------------------------
>> *From:* Matthew Knepley [knepley at gmail.com]
>> *Sent:* Friday, September 09, 2016 12:21 PM
>> *To:* Morten Nobel-Jørgensen
>> *Cc:* PETSc ‎[petsc-users at mcs.anl.gov]‎
>> *Subject:* Re: [petsc-users] DMPlex problem
>>
>> On Fri, Sep 9, 2016 at 4:04 AM, Morten Nobel-Jørgensen <mono at dtu.dk>
>> wrote:
>>
>>> Dear PETSc developers and users,
>>>
>>> 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?
>>>
>>
>> This is my fault. You have not gotten a response because everyone else
>> was waiting for me, and I have been
>> slow because I just moved houses at the same time as term started here.
>> Sorry about that.
>>
>> The example ran for me and I saw your problem. The local-tp-global map is
>> missing for some reason.
>> I am tracking it down now. It should be made by DMCreateMatrix(), so this
>> is mysterious. I hope to have
>> this fixed by early next week.
>>
>>   Thanks,
>>
>>     Matt
>>
>>
>>> We have created a small example problem that demonstrates the error in
>>> the matrix assembly.
>>>
>>> Thanks,
>>> Morten
>>>
>>>
>>>
>>
>>
>> --
>> What most experimenters take for granted before they begin their
>> experiments is infinitely more interesting than any results to which their
>> experiments lead.
>> -- Norbert Wiener
>>
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160923/aadf6f6f/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex18.c
Type: application/octet-stream
Size: 4946 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160923/aadf6f6f/attachment-0001.obj>