def test_constraintSub(self):
m = ri.create_model()
m.pert_a = Param(initialize=0.01)
m.pert_b = Param(initialize=1.01)
m_sipopt = sensitivity_calculation('sipopt', m, [m.a, m.b],
[m.pert_a, m.pert_b])
# verify substitutions in equality constraint
self.assertTrue(m_sipopt.C_equal.lower.ctype is Param
and m_sipopt.C_equal.upper.ctype is Param)
self.assertFalse(m_sipopt.C_equal.active)
self.assertTrue(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[3].lower == 0.0
and m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[3].upper == 0.0
and len(
list(
identify_variables(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[3].body)))
== 2)
# verify substitutions in one-sided bounded constraint
self.assertTrue(m_sipopt.C_singleBnd.lower is None
and m_sipopt.C_singleBnd.upper.ctype is Param)
self.assertFalse(m_sipopt.C_singleBnd.active)
self.assertTrue(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[4].lower is None
and m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[4].upper == 0.0
and len(
list(
identify_variables(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[4].body)))
== 2)
# verify substitutions in ranged inequality constraint
self.assertTrue(m_sipopt.C_rangedIn.lower.ctype is Param
and m_sipopt.C_rangedIn.upper.ctype is Param)
self.assertFalse(m_sipopt.C_rangedIn.active)
self.assertTrue(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[1].lower is None
and m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[1].upper == 0.0
and len(
list(
identify_variables(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[1].body)))
== 2)
self.assertTrue(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[2].lower is None
and m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[2].upper == 0.0
and len(
list(
identify_variables(
m_sipopt._SENSITIVITY_TOOLBOX_DATA.constList[2].body)))
== 2)
def test_constraintSub(self):
m = ri.create_model()
m.pert_a = Param(initialize=0.01)
m.pert_b = Param(initialize=1.01)
m_sipopt = sipopt(m,[m.a,m.b], [m.pert_a,m.pert_b])
#verify substitutions in equality constraint
self.assertTrue(m_sipopt.C_equal.lower.type() is Param and
m_sipopt.C_equal.upper.type() is Param)
self.assertFalse(m_sipopt.C_equal.active)
self.assertTrue(m_sipopt._sipopt_data.constList[3].lower == 0.0 and
m_sipopt._sipopt_data.constList[3].upper == 0.0 and
len(list(identify_variables(
m_sipopt._sipopt_data.constList[3].body))) == 2)
#verify substitutions in one-sided bounded constraint
self.assertTrue(m_sipopt.C_singleBnd.lower is None and
m_sipopt.C_singleBnd.upper.type() is Param)
self.assertFalse(m_sipopt.C_singleBnd.active)
self.assertTrue(m_sipopt._sipopt_data.constList[4].lower is None and
m_sipopt._sipopt_data.constList[4].upper == 0.0 and
len(list(identify_variables(
m_sipopt._sipopt_data.constList[4].body))) == 2)
#verify substitutions in ranged inequality constraint
self.assertTrue(m_sipopt.C_rangedIn.lower.type() is Param and
m_sipopt.C_rangedIn.upper.type() is Param)
self.assertFalse(m_sipopt.C_rangedIn.active)
self.assertTrue(m_sipopt._sipopt_data.constList[1].lower is None and
m_sipopt._sipopt_data.constList[1].upper == 0.0 and
len(list(identify_variables(
m_sipopt._sipopt_data.constList[1].body))) == 2)
self.assertTrue(m_sipopt._sipopt_data.constList[2].lower is None and
m_sipopt._sipopt_data.constList[2].upper == 0.0 and
len(list(identify_variables(
m_sipopt._sipopt_data.constList[2].body))) == 2)