def _convert_to_const(val):
if val.__class__ in native_numeric_types:
return val
elif val.is_expression_type():
return EXPR.evaluate_expression(val)
else:
return value(val)
def evalHessian(hessian):
n = len(hessian)
H = np.zeros((n, n))
for i, row in enumerate(hessian):
for j, expression in enumerate(row):
H[i][j] = evaluate_expression(expression)
return H
def _processUnhashableIndex(self, idx):
"""Process a call to __getitem__ with unhashable elements
There are three basic ways to get here:
1) the index contains one or more slices or ellipsis
2) the index contains an unhashable type (e.g., a Pyomo
(Simple)Component
3) the index contains an IndexTemplate
"""
from pyomo.core.expr import current as EXPR
#
# Iterate through the index and look for slices and constant
# components
#
fixed = {}
sliced = {}
ellipsis = None
_found_numeric = False
#
# Setup the slice template (in fixed)
#
if type(idx) is tuple:
# We would normally do "flatten()" here, but the current
# (10/2016) implementation of flatten() is too aggressive:
# it will attempt to expand *any* iterable, including
# SimpleParam.
idx = pyutilib.misc.flatten_tuple(idx)
elif type(idx) is list:
idx = pyutilib.misc.flatten_tuple(tuple(idx))
else:
idx = (idx,)
for i,val in enumerate(idx):
if type(val) is slice:
if val.start is not None or val.stop is not None:
raise IndexError(
"Indexed components can only be indexed with simple "
"slices: start and stop values are not allowed.")
if val.step is not None:
logger.warning(
"DEPRECATION WARNING: The special wildcard slice "
"(::0) is deprecated. Please use an ellipsis (...) "
"to indicate '0 or more' indices")
val = Ellipsis
else:
if ellipsis is None:
sliced[i] = val
else:
sliced[i-len(idx)] = val
continue
if val is Ellipsis:
if ellipsis is not None:
raise IndexError(
"Indexed components can only be indexed with simple "
"slices: the Pyomo wildcard slice (Ellipsis; "
"e.g., '...') can only appear once")
ellipsis = i
continue
if hasattr(val, 'is_expression_type'):
_num_val = val
# Attempt to retrieve the numeric value .. if this
# is a template expression generation, then it
# should raise a TemplateExpressionError
try:
val = EXPR.evaluate_expression(val, constant=True)
_found_numeric = True
except TemplateExpressionError:
#
# The index is a template expression, so return the
# templatized expression.
#
from pyomo.core.expr import current as EXPR
return EXPR.GetItemExpression(tuple(idx), self)
except EXPR.NonConstantExpressionError:
#
# The expression contains an unfixed variable
#
raise RuntimeError(
"""Error retrieving the value of an indexed item %s:
index %s is not a constant value. This is likely not what you meant to
do, as if you later change the fixed value of the object this lookup
will not change. If you understand the implications of using
non-constant values, you can get the current value of the object using
the value() function.""" % ( self.name, i ))
except EXPR.FixedExpressionError:
#
# The expression contains a fixed variable
#
raise RuntimeError(
"""Error retrieving the value of an indexed item %s:
index %s is a fixed but not constant value. This is likely not what you
meant to do, as if you later change the fixed value of the object this
lookup will not change. If you understand the implications of using
fixed but not constant values, you can get the current value using the
value() function.""" % ( self.name, i ))
#
# There are other ways we could get an exception such as
# evaluating a Param / Var that is not initialized.
# These exceptions will continue up the call stack.
#
# verify that the value is hashable
hash(val)
if ellipsis is None:
fixed[i] = val
else:
fixed[i - len(idx)] = val
if sliced or ellipsis is not None:
return _IndexedComponent_slice(self, fixed, sliced, ellipsis)
elif _found_numeric:
if len(idx) == 1:
return fixed[0]
else:
return tuple( fixed[i] for i in range(len(idx)) )
else:
raise DeveloperError(
"Unknown problem encountered when trying to retrieve "
"index for component %s" % (self.name,) )
def _processUnhashableIndex(self, idx, _exception=None):
"""Process a call to __getitem__ with unhashable elements
There are three basic ways to get here:
1) the index contains one or more slices or ellipsis
2) the index contains an unhashable type (e.g., a Pyomo
(Simple)Component
3) the index contains an IndexTemplate
"""
from pyomo.core.expr import current as EXPR
#
# Iterate through the index and look for slices and constant
# components
#
fixed = {}
sliced = {}
ellipsis = None
_found_numeric = False
#
# Setup the slice template (in fixed)
#
if type(idx) is tuple:
# We would normally do "flatten()" here, but the current
# (10/2016) implementation of flatten() is too aggressive:
# it will attempt to expand *any* iterable, including
# SimpleParam.
idx = pyutilib.misc.flatten_tuple(idx)
elif type(idx) is list:
idx = pyutilib.misc.flatten_tuple(tuple(idx))
else:
idx = (idx, )
for i, val in enumerate(idx):
if type(val) is slice:
if val.start is not None or val.stop is not None:
raise IndexError(
"Indexed components can only be indexed with simple "
"slices: start and stop values are not allowed.")
if val.step is not None:
logger.warning(
"DEPRECATION WARNING: The special wildcard slice "
"(::0) is deprecated. Please use an ellipsis (...) "
"to indicate '0 or more' indices")
val = Ellipsis
else:
if ellipsis is None:
sliced[i] = val
else:
sliced[i - len(idx)] = val
continue
if val is Ellipsis:
if ellipsis is not None:
raise IndexError(
"Indexed components can only be indexed with simple "
"slices: the Pyomo wildcard slice (Ellipsis; "
"e.g., '...') can only appear once")
ellipsis = i
continue
if hasattr(val, 'as_numeric'):
try:
val = EXPR.evaluate_expression(val, constant=True)
_found_numeric = True
except TemplateExpressionError:
#
# The index is a template expression, so return the
# templatized expression.
#
from pyomo.core.expr import current as EXPR
return EXPR.GetItemExpression(tuple(idx), self)
except EXPR.NonConstantExpressionError:
#
# The expression contains an unfixed variable
#
raise RuntimeError(
"""Error retrieving the value of an indexed item %s:
index %s is not a constant value. This is likely not what you meant to
do, as if you later change the fixed value of the object this lookup
will not change. If you understand the implications of using
non-constant values, you can get the current value of the object using
the value() function.""" % (self.name, i))
except EXPR.FixedExpressionError:
#
# The expression contains a fixed variable
#
raise RuntimeError(
"""Error retrieving the value of an indexed item %s:
index %s is a fixed but not constant value. This is likely not what you
meant to do, as if you later change the fixed value of the object this
lookup will not change. If you understand the implications of using
fixed but not constant values, you can get the current value using the
value() function.""" % (self.name, i))
#
# There are other ways we could get an exception such as
# evaluating a Param / Var that is not initialized.
# These exceptions will continue up the call stack.
#
# verify that the value is hashable
hash(val)
if ellipsis is None:
fixed[i] = val
else:
fixed[i - len(idx)] = val
if sliced or ellipsis is not None:
return _IndexedComponent_slicer(self, fixed, sliced, ellipsis)
elif _found_numeric:
if len(idx) == 1:
return fixed[0]
else:
return tuple(fixed[i] for i in range(len(idx)))
elif _exception is not None:
raise
else:
raise DeveloperError(
"Unknown problem encountered when trying to retrieve "
"index for component %s" % (self.name, ))
def _processUnhashableIndex(self, idx):
"""Process a call to __getitem__ with unhashable elements
There are three basic ways to get here:
1) the index contains one or more slices or ellipsis
2) the index contains an unhashable type (e.g., a Pyomo
(Simple)Component
3) the index contains an IndexTemplate
"""
from pyomo.core.expr import current as EXPR
#
# Iterate through the index and look for slices and constant
# components
#
fixed = {}
sliced = {}
ellipsis = None
_found_numeric = False
#
# Setup the slice template (in fixed)
#
if normalize_index.flatten:
idx = normalize_index(idx)
if idx.__class__ is not tuple:
idx = (idx, )
for i, val in enumerate(idx):
if type(val) is slice:
if val.start is not None or val.stop is not None:
raise IndexError(
"Indexed components can only be indexed with simple "
"slices: start and stop values are not allowed.")
if val.step is not None:
logger.warning(
"DEPRECATION WARNING: The special wildcard slice "
"(::0) is deprecated. Please use an ellipsis (...) "
"to indicate '0 or more' indices")
val = Ellipsis
else:
if ellipsis is None:
sliced[i] = val
else:
sliced[i - len(idx)] = val
continue
if val is Ellipsis:
if ellipsis is not None:
raise IndexError(
"Indexed components can only be indexed with simple "
"slices: the Pyomo wildcard slice (Ellipsis; "
"e.g., '...') can only appear once")
ellipsis = i
continue
if hasattr(val, 'is_expression_type'):
_num_val = val
# Attempt to retrieve the numeric value .. if this
# is a template expression generation, then it
# should raise a TemplateExpressionError
try:
val = EXPR.evaluate_expression(val, constant=True)
_found_numeric = True
except TemplateExpressionError:
#
# The index is a template expression, so return the
# templatized expression.
#
from pyomo.core.expr import current as EXPR
return EXPR.GetItemExpression((self, ) + tuple(idx))
except EXPR.NonConstantExpressionError:
#
# The expression contains an unfixed variable
#
raise RuntimeError(
"""Error retrieving the value of an indexed item %s:
index %s is not a constant value. This is likely not what you meant to
do, as if you later change the fixed value of the object this lookup
will not change. If you understand the implications of using
non-constant values, you can get the current value of the object using
the value() function.""" % (self.name, i))
except EXPR.FixedExpressionError:
#
# The expression contains a fixed variable
#
raise RuntimeError(
"""Error retrieving the value of an indexed item %s:
index %s is a fixed but not constant value. This is likely not what you
meant to do, as if you later change the fixed value of the object this
lookup will not change. If you understand the implications of using
fixed but not constant values, you can get the current value using the
value() function.""" % (self.name, i))
#
# There are other ways we could get an exception such as
# evaluating a Param / Var that is not initialized.
# These exceptions will continue up the call stack.
#
# verify that the value is hashable
hash(val)
if ellipsis is None:
fixed[i] = val
else:
fixed[i - len(idx)] = val
if sliced or ellipsis is not None:
slice_dim = len(idx)
if ellipsis is not None:
slice_dim -= 1
if normalize_index.flatten:
set_dim = self.dim()
elif self._implicit_subsets is None:
# Scalar component.
set_dim = 0
else:
set_dim = len(self._implicit_subsets)
structurally_valid = False
if slice_dim == set_dim or set_dim is None:
structurally_valid = True
elif ellipsis is not None and slice_dim < set_dim:
structurally_valid = True
elif set_dim == 0 and idx == (slice(None), ):
# If dim == 0 and idx is slice(None), the component was
# a scalar passed a single slice. Since scalar components
# can be accessed with a "1-dimensional" index of None,
# this behavior is allowed.
#
# Note that x[...] is caught above, as slice_dim will be
# 0 in that case
structurally_valid = True
if not structurally_valid:
raise IndexError(
"Index %s contains an invalid number of entries for "
"component %s. Expected %s, got %s." %
(idx, self.name, set_dim, slice_dim))
return IndexedComponent_slice(self, fixed, sliced, ellipsis)
elif _found_numeric:
if len(idx) == 1:
return fixed[0]
else:
return tuple(fixed[i] for i in range(len(idx)))
else:
raise DeveloperError(
"Unknown problem encountered when trying to retrieve "
"index for component %s" % (self.name, ))
def get_objective(instance, name="objective"):
return evaluate_expression(instance.objective)
