Module solvers 

GPU-ready iterative sparse linear solvers.

Implements Conjugate Gradient (CG) for symmetric positive definite systems and BiCGSTAB for general non-symmetric systems. Both use a BLAS-friendly memory layout and an optional Jacobi (diagonal) preconditioner. The CPU computation paths mirror what GPU compute shaders would execute.

Structs

GpuSolverConfig

Configuration for GPU iterative solvers.

SolverResult

Result returned by iterative solvers.

Enums

GpuSolverBackend

GPU compute backend used by the solvers.

Functions

bicgstab_csr

Solve a general (possibly non-symmetric) system A * x = b using Bi-Conjugate Gradient Stabilised (BiCGSTAB).

cg_csr

Solve the symmetric positive definite system A * x = b using the Preconditioned Conjugate Gradient method.