Added threads options for easy hybrid MPI/openmp
Some documentation, bug fix.