Re: [petsc-users] Question about SuperLU

2022-06-09 Thread Xiaoye S. Li
Are you using serial SuperLU, or distributed-memory SuperLU_DIST? What are the algorithm options are you using? Sherry Li On Thu, Jun 9, 2022 at 2:20 PM Jorti, Zakariae via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hi, > > I am solving non-linear problem that has 5 unknowns {ni, T, E, B,

Re: [petsc-users] VecConcatenate in petsc4py

2022-06-09 Thread Jed Brown
You don't want to create a new vector here, but read from (and write to) multiple parts of the same vector. You can use PETSc interfaces for subvectors, or do it with NumPy slices (perhaps more natural and ergonomic, depending on how your code is written). Armand Touminet via petsc-users

[petsc-users] Question about SuperLU

2022-06-09 Thread Jorti, Zakariae via petsc-users
Hi, I am solving non-linear problem that has 5 unknowns {ni, T, E, B, V}, and for the preconditioning part, I am using a FieldSplit preconditioner. At the last fieldsplit/level, we are left with a {B,V} block that tried to precondition in 2 different ways: a) SuperLU:

[petsc-users] VecConcatenate in petsc4py

2022-06-09 Thread Armand Touminet via petsc-users
Dear Petsc team, I'm trying to implement PDE constrained optimization using TAO from the petsc4py interface. Since my problem has multiple parameter fields to optimize, I need to combine them into a single Vec object to supply to TAO. I've found the VecConcatenate in the C documentation, which