Re: [petsc-users] Python PETSc performance vs scipy ZVODE

2023-08-14 Thread Zhang, Hong via petsc-users
PETSs is not necessarily faster than scipy for your problem when executed in serial. But you get benefits when running in parallel. The PETSc code you wrote uses float64 while your scipy code uses complex128, so the comparison may not be fair. In addition, using the RHS Jacobian does not necess

Re: [petsc-users] DM/DS crash after update to 3.19 (Fortran)

2023-08-14 Thread Matthew Knepley
On Tue, Aug 15, 2023 at 2:06 AM Martin Diehl wrote: > Dear PETSc team, > > my simulation crashes after updating from 3.18.5 to 3.19.4. > > The error message is attached, so is the main code. The mesh (variable > named geomMesh) is read with DMPlexCreateFromFile in a different part > of the code).

Re: [petsc-users] CUDA error trying to run a job with two mpi processes and 1 GPU

2023-08-14 Thread Junchao Zhang
I don't see a problem in the matrix assembly. If you point me to your repo and show me how to build it, I can try to reproduce. --Junchao Zhang On Mon, Aug 14, 2023 at 2:53 PM Vanella, Marcos (Fed) < marcos.vane...@nist.gov> wrote: > Hi Junchao, I've tried for my case using the -ksp_type gmres

Re: [petsc-users] performance regression with GAMG

2023-08-14 Thread Mark Adams
On Mon, Aug 14, 2023 at 11:03 AM Stephan Kramer wrote: > Many thanks for looking into this, Mark > > My 3D tests were not that different and I see you lowered the threshold. > > Note, you can set the threshold to zero, but your test is running so much > > differently than mine there is something

Re: [petsc-users] CUDA error trying to run a job with two mpi processes and 1 GPU

2023-08-14 Thread Junchao Zhang
Yeah, it looks like ex60 was run correctly. Double check your code again and if you still run into errors, we can try to reproduce on our end. Thanks. --Junchao Zhang On Mon, Aug 14, 2023 at 1:05 PM Vanella, Marcos (Fed) < marcos.vane...@nist.gov> wrote: > Hi Junchao, I compiled and run ex60 th

[petsc-users] DM/DS crash after update to 3.19 (Fortran)

2023-08-14 Thread Martin Diehl
Dear PETSc team, my simulation crashes after updating from 3.18.5 to 3.19.4. The error message is attached, so is the main code. The mesh (variable named geomMesh) is read with DMPlexCreateFromFile in a different part of the code).  I did not start serious debugging yet in the hope that you can p

Re: [petsc-users] performance regression with GAMG

2023-08-14 Thread Stephan Kramer
Many thanks for looking into this, Mark My 3D tests were not that different and I see you lowered the threshold. Note, you can set the threshold to zero, but your test is running so much differently than mine there is something else going on. Note, the new, bad, coarsening rate of 30:1 is what we

Re: [petsc-users] eigenvalue problem involving inverse of a matrix

2023-08-14 Thread maitri ksh
got it, thanks Pierre & Jose. On Mon, Aug 14, 2023 at 12:50 PM Jose E. Roman wrote: > See for instance ex3.c and ex9.c > https://slepc.upv.es/documentation/current/src/eps/tutorials/index.html > > Jose > > > > El 14 ago 2023, a las 10:45, Pierre Jolivet > escribió: > > > > > > > >> On 14 Aug 20

Re: [petsc-users] eigenvalue problem involving inverse of a matrix

2023-08-14 Thread Jose E. Roman
See for instance ex3.c and ex9.c https://slepc.upv.es/documentation/current/src/eps/tutorials/index.html Jose > El 14 ago 2023, a las 10:45, Pierre Jolivet escribió: > > > >> On 14 Aug 2023, at 10:39 AM, maitri ksh wrote: >> >>  >> Hi, >> I need to solve an eigenvalue problem Ax=lmbda*x

Re: [petsc-users] eigenvalue problem involving inverse of a matrix

2023-08-14 Thread Pierre Jolivet
> On 14 Aug 2023, at 10:39 AM, maitri ksh wrote: > >  > Hi, > I need to solve an eigenvalue problem Ax=lmbda*x, where A=(B^-H)*Q*B^-1 is a > hermitian matrix, 'B^-H' refers to the hermitian of the inverse of the matrix > B. Theoretically it would take around 1.8TB to explicitly compute the

[petsc-users] eigenvalue problem involving inverse of a matrix

2023-08-14 Thread maitri ksh
Hi, I need to solve an eigenvalue problem *Ax=lmbda*x*, where A=(B^-H)*Q*B^-1 is a hermitian matrix, 'B^-H' refers to the hermitian of the inverse of the matrix B. Theoretically it would take around 1.8TB to explicitly compute the matrix B^-1 . A feasible way to solve this eigenvalue problem would