[petsc-users] Question about how to solve the DAE step by step using PETSC4PY.

Good afternoon, Thanks for all your help. I got a question about how to solve the DAE step by step using PETSC4PY. I got two DAE systems, let's say f1(x1, u1) and f2(x2, u2). The simplified algorithm I need to implement is as follows: Step 1: solve f1 for 1 step. Step 2: use part of x1 as input

Re: [petsc-users] Configuring with CMake

I didn't get around to writing a reproducer, but I did send an extensive complaint to the CMake user group. I haven't heard anything back yet. I finally got this to work by adding this to the CMakeLists.txt file if (NOT BUILD_SHARED_LIBS) target_link_libraries(gridpack_math PUBLIC

Re: [petsc-users] Question on MATMFFD_WP

This uses a PETSc "trick". When a norm is computed on a vector it is "stashed" in the object and retrieved quickly if requested again. https://petsc.org/main/src/vec/vec/interface/rvector.c.html#VecNorm Because PETSc

Re: [petsc-users] DMSwarm

Hi, On Wed 23. Mar 2022 at 18:52, Matthew Knepley wrote: > On Wed, Mar 23, 2022 at 11:09 AM Joauma Marichal < > joauma.maric...@uclouvain.be> wrote: > >> Hello, >> >> I sent an email last week about an issue I had with DMSwarm but did not >> get an answer yet. If there is any other information

[petsc-users] Question on MATMFFD_WP

Hello PETSc team and users, I have a question regarding MATMFFD_WP : the documentation states that 1) || U || does not change between linear iterations and 2) in GMRES || a || = 1 except at restart. In src/mat/impls/mffd/wp.c, in MatMFFDCompute_WP, I see that the computation of || U || is