Xiangdong, cuSPARSE does support block compressed sparse row (BAIJ) format, but we don't currently support that cuSPARSE functionality in PETSc. It should be easy to add, but we are currently refactoring the way we interface with third party GPU libraries such as cuSPARSE, and it would probably make more sense to add this support after that refactor is done. Do you need this right away, or could it wait maybe a few weeks until this is completed?
Best regards, Richard On Fri, Oct 25, 2019 at 1:50 PM Smith, Barry F. via petsc-users <petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>> wrote: 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@mcs.anl.gov<mailto:petsc-users@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 > <epsco...@gmail.com<mailto:epsco...@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 >
