Re: [petsc-dev] [petsc-users] CUDA MatSetValues test

2021-06-02 Thread Stefano Zampini
Barry If I remember correctly, the ctable business was a limitation of Mark's initial implementation, which I removed here https://gitlab.com/petsc/petsc/-/merge_requests/3411/diffs?commit_id=36173707431f5583889720868beae9046e6e1fd2. Here colmap was "int" already. We made a mistake letting this

Re: [petsc-dev] [petsc-users] CUDA MatSetValues test

2021-06-02 Thread Barry Smith
Switching the thread to petsc-dev since users don't need this. Had a nasty failure on one the CI, turned out the GPU colmap was declared as PetscInt* but it was copied down from the CPU into the GPU memory as int. The power of void* interface. Thus ex5cu.cu which was never CI tested was

Re: [petsc-dev] I think main is broken

2021-06-02 Thread Satish Balay via petsc-dev
Fixes at https://gitlab.com/petsc/petsc/-/merge_requests/4059 Satish On Wed, 2 Jun 2021, Matthew Knepley wrote: > I just rebased on main and get this: > > config/BuildSystem/config/packages/kokkos-kernels.py:7: self.gitcommit = > '54e4213ca028163a76639d23b1204b213203bc79' # origin/develop of

Re: [petsc-dev] I think main is broken

2021-06-02 Thread Matthew Knepley
Thanks! Matt On Wed, Jun 2, 2021 at 11:14 AM Satish Balay wrote: > Fixes at https://gitlab.com/petsc/petsc/-/merge_requests/4059 > > Satish > > On Wed, 2 Jun 2021, Matthew Knepley wrote: > > > I just rebased on main and get this: > > > >

[petsc-dev] I think main is broken

2021-06-02 Thread Matthew Knepley
I just rebased on main and get this: config/BuildSystem/config/packages/kokkos-kernels.py:7: self.gitcommit = '54e4213ca028163a76639d23b1204b213203bc79' # origin/develop of 2021-04-28 569 config/BuildSystem/config/packages/kokkos.py:7:

[petsc-dev] Petsc Parallel Matrix Vector Multiplication for Skew-Symmetric Matrix

2021-06-02 Thread Yağmur Aksan
Hello, I need a function that parallelly multiplies skew-symmetric matrix and vector. Is there any function that Petsc provides? To do so, I think of using matMult function in Petsc and changing it to work with skew-symmetric matrices. But, I haven't found the implementation of MatMult function.

Re: [petsc-dev] [petsc-users] CUDA MatSetValues test

2021-06-02 Thread Barry Smith
> On Jun 2, 2021, at 2:06 AM, Stefano Zampini wrote: > > Barry > > If I remember correctly, the ctable business was a limitation of Mark's > initial implementation, which I removed here >

Re: [petsc-dev] Petsc Parallel Matrix Vector Multiplication for Skew-Symmetric Matrix

2021-06-02 Thread Matthew Knepley
On Wed, Jun 2, 2021 at 3:23 PM Yağmur Aksan wrote: > Hello, > > I need a function that parallelly multiplies skew-symmetric matrix and > vector. Is there any function that Petsc provides? > To do so, I think of using matMult function in Petsc and changing it to > work with skew-symmetric

Re: [petsc-dev] Petsc Parallel Matrix Vector Multiplication for Skew-Symmetric Matrix

2021-06-02 Thread Barry Smith
The easiest way to find implementations on the web is to type the function name MatMult() in google and look at the bottom of the page where it has links to all the implementations. Of course you have to guess which implementation you need based on the shortened names. Here is the one you

Re: [petsc-dev] Petsc Parallel Matrix Vector Multiplication for Skew-Symmetric Matrix

2021-06-02 Thread Mark Adams
I guess we still have a MATSBAIJ. In theory you could modify the MatVec to do skew-symmetric, easily, but this is such a niche operation (never seen it requested in 20+ years) that I would not want to put it in the repo. This would save <2x memory and <2x time. If you have a blocked system then it