[needs testing] make distutils compile Cython extensions in parallel

[craigcitro]

Currently, we've got excellent code to run Cython in parallel as part of the Sage build process. However, once we ask distutils to begin compiling that code, everything is done in serial, because distutils works solely in serial.

The Code

The attached file changes that. This (somewhat brutally) hacks distutils to dispatch the calls to build C/C++ extensions in parallel, using pyprocessing. (I totally jacked William's code from #3765 for this.) Here's how to put this code in place:

• download the attached file (build_ext.py)
• replace $SAGE_ROOT/local/lib/python2.5/distutils/command/build_ext.py with the new version.

Now, if you want to test the new code, do the following:

• set the environment variable SAGE_PARALLEL_DIST to something. (The code just checks to see if the variable is defined at all.)
• set the environment variable MAKE to MAKE -j2, where 2 is replaced by the number of simultaneous build processes you want.
• build.

Notes

If you want to test this, don't go around touching the .pyx files in the Sage library, since Cython is much slower than gcc. Instead, simply go around touching the .c and .cpp files in $SAGE_ROOT/devel/sage/sage. One of the cool features we added with the new build system is that these files get recompiled when they change.

There is now a line that prints as part of the build process that looks like:

Total time spent compiling C/C++ extensions:  5.2876701355 seconds.

So try touching a bunch of files in the Sage library, and seeing what kind of speedups you get.

There are two caveats I want to offer with this code:

• Michael points out that numpy does a lot with distutils. I could very well have broken their use of distutils.
• I don't do anything involving dependency tracking between extensions. In particular, if there are extensions that have to be built in a certain order, this code could break it. (This code still compiles all source files within an extension in the usual way.) Neither Michael nor I could think of a situation where this would break anything in Sage, but who knows ...

The Plan

So while this code is cool, and will definitely save a ton of CPU time on sage.math (cough mabshoff cough), the plan is not to maintain it as a part of the Sage Python spkg indefinitely. Instead, if it seems to work well, then we should try to clean this code up a bit more and upstream it, since pyprocessing is standard in Python 2.6.

Component: build

Author: Craig Citro

Merged: 4.0.1.rc3

Issue created by migration from https://trac.sagemath.org/ticket/4652

[williamstein]

comment:1

Attachment: build_ext.py.gz

Just out of curiosity, did you try to get this to work without modifying build_ext.py at all, but by simply changing stuff in that module at runtime, like I illustrate in #3765?

Here's the diff of what you did:

--- build_ext.py.orig	2008-11-29 10:32:57.000000000 -0800
+++ build_ext.py	2008-11-29 03:38:35.000000000 -0800
@@ -409,13 +409,57 @@
     # get_outputs ()
 
     def build_extensions(self):
+
         # First, sanity-check the 'extensions' list
         self.check_extensions_list(self.extensions)
 
-        for ext in self.extensions:
-            self.build_extension(ext)
+        if not os.environ.has_key('MAKE'):
+            ncpus = 1
+        else:
+            MAKE = os.environ['MAKE']
+            z = [w[2:] for w in MAKE.split() if w.startswith('-j')]
+            if len(z) == 0:  # no command line option
+                ncpus = 1
+            else:
+                # Determine number of cpus from command line argument.
+                # Also, use the OS to cap the number of cpus, in case
+                # user annoyingly makes a typo and asks to use 10000
+                # cpus at once.
+                try:
+                    ncpus = int(z[0])
+                    n = 2*number_of_cpus()
+                    if n:  # prevent dumb typos.
+                        ncpus = min(ncpus, n)
+                except ValueError:
+                    ncpus = 1
+
+        import time
+        t = time.time()
+        
+        # See if we're trying out the experimental parallel build
+        # code.
+        if ncpus > 1 and os.environ.has_key('SAGE_PARALLEL_DIST'):
+            extensions_to_compile = []
+            for ext in self.extensions:
+                need_to_compile, p = self.prepare_extension(ext)
+                if need_to_compile:
+                    extensions_to_compile.append(p)
+
+            if extensions_to_compile:
+               from processing import Pool
+               p = Pool(min(ncpus, len(extensions_to_compile)))
+               for r in p.imap(self.build_extension, extensions_to_compile):
+                   pass
+            
+        else:
+            for ext in self.extensions:
+                need_to_compile, p = self.prepare_extension(ext)
+                if need_to_compile:
+                    self.build_extension(p)
+
+        print "Total time spent compiling C/C++ extensions: ", time.time() - t, "seconds."
 
-    def build_extension(self, ext):
+    def prepare_extension(self, ext):
         sources = ext.sources
         if sources is None or type(sources) not in (ListType, TupleType):
             raise DistutilsSetupError, \
@@ -443,9 +487,16 @@
         depends = sources + ext.depends
         if not (self.force or newer_group(depends, ext_filename, 'newer')):
             log.debug("skipping '%s' extension (up-to-date)", ext.name)
-            return
+            need_to_compile = False
         else:
             log.info("building '%s' extension", ext.name)
+            need_to_compile = True
+
+        return need_to_compile, (sources, ext, ext_filename)
+
+    def build_extension(self, p):
+
+        sources, ext, ext_filename = p
 
         # First, scan the sources for SWIG definition files (.i), run
         # SWIG on 'em to create .c files, and modify the sources list
@@ -715,3 +766,23 @@
                 return ext.libraries
 
 # class build_ext
+
+
+
+def number_of_cpus():
+    """
+    Try to determine the number of cpu's on this system.
+    If successful return that number.  Otherwise return 0
+    to indicate failure.
+
+    OUTPUT:
+        int
+    """
+    if hasattr(os, "sysconf") and os.sysconf_names.has_key("SC_NPROCESSORS_ONLN"): # Linux and Unix
+        n = os.sysconf("SC_NPROCESSORS_ONLN")
+        if isinstance(n, int) and n > 0:
+            return n
+    try:
+        return int(os.popen2("sysctl -n hw.ncpu")[1].read().strip())
+    except:
+        return 0

[williamstein]

comment:2

POSITIVE TESTING: I just want to comment that I used this a lot all day today on both OS X and Linux, and didn't have a single problem. I like it!

[ncalexan]

comment:3

This has been in my tree for more than a week on a Mac OS X 10.5 machine. I have been doing some fairly heavy development, although most has not been in Cython. No problems to report; I say merge.

sage: version()
'Sage Version 3.2.1.alpha2, Release Date: 2008-11-26'

[sagetrac-mabshoff]

comment:4

There is a slight chicken and egg problem here since we build pyprocessing after python, so if someone sets the wrong environment variable the build of Python itself will blow up.

To fix this:

• make sure to unset SAGE_PARALLEL_DIST for the python.spkg as well as the pyprocessing.spkg
• make all the other python based spkgs depend on pyprocessing, so it is immediately build after Python

Then this should work out of the box, even though we should still test this with a couple rounds of builds of Sage with a high parallel level to make sure we do not hit any race conditions.

Cheers,

Michael

[sagetrac-mabshoff]

comment:5

After some chit chat with Craig we agreed that this should go into the enf of the pyprocessing spkg-install. That way everything is in place even for a fresh build.

Cheers,

Michael

[sagetrac-mabshoff]

comment:6

The spkg at

http://sage.math.washington.edu/home/mabshoff/release-cycles-3.2.2/alpha2/pyprocessing-0.52.p0.spkg

contains the fix. I tested it in various configs and it works for me :)

Cheers,

Michael