def validate_nonanticipativity(options, scenario_tree, scenario_solutions,
node_solutions):
# Check that variables with non-anticipativity constraints have
# converged
for node_name, node in scenario_tree['nodes'].items():
assert node_name == node['name']
node_sol = node_solutions[node_name]
for scenario_name in node['scenarios']:
scenario_sol = scenario_solutions[scenario_name]
# variables, objective, stage costs
for variable_name, variable_node_sol in node_sol[
'variables'].items():
variable_scenario_sol = scenario_sol['variables'][
variable_name]
try:
assert_float_equals(variable_node_sol['solution'],
variable_scenario_sol['value'],
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating nonanticipativity within "
"node %s for variable %s in scenario %s. "
"Variable value in scenario solution does not "
"match node solution.\n" %
(node_name, variable_name, scenario_name))
msg += ("Node value: %r\n" %
(variable_node_sol['solution']))
msg += ("Scenario value: %r\n" %
(variable_scenario_sol['value']))
_update_exception_message(msg, e)
raise
try:
assert_value_equals(variable_node_sol['fixed'],
variable_scenario_sol['fixed'])
except AssertionError as e:
msg = ""
msg += (
"Problem validating nonanticipativity within "
"node %s for variable %s in scenario %s. "
"Variable fixed status in scenario solution does not "
"match node fixed status.\n" %
(node_name, variable_name, scenario_name))
msg += ("Node fixed status: %s\n" %
(variable_node_sol['fixed']))
msg += ("Scenario fixed status: %s\n" %
(variable_scenario_sol['fixed']))
_update_exception_message(msg, e)
raise
def validate_nonanticipativity(options,
scenario_tree,
scenario_solutions,
node_solutions):
# Check that variables with non-anticipativity constraints have
# converged
for node_name, node in scenario_tree['nodes'].items():
assert node_name == node['name']
node_sol = node_solutions[node_name]
for scenario_name in node['scenarios']:
scenario_sol = scenario_solutions[scenario_name]
# variables, objective, stage costs
for variable_name, variable_node_sol in node_sol['variables'].items():
variable_scenario_sol = scenario_sol['variables'][variable_name]
try:
assert_float_equals(variable_node_sol['solution'],
variable_scenario_sol['value'],
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating nonanticipativity within "
"node %s for variable %s in scenario %s. "
"Variable value in scenario solution does not "
"match node solution.\n"
% (node_name, variable_name, scenario_name))
msg += ("Node value: %r\n"
% (variable_node_sol['solution']))
msg += ("Scenario value: %r\n"
% (variable_scenario_sol['value']))
_update_exception_message(msg, e)
raise
try:
assert_value_equals(variable_node_sol['fixed'],
variable_scenario_sol['fixed'])
except AssertionError as e:
msg = ""
msg += ("Problem validating nonanticipativity within "
"node %s for variable %s in scenario %s. "
"Variable fixed status in scenario solution does not "
"match node fixed status.\n"
% (node_name, variable_name, scenario_name))
msg += ("Node fixed status: %s\n"
% (variable_node_sol['fixed']))
msg += ("Scenario fixed status: %s\n"
% (variable_scenario_sol['fixed']))
_update_exception_message(msg, e)
raise
def validate_expected_node_costs(options, scenario_tree, scenario_solutions,
node_solutions):
# Check that reported expected node costs make sense
# and that the expected node cost of the root node
# equals the ef objective function
# sorted by time
sorted_stages = \
sorted(scenario_tree['stages'].values(),key=lambda x: x['order'])
# leaf to root
reverse_sorted_stages = reversed(sorted_stages)
node_expected_costs = \
dict((node_name,0.0) for node_name in scenario_tree['nodes'])
for stage in reverse_sorted_stages:
assert stage is scenario_tree['stages'][stage['name']]
stage_name = stage['name']
for node_name in stage['nodes']:
node = scenario_tree['nodes'][node_name]
assert node['name'] == node_name
node_scenario_names = node['scenarios']
single_scenario_name = node_scenario_names[0]
single_scenario_stage_cost = \
scenario_solutions[single_scenario_name]\
['stage costs'][stage_name]
node_cost = 0.0
for scenario_name in node_scenario_names:
scenario = scenario_tree['scenarios'][scenario_name]
scenario_stage_cost = \
scenario_solutions[scenario_name]['stage costs'][stage_name]
node_cost += \
scenario['probability'] * \
scenario_solutions[scenario_name]['stage costs'][stage_name]
try:
assert_float_equals(single_scenario_stage_cost,
scenario_stage_cost,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating stage costs "
"for stage %s of scenarios within node %s. "
"Not all scenarios report the same stage cost.\n" %
(stage_name, node_name))
msg += ("Scenario %s stage cost: %r\n" %
(scenario_name, scenario_stage_cost))
msg += ("Scenario %s stage cost: %r\n" %
(single_scenario_name, single_scenario_stage_cost))
_update_exception_message(msg, e)
raise
node_cost /= node['probability']
node_expected_costs[node_name] += node_cost
for child_node_name in node['children']:
child_node = scenario_tree['nodes'][child_node_name]
node_expected_costs[node_name] += \
child_node['conditional probability'] * \
node_expected_costs[child_node_name]
del child_node
del node
del node_scenario_names
del node_cost
del single_scenario_stage_cost
del stage_name
for node_name, expected_cost in node_expected_costs.items():
node_sol = node_solutions[node_name]
try:
assert_float_equals(expected_cost, node_sol['expected cost'],
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating node expected cost for node %s\n" %
(node_name))
msg += ("Node expected cost in solution: %r\n" %
(node_sol['expected cost']))
msg += ("Computed average from stage costs: %r\n" %
(expected_cost))
_update_exception_message(msg, e)
raise
del node_sol
# get the root node
assert len(scenario_tree['stages'][sorted_stages[0]['name']]['nodes']) == 1
root_node_name = scenario_tree['stages'][sorted_stages[0]
['name']]['nodes'][0]
root_node_expected_cost = node_expected_costs[root_node_name]
expected_scenario_objective = 0.0
expected_scenario_cost1 = 0.0
expected_scenario_cost2 = 0.0
for scenario_name, scenario in scenario_tree['scenarios'].items():
scenario_sol = scenario_solutions[scenario_name]
expected_scenario_objective += \
scenario['probability']*scenario_sol['objective']
expected_scenario_cost1 += \
scenario['probability']*scenario_sol['cost']
expected_scenario_cost2 += \
scenario['probability'] * \
sum(stage_cost for stage_cost in scenario_sol['stage costs'].values())
try:
assert_float_equals(root_node_expected_cost, expected_scenario_cost1,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating root node expected cost\n")
msg += ("Root node expected cost in solution: %r\n" %
(root_node_expected_cost))
msg += ("Computed average from scenario costs in solution: %r\n" %
(expected_scenario_cost1))
_update_exception_message(msg, e)
raise
try:
assert_float_equals(root_node_expected_cost, expected_scenario_cost2,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating root node expected cost\n")
msg += ("Root node expected cost in solution: %r\n" %
(root_node_expected_cost))
msg += (
"Computed average from scenario stage costs in solution: %r\n" %
(expected_scenario_cost2))
_update_exception_message(msg, e)
raise
try:
assert_float_equals(root_node_expected_cost,
expected_scenario_objective,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating root node expected cost\n")
msg += ("Root node expected cost in solution: %r\n" %
(root_node_expected_cost))
msg += ("Computed average from scenario objectives in solution: %r\n" %
(expected_scenario_objective))
_update_exception_message(msg, e)
raise
def validate_expected_node_costs(options,
scenario_tree,
scenario_solutions,
node_solutions):
# Check that reported expected node costs make sense
# and that the expected node cost of the root node
# equals the ef objective function
# sorted by time
sorted_stages = \
sorted(scenario_tree['stages'].values(),key=lambda x: x['order'])
# leaf to root
reverse_sorted_stages = reversed(sorted_stages)
node_expected_costs = \
dict((node_name,0.0) for node_name in scenario_tree['nodes'])
for stage in reverse_sorted_stages:
assert stage is scenario_tree['stages'][stage['name']]
stage_name = stage['name']
for node_name in stage['nodes']:
node = scenario_tree['nodes'][node_name]
assert node['name'] == node_name
node_scenario_names = node['scenarios']
single_scenario_name = node_scenario_names[0]
single_scenario_stage_cost = \
scenario_solutions[single_scenario_name]\
['stage costs'][stage_name]
node_cost = 0.0
for scenario_name in node_scenario_names:
scenario = scenario_tree['scenarios'][scenario_name]
scenario_stage_cost = \
scenario_solutions[scenario_name]['stage costs'][stage_name]
node_cost += \
scenario['probability'] * \
scenario_solutions[scenario_name]['stage costs'][stage_name]
try:
assert_float_equals(single_scenario_stage_cost,
scenario_stage_cost,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating stage costs "
"for stage %s of scenarios within node %s. "
"Not all scenarios report the same stage cost.\n"
% (stage_name, node_name))
msg += ("Scenario %s stage cost: %r\n"
% (scenario_name, scenario_stage_cost))
msg += ("Scenario %s stage cost: %r\n"
% (single_scenario_name, single_scenario_stage_cost))
_update_exception_message(msg, e)
raise
node_cost /= node['probability']
node_expected_costs[node_name] += node_cost
for child_node_name in node['children']:
child_node = scenario_tree['nodes'][child_node_name]
node_expected_costs[node_name] += \
child_node['conditional probability'] * \
node_expected_costs[child_node_name]
del child_node
del node
del node_scenario_names
del node_cost
del single_scenario_stage_cost
del stage_name
for node_name, expected_cost in node_expected_costs.items():
node_sol = node_solutions[node_name]
try:
assert_float_equals(expected_cost,
node_sol['expected cost'],
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating node expected cost for node %s\n"
% (node_name))
msg += ("Node expected cost in solution: %r\n"
% (node_sol['expected cost']))
msg += ("Computed average from stage costs: %r\n"
% (expected_cost))
_update_exception_message(msg, e)
raise
del node_sol
# get the root node
assert len(scenario_tree['stages'][sorted_stages[0]['name']]['nodes']) == 1
root_node_name = scenario_tree['stages'][sorted_stages[0]['name']]['nodes'][0]
root_node_expected_cost = node_expected_costs[root_node_name]
expected_scenario_objective = 0.0
expected_scenario_cost1 = 0.0
expected_scenario_cost2 = 0.0
for scenario_name, scenario in scenario_tree['scenarios'].items():
scenario_sol = scenario_solutions[scenario_name]
expected_scenario_objective += \
scenario['probability']*scenario_sol['objective']
expected_scenario_cost1 += \
scenario['probability']*scenario_sol['cost']
expected_scenario_cost2 += \
scenario['probability'] * \
sum(stage_cost for stage_cost in scenario_sol['stage costs'].values())
try:
assert_float_equals(root_node_expected_cost,
expected_scenario_cost1,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating root node expected cost\n")
msg += ("Root node expected cost in solution: %r\n"
% (root_node_expected_cost))
msg += ("Computed average from scenario costs in solution: %r\n"
% (expected_scenario_cost1))
_update_exception_message(msg, e)
raise
try:
assert_float_equals(root_node_expected_cost,
expected_scenario_cost2,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating root node expected cost\n")
msg += ("Root node expected cost in solution: %r\n"
% (root_node_expected_cost))
msg += ("Computed average from scenario stage costs in solution: %r\n"
% (expected_scenario_cost2))
_update_exception_message(msg, e)
raise
try:
assert_float_equals(root_node_expected_cost,
expected_scenario_objective,
options.absolute_tolerance)
except AssertionError as e:
msg = ""
msg += ("Problem validating root node expected cost\n")
msg += ("Root node expected cost in solution: %r\n"
% (root_node_expected_cost))
msg += ("Computed average from scenario objectives in solution: %r\n"
% (expected_scenario_objective))
_update_exception_message(msg, e)
raise