Update PTDF-based DC powerflow formulation to handle multiple islands

[DillonJ]

Currently, SpineOpt will fail due to an uninvertible admittance matrix if there are mutliple island. We already have code to determine the number of islands and their constituent nodes... we would just need to additonally dimension the ptdf and lodf arrays on network ID.

This would be a nice consolidation of existing functionality.