Re: [petsc-users] Advice on small block matrix vector multiplication

2023-06-18 Thread Matthew Knepley
On Sun, Jun 18, 2023 at 2:49 PM Duan Junming wrote: > *From:* knep...@gmail.com > *Sent:* Sunday, June 18, 2023 20:35 > *To:* Duan Junming > *Cc:* petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] Advice on small block matrix vector > multiplication > > > On S

Re: [petsc-users] Advice on small block matrix vector multiplication

2023-06-18 Thread Duan Junming via petsc-users
From: knep...@gmail.com Sent: Sunday, June 18, 2023 20:35 To: Duan Junming Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Advice on small block matrix vector multiplication On Sun, Jun 18, 2023 at 2:13 PM Duan Junming via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Dear a

Re: [petsc-users] Advice on small block matrix vector multiplication

2023-06-18 Thread Matthew Knepley
On Sun, Jun 18, 2023 at 2:13 PM Duan Junming via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear all, > > I am using DMPlex to manage the unknowns, two fields, one for pressure, > and one for velocities with two/three components, defined in each cell. > They're represented by polynomials, wit

[petsc-users] Advice on small block matrix vector multiplication

2023-06-18 Thread Duan Junming via petsc-users
Dear all, I am using DMPlex to manage the unknowns, two fields, one for pressure, and one for velocities with two/three components, defined in each cell. They're represented by polynomials, with N (10~50) dofs for each component. I have an operator which can be written in a matrix form (N-by-N,