Re: [petsc-users] Question on ISLocalToGlobalMappingGetIndices Fortran Interface

2023-05-08 Thread Danyang Su
Thanks, Mark. Yes, it actually works when I update to ISLocalToGlobalMappingGetIndicesF90. I made a mistake reporting this does not work. Danyang From: Mark Adams Date: Monday, May 8, 2023 at 7:22 PM To: Danyang Su Cc: petsc-users Subject: Re: [petsc-users] Question on ISLocalToGlobalM

Re: [petsc-users] Question on ISLocalToGlobalMappingGetIndices Fortran Interface

2023-05-08 Thread Mark Adams
On Mon, May 8, 2023 at 6:50 PM Danyang Su wrote: > Dear PETSc-Users, > > > > Is there any changes in ISLocalToGlobalMappingGetIndices function after > PETSc 3.17? > > > > In the previous PETSc version (<= 3.17), the function > ‘ISLocalToGlobalMappingGetIndices(ltogm,ltog,idltog,ierr)’ works fine,

Re: [petsc-users] Help with KSPSetConvergenceTest

2023-05-08 Thread Barry Smith
See https://gitlab.com/petsc/petsc/-/merge_requests/6436 with commit barry/2023-05-08/add-ksp-min-its > On May 7, 2023, at 2:22 PM, Edoardo alinovi wrote: > > Hello Barry, > > Mega! Thank you Berry much for providing me with a working example! I ended > up in writing this: > >

[petsc-users] Question on ISLocalToGlobalMappingGetIndices Fortran Interface

2023-05-08 Thread Danyang Su
Dear PETSc-Users, Is there any changes in ISLocalToGlobalMappingGetIndices function after PETSc 3.17? In the previous PETSc version (<= 3.17), the function ‘ISLocalToGlobalMappingGetIndices(ltogm,ltog,idltog,ierr)’ works fine, even though the value of idltog looks out of bound (-1147265

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

2023-05-08 Thread Jed Brown
Sebastian Blauth writes: > Hello everyone, > > I wanted to briefly follow up on my question (see my last reply). > Does anyone know / have an idea why the LSC preconditioner in PETSc does > not seem to scale well with the problem size (the outer fgmres solver I > am using nearly scale nearly li

Re: [petsc-users] Petsc ObjectStateIncrease without proivate header

2023-05-08 Thread Stefano Zampini
You can achieve the same effect by calling MatAssemblyBegin/End Il giorno lun 8 mag 2023 alle ore 15:54 Pichler, Franz < franz.pich...@v2c2.at> ha scritto: > Hello, > i am using petsc in a single cpu setup where I have preassembled crs > matrices that I wrap via PetSC’s > > MatCreateSeqAIJWithArr

[petsc-users] Petsc ObjectStateIncrease without proivate header

2023-05-08 Thread Pichler, Franz
Hello, i am using petsc in a single cpu setup where I have preassembled crs matrices that I wrap via PetSC's MatCreateSeqAIJWithArrays Functionality. Now I manipulate the values of these matrices (wohtout changing the sparsity) without using petsc, When I now want to solve again I have to call

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

2023-05-08 Thread Matthew Knepley
On Mon, May 8, 2023 at 7:27 AM Mark Adams wrote: > > > On Mon, May 8, 2023 at 2:32 AM Sebastian Blauth < > sebastian.bla...@itwm.fraunhofer.de> wrote: > >> Hello everyone, >> >> I wanted to briefly follow up on my question (see my last reply). >> Does anyone know / have an idea why the LSC precon

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

2023-05-08 Thread Mark Adams
On Mon, May 8, 2023 at 2:32 AM Sebastian Blauth < sebastian.bla...@itwm.fraunhofer.de> wrote: > Hello everyone, > > I wanted to briefly follow up on my question (see my last reply). > Does anyone know / have an idea why the LSC preconditioner in PETSc does > not seem to scale well with the problem