Update to library to use generic factorization method