On Tue, Mar 12, 2024 at 10:54 PM Adrian Croucher
wrote:
> hi Matt,
> On 13/03/24 12:46 pm, Matthew Knepley wrote:
>
>
> Thanks for finding this. This was indeed rewritten because the old
> algorithm was too expensive at large sizes/parallelism. The function where
> this breaks in DMPlexStratifyMi
hi Matt,
On 13/03/24 12:46 pm, Matthew Knepley wrote:
Thanks for finding this. This was indeed rewritten because the old
algorithm was too expensive at large sizes/parallelism. The function
where this breaks in DMPlexStratifyMigrationSF() in plexdistribute.c.
The problem is here:
https://u
On Tue, Mar 12, 2024 at 6:49 PM Adrian Croucher
wrote:
> hi Matt,
>
> Thanks for that. I think that is a useful change, even though it doesn't
> actually fix my problem (i.e. my code still gets an error). I think now
> that although the error message I'm getting is similar to the one that used
>
hi Matt,
Thanks for that. I think that is a useful change, even though it doesn't
actually fix my problem (i.e. my code still gets an error). I think now
that although the error message I'm getting is similar to the one that
used to result from just a simple repeated stratification, the
under
On Wed, Mar 6, 2024 at 5:44 PM Adrian Croucher
wrote:
> hi Matt,
> On 7/03/24 2:55 am, Matthew Knepley wrote:
>
> Let's do this. Can you send me a dual porosity mesh, and I will figure out
> why the stratification
> algorithm does not work? I think you can just output the Plex as binary or
> HDF5
On Thu, Mar 7, 2024 at 12:08 AM Adrian Croucher
wrote:
> hi Matt,
>
> With a bit more experimenting, I find I get the same sort of error in
> DMPlexStratify() if I run a very simple program that just reads in a
> DMPlex from file and then stratifies it (no dual porosity, no
> distribution etc.),
hi Matt, With a bit more experimenting, I find I get the same sort of error in DMPlexStratify() if I run a very simple program that just reads in a DMPlex from file and then stratifies it (no dual porosity, no distribution etc. ), i. e. : call
ZjQcmQRYFpfptBannerStart
On Tue, Mar 5, 2024 at 10:43 PM Adrian Croucher
wrote:
> hi Matt,
>
> I've made a bit of progress with this, using a debugger.
>
> I wasn't setting the overlap for my dual-porosity DM, so it was
> defaulting to zero, and when it got to the DMPlexCheckPointSF() call it
> was doing the check for ce
hi Matt, I've made a bit of progress with this, using a debugger. I wasn't setting the overlap for my dual-porosity DM, so it was defaulting to zero, and when it got to the DMPlexCheckPointSF() call it was doing the check for cells in the overlap
ZjQcmQRYFpfptBannerStart
hi Matt,
On 5/03/24 1:40 am, Matthew Knepley wrote:
Let's start with the distribution. I think the problem may be here.
Even with dual porosity cells, I don't think we should have cells in
the pointSF, just the internal face between cells. Do you agree?
You mean the re-distribution? If you ar
On Sun, Mar 3, 2024 at 9:48 PM Adrian Croucher
wrote:
> hi Matt,
> On 4/03/24 3:25 pm, Matthew Knepley wrote:
>
>
> So here I am guessing that the interior cells (rock matrix) mess with the
> stratum consistency
> (ordering of dimension by stratum in the DAG). So some questions:
>
> 1) Are rock m
hi Matt,
On 4/03/24 3:25 pm, Matthew Knepley wrote:
So here I am guessing that the interior cells (rock matrix) mess with
the stratum consistency
(ordering of dimension by stratum in the DAG). So some questions:
1) Are rock matrix cells always associated with only one regular cell?
Yes. In t
On Sun, Mar 3, 2024 at 9:18 PM Adrian Croucher
wrote:
> On 4/03/24 2:36 pm, Adrian Croucher wrote:
> > I just re-tested my code with this "no cells in the interface" check
> > commented out of the PETSc code, and it appears to run fine.
>
I put in that check to look for improper handling of over
On 4/03/24 2: 36 pm, Adrian Croucher wrote: > I just re-tested my code with this "no cells in the interface" check > commented out of the PETSc code, and it appears to run fine. Whoops, no it doesn't. . . I checked more closely and there
ZjQcmQRYFpfptBannerStart
This Me
hi again, On 1/03/24 3: 41 pm, Adrian Croucher wrote: > > I am still getting problems with redistribution of MINC (dual > porosity) meshes. The function DMPlexCheckPointSF() is now raising > errors like this: > > [0]PETSC ERROR:
ZjQcmQRYFpfptBannerStart
This Message Is
hi Matt,
On 20/01/24 6:03 am, Matthew Knepley wrote:
I put in a test for redistribution with overlap, Plex test ex12, test
7. You can run it using
make -f ./gmakefile test search="dm_impls_plex_tests-ex12_7"
and see the output in src/dm/impls/plex/tests/output/ex12_7.out. The
second DMPl
16 matches
Mail list logo
