Solving embarassingly parallel problems using Python multiprocessing

Question

How does one use multiprocessing to tackle embarrassingly parallel problems?

Embarassingly parallel problems typically consist of three basic parts:

1. Read input data (from a file, database, tcp connection, etc.).
2. Run calculations on the input data, where each calculation is independent of any other calculation.
3. Write results of calculations (to a file, database, tcp connection, etc.).

We can parallelize the program in two dimensions:

• Part 2 can run on multiple cores, since each calculation is independent; order of processing doesn't matter.
• Each part can run independently. Part 1 can place data on an input queue, part 2 can pull data off the input queue and put results onto an output queue, and part 3 can pull results off the output queue and write them out.

This seems a most basic pattern in concurrent programming, but I am still lost in trying to solve it, so let's write a canonical example to illustrate how this is done using multiprocessing.

Here is the example problem: Given a CSV file with rows of integers as input, compute their sums. Separate the problem into three parts, which can all run in parallel:

1. Process the input file into raw data (lists/iterables of integers)
2. Calculate the sums of the data, in parallel
3. Output the sums

Below is traditional, single-process bound Python program which solves these three tasks:

#!/usr/bin/env python # -*- coding: UTF-8 -*- # basicsums.py """A program that reads integer values from a CSV file and writes out their sums to another CSV file. """  import csv import optparse import sys  def make_cli_parser():     """Make the command line interface parser."""     usage = "\n\n".join(["python %prog INPUT_CSV OUTPUT_CSV",             __doc__,             """ ARGUMENTS:     INPUT_CSV: an input CSV file with rows of numbers     OUTPUT_CSV: an output file that will contain the sums\ """])     cli_parser = optparse.OptionParser(usage)     return cli_parser   def parse_input_csv(csvfile):     """Parses the input CSV and yields tuples with the index of the row     as the first element, and the integers of the row as the second     element.      The index is zero-index based.      :Parameters:     - `csvfile`: a `csv.reader` instance      """     for i, row in enumerate(csvfile):         row = [int(entry) for entry in row]         yield i, row   def sum_rows(rows):     """Yields a tuple with the index of each input list of integers     as the first element, and the sum of the list of integers as the     second element.      The index is zero-index based.      :Parameters:     - `rows`: an iterable of tuples, with the index of the original row       as the first element, and a list of integers as the second element      """     for i, row in rows:         yield i, sum(row)   def write_results(csvfile, results):     """Writes a series of results to an outfile, where the first column     is the index of the original row of data, and the second column is     the result of the calculation.      The index is zero-index based.      :Parameters:     - `csvfile`: a `csv.writer` instance to which to write results     - `results`: an iterable of tuples, with the index (zero-based) of       the original row as the first element, and the calculated result       from that row as the second element      """     for result_row in results:         csvfile.writerow(result_row)   def main(argv):     cli_parser = make_cli_parser()     opts, args = cli_parser.parse_args(argv)     if len(args) != 2:         cli_parser.error("Please provide an input file and output file.")     infile = open(args[0])     in_csvfile = csv.reader(infile)     outfile = open(args[1], 'w')     out_csvfile = csv.writer(outfile)     # gets an iterable of rows that's not yet evaluated     input_rows = parse_input_csv(in_csvfile)     # sends the rows iterable to sum_rows() for results iterable, but     # still not evaluated     result_rows = sum_rows(input_rows)     # finally evaluation takes place as a chain in write_results()     write_results(out_csvfile, result_rows)     infile.close()     outfile.close()   if __name__ == '__main__':     main(sys.argv[1:]) 

Let's take this program and rewrite it to use multiprocessing to parallelize the three parts outlined above. Below is a skeleton of this new, parallelized program, that needs to be fleshed out to address the parts in the comments:

#!/usr/bin/env python # -*- coding: UTF-8 -*- # multiproc_sums.py """A program that reads integer values from a CSV file and writes out their sums to another CSV file, using multiple processes if desired. """  import csv import multiprocessing import optparse import sys  NUM_PROCS = multiprocessing.cpu_count()  def make_cli_parser():     """Make the command line interface parser."""     usage = "\n\n".join(["python %prog INPUT_CSV OUTPUT_CSV",             __doc__,             """ ARGUMENTS:     INPUT_CSV: an input CSV file with rows of numbers     OUTPUT_CSV: an output file that will contain the sums\ """])     cli_parser = optparse.OptionParser(usage)     cli_parser.add_option('-n', '--numprocs', type='int',             default=NUM_PROCS,             help="Number of processes to launch [DEFAULT: %default]")     return cli_parser   def main(argv):     cli_parser = make_cli_parser()     opts, args = cli_parser.parse_args(argv)     if len(args) != 2:         cli_parser.error("Please provide an input file and output file.")     infile = open(args[0])     in_csvfile = csv.reader(infile)     outfile = open(args[1], 'w')     out_csvfile = csv.writer(outfile)      # Parse the input file and add the parsed data to a queue for     # processing, possibly chunking to decrease communication between     # processes.      # Process the parsed data as soon as any (chunks) appear on the     # queue, using as many processes as allotted by the user     # (opts.numprocs); place results on a queue for output.     #     # Terminate processes when the parser stops putting data in the     # input queue.      # Write the results to disk as soon as they appear on the output     # queue.      # Ensure all child processes have terminated.      # Clean up files.     infile.close()     outfile.close()   if __name__ == '__main__':     main(sys.argv[1:]) 

