parallel_dirty_map.py has option for greater threading.  Each thread acts on views of node's covariance chunk.