Add implementation of ILUT(p, tau) factorization

This allows us to study how ILUT(p, tau) factorization performs. In
the next step we might want to see how SPIKE can be combined with
ILUT factorization.