For the biggest potential improvement, see this issue: #23 all pairs matrix construction functions. matrix free functions for low memory usage settings. block-wise functions that are especially ...