pyomo__pyomo-3017

{
  "status": "error",
  "error_type": "CommandError",
  "error_message": "Traceback (most recent call last):\n  File \"<stdin>\", line 6, in <module>\n  File \"/workspace/pyomo/core/base/units_container.py\", line 1532, in __getattribute__\n    return getattr(pint_module, attr)\n  File \"/workspace/pyomo/common/dependencies.py\", line 75, in __getattr__\n    raise DeferredImportError(self._moduleunavailable_message())\npyomo.common.errors.DeferredImportError: The \"pint\" package failed to import. This package is necessary to use Pyomo units. (import raised ModuleNotFoundError: No module named 'pint')\n",
  "details": {
    "stdout": "",
    "stderr": "Traceback (most recent call last):\n  File \"<stdin>\", line 6, in <module>\n  File \"/workspace/pyomo/core/base/units_container.py\", line 1532, in __getattribute__\n    return getattr(pint_module, attr)\n  File \"/workspace/pyomo/common/dependencies.py\", line 75, in __getattr__\n    raise DeferredImportError(self._moduleunavailable_message())\npyomo.common.errors.DeferredImportError: The \"pint\" package failed to import. This package is necessary to use Pyomo units. (import raised ModuleNotFoundError: No module named 'pint')\n",
    "exit_code": 1
  }
}

{
  "status": "ok",
  "stdout": "",
  "stderr": "Traceback (most recent call last):\n  File \"/workspace/pyomo/common/collections/component_map.py\", line 71, in __getitem__\n    return self._dict[id(obj)][1]\nKeyError: 129418845037200\n\nDuring handling of the above exception, another exception occurred:\n\nTraceback (most recent call last):\n  File \"<stdin>\", line 17, in <module>\n  File \"/workspace/pyomo/core/expr/calculus/derivatives.py\", line 96, in differentiate\n    res = reverse_sd(expr=expr)\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 484, in reverse_sd\n    return _reverse_diff_helper(expr, False)\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 442, in _reverse_diff_helper\n    _diff_GeneralExpression(e, val_dict, der_dict)\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 331, in _diff_GeneralExpression\n    der_dict[node.expr] += der_dict[node]\n  File \"/workspace/pyomo/common/collections/component_map.py\", line 73, in __getitem__\n    raise KeyError(\"Component with id '%s': %s\" % (id(obj), str(obj)))\nKeyError: \"Component with id '129418845037200': 2.0\"\n",
  "exit_code": 0
}

{
  "status": "ok",
  "filepath": "pyomo/core/expr/calculus/diff_with_pyomo.py",
  "content": "    val_dict: ComponentMap\n    der_dict: ComponentMap\n    \"\"\"\n    if node.getname() in _unary_map:\n        _unary_map[node.getname()](node, val_dict, der_dict)\n    else:\n        raise DifferentiationException(\n            'Unsupported expression type for differentiation: {0}'.format(type(node))\n        )\n\n\ndef _diff_GeneralExpression(node, val_dict, der_dict):\n    \"\"\"\n    Reverse automatic differentiation for named expressions.\n\n    Parameters\n    ----------\n    node: The named expression\n    val_dict: ComponentMap\n    der_dict: ComponentMap\n    \"\"\"\n    der_dict[node.expr] += der_dict[node]\n\n\ndef _diff_ExternalFunctionExpression(node, val_dict, der_dict):\n    \"\"\"\n\n    Parameters\n    ----------\n    node: pyomo.core.expr.numeric_expr.ExternalFunctionExpression\n    val_dict: ComponentMap\n    der_dict: ComponentMap\n    \"\"\"\n    der = der_dict[node]\n    vals = tuple(val_dict[i] for i in node.args)\n    derivs = node._fcn.evaluate_fgh(vals, fgh=1)[1]\n    for ndx, arg in enumerate(node.args):\n        der_dict[arg] += der * derivs[ndx]\n\n\n_diff_map = dict()",
  "start_line": 310,
  "end_line": 350,
  "total_lines": 484,
  "is_truncated": false
}

