Bug fixes and efficiency gains

1. Better exception handling in equation_solve(). Removes trivially erroneous, all-zero columns.
2. Efficiency gains in indicator and multiplier creation.
3. New function to create backward linkages.