Re: [petsc-users] BAIJCUSPARSE?

We will let you know when this is ready, Xiangdong. Let me address a part of your original question that I don't think anyone else noticed: In my current code, the Jacobian matrix preallocated and assembled as BAIJ format. Do I have to rewrite this part of code to preallocate and assemble the

Re: [petsc-users] PETSc 3.12 with .f90 files

On Tue, 29 Oct 2019, Matthew Knepley via petsc-users wrote: > On Tue, Oct 29, 2019 at 2:35 PM Smith, Barry F. wrote: > > > > >The problem is that this change DOES use the preprocessor on the f90 > > file, does it not? We need a rule that does not use the preprocessor. > > > > This change

Re: [petsc-users] PETSc 3.12 with .f90 files

The problem is that this change DOES use the preprocessor on the f90 file, does it not? We need a rule that does not use the preprocessor. Barry > On Oct 29, 2019, at 10:50 AM, Matthew Knepley via petsc-users > wrote: > > On Tue, Oct 29, 2019 at 11:38 AM Randall Mackie wrote: > Hi

Re: [petsc-users] PETSc 3.12 with .f90 files

On Tue, 29 Oct 2019, Randall Mackie via petsc-users wrote: > Dear PETSc users: > > In our code, we have one or two small .f90 files that are part of the > software, and they have always compiled without any issues with previous > versions of PETSc, using standard PETSc make files. > >

Re: [petsc-users] PETSc 3.12 with .f90 files

On Tue, Oct 29, 2019 at 10:54 AM Randall Mackie via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear PETSc users: > > In our code, we have one or two small .f90 files that are part of the > software, and they have always compiled without any issues with previous > versions of PETSc, using

[petsc-users] PETSc 3.12 with .f90 files

Dear PETSc users: In our code, we have one or two small .f90 files that are part of the software, and they have always compiled without any issues with previous versions of PETSc, using standard PETSc make files. However, starting with PETSc 3.12, they no longer compile. Was there some

Re: [petsc-users] VI: RS vs SS

Hi Todd, On October 29, 2019 at 7:15:50 AM, Munson, Todd (tmun...@mcs.anl.gov) wrote: Hi, Is the matrix for the linear PDE symmetric? If so, then the VI is equivalent to finding the stationary points of a bound-constrained quadratic program and you may want to use

Re: [petsc-users] calculating eigenvalues for the Stokes equations in 3D

Patrick, thank you very much for your comment! I'll try to use the example provided by Matt as a starting point. I've never used SLEPc before, so my knowledge of the subject is a bit insufficient, to say the least. 29.10.2019, 16:11, "Patrick Sanan" :"Solve the forward problem" means "solve the

Re: [petsc-users] calculating eigenvalues for the Stokes equations in 3D

Thank you, Matthew! I'll give it a try. 29.10.2019, 13:30, "Matthew Knepley" :On Tue, Oct 29, 2019 at 5:06 AM Dave May via petsc-users wrote:  On Tue, 29 Oct 2019 at 08:59, wrote:Hi, Dave! Thank you for your assistance. The problem is that I don't

Re: [petsc-users] calculating eigenvalues for the Stokes equations in 3D

Sorry, perhaps I stated my problem a bit inaccurately and that may have caused misunderstanding. In fact, I need to solve the forward problem, i.e. I have to find eigenvalues for the three dimensional Stokes equations, but I don't know how to do it in SLEPc. 29.10.2019, 12:05, "Dave May" :  On

Re: [petsc-users] VI: RS vs SS

Hi, Is the matrix for the linear PDE symmetric? If so, then the VI is equivalent to finding the stationary points of a bound-constrained quadratic program and you may want to use the TAO Newton Trust-Region or Line-Search methods for bound-constrained optimization problems. Alp: are there

Re: [petsc-users] calculating eigenvalues for the Stokes equations in 3D

On Tue, Oct 29, 2019 at 5:06 AM Dave May via petsc-users < petsc-users@mcs.anl.gov> wrote: > > > On Tue, 29 Oct 2019 at 08:59, wrote: > >> Hi, Dave! Thank you for your assistance. The problem is that I don't have >> the matrix representation for my problem. >> > > You will have to explain in

Re: [petsc-users] calculating eigenvalues for the Stokes equations in 3D

Hi, Dave! Thank you for your assistance. The problem is that I don't have the matrix representation for my problem.  I've searched on the internet for a solution, but found only some examples of solving 2D eigenvalue problems for the Stokes equations. Maybe you could help me with the matrix

[petsc-users] calculating eigenvalues for the Stokes equations in 3D

Hello!In my project I have to calculate eigenvalues for the three dimensional Stokes equations with Dirichlet boundary conditions. I would like to use SLEPc iterative solver for this purpose, but I'm not quite sure how I should set up the matrices for the solver. I would be grateful if anyone