Re: [petsc-users] Scalable Solver for Incompressible Flow

Alexander Lindsay writes: > This has been a great discussion to follow. Regarding > >> when time stepping, you have enough mass matrix that cheaper preconditioners >> are good enough > > I'm curious what some algebraic recommendations might be for high Re in > transients. What mesh aspect

Re: [petsc-users] Initializing kokkos before petsc causes a problem

Hi, Philip, Thanks for reporting. I will have a look at the issue. --Junchao Zhang On Wed, Jun 7, 2023 at 9:30 AM Fackler, Philip via petsc-users < petsc-users@mcs.anl.gov> wrote: > I'm encountering a problem in xolotl. We initialize kokkos before > initializing petsc. Therefore... > > The

Re: [petsc-users] Scalable Solver for Incompressible Flow

This has been a great discussion to follow. Regarding > when time stepping, you have enough mass matrix that cheaper preconditioners are good enough I'm curious what some algebraic recommendations might be for high Re in transients. I've found one-level DD to be ineffective when applied

[petsc-users] Interpolation Between DMSTAG Objects

Hello, I am new to PETSc so apologies in advance if there is an easy answer to this question I've overlooked. I have a problem in which the computational domain is divided into two overlapping regions (overset grids). I would like to discretize each region as a separate DMSTAG object. What I do

[petsc-users] PETSc :: FEM Help

Good Morning, I'm trying to verify that the CAD -> PETSc/DMPlex methods I've developed can be used for FEM analyses using PETSc. Attached is my current attempt where I import a CAD STEP file to create a volumetric tetrahedral discretization (DMPlex), designate boundary condition points using

[petsc-users] Initializing kokkos before petsc causes a problem

I'm encountering a problem in xolotl. We initialize kokkos before initializing petsc. Therefore... The pointer referenced here: https://gitlab.com/petsc/petsc/-/blob/main/src/vec/is/sf/impls/basic/kokkos/sfkok.kokkos.cxx#L363

Re: [petsc-users] Petsc using VecDuplicate in solution process

Can you please present the all output that callgrind is outputing to you that provides this indication. > On Jun 7, 2023, at 4:11 AM, Pichler, Franz wrote: > > Hello thanks for the reply, > > I created a working minimal example (as minimal as I can think of?) that I > include here, even

Re: [petsc-users] Petsc using VecDuplicate in solution process

Hello thanks for the reply, I created a working minimal example (as minimal as I can think of?) that I include here, even though I am not sure which is the best format to do this, I just add some plain text: