Currently we don't support Kokkos dense matrix and its solvers. You can use MATSEQDENSECUDA/HIP
--Junchao Zhang On Wed, Jul 24, 2024 at 2:08 PM Barry Smith <bsm...@petsc.dev> wrote: > > For one MPI rank, it looks like you can use -pc_type cholesky > -pc_factor_mat_solver_type cupm though it is not documented in > https://urldefense.us/v3/__https://petsc.org/release/overview/linear_solve_table/*direct-solvers__;Iw!!G_uCfscf7eWS!fvASLTU48U_NIIf2O2CcYRSSki2GbUmrm4zw6CBOWx1rIY8-CnqmhboVIA5-aey5_QOPCQhaI2nbv4FJhCDdQMiJrHPc$ > > > Of if you also ./configure --download-kokkos --download-kokkos-kernels > you can use -pc_factor_mat_solver_type kokkos if you also this may also > work for multiple GPUs but that is not documented in the table either > (Junchao) Nor are sparse Kokkos or CUDA stuff documented (if they exist) in > the table. > > > Barry > > > > On Jul 24, 2024, at 2:44 PM, Sreeram R Venkat <srven...@utexas.edu> wrote: > > This Message Is From an External Sender > This message came from outside your organization. > I have an SPD dense matrix of size NxN, where N can range from 10^4-10^5. > Are there any Cholesky factorization/solve routines for it in PETSc (or in > any of the external libraries)? If possible, I want to use GPU acceleration > with 1 or more GPUs. The matrix type can be MATSEQDENSE/MATMPIDENSE or > MATSEQDENSECUDA/MATMPIDENSECUDA accordingly. If it is possible to do the > factorization beforehand and store it to do the triangular solves later, > that would be great. > > Thanks, > Sreeram > > >
