Added code for local/global event chain computation.  Currently not working beyond my base test case.
the chain computation is correct; reformatting and storing into an LP matrix for summation remains