Simplifying LowerSolve and LowerMultiply by forcing the frontal trees to explicitly store ones on the diagonal when the triangular matrices are unit-diagonal (and adding the appropriate extra member variables for storing the diagonal explicitly). Only the workspaces are now mutable.