Re: [petsc-users] PETSc User Meeting 2017, June 14-16 in Boulder, Colorado

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 11:26 PM, Jed Brown wrote: > The program is up on the website: > > https://www.mcs.anl.gov/petsc/meetings/2017/ Put Toby on the oanel. Matt > > If you haven't registered yet, we can still accommodate you, but please > register soon. If you

Re: [petsc-users] PETSc User Meeting 2017, June 14-16 in Boulder, Colorado

2017-05-25 Thread Jed Brown
The program is up on the website: https://www.mcs.anl.gov/petsc/meetings/2017/ If you haven't registered yet, we can still accommodate you, but please register soon. If you haven't booked lodging, please do that soon -- the on-campus lodging option will close on *Tuesday, May 30*.

Re: [petsc-users] PCFactorSetShiftType does not work in code but -pc_factor_set_shift_type works

2017-05-25 Thread Danyang Su
Hi Hong, It works like a charm. I really appreciate your help. Regards, Danyang On 17-05-25 07:49 AM, Hong wrote: Danyang: You must access inner pc, then set shift. See petsc/src/ksp/ksp/examples/tutorials/ex7.c For example, I add following to petsc/src/ksp/ksp/examples/tutorials/ex2.c,

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 2:22 PM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 25 May 2017, at 20:03, Matthew Knepley wrote: > > > > > > Hmm, I thought I made adjacency per field. I have to look. That way, no > problem with the Stokes example. DG is

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Lawrence Mitchell
> On 25 May 2017, at 20:03, Matthew Knepley wrote: > > > Hmm, I thought I made adjacency per field. I have to look. That way, no > problem with the Stokes example. DG is still weird. You might, we don't right now. We just make the topological adjacency that is "large

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 1:58 PM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 25 May 2017, at 19:46, Matthew Knepley wrote: > > > > Sounds like DG. I will get out my dead chicken for the incantation > > Actually no! Mixed H(div)-L2 for Stokes. Which

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Lawrence Mitchell
> On 25 May 2017, at 19:46, Matthew Knepley wrote: > > Sounds like DG. I will get out my dead chicken for the incantation Actually no! Mixed H(div)-L2 for Stokes. Which has facet integrals for partially discontinuous fields. If you do redundant compute for such terms,

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 1:38 PM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 25 May 2017, at 19:23, Matthew Knepley wrote: > > > > Ok, let me clarify. > > > > Given shared facets, I'd like closure(support(facet)) this is a subset > of the fem

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Lawrence Mitchell
> On 25 May 2017, at 19:23, Matthew Knepley wrote: > > Ok, let me clarify. > > Given shared facets, I'd like closure(support(facet)) this is a subset of the > fem adjacency. "Add in the cell and its closure from the remote rank". This > doesn't include remote cells I can

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 1:10 PM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 25 May 2017, at 18:05, Matthew Knepley wrote: > > If you want that, is there a reason you cannot use the FEM style > FALSE+TRUE? > If you already want the closure, usually

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Lawrence Mitchell
> On 25 May 2017, at 18:05, Matthew Knepley wrote: > > If you want that, is there a reason you cannot use the FEM style FALSE+TRUE? > If you already want the closure, usually the star is not really adding > anything new. Ok, let me clarify. Given shared facets, I'd like

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 11:27 AM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > On 25/05/17 16:25, Matthew Knepley wrote: > > On Thu, May 25, 2017 at 9:27 AM, Lawrence Mitchell > > > > wrote: > > > >

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Lawrence Mitchell
On 25/05/17 16:25, Matthew Knepley wrote: > On Thu, May 25, 2017 at 9:27 AM, Lawrence Mitchell > > wrote: > > Dear petsc-users, > > I am trying to distribute a triangle mesh with a cell halo defined by >

Re: [petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Matthew Knepley
On Thu, May 25, 2017 at 9:27 AM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > Dear petsc-users, > > I am trying to distribute a triangle mesh with a cell halo defined by > FVM adjacency (i.e. if I have a facet in the initial (0-overlap) > distribution, I want the cell on the

Re: [petsc-users] PCFactorSetShiftType does not work in code but -pc_factor_set_shift_type works

2017-05-25 Thread Hong
Danyang: You must access inner pc, then set shift. See petsc/src/ksp/ksp/examples/tutorials/ex7.c For example, I add following to petsc/src/ksp/ksp/examples/tutorials/ex2.c, line 191: PetscBool isbjacobi; PCpc; ierr = KSPGetPC(ksp,);CHKERRQ(ierr); ierr =

[petsc-users] DMPlex distribution with FVM adjacency

2017-05-25 Thread Lawrence Mitchell
Dear petsc-users, I am trying to distribute a triangle mesh with a cell halo defined by FVM adjacency (i.e. if I have a facet in the initial (0-overlap) distribution, I want the cell on the other side of it). Reading the documentation, I think I do: DMPlexSetAdjacencyUseCone(PETSC_TRUE)

Re: [petsc-users] PETSC OO C guide/standard?

2017-05-25 Thread John Chludzinski
Thanks. C++ has now become the apotheosis of "no value-added complexity". Even Bjarne Stroustrup admits to understanding only a small fraction of the whole. On Wed, May 24, 2017 at 9:53 AM, Matthew Knepley wrote: > On Wed, May 24, 2017 at 8:50 AM, John Chludzinski

[petsc-users] Postdoc Position in the Computational Seismology Group at ETH Zurich.

2017-05-25 Thread Michael Afanasiev
Hi everyone, The computational seismology group at ETH Zurich is looking for a postdoc to work with us on Salvus (www.salvus.io ) - a spectral-element software package for full-waveform modelling and inversion. The exact focus of the job is tied to the applicant's

[petsc-users] PCFactorSetShiftType does not work in code but -pc_factor_set_shift_type works

2017-05-25 Thread Danyang Su
Dear Hong and Barry, I have implemented this option in the code, as we also need to use configuration from file for convenience. When I run the code using options, it works fine, however, when I run the code using configuration file, it does not work. The code has two set of equations, flow