 Cynopolis
|
ac4bef88fc
|
Fix SVD pipeline bugs; extract and unit-test bidiagonal block solvers
- Snapshot bidiagonal block before solving (W overwrite corrupted Ublock)
- Wide-matrix QL row extent fix (rowsQL = n when transposed)
- Vt output copy bound j<columns (was OOB write clobbering caller sigma)
- Rank-deficient blocks: residual singular values sigma_i = ||B*v_i||
instead of sqrt(eigs of B*B) (condition-number squaring); Gram-Schmidt
orthonormalization of U columns with rank-deficiency completion
- JacobiEigenSymmetric: return signed eigenvalues (removed fabsf)
- Fix 28 WithinRel(0.0f,.) matcher misuse -> WithinAbs in matrix-tests
- New unit tests: SolveBidiagonalBlock2x2 (7 scipy-referenced cases),
JacobiEigenSymmetric (3 cases), DeflateBidiagonal/BidiagonalIsDiagonal
|
2026-08-19 14:05:44 -04:00 |
|