Build a sorted list of unique basis functions based on columns, where each basis function is a list
make_basis_list(X_sub, cols, order_map)
Note that sorting of columns is performed such that the basis order equals cols.length() and each basis function is a list(cols, cutoffs).