Adding the modification to Druinsky-Toledo matrices alluded to in the corresponding paper and adding python wrappers for control solvers