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 > <[email protected]> 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 <[email protected]> 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 >
