Open tbetcke opened 2 years ago
Once we have a triangular solver, I might be motivated to implement a dense LDL (https://en.wikipedia.org/wiki/Cholesky_decomposition#LDL_decomposition) decomposition, to solve linear systems Ax=b, with A being positive semi-definite.
Sounds good. My goal for Householder is that eventually we have modern task based algorithms for all decompositions. For LDL^T see for example (https://hal.inria.fr/hal-00809663/document). However,
Triangular solves are a core task of many more advanced matrix algorithms. We want to implement.
The multithreaded version also requires serial triangular solves. So it is natural that both needs to be done. A precondition for this issue is a decision on the parallelization framework in #11 .