Various small fixes (#1443)

* auto-rebuild if settings indicates old gpkit version
* fix ohm printing for kiloohms, etc.
* fix constraint sensitivities for ccp_localsolve

gpkit/globals.py

@@ -26,7 +26,9 @@ def load_settings(path=None, firstattempt=True):
                         value.remove("mosek")
     except IOError:
         settings_ = {"installed_solvers": [""]}
-    if settings_["installed_solvers"] == [""]:
+    if (settings_["installed_solvers"] == [""]
+            or ("mosek" in settings_["installed_solvers"]
+                and "mosek_version" not in settings_)):
         if firstattempt:
             print("Found no installed solvers, beginning a build.")
             build()

gpkit/nomials/math.py

@@ -375,6 +375,7 @@ def mono_approximation(self, x0):
 class ScalarSingleEquationConstraint(SingleEquationConstraint):
     "A SingleEquationConstraint with scalar left and right sides."
     nomials = []
+    sgp_parent = None
 
     def __init__(self, left, oper, right):
         lr = [left, right]
@@ -394,14 +395,18 @@ def __init__(self, left, oper, right):
     def relaxed(self, relaxvar):
         "Returns the relaxation of the constraint in a list."
         if self.oper == ">=":
-            return [relaxvar*self.left >= self.right]
+            relaxed = [relaxvar*self.left >= self.right]
         elif self.oper == "<=":
-            return [self.left <= relaxvar*self.right]
+            relaxed = [self.left <= relaxvar*self.right]
         elif self.oper == "=":
-            return [self.left <= relaxvar*self.right,
-                    relaxvar*self.left >= self.right]
-        raise ValueError(
-            "Constraint %s had unknown operator %s." % self.oper, self)
+            relaxed = [self.left <= relaxvar*self.right,
+                       relaxvar*self.left >= self.right]
+        else:
+            raise ValueError(
+                "Constraint %s had unknown operator %s." % self.oper, self)
+        for constr in relaxed:
+            constr.sgp_parent = self
+        return relaxed
 
 
 # pylint: disable=too-many-instance-attributes, invalid-unary-operand-type
@@ -412,7 +417,6 @@ class PosynomialInequality(ScalarSingleEquationConstraint):
     """
 
     feastol = 1e-3
-    sgp_parent = None
     relax_sensitivity = None
     # NOTE: follows .check_result's max default, but 1e-3 seems a bit lax...
 
@@ -519,6 +523,8 @@ def sens_from_dual(self, la, nu, result):  # pylint: disable=unused-argument
         self.relax_sensitivity = la
         if self.sgp_parent:
             self.sgp_parent.relax_sensitivity = la
+            if getattr(self.sgp_parent, "sgp_parent", None):
+                self.sgp_parent.sgp_parent.relax_sensitivity = la
         nu, = nu
         presub, = self.unsubbed
         if hasattr(self, "pmap"):
@@ -600,6 +606,9 @@ def sens_from_dual(self, la, nu, result):
         self.relax_sensitivity = la[0] - la[1]
         if self.sgp_parent:
             self.sgp_parent.relax_sensitivity = self.relax_sensitivity
+            if getattr(self.sgp_parent, "sgp_parent", None):
+                self.sgp_parent.sgp_parent.relax_sensitivity = \
+                    self.relax_sensitivity
         var_senss = {}
         for var in self.varkeys:
             for i, m in enumerate(self.unsubbed):

gpkit/repr_conventions.py

@@ -34,7 +34,7 @@ def unitstr(units, into="%s", options=UNIT_FORMATTING, dimless=""):
     if not isinstance(units, Quantity):
         return dimless
     if options == ":~" and "ohm" in str(units.units):
-        rawstr = "ohm"  # otherwise it'll be a capital Omega
+        rawstr = str(units.units)  # otherwise it'll be a capital Omega
     else:
         rawstr = ("{%s}" % options).format(units.units)
     units = rawstr.replace(" ", "").replace("dimensionless", dimless)

gpkit/solution_array.py

@@ -70,9 +70,14 @@ def tight_table(self, _, ntightconstrs=5, tight_senss=1e-2, **kwargs):
     if not self.model:
         return []
     title = "Tightest Constraints"
-    data = [((-float("%+6.2g" % c.relax_sensitivity), str(c)),
-             "%+6.2g" % c.relax_sensitivity, id(c), c)
-            for c in self.model.flat() if c.relax_sensitivity >= tight_senss]
+    try:
+        data = [((-float("%+6.2g" % c.relax_sensitivity), str(c)),
+                 "%+6.2g" % c.relax_sensitivity, id(c), c)
+                for c in self.model.flat()
+                if c.relax_sensitivity >= tight_senss]
+    except AttributeError:
+        print("Constraint %s had no `relax_sensitivity` attribute." % c)
+        return []
     if not data:
         lines = ["No constraints had a sensitivity above %+5.1g."
                  % tight_senss]