[petsc-users] Refine DMPlex with a Refinement Function

[Matthew Knepley]

On Fri, Aug 2, 2019 at 9:10 AM Daniel Mckinnell via petsc-users <
petsc-users at mcs.anl.gov> wrote:

> Thank you a lot, that works well for Simplex cells when running on one
> processor. Do you know anything that works for Tensor cells and/or when
> running in parallel?
>

For simplex in parallel, we use Pragmatic, and for tensor cells in parallel
we use p4est. Both of them operate through the
very new DMAdapt() interface. We have some initial experience, but that
part of PETSc is still unstable.

Note that you can regularly refine simplices or hexes in parallel using
just DMRefine().

  Thanks,

    Matt


> Thanks again for all of your help.
> Daniel
>
>
> -----Original Message-----
> From: Lawrence Mitchell <wence at gmx.li>
> To: Daniel Mckinnell <d_mckinnell at aol.co.uk>
> CC: petsc-users <petsc-users at mcs.anl.gov>
> Sent: Thu, 1 Aug 2019 17:36
> Subject: Re: [petsc-users] Refine DMPlex with a Refinement Function
>
>
>
> On Thu, 1 Aug 2019 at 16:59, Daniel Mckinnell via petsc-users <
> petsc-users at mcs.anl.gov> wrote:
>
> Hi,
>
> I have been having some trouble trying to refine a DMPlex object using a
> Refinement Function. I have been working with reference to the code
> discussed here:
> https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2019-April/038341.html
> including the code from the GitHub directory mentioned. The error I have
> been getting is as follows:
>
> [0]PETSC ERROR: Argument out of range
> [0]PETSC ERROR: No grid refiner of dimension 2 registered
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.11.2, unknown
> [0]PETSC ERROR: ./fpf on a arch-linux2-c-debug named cromars by daniel Thu
> Aug  1 15:40:18 2019
> [0]PETSC ERROR: Configure options −−download−f2cblaslapack=yes
> −−with−debugging=1 −−download−metis=yes −−download−parmetis=yes
> −−with−fortran−bindings=0 −−with−python=0 -download-openmpi=1
> −−with−c−support −−with−clanguage=cxx
> [0]PETSC ERROR: #1 DMPlexRefine_Internal() line 215 in
> /home/daniel/petsc/src/dm/impls/plex/plexadapt.c
> [0]PETSC ERROR: #2 DMRefine_Plex() line 10381 in
> /home/daniel/petsc/src/dm/impls/plex/plexrefine.c
> [0]PETSC ERROR: #3 DMRefine() line 1881 in
> /home/daniel/petsc/src/dm/interface/dm
>
>
> I suspect you need to configure with --download-triangle? When you mark
> cells for refinement like this, petsc uses and external package to do the
> adaptation.
>
> I'm not sure what options are available as packages, but triangle I think
> is one.
>
> Cheers,
>
> Lawrence
>
>
>

-- 
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

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190802/9dd20e6f/attachment.html>