Block-diag-ilu

======
- solve may now return non-zero exit (in C++, raises error in Python)
- check for NaN in solve
- check for zero diagonal elements in U

======
- LU bindings in Python

====
- Initial release with a minimal Python interface