May 2023 - petsc-users - mailman3.cels.anl.gov

Final reminder to register for the PETSc user conference June 5-7 if you plan to attend
by Barry Smith 09 May '23

09 May '23

One final reminder to please register for the PETSc user conference at https://www.eventbrite.com/e/petsc-2023-user-meeting-tickets-494165441137 now if you plan to attend. We need numbers for ordering food etc. The preliminary agenda can be found at https://petsc.org/release/community/meetings/2023/#meeting Thanks Barry

1 0

Understanding index sets for PCGASM
by LEONARDO MUTTI 09 May '23

09 May '23

Hello. I am having a hard time understanding the index sets to feed PCGASMSetSubdomains, and I am working in Fortran (as a PETSc novice). To get more intuition on how the IS objects behave I tried the following minimal (non) working example, which should tile a 16x16 matrix into 16 square, non-overlapping submatrices: #include <petsc/finclude/petscmat.h> #include <petsc/finclude/petscksp.h> #include <petsc/finclude/petscpc.h> USE petscmat USE petscksp USE petscpc Mat :: A PetscInt :: M, NSubx, dof, overlap, NSub INTEGER :: I,J PetscErrorCode :: ierr PetscScalar :: v KSP :: ksp PC :: pc IS :: subdomains_IS, inflated_IS call PetscInitialize(PETSC_NULL_CHARACTER , ierr) !-----Create a dummy matrix M = 16 call MatCreateAIJ(MPI_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, & M, M, & PETSC_DEFAULT_INTEGER, PETSC_NULL_INTEGER, & PETSC_DEFAULT_INTEGER, PETSC_NULL_INTEGER, & A, ierr) DO I=1,M DO J=1,M v = I*J CALL MatSetValue (A,I-1,J-1,v, & INSERT_VALUES , ierr) END DO END DO call MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY , ierr) call MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY , ierr) !-----Create KSP and PC call KSPCreate(PETSC_COMM_WORLD,ksp, ierr) call KSPSetOperators(ksp,A,A, ierr) call KSPSetType(ksp,"bcgs",ierr) call KSPGetPC(ksp,pc,ierr) call KSPSetUp(ksp, ierr) call PCSetType(pc,PCGASM, ierr) call PCSetUp(pc , ierr) !-----GASM setup NSubx = 4 dof = 1 overlap = 0 call PCGASMCreateSubdomains2D(pc, & M, M, & NSubx, NSubx, & dof, overlap, & NSub, subdomains_IS, inflated_IS, ierr) call ISView(subdomains_IS, PETSC_VIEWER_STDOUT_WORLD, ierr) call KSPDestroy(ksp, ierr) call PetscFinalize(ierr) Running this on one processor, I get NSub = 4. If PCASM and PCASMCreateSubdomains2D are used instead, I get NSub = 16 as expected. Moreover, I get in the end "forrtl: severe (157): Program Exception - access violation". So: 1) why do I get two different results with ASM, and GASM? 2) why do I get access violation and how can I solve this? In fact, in C, subdomains_IS, inflated_IS should pointers to IS objects. As I see on the Fortran interface, the arguments to PCGASMCreateSubdomains2D are IS objects: subroutine PCGASMCreateSubdomains2D(a,b,c,d,e,f,g,h,i,j,z) import tPC,tIS PC a ! PC PetscInt b ! PetscInt PetscInt c ! PetscInt PetscInt d ! PetscInt PetscInt e ! PetscInt PetscInt f ! PetscInt PetscInt g ! PetscInt PetscInt h ! PetscInt IS i ! IS IS j ! IS PetscErrorCode z end subroutine PCGASMCreateSubdomains2D Thus: 3) what should be inside e.g., subdomains_IS? I expect it to contain, for every created subdomain, the list of rows and columns defining the subblock in the matrix, am I right? Context: I have a block-tridiagonal system arising from space-time finite elements, and I want to solve it with GMRES+PCGASM preconditioner, where each overlapping submatrix is on the diagonal and of size 3x3 blocks (and spanning multiple processes). This is PETSc 3.17.1 on Windows. Thanks in advance, Leonardo