These pieces of code, as well as another piece of code that can generate example CSV files for testing purposes, can be found on github.

I would appreciate any insight here as to how you concurrency gurus would approach this problem.

---

Here are some questions I had when thinking about this problem. Bonus points for addressing any/all:

• Should I have child processes for reading in the data and placing it into the queue, or can the main process do this without blocking until all input is read?
• Likewise, should I have a child process for writing the results out from the processed queue, or can the main process do this without having to wait for all the results?
• Should I use a processes pool for the sum operations?
  • If yes, what method do I call on the pool to get it to start processing the results coming into the input queue, without blocking the input and output processes, too? apply_async()? map_async()? imap()? imap_unordered()?
• Suppose we didn't need to siphon off the input and output queues as data entered them, but could wait until all input was parsed and all results were calculated (e.g., because we know all the input and output will fit in system memory). Should we change the algorithm in any way (e.g., not run any processes concurrently with I/O)?

Answer 1

My solution has an extra bell and whistle to make sure that the order of the output has the same as the order of the input. I use multiprocessing.queue's to send data between processes, sending stop messages so each process knows to quit checking the queues. I think the comments in the source should make it clear what's going on but if not let me know.

#!/usr/bin/env python # -*- coding: UTF-8 -*- # multiproc_sums.py """A program that reads integer values from a CSV file and writes out their sums to another CSV file, using multiple processes if desired. """  import csv import multiprocessing import optparse import sys  NUM_PROCS = multiprocessing.cpu_count()  def make_cli_parser():     """Make the command line interface parser."""     usage = "\n\n".join(["python %prog INPUT_CSV OUTPUT_CSV",             __doc__,             """ ARGUMENTS:     INPUT_CSV: an input CSV file with rows of numbers     OUTPUT_CSV: an output file that will contain the sums\ """])     cli_parser = optparse.OptionParser(usage)     cli_parser.add_option('-n', '--numprocs', type='int',             default=NUM_PROCS,             help="Number of processes to launch [DEFAULT: %default]")     return cli_parser  class CSVWorker(object):     def __init__(self, numprocs, infile, outfile):         self.numprocs = numprocs         self.infile = open(infile)         self.outfile = outfile         self.in_csvfile = csv.reader(self.infile)         self.inq = multiprocessing.Queue()         self.outq = multiprocessing.Queue()          self.pin = multiprocessing.Process(target=self.parse_input_csv, args=())         self.pout = multiprocessing.Process(target=self.write_output_csv, args=())         self.ps = [ multiprocessing.Process(target=self.sum_row, args=())                         for i in range(self.numprocs)]          self.pin.start()         self.pout.start()         for p in self.ps:             p.start()          self.pin.join()         i = 0         for p in self.ps:             p.join()             print "Done", i             i += 1          self.pout.join()         self.infile.close()      def parse_input_csv(self):             """Parses the input CSV and yields tuples with the index of the row             as the first element, and the integers of the row as the second             element.              The index is zero-index based.              The data is then sent over inqueue for the workers to do their             thing.  At the end the input process sends a 'STOP' message for each             worker.             """             for i, row in enumerate(self.in_csvfile):                 row = [ int(entry) for entry in row ]                 self.inq.put( (i, row) )              for i in range(self.numprocs):                 self.inq.put("STOP")      def sum_row(self):         """         Workers. Consume inq and produce answers on outq         """         tot = 0         for i, row in iter(self.inq.get, "STOP"):                 self.outq.put( (i, sum(row)) )         self.outq.put("STOP")      def write_output_csv(self):         """         Open outgoing csv file then start reading outq for answers         Since I chose to make sure output was synchronized to the input there         is some extra goodies to do that.          Obviously your input has the original row number so this is not         required.         """         cur = 0         stop = 0         buffer = {}         # For some reason csv.writer works badly across processes so open/close         # and use it all in the same process or else you'll have the last         # several rows missing         outfile = open(self.outfile, "w")         self.out_csvfile = csv.writer(outfile)          #Keep running until we see numprocs STOP messages         for works in range(self.numprocs):             for i, val in iter(self.outq.get, "STOP"):                 # verify rows are in order, if not save in buffer                 if i != cur:                     buffer[i] = val                 else:                     #if yes are write it out and make sure no waiting rows exist                     self.out_csvfile.writerow( [i, val] )                     cur += 1                     while cur in buffer:                         self.out_csvfile.writerow([ cur, buffer[cur] ])                         del buffer[cur]                         cur += 1          outfile.close()  def main(argv):     cli_parser = make_cli_parser()     opts, args = cli_parser.parse_args(argv)     if len(args) != 2:         cli_parser.error("Please provide an input file and output file.")      c = CSVWorker(opts.numprocs, args[0], args[1])  if __name__ == '__main__':     main(sys.argv[1:]) 

Answer 2

Coming late to the party...

joblib has a layer on top of multiprocessing to help making parallel for loops. It gives you facilities like a lazy dispatching of jobs, and better error reporting in addition to its very simple syntax.

As a disclaimer, I am the original author of joblib.