Re: [petsc-users] Is it possible to keep track of original elements # after a call to DMPlexDistribute ?

Hi Matthew, I have joined the ex44.c example so you can see that we see. The problem is about the "f" field number which I can't see any clue in the documentation to what to put there... We are missing some information maybe written elsewhere? I have tried 3 different calls (lines 180-182)

Re: [petsc-users] Error running make on MUMPS

You are using petsc-3.15.1 - and likely the mumps build change between then and current 3.16. Can you use latest PETSc release? If not - Suggest removing --download-mumps=https://bitbucket.org/petsc/pkg-mumps.git --download-mumps-commit=v5.4.1-p1 options [and PETSC_ARCH] and going back to

Re: [petsc-users] Error running make on MUMPS

Hello April, You’ll need PETSc 3.16.0 to use --download-mumps=https://bitbucket.org/petsc/pkg-mumps.git --download-mumps-commit=v5.4.1-p1 The MUMPS Makefile rule changed between 3.15.X (alllib) and 3.16.Y (all). Thanks, Pierre > On 10 Nov 2021, at 8:00 PM, Novak, April via petsc-users >

Re: [petsc-users] Error running make on MUMPS

Thanks, Barry, We will try this, and report back Fande On Tue, Nov 9, 2021 at 5:41 PM Barry Smith wrote: > >This version of MUMPS has a bug in its build system; it does not have > all the dependencies on Fortran modules properly listed so Fortran files > can get compiled too early

Re: [petsc-users] How to compute the condition number of SchurComplementMat preconditioned with PCSHELL.

> P is a diffusion matrix, which itself is inverted by KSPCG) This worries me. Unless solved to full precision the action of solving with CG is not a linear operator in the input variable b, this means that the action of your Schur complement is not a linear operator and so iterative

Re: [petsc-users] How to compute the condition number of SchurComplementMat preconditioned with PCSHELL.

Great thanks, Matt! The second option is what I was looking for. Best Regards, Vlad От: Matthew Knepley Отправлено: 10 ноября 2021 г. 16:45:00 Кому: Vladislav Pimanov Копия: petsc-users@mcs.anl.gov Тема: Re: [petsc-users] How to compute the condition number

Re: [petsc-users] How to compute the condition number of SchurComplementMat preconditioned with PCSHELL.

On Wed, Nov 10, 2021 at 8:42 AM Vladislav Pimanov < vladislav.pima...@skoltech.ru> wrote: > Dear PETSc community, > > > I wonder if you could give me a hint on how to compute the condition > number of a preconditioned matrix in a proper way. > > I have a *MatSchurComplement* matrix S and a

[petsc-users] How to compute the condition number of SchurComplementMat preconditioned with PCSHELL.

Dear PETSc community, I wonder if you could give me a hint on how to compute the condition number of a preconditioned matrix in a proper way. I have a MatSchurComplement matrix S and a preconditioner P of the type PCSHELL (P is a diffusion matrix, which itself is inverted by KSPCG). I tried

Re: [petsc-users] dmplex extrude in parallel

On Wed, Nov 10, 2021 at 8:26 AM Bhargav Subramanya < bhargav.subrama...@kaust.edu.sa> wrote: > Dear Matt, > > Thanks a lot for the reply. I am now able to generate the prismatic mesh > properly. > Cool. > In the case of mpiexec -n 2 ./meshtest -dm_plex_shape sphere > -dm_refine_pre 3

Re: [petsc-users] dmplex extrude in parallel

Dear Matt, Thanks a lot for the reply. I am now able to generate the prismatic mesh properly. In the case of mpiexec -n 2 ./meshtest -dm_plex_shape sphere -dm_refine_pre 3 -dm_distribute -dm_refine 2 and DMExtrude(), where I am ordering the extruded cells on the layers first; Is the mesh

Re: [petsc-users] Strange behavior of TS after setting hand-coded Jacobian

Thanks for your valuable advice. Your patience and recognition is highly appreciated. I will do some reading and think it over. Some tests will be performed in recent days. Regards. Yu Barry Smith 于2021年11月7日周日 下午11:15写道： > > >I think you can use the following approach. Instead of making