{
  "status": "ok",
  "filepath": "pyomo/common/collections/component_map.py",
  "content": "    def __str__(self):\n        \"\"\"String representation of the mapping.\"\"\"\n        tmp = {str(c) + \" (id=\" + str(id(c)) + \")\": v for c, v in self.items()}\n        return \"ComponentMap(\" + str(tmp) + \")\"\n\n    #\n    # Implement MutableMapping abstract methods\n    #\n\n    def __getitem__(self, obj):\n        try:\n            return self._dict[id(obj)][1]\n        except KeyError:\n            raise KeyError(\"Component with id '%s': %s\" % (id(obj), str(obj)))\n\n    def __setitem__(self, obj, val):\n        self._dict[id(obj)] = (obj, val)\n\n    def __delitem__(self, obj):\n        try:\n            del self._dict[id(obj)]",
  "start_line": 60,
  "end_line": 80,
  "total_lines": 135,
  "is_truncated": false
}

{
  "status": "ok",
  "stdout": "node: product, node.expr: x*2, type(node.expr): <class 'pyomo.core.expr.numeric_expr.ProductExpression'>\nnode: y, node.expr: 2.0, type(node.expr): <class 'pyomo.core.expr.numvalue.NumericConstant'>\n",
  "stderr": "Traceback (most recent call last):\n  File \"/workspace/pyomo/common/collections/component_map.py\", line 71, in __getitem__\n    return self._dict[id(obj)][1]\nKeyError: 128142104586096\n\nDuring handling of the above exception, another exception occurred:\n\nTraceback (most recent call last):\n  File \"<stdin>\", line 25, in <module>\n  File \"/workspace/pyomo/core/expr/calculus/derivatives.py\", line 96, in differentiate\n    res = reverse_sd(expr=expr)\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 484, in reverse_sd\n    return _reverse_diff_helper(expr, False)\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 442, in _reverse_diff_helper\n    _diff_GeneralExpression(e, val_dict, der_dict)\n  File \"<stdin>\", line 20, in new_diff_general_expression\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 331, in _diff_GeneralExpression\n    der_dict[node.expr] += der_dict[node]\n  File \"/workspace/pyomo/common/collections/component_map.py\", line 73, in __getitem__\n    raise KeyError(\"Component with id '%s': %s\" % (id(obj), str(obj)))\nKeyError: \"Component with id '128142104586096': 2.0\"\n",
  "exit_code": 0
}

{
  "status": "ok",
  "filepath": "pyomo/core/expr/numvalue.py",
  "content": "#  ___________________________________________________________________________\n#\n#  Pyomo: Python Optimization Modeling Objects\n#  Copyright (c) 2008-2022\n#  National Technology and Engineering Solutions of Sandia, LLC\n#  Under the terms of Contract DE-NA0003525 with National Technology and\n#  Engineering Solutions of Sandia, LLC, the U.S. Government retains certain\n#  rights in this software.\n#  This software is distributed under the 3-clause BSD License.\n#  ___________________________________________________________________________\n\n__all__ = (\n    'value',\n    'is_constant',\n    'is_fixed',\n    'is_variable_type',\n    'is_potentially_variable',\n    'NumericValue',\n    'ZeroConstant',\n    'native_numeric_types',\n    'native_types',\n    'nonpyomo_leaf_types',\n    'polynomial_degree',\n)\n\nimport collections\nimport sys\nimport logging\n\nfrom pyomo.common.deprecation import (\n    deprecated,\n    deprecation_warning,\n    relocated_module_attribute,\n)\nfrom pyomo.core.expr.expr_common import ExpressionType\nfrom pyomo.core.expr.numeric_expr import NumericValue\nimport pyomo.common.numeric_types as _numeric_types\n\n# TODO: update Pyomo to import these objects from common.numeric_types\n#   (and not from here)\nfrom pyomo.common.numeric_types import (\n    nonpyomo_leaf_types,\n    native_types,\n    native_numeric_types,\n    native_integer_types,\n    native_logical_types,\n    pyomo_constant_types,\n    check_if_numeric_type,\n    value,\n)\nfrom pyomo.core.pyomoobject import PyomoObject\n\nrelocated_module_attribute(\n    'native_boolean_types',\n    'pyomo.common.numeric_types._native_boolean_types',\n    version='6.6.0',\n    f_globals=globals(),\n    msg=\"The native_boolean_types set will be removed in the future: the set \"\n    \"contains types that were convertible to bool, and not types that should \"\n    \"be treated as if they were bool (as was the case for the other \"\n    \"native_*_types sets).  Users likely should use native_logical_types.\",\n)\nrelocated_module_attribute(\n    'RegisterNumericType',\n    'pyomo.common.numeric_types.RegisterNumericType',\n    version='6.6.0',\n    f_globals=globals(),\n)\nrelocated_module_attribute(\n    'RegisterIntegerType',\n    'pyomo.common.numeric_types.RegisterIntegerType',\n    version='6.6.0',\n    f_globals=globals(),\n)\nrelocated_module_attribute(\n    'RegisterBooleanType',\n    'pyomo.common.numeric_types.RegisterBooleanType',\n    version='6.6.0',\n    f_globals=globals(),\n)\nrelocated_module_attribute(\n    'NumericValue',\n    'pyomo.core.expr.numeric_expr.NumericValue',\n    version='6.6.2',\n    f_globals=globals(),\n)\nrelocated_module_attribute(\n    'NumericNDArray',\n    'pyomo.core.expr.numeric_expr.NumericNDArray',\n    version='6.6.2',\n    f_globals=globals(),\n)\n\nlogger = logging.getLogger('pyomo.core')\n\n\n##------------------------------------------------------------------------\n##\n## Standard types of expressions\n##\n##------------------------------------------------------------------------\n\n\nclass NonNumericValue(object):\n    \"\"\"An object that contains a non-numeric value\n\n    Constructor Arguments:\n        value           The initial value.\n    \"\"\"\n\n    __slots__ = ('value',)\n\n    def __init__(self, value):\n        self.value = value\n\n    def __str__(self):\n        return str(self.value)\n\n\nnonpyomo_leaf_types.add(NonNumericValue)\n\n\ndef is_constant(obj):\n    \"\"\"\n    A utility function that returns a boolean that indicates\n    whether the object is a constant.\n    \"\"\"\n    # JDS: NB: I am not sure why we allow str to be a constant, but\n    # since we have historically done so, we check for type membership\n    # in native_types and not in native_numeric_types.\n    #\n    if obj.__class__ in native_types:\n        return True\n    try:\n        return obj.is_constant()\n    except AttributeError:\n        pass\n    # Now we need to confirm that we have an unknown numeric type\n    #\n    # As this branch is only hit for previously unknown (to Pyomo)\n    # types that behave reasonably like numbers, we know they *must*\n    # be constant.\n    if check_if_numeric_type(obj):\n        return True\n    else:\n        raise TypeError(\n            \"Cannot assess properties of object with unknown type: %s\"\n            % (type(obj).__name__,)\n        )\n",
  "start_line": 1,
  "end_line": 150,
  "total_lines": 433,
  "is_truncated": true
}

{
  "status": "ok",
  "filepath": "pyomo/core/expr/calculus/diff_with_pyomo.py",
  "content": "#  ___________________________________________________________________________\n#\n#  Pyomo: Python Optimization Modeling Objects\n#  Copyright (c) 2008-2022\n#  National Technology and Engineering Solutions of Sandia, LLC\n#  Under the terms of Contract DE-NA0003525 with National Technology and\n#  Engineering Solutions of Sandia, LLC, the U.S. Government retains certain\n#  rights in this software.\n#  This software is distributed under the 3-clause BSD License.\n#  ___________________________________________________________________________\n\nfrom pyomo.common.collections import ComponentMap, ComponentSet\nimport pyomo.core.expr as _expr\nfrom pyomo.core.expr.visitor import ExpressionValueVisitor, nonpyomo_leaf_types\nfrom pyomo.core.expr.numvalue import value, is_constant\nfrom pyomo.core.expr import exp, log, sin, cos\nimport math\n\n\n\"\"\"\nThe purpose of this file is to perform symbolic differentiation and \nfirst order automatic differentiation directly with pyomo \nexpressions. This is certainly not as efficient as doing AD in C or \nC++, but it avoids the translation from pyomo expressions to a form \nwhere AD can be performed efficiently. The only functions that are \nmeant to be used by users are reverse_ad and reverse_sd. First, \nvalues are propagated from the leaves to each node in the tree with \nthe LeafToRoot visitors. Then derivative values are propagated from \nthe root to the leaves with the RootToLeaf visitors.\n\"\"\"\n\n\nclass DifferentiationException(Exception):\n    pass\n\n\ndef _diff_ProductExpression(node, val_dict, der_dict):\n    \"\"\"\n\n    Parameters\n    ----------\n    node: pyomo.core.expr.numeric_expr.ProductExpression\n    val_dict: ComponentMap\n    der_dict: ComponentMap\n    \"\"\"\n    assert len(node.args) == 2\n    arg1, arg2 = node.args\n    der = der_dict[node]\n    der_dict[arg1] += der * val_dict[arg2]\n    der_dict[arg2] += der * val_dict[arg1]",
  "start_line": 1,
  "end_line": 50,
  "total_lines": 484,
  "is_truncated": false
}

