Constants about gpkit

[bqpd]

Ooh, neat. How about a new 'sub' function in Monomial and Posynomial?

(moving 'matrix of exponents associated with each constant in each constraint' to the GP class thread)

L , D, rho, C_L, C_D, S = gpkit.monify('L D rho C_L C_D S')
slf = [ L >= 0.5*rho*C_L*S*V**2,
        D >= 0.5*rho*C_D*S*V**2 ]

constants = {'rho': 1.2, 'V': 50}

# apply constants
cost = (0.5*rho*C_D*S*V**2).sub(constants)
constraints = [ p.sub(constants) for p in slf+other_constraints ]
sol = gpkit.minimize(cost, constraints)

[whoburg]

Finally getting around to seeing how this works -- this is pretty cool.

[bqpd]

Still a bit hacky, since it just replaces all of the constants each time, instead of updating them efficiently.