07/04/03: Initial version 10/28/09: Replaced GMP subset by MPI