Implemented working beta doubling and a variable, nnmax, to keep track of the max number of operators.
Problems: cutoff is routinely way to high >2 x cutoff, which really slows the sim down. Don’t  know how to fix it.