Re: [petsc-users] Tough to reproduce petsctablefind error

2020-09-25 Thread Barry Smith
I May 2019 Lisandro changed the versions of Metis and ParMetis that PETSc uses to use a portable machine independent random number generator so if you are having PETSc install Metis then its random number generator should generate the exact same random numbers on repeated identical runs on

Re: [petsc-users] Tough to reproduce petsctablefind error

2020-09-25 Thread Matthew Knepley
On Fri, Sep 25, 2020 at 1:29 PM Mark McClure wrote: > Hello, > > I work with Chris, and have a few comments, hopefully helpful. Thank you > all, for your help. > > Our program is unfortunately behaving a little bit nondeterministically. I > am not sure why because for the OpenMP parts, I test it

Re: [petsc-users] Tough to reproduce petsctablefind error

2020-09-25 Thread Mark McClure
Hello, I work with Chris, and have a few comments, hopefully helpful. Thank you all, for your help. Our program is unfortunately behaving a little bit nondeterministically. I am not sure why because for the OpenMP parts, I test it for race conditions using Intel Inspector and see none. We are

Re: [petsc-users] Finding which cell an arbitrary point belongs to in DMPlex

2020-09-25 Thread Matthew Knepley
On Fri, Sep 25, 2020 at 7:44 AM Jeremy Theler wrote: > On Thu, 2020-09-24 at 18:30 -0400, Matthew Knepley wrote: > > > There is also DMPlexFindVertices() which finds the nearest vertex > > > to > > > > the given coords in the given radius. > > > > > > At first I had understood that this function

Re: [petsc-users] Finding which cell an arbitrary point belongs to in DMPlex

2020-09-25 Thread Jeremy Theler
On Thu, 2020-09-24 at 18:30 -0400, Matthew Knepley wrote: > > There is also DMPlexFindVertices() which finds the nearest vertex > > to > > > the given coords in the given radius. > > > > At first I had understood that this function performed a nearest- > > neighbor search but after a closer look

Re: [petsc-users] Compute C*Ct using MatCreateTranspose for Ct

2020-09-25 Thread Olivier Jamond
@Mark: Sorry, I totally misunderstood your answer, so my answer to your answer about petsc versions is non-sense... @Hong: MatTransposeMatMult gives (Ct*C), no? And I feel that MatCreateTranspose actually works for mpiaij. On 24/09/2020 16:16, Zhang, Hong via petsc-users wrote: Olivier

Re: [petsc-users] something wrong with digest?

2020-09-25 Thread Klaij, Christiaan
That could be the reason, there were some lengthy emails indeed. But the attachments are removed in Digest. Chris dr. ir. Christiaan Klaij | Senior Researcher | Research & Development MARIN | T +31 317 49 33 44 | mailto:c.kl...@marin.nl | http://www.marin.nl

Re: [petsc-users] something wrong with digest?

2020-09-25 Thread Satish Balay via petsc-users
I see one of the mailman settings is: How big in Kb should a digest be before it gets sent out? 0 implies no maximum size. (Edit digest_size_threshhold) 40 I don't think we ever changed this value [so likely a default] And I see a bunch of e-mails on the list exceeding this number. Perhaps

[petsc-users] something wrong with digest?

2020-09-25 Thread Klaij, Christiaan
Today I got more than 20 petsc-user Digest emails in my inbox (Vol 141, Issue 78 to 101), many containing only one message and being sent within a few minutes of each other (10:39, 10:41 10:44, ...). Is this really how Digest is supposed to work? Chris dr. ir. Christiaan Klaij | Senior