linbox icon indicating copy to clipboard operation
linbox copied to clipboard

Do not compute inverse for Blas Dixon

Open jgdumas opened this issue 6 years ago • 1 comments

No need to compute the inverse mod p to lift in Dixon. The PLUQ factorization is sufficient to solve the system at each iteration.

jgdumas avatar Dec 07 '18 10:12 jgdumas

It may depend on the dimension/magnitude: although the overall arithmetic count is always in favor of the PLUQ variant, solving with a PLUQ at each iteration may be expensive: poorer cache behavious of the permutations and 2 trsv are less efficient than 1 fgemm (due to more modular reductions). Some benchmarking is required for a comparison.

ClementPernet avatar Dec 07 '18 10:12 ClementPernet