Reducible polynomials in nf

Currently ``nf`` allows constructing number fields with a reducible defining polynomial, and operations on number field elements don't perform error handling, e.g. ``nf_elem_inv`` silently returns a wrong result.

This all ought to be done more robustly and explicitly documented.