Speedup

Some core cobrapy functions are too slow, or scale poorly with large model size. @jlerman44 pointed out some performance issues with addiing and removing reactions in large, so those functions are the first targets for optimization.