Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Barry Smith
, 2024 4:43 PM > To: Yesypenko, Anna mailto:a...@oden.utexas.edu>> > Cc: petsc-users@mcs.anl.gov <mailto:petsc-users@mcs.anl.gov> > mailto:petsc-users@mcs.anl.gov>>; Victor Eijkhout > mailto:eijkh...@tacc.utexas.edu>> > Subject: Re: [petsc-users] HashMap

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Yesypenko, Anna
Hi Barry, I'm using version 3.20.3. The tacc system is lonestar6. Best, Anna From: Barry Smith Sent: Thursday, January 18, 2024 4:43 PM To: Yesypenko, Anna Cc: petsc-users@mcs.anl.gov ; Victor Eijkhout Subject: Re: [petsc-users] HashMap Error when populating

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Barry Smith
-1, 2) >>> >>> for index in range(1, n - 1): >>> A.setValue(index, index - 1, -1) >>> A.setValue(index, index, 2) >>> A.setValue(index, index + 1, -1) >>> A.assemble() >>> ``` >>> If it means anything to you,

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Yesypenko, Anna
From: Barry Smith Sent: Thursday, January 18, 2024 3:38 PM To: Yesypenko, Anna Cc: petsc-users@mcs.anl.gov ; Victor Eijkhout Subject: Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix It is using the hash map system for inserting values which only inserts on the CPU

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Barry Smith
es. >> >> Thank you for your help and suggestions! >> Anna >> >> From: Barry Smith mailto:bsm...@petsc.dev>> >> Sent: Thursday, January 18, 2024 2:35 PM >> To: Yesypenko, Anna mailto:a...@oden.utexas.edu>> >> Cc: petsc-users@mcs.anl.gov <

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Matthew Knepley
Sent:* Thursday, January 18, 2024 2:35 PM > *To:* Yesypenko, Anna > *Cc:* petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] HashMap Error when populating AIJCUSPARSE > matrix > > >Do you ever get a problem with 'aij` ? Can you run in a loop with > 'aij' to confirm it d

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Yesypenko, Anna
From: Barry Smith Sent: Thursday, January 18, 2024 2:35 PM To: Yesypenko, Anna Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix Do you ever get a problem with 'aij` ? Can you run in a loop with 'aij

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Barry Smith
It appears to be crashing in kh_resize() in khash.h on a memory allocation failure when it tries to get additional memory for storing the matrix. This code seems to be only using the CPU memory so it should also fail in a similar way with 'aij'. But the matrix is not large and so I

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Barry Smith
Do you ever get a problem with 'aij` ? Can you run in a loop with 'aij' to confirm it doesn't fail then? Barry > On Jan 17, 2024, at 4:51 PM, Yesypenko, Anna wrote: > > Dear Petsc users/developers, > > I'm experiencing a bug when using petsc4py with GPU support. It may be my

Re: [petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Matthew Knepley
On Thu, Jan 18, 2024 at 9:04 AM Yesypenko, Anna wrote: > Dear Petsc users/developers, > > I'm experiencing a bug when using petsc4py with GPU support. It may be my > mistake in how I set up a AIJCUSPARSE matrix. > For larger matrices, I sometimes encounter a error in assigning matrix > values;

[petsc-users] HashMap Error when populating AIJCUSPARSE matrix

2024-01-18 Thread Yesypenko, Anna
Dear Petsc users/developers, I'm experiencing a bug when using petsc4py with GPU support. It may be my mistake in how I set up a AIJCUSPARSE matrix. For larger matrices, I sometimes encounter a error in assigning matrix values; the error is thrown in PetscHMapIJVQuerySet(). Here is a minimum