{
  "status": "ok",
  "stdout": "Differentiation successful!\ny : Size=1, Index=None\n    Key  : Expression\n    None :        2.0\n",
  "stderr": "",
  "exit_code": 0
}

{
  "status": "ok",
  "filepath": "pyomo/core/expr/calculus/diff_with_pyomo.py",
  "content": "_diff_map = dict()\n_diff_map[_expr.ProductExpression] = _diff_ProductExpression\n_diff_map[_expr.DivisionExpression] = _diff_DivisionExpression\n_diff_map[_expr.PowExpression] = _diff_PowExpression\n_diff_map[_expr.SumExpression] = _diff_SumExpression\n_diff_map[_expr.MonomialTermExpression] = _diff_ProductExpression\n_diff_map[_expr.NegationExpression] = _diff_NegationExpression\n_diff_map[_expr.UnaryFunctionExpression] = _diff_UnaryFunctionExpression\n_diff_map[_expr.ExternalFunctionExpression] = _diff_ExternalFunctionExpression\n_diff_map[_expr.LinearExpression] = _diff_SumExpression\n_diff_map[_expr.AbsExpression] = _diff_abs\n\n_diff_map[_expr.NPV_ProductExpression] = _diff_ProductExpression\n_diff_map[_expr.NPV_DivisionExpression] = _diff_DivisionExpression\n_diff_map[_expr.NPV_PowExpression] = _diff_PowExpression\n_diff_map[_expr.NPV_SumExpression] = _diff_SumExpression\n_diff_map[_expr.NPV_NegationExpression] = _diff_NegationExpression\n_diff_map[_expr.NPV_UnaryFunctionExpression] = _diff_UnaryFunctionExpression\n_diff_map[_expr.NPV_ExternalFunctionExpression] = _diff_ExternalFunctionExpression\n_diff_map[_expr.NPV_AbsExpression] = _diff_abs\n\n\ndef _symbolic_value(x):\n    return x\n\n\ndef _numeric_apply_operation(node, values):\n    return node._apply_operation(values)\n\n\ndef _symbolic_apply_operation(node, values):\n    return node\n\n\nclass _LeafToRootVisitor(ExpressionValueVisitor):\n    def __init__(self, val_dict, der_dict, expr_list, numeric=True):\n        \"\"\"\n        Parameters\n        ----------\n        val_dict: ComponentMap\n        der_dict: ComponentMap\n        \"\"\"\n        self.val_dict = val_dict\n        self.der_dict = der_dict\n        self.expr_list = expr_list\n        assert len(self.expr_list) == 0\n        assert len(self.val_dict) == 0\n        assert len(self.der_dict) == 0\n        if numeric:\n            self.value_func = value\n            self.operation_func = _numeric_apply_operation\n        else:\n            self.value_func = _symbolic_value\n            self.operation_func = _symbolic_apply_operation\n\n    def visit(self, node, values):\n        self.val_dict[node] = self.operation_func(node, values)\n        self.der_dict[node] = 0\n        self.expr_list.append(node)\n        return self.val_dict[node]\n\n    def visiting_potential_leaf(self, node):\n        if node in self.val_dict:\n            return True, self.val_dict[node]\n\n        if node.__class__ in nonpyomo_leaf_types:\n            self.val_dict[node] = node\n            self.der_dict[node] = 0\n            return True, node\n\n        if not node.is_expression_type():\n            val = self.value_func(node)\n            self.val_dict[node] = val\n            self.der_dict[node] = 0\n            return True, val\n\n        return False, None\n\n\ndef _reverse_diff_helper(expr, numeric=True):\n    val_dict = ComponentMap()\n    der_dict = ComponentMap()\n    expr_list = list()\n\n    visitorA = _LeafToRootVisitor(val_dict, der_dict, expr_list, numeric=numeric)\n    visitorA.dfs_postorder_stack(expr)\n\n    der_dict[expr] = 1\n    for e in reversed(expr_list):\n        if e.__class__ in _diff_map:\n            _diff_map[e.__class__](e, val_dict, der_dict)\n        elif e.is_named_expression_type():\n            _diff_GeneralExpression(e, val_dict, der_dict)\n        else:\n            raise DifferentiationException(\n                'Unsupported expression type for differentiation: {0}'.format(type(e))\n            )\n\n    return der_dict\n\n",
  "start_line": 350,
  "end_line": 450,
  "total_lines": 484,
  "is_truncated": false
}

