def solve(casefile):
""" given a matpower casefile, solves the power flow using the Jabr method
and returns a dictionary mapping bus number to
(voltage magnitude, voltage angle) tuples. angles are in radians
"""
case = load_case(casefile)
u, R, I = build_gurobi_model(case)
V, theta = recover_original_variables(u, I)
i2e = case.i2e
answer = {i2e[i]: (v, t) for (i, (v, t)) in enumerate(zip(V, theta))}
return answer
def solve(casefile):
""" given a matpower casefile, solves the power flow using the Jabr method
and returns a dictionary mapping bus number to
(voltage magnitude, voltage angle) tuples. angles are in radians
"""
case = load_case(casefile)
u, R, I = build_gurobi_model(case)
V, theta = recover_original_variables(u, I)
i2e = case.i2e
answer = {i2e[i]: (v, t) for (i, (v, t)) in enumerate(zip(V, theta))}
return answer
def case118():
c = CASE_DIRECTORY + 'case118_v2.m'
return load_case(c)
def case14():
c = CASE_DIRECTORY + 'case14_tree.m'
return load_case(c)
def case5q():
c = CASE_DIRECTORY + 'case5_with_q.m'
return load_case(c)
def case5():
c = CASE_DIRECTORY + 'case5_renumber_tree.m'
return load_case(c)
