Source code for presentation made to SIAM CSE 2013.