Open
Description
Summary
Write a different version of this test within the ReSolve framework. Implement the necessary classes and functions on cpu/cuda/hip (or copy implementations as possible). Note that it is crucial to separate numeric and symbolic factorizations, because symbolic factorization is only called once.
Rationale
Cholesky factorization and solve are essential for HyKKT to work. A simpler version could be used for a single matrix in a series for other solvers for problems that are fully SPD.
Description
- Any originally cpu functions need only have cpu implementations
- Any originally gpu functions should have cpu/cuda/hip support.
- Write tests as within the rest of ReSolve.
- Write a test that generates an SPD matrix of arbitrary (square) size, multiplies by a known solution x (with current matvec in ReSolve) to get the right hand side b, and then solves and verifies that the solution is correct. The simplest solution would be to test
$||Ax-b||/||b||$
Additional information
No response