{
  "status": "ok",
  "stdout": "der_dict after visitor: dict_keys([130929167799280, 130929185784080, 130929167866416, 130929167822528, 130929167875136])\nProcessing e: product, type: <class 'pyomo.core.base.expression.ScalarExpression'>\n  Calling _diff_GeneralExpression for product, e.expr: x*2\nProcessing e: x*2, type: <class 'pyomo.core.expr.numeric_expr.ProductExpression'>\nProcessing e: y, type: <class 'pyomo.core.base.expression.ScalarExpression'>\n  Calling _diff_GeneralExpression for y, e.expr: 2.0\n",
  "stderr": "Traceback (most recent call last):\n  File \"/workspace/pyomo/common/collections/component_map.py\", line 71, in __getitem__\n    return self._dict[id(obj)][1]\nKeyError: 130929167822768\n\nDuring handling of the above exception, another exception occurred:\n\nTraceback (most recent call last):\n  File \"<stdin>\", line 45, in <module>\n  File \"/workspace/pyomo/core/expr/calculus/derivatives.py\", line 96, in differentiate\n    res = reverse_sd(expr=expr)\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 484, in reverse_sd\n    return _reverse_diff_helper(expr, False)\n  File \"<stdin>\", line 34, in patched_reverse_diff_helper\n  File \"/workspace/pyomo/core/expr/calculus/diff_with_pyomo.py\", line 331, in _diff_GeneralExpression\n    der_dict[node.expr] += der_dict[node]\n  File \"/workspace/pyomo/common/collections/component_map.py\", line 73, in __getitem__\n    raise KeyError(\"Component with id '%s': %s\" % (id(obj), str(obj)))\nKeyError: \"Component with id '130929167822768': 2.0\"\n",
  "exit_code": 0
}

{
  "status": "ok",
  "filepath": "pyomo/core/expr/calculus/diff_with_pyomo.py",
  "occurrences": 1,
  "strategy": "unescaped",
  "diff": "--- a/pyomo/core/expr/calculus/diff_with_pyomo.py\n+++ b/pyomo/core/expr/calculus/diff_with_pyomo.py\n@@ -328,7 +328,8 @@\n     val_dict: ComponentMap\n     der_dict: ComponentMap\n     \"\"\"\n-    der_dict[node.expr] += der_dict[node]\n+    if not is_constant(node.expr):\n+        der_dict[node.expr] += der_dict[node]\n \n \n def _diff_ExternalFunctionExpression(node, val_dict, der_dict):\n",
  "is_truncated": false
}

{
  "status": "ok",
  "stdout": "Differentiation successful!\ny : Size=1, Index=None\n    Key  : Expression\n    None :        2.0\n",
  "stderr": "",
  "exit_code": 0
}