[petsc-users] BAIJCUSPARSE?

[Smith, Barry F.]

  You would need to investigate if the Nvidia cuSPARSE package supports such a format. If it does then it would be reasonably straightforward for you to hook up the required interface from PETSc. If it does not then it is a massive job to provide such code and you should see if any open source packages provide such CUDA support and then you could hook PETSc up to use that.

  Barry


> On Oct 25, 2019, at 3:43 PM, Xiangdong via petsc-users <petsc-users at mcs.anl.gov> wrote:
> 
> Can anyone comment on the PETSc's GPU version of Block CSR, say BAIJCUSPARSE? Does it make sense to have such format on GPU? Is it under development?
> 
> Thank you.
> 
> Xiangdong
> 
> On Wed, Oct 23, 2019 at 11:36 AM Xiangdong <epscodes at gmail.com> wrote:
> Hello everyone,
> 
> I am wondering whether there is a format BAIJCUSPARSE for Block CSR on GPU.
> 
> 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 matrix as AIJ in order to use aijcusparse?
> 
> Thank you.
> 
> Xiangdong
>