dallan-keylogic

It isn't clear to me that there's anything there to "fix". The Jacobian is singular to machine precision. Something is clearly wrong with the example numerically. I can go through...

Okay, if we tighten `dt.config.parallel_component_tolerance` from `1e-8` to `1e-15`, we get the expeced output: ``` dt.config.parallel_component_tolerance = 1e-15 dt.display_near_parallel_constraints() ``` ``` ==================================================================================== The following pairs of constraints are nearly parallel:...

> Maybe we should use sqrt(unorm*vnorm - uvnorm) as the measure? Is the goal here to make the measure of colinearity extensive? Then we run into problems if we have...

One of the tests we're now failing is this one: ``` @pytest.mark.component def test_display_near_parallel_constraints(self): model = ConcreteModel() model.v1 = Var(initialize=1e-8) model.v2 = Var() model.v3 = Var() model.c1 = Constraint(expr=model.v1 ==...

@michaelbynum , they're getting used as part of building up an expression for excess molar enthalpy in an activity coefficient model. Initially, I was using AD to build up the...

This is in the context of building up a big Pyomo expression, so it has to be symbolic.

Somewhat tangential to this issue, I have another entry in the "zero should be treated as being equivalent with any set of units" file: ``` pyunits.convert( N_A * (q_e)**2/(8 *...