Re: [petsc-users] Fortran preprocessor not work in pets-dev

2023-05-06 Thread Satish Balay via petsc-users
On Sat, 6 May 2023, Danyang Su wrote: > 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 form

Re: [petsc-users] Step size setting in TS

2023-05-06 Thread Matthew Knepley
On Sat, May 6, 2023 at 7:25 PM Jorti, Zakariae via petsc-users < petsc-users@mcs.anl.gov> wrote: > 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.

[petsc-users] Step size setting in TS

2023-05-06 Thread Jorti, Zakariae via petsc-users
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 th

Re: [petsc-users] question about leap-frog for wave equation in petsc

2023-05-06 Thread Matthew Knepley
On Sat, May 6, 2023 at 11:47 AM Huidong Yang wrote: > Hi Petsc developer. > > may I ask if there is any available implementations in petsc > using leap-frog scheme? > I don't think we have leapfrog, but we do have Stormer-Verlet, which is also a 2nd order symplectic method. Thanks, Matt

[petsc-users] question about leap-frog for wave equation in petsc

2023-05-06 Thread Huidong Yang
Hi Petsc developer. may I ask if there is any available implementations in petsc using leap-frog scheme? Thanks.

Re: [petsc-users] Is it necessary to call MatAssembly*() with MatSetValue()

2023-05-06 Thread Jose E. Roman
Fixed in https://gitlab.com/petsc/petsc/-/merge_requests/6423 Jose > El 6 may 2023, a las 9:16, gaochenyi14 escribió: > > 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

[petsc-users] Is it necessary to call MatAssembly*() with MatSetValue()

2023-05-06 Thread gaochenyi14
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 c