3 13

Step size setting in TS
by Jorti, Zakariae 09 May '23

09 May '23

Hello, I have a time-dependent model that I solve using TSSolve. And I am trying to adaptively change the step size (dt). I found that there are some TSAdapt schemes already available. I have tried TSADAPTBASIC and TSADAPTCFL. The former runs without any problems, whereas the latter yields the following error: " [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: No support for this operation for this object type [0]PETSC ERROR: Step rejection not implemented. The CFL implementation is incomplete/unusable " How/where the step rejection should be implemented? Are there any examples available? Besides, I was also attempting to change the step size through the monitor: " SNES snes; TSGetSNES(ts, & snes); SNESConvergedReason reason = SNES_CONVERGED_ITERATING; PetscCall(SNESGetConvergedReason(snes, &reason)); if(reason < 0){ TSSetTimeStep(ts,28618.7); } else{ TSSetTimeStep(ts,57237.4); } " But when I try this solution, the TSStep seems to diverge even though the SNES solver converges (see log file below). Am I doing something wrong here by changing the value of the step size inside the monitor? Thank you. Best, Zakariae ----------------------------------------------------------------- Timestep 0: step size = 57237.4, time = 0., 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 8.393431982129e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.985106864162e-01 1 KSP Residual norm 2.509248490095e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 1.985106864162e-01 0 KSP Residual norm 1.985106864162e-01 1 KSP Residual norm 2.739697200542e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.885623182909e-01 1 KSP Residual norm 8.784869181255e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.664114929279e-01 1 KSP Residual norm 2.277701872167e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.664114929279e-01 0 KSP Residual norm 1.664114929279e-01 1 KSP Residual norm 2.062369048201e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 6.155996487462e-02 1 KSP Residual norm 8.830688649637e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 3 SNES Function norm 6.155996487462e-02 0 KSP Residual norm 6.155996487462e-02 1 KSP Residual norm 6.231091742747e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.828806540636e-02 1 KSP Residual norm 3.623270770725e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 4 SNES Function norm 1.828806540636e-02 0 KSP Residual norm 1.828806540636e-02 1 KSP Residual norm 7.390325185746e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.659571588896e-02 1 KSP Residual norm 5.280284947190e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 5 SNES Function norm 1.659571588896e-02 0 KSP Residual norm 1.659571588896e-02 1 KSP Residual norm 8.311611006756e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.475636107890e-02 1 KSP Residual norm 1.890387034462e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.635129206894e-02 1 KSP Residual norm 8.192897103059e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 6 SNES Function norm 1.635129206894e-02 0 KSP Residual norm 1.635129206894e-02 1 KSP Residual norm 9.049749505225e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.634946381608e-02 1 KSP Residual norm 9.038879304284e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.635085731721e-02 1 KSP Residual norm 9.050535714933e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 7 SNES Function norm 1.635085731721e-02 0 KSP Residual norm 1.635085731721e-02 1 KSP Residual norm 1.248783858778e-03 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.635003937816e-02 1 KSP Residual norm 1.248746884526e-03 Linear solve converged due to CONVERGED_RTOL iterations 1 8 SNES Function norm 1.635003937816e-02 0 KSP Residual norm 1.635003937816e-02 1 KSP Residual norm 1.631242803481e-02 2 KSP Residual norm 1.551854905139e-02 3 KSP Residual norm 7.586491369138e-03 4 KSP Residual norm 2.304721992110e-04 Linear solve converged due to CONVERGED_RTOL iterations 4 0 KSP Residual norm 1.635003936717e-02 1 KSP Residual norm 1.631242801091e-02 2 KSP Residual norm 1.551784322549e-02 3 KSP Residual norm 7.706018759197e-03 4 KSP Residual norm 2.652435205967e-04 Linear solve converged due to CONVERGED_RTOL iterations 4 9 SNES Function norm 1.282692400720e+06 Nonlinear solve did not converge due to DIVERGED_DTOL iterations 9 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 5.698410682663e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 8.618385788411e-02 1 KSP Residual norm 9.179214337527e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 8.618385788411e-02 0 KSP Residual norm 8.618385788411e-02 1 KSP Residual norm 9.659611612697e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.763975792565e-02 1 KSP Residual norm 2.674728432525e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 2.763975792565e-02 0 KSP Residual norm 2.763975792565e-02 1 KSP Residual norm 1.836844517206e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.773540750989e-03 1 KSP Residual norm 1.340353004310e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 3 SNES Function norm 2.773540750989e-03 0 KSP Residual norm 2.773540750989e-03 1 KSP Residual norm 5.382019472069e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 3.441017533038e-04 1 KSP Residual norm 5.218147408275e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 4 SNES Function norm 3.441017533038e-04 0 KSP Residual norm 3.441017533038e-04 1 KSP Residual norm 4.811955757758e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 3.991606932979e-06 1 KSP Residual norm 8.670005394864e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 5 SNES Function norm 3.991606932979e-06 0 KSP Residual norm 3.991606932979e-06 1 KSP Residual norm 8.920137762569e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 7.982782045945e-10 1 KSP Residual norm 3.185856555648e-12 Linear solve converged due to CONVERGED_RTOL iterations 1 6 SNES Function norm 7.982782045945e-10 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 6 0 SNES Function norm 3.475518085003e-04 0 KSP Residual norm 3.475518085003e-04 1 KSP Residual norm 2.328087365783e-04 2 KSP Residual norm 5.062641920970e-09 Linear solve converged due to CONVERGED_RTOL iterations 2 0 KSP Residual norm 3.256441726122e-02 1 KSP Residual norm 8.957860841275e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 3.256441726122e-02 0 KSP Residual norm 3.256441726122e-02 1 KSP Residual norm 5.886144296665e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 4.876417342600e-03 1 KSP Residual norm 5.064113356579e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 4.876417342600e-03 0 KSP Residual norm 4.876417342600e-03 1 KSP Residual norm 7.010281488349e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 6.144574635801e-04 1 KSP Residual norm 9.768499579865e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 3 SNES Function norm 6.144574635801e-04 0 KSP Residual norm 6.144574635801e-04 1 KSP Residual norm 8.324370478600e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 8.721913471286e-06 1 KSP Residual norm 1.370545428753e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 4 SNES Function norm 8.721913471286e-06 0 KSP Residual norm 8.721913471286e-06 1 KSP Residual norm 1.353680905072e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 3.399977638261e-09 1 KSP Residual norm 5.989016835446e-12 Linear solve converged due to CONVERGED_RTOL iterations 1 5 SNES Function norm 3.399977638261e-09 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 5 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 1.733028652796e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 7.490800616133e-03 1 KSP Residual norm 8.231034393757e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 7.490800616133e-03 0 KSP Residual norm 7.490800616133e-03 1 KSP Residual norm 8.349300836114e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.394207859364e-04 1 KSP Residual norm 5.991798196295e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.394207859364e-04 0 KSP Residual norm 1.394207859364e-04 1 KSP Residual norm 5.941984242417e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.985895331367e-08 1 KSP Residual norm 1.522240176483e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 3 SNES Function norm 2.985895331367e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 3 0 SNES Function norm 1.026110861069e-03 0 KSP Residual norm 1.026110861069e-03 1 KSP Residual norm 1.768297003974e-04 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 7.742520550463e-03 1 KSP Residual norm 4.379214132544e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 7.742520550463e-03 0 KSP Residual norm 7.742520550463e-03 1 KSP Residual norm 4.526371657285e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.574467856482e-05 1 KSP Residual norm 2.285626422848e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 2.574467856482e-05 0 KSP Residual norm 2.574467856482e-05 1 KSP Residual norm 2.274220735999e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 4.797870687759e-09 1 KSP Residual norm 2.841037171618e-12 Linear solve converged due to CONVERGED_RTOL iterations 1 3 SNES Function norm 4.797870687759e-09 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 3 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 5.979305049211e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 8.195748548281e-04 1 KSP Residual norm 1.017772307359e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 8.195748548281e-04 0 KSP Residual norm 8.195748548281e-04 1 KSP Residual norm 1.026961867317e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.717269750886e-06 1 KSP Residual norm 4.675717893044e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.717269750886e-06 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 1.325295461903e-03 0 KSP Residual norm 1.325295461903e-03 1 KSP Residual norm 7.421950279125e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.330639207844e-03 1 KSP Residual norm 3.764830059967e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 1.330639207844e-03 0 KSP Residual norm 1.330639207844e-03 1 KSP Residual norm 3.764587772248e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.146996236879e-07 1 KSP Residual norm 1.459926588396e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 2.146996236879e-07 0 KSP Residual norm 2.146996236879e-07 1 KSP Residual norm 1.470407475325e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.171052440924e-11 1 KSP Residual norm 1.397958787200e-15 Linear solve converged due to CONVERGED_RTOL iterations 1 3 SNES Function norm 2.171052440924e-11 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 3 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 3.067250018894e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.608120509564e-04 1 KSP Residual norm 2.171575389779e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 1.608120509564e-04 0 KSP Residual norm 1.608120509564e-04 1 KSP Residual norm 2.240227170220e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 9.106906778906e-08 1 KSP Residual norm 1.799636201836e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 9.106906778906e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 1.418151440967e-03 0 KSP Residual norm 1.418151440967e-03 1 KSP Residual norm 3.515151918636e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.938727941997e-04 1 KSP Residual norm 6.040565904176e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 2.938727941997e-04 0 KSP Residual norm 2.938727941997e-04 1 KSP Residual norm 6.089297895562e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.093551284505e-08 1 KSP Residual norm 2.891462873664e-12 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.093551284505e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 2.050867218313e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 3.054570961709e-05 1 KSP Residual norm 4.596885059357e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 3.054570961709e-05 0 KSP Residual norm 3.054570961709e-05 1 KSP Residual norm 4.828909801995e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.248104282690e-08 1 KSP Residual norm 8.831951540120e-13 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.248104282690e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 1.460190420374e-03 0 KSP Residual norm 1.460190420374e-03 1 KSP Residual norm 1.580519450452e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 5.884986293527e-05 1 KSP Residual norm 9.020937669734e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 5.884986293527e-05 0 KSP Residual norm 5.884986293527e-05 1 KSP Residual norm 9.054367061995e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.698647259662e-09 1 KSP Residual norm 1.767366286599e-13 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.698647259662e-09 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 1.834492501580e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 5.663374849972e-06 1 KSP Residual norm 9.752865227095e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 5.663374849972e-06 0 KSP Residual norm 5.663374849972e-06 1 KSP Residual norm 1.014567990183e-08 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.600155584863e-09 1 KSP Residual norm 8.580477366667e-14 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 2.600155584863e-09 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 1.479145244919e-03 0 KSP Residual norm 1.479145244919e-03 1 KSP Residual norm 6.804553002336e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.118291636442e-05 1 KSP Residual norm 1.397149309582e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 1.118291636442e-05 0 KSP Residual norm 1.118291636442e-05 1 KSP Residual norm 1.386385106529e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.800461695540e-10 1 KSP Residual norm 3.217319587127e-14 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 2.800461695540e-10 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 1.776151497065e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.076105708636e-06 1 KSP Residual norm 2.155085223926e-09 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 1.076105708636e-06 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 0 SNES Function norm 1.489507915803e-03 0 KSP Residual norm 1.489507915803e-03 1 KSP Residual norm 2.921538148329e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.172405782412e-06 1 KSP Residual norm 2.719699969205e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 2.172405782412e-06 0 KSP Residual norm 2.172405782412e-06 1 KSP Residual norm 2.731813669182e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 5.323734019212e-11 1 KSP Residual norm 6.095298326433e-15 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 5.323734019212e-11 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 1.746846293829e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.376406135686e-07 1 KSP Residual norm 5.310287317591e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 2.376406135686e-07 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 0 SNES Function norm 1.497239199277e-03 0 KSP Residual norm 1.497239199277e-03 1 KSP Residual norm 1.251830706962e-06 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 5.023668872625e-07 1 KSP Residual norm 7.656373075632e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 5.023668872625e-07 0 KSP Residual norm 5.023668872625e-07 1 KSP Residual norm 7.473091471413e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.003331377435e-11 1 KSP Residual norm 1.347408932811e-15 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 2.003331377435e-11 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 1.757419385096e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 7.310279661922e-08 1 KSP Residual norm 1.417918576189e-10 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 7.310279661922e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 0 SNES Function norm 1.503392083665e-03 0 KSP Residual norm 1.503392083665e-03 1 KSP Residual norm 5.327548050375e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.611770751552e-07 1 KSP Residual norm 2.611036263545e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 1.611770751552e-07 0 KSP Residual norm 1.611770751552e-07 1 KSP Residual norm 2.607692240234e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 1.627978707015e-11 1 KSP Residual norm 4.254078683385e-16 Linear solve converged due to CONVERGED_RTOL iterations 1 2 SNES Function norm 1.627978707015e-11 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2 0 SNES Function norm 2.854104379157e-02 0 KSP Residual norm 2.854104379157e-02 1 KSP Residual norm 1.747957348554e-05 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 2.999731614900e-08 1 KSP Residual norm 3.643836796648e-11 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 2.999731614900e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 0 SNES Function norm 1.507740473908e-03 0 KSP Residual norm 1.507740473908e-03 1 KSP Residual norm 2.263876724496e-07 Linear solve converged due to CONVERGED_RTOL iterations 1 0 KSP Residual norm 6.390049329095e-08 1 KSP Residual norm 9.745723464126e-12 Linear solve converged due to CONVERGED_RTOL iterations 1 1 SNES Function norm 6.390049329095e-08 Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1 [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: TSStep has failed due to DIVERGED_STEP_REJECTED

4 4

Bug report LMVM matrix class
by Stephan Köhler 09 May '23

09 May '23

Dear PETSc/Tao team, it seems to be that there is a bug in the LMVM matrix class: The function MatMultAdd_LMVM, see, e.g., https://petsc.org/release/src/ksp/ksp/utils/lmvm/lmvmimpl.c.html at line 114, if the vectors Y and Z are the same, then the result is wrong, since the first MatMult overwrites also the value in Y. Best regards Stephan Köhler -- Stephan Köhler TU Bergakademie Freiberg Institut für numerische Mathematik und Optimierung Akademiestraße 6 09599 Freiberg Gebäudeteil Mittelbau, Zimmer 2.07 Telefon: +49 (0)3731 39-3173 (Büro)

2 1

Question on ISLocalToGlobalMappingGetIndices Fortran Interface
by Danyang Su 08 May '23

08 May '23

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 (-11472655627481), https://www.mcs.anl.gov/petsc/petsc-3.14/src/ksp/ksp/tutorials/ex14f.F90.ht…. The value of idltog is not clear. In the latest PETSc version, this function can be called, but due to the extreme value of idltog, the code fails. I also tried to use ‘ISLocalToGlobalMappingGetIndicesF90(ltogm,ltog,ierr)’ () but no success. #if (PETSC_VERSION_MAJOR == 3 && PETSC_VERSION_MINOR <= 4) call DMDAGetGlobalIndicesF90(dmda_flow%da,PETSC_NULL_INTEGER, & idx,ierr) CHKERRQ(ierr) #else call DMGetLocalToGlobalMapping(dmda_flow%da,ltogm,ierr) CHKERRQ(ierr) call ISLocalToGlobalMappingGetIndices(ltogm,ltog,idltog,ierr) CHKERRQ(ierr) #endif dof = dmda_flow%dof do ivol = 1, nngl #if (PETSC_VERSION_MAJOR == 3 && PETSC_VERSION_MINOR <= 4) node_idx_lg2pg(ivol) = (idx(ivol*dof)+1)/dof #else node_idx_lg2pg(ivol) = (ltog(ivol*dof + idltog)+1)/dof #endif end do Any suggestions on that. Thanks, Danyang

2 2

Fortran preprocessor not work in pets-dev
by Danyang Su 08 May '23

08 May '23

Hi All, My code has some FPP. It works fine in PETSc 3.18 and earlier version, but stops working in the latest PETSc-Dev. For example the following FPP STANDARD_FORTRAN is not recognized. #ifdef STANDARD_FORTRAN 1 format(15x,1000a15) 2 format(1pe15.6e3,1000(1pe15.6e3)) #else 1 format(15x,<ncol>a15) 2 format(1pe15.6e3,<ncol>(1pe15.6e3)) #endif In the makefile, I define the preprocessor as PPFLAGS. PPFLAGS := -DLINUX -DRELEASE -DRELEASE_X64 -DSTANDARD_FORTRAN … exe: $(OBJS) chkopts -${FLINKER} $(FFLAGS) $(FPPFLAGS) $(CPPFLAGS) -o $(EXENAME) $(OBJS) ${PETSC_LIB} ${LIS_LIB} ${DLIB} ${SLIB} Any idea on this problem? All the best,

2 4

question about leap-frog for wave equation in petsc
by Huidong Yang 06 May '23

06 May '23

Hi Petsc developer. may I ask if there is any available implementations in petsc using leap-frog scheme? Thanks.

2 1

Is it necessary to call MatAssembly*() with MatSetValue()
by gaochenyi14 06 May '23

06 May '23

Hi, By `find` and `grep`, I find that in many PETSc examples `MatAssembly*()` are called after `MatSetValue()`. But in the C/Fortran API manual, the man page of `MatAssembly*()` does not say it is necessary for `MatSetValue()`. And the man page of `MatSetValue()` does not say it is a must to call `MatAssembly*()` afterwards. Is the man page of `MatSetValue()` up-to-date, or it lacks the reminder for the necessity of calling `MatAssembly*()`? Best regards, C.-Y. GAO

2 1

sources of floating point randomness in JFNK in serial
by Mark Lohry 05 May '23

05 May '23

I'm running multiple iterations of newtonls with an MFFD/JFNK nonlinear solver where I give it the sparsity. PC asm, KSP gmres, with SNESSetLagJacobian -2 (compute once and then frozen jacobian). I'm seeing slight (<1%) but nonzero differences in residuals from run to run. I'm wondering where randomness might enter here -- does the jacobian coloring use a random seed?

5 25

parallel computing error
by 권승리 / 학생 / 항공우주공학과 05 May '23

05 May '23

Dear developers I'm trying to use parallel computing and I ran the command 'mpirun -np 4 ./app' In this case, there are two problems. *First,* I encountered error message /// [0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [1]PETSC ERROR: Invalid argument [1]PETSC ERROR: Comm must be of size 1 /// The code on the error position is MatCreateSeqDense(PETSC_COMM_SELF, nns, ns, NULL, &Kns)); Could "MatCreateSeqDense" not be used in parallel computing? *Second*, the same error message is repeated as many times as the number of cores. if I use command -np 4, then the error message is repeated 4 times. Could you recommend some advice related to this? Best, Seung Lee Kwon -- Seung Lee Kwon, Ph.D.Candidate Aerospace Structures and Materials Laboratory Department of Mechanical and Aerospace Engineering Seoul National University Building 300 Rm 503, Gwanak-ro 1, Gwanak-gu, Seoul, South Korea, 08826 E-mail : ksl7912(a)snu.ac.kr Office : +82-2-880-7389 C. P : +82-10-4695-1062

4 11