On Fri, Sep 23, 2016 at 7:45 AM, Matthew Knepley <[email protected]> wrote:
> On Fri, Sep 23, 2016 at 3:48 AM, Morten Nobel-Jørgensen <[email protected]> > 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 [[email protected]] >> *Sent:* Friday, September 09, 2016 12:21 PM >> *To:* Morten Nobel-Jørgensen >> *Cc:* PETSc [[email protected]] >> *Subject:* Re: [petsc-users] DMPlex problem >> >> On Fri, Sep 9, 2016 at 4:04 AM, Morten Nobel-Jørgensen <[email protected]> >> 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
ex18.c
Description: Binary data
