dgels_ instead of dgelsd_, -Q (alt. epsilon, value is uniformly distributed across all accessible positions rather than added to each accessible position), additional error handling and better inference of BP distance when -X -N are provided (handles case when (0, d_0) XOR (d_0, 0) are not present).