We present a collection of public-domain \mbox{Fortran 77} routines for the solution of systems of linear equations using a variety of iterative methods. The routines implement methods which have been modified for their efficient use on parallel architectures with either shared- or distributed-memory. PIM was designed to be portable across different machines. Results are presented for a variety of parallel computers.