Modifying Bunch-Kaufman to return the subdiagonal entries of D so that only L's main diagonal need be implicitly modified, and returning pivot vector in standard form which does not use negative entries to signal 2x2 pivots. QuasiDiagonalSolve was then added, as well as ldl::SolveAfter for BunchKaufman.