Current gas flow implementation using first order Taylor approximation and big M could be improved. (Related to #1001)