OOM for large order and small cluster size. 

Can we avoid allocating interpolation points when they will fail the size check?  Or is that too hairy?  Maybe we add a compatibility check of the input parameters, and increase the cluster size within the code if the user uses too small a value.