Added pair.code function. Given a pair object return the PC (number from 1 to 16)
Implementation of xCM
New helper functions. Change how we handle IR constraints in the LP model
Cleaned up rke2.R  -- Added initialization for all mechanisms to save time. Started implementing Bonus (incomplete)