Skip to content

Commit

Permalink
Ajout du isclose et modification des tests pour les electrolyzer n_in…
Browse files Browse the repository at this point in the history 
…puts
  • Loading branch information
@Yann-Temudjin
Yann-Temudjin committed Mar 13, 2024
1 parent 18a4b1f commit d0e8e26
Show file tree
Hide file tree
Showing 2 changed files with 119 additions and 51 deletions.
148 changes: 113 additions & 35 deletions ...s/andromede/test_electrolyzer_n_inputs.py → ...s/unittests/test_electrolyzer_n_inputs.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -10,14 +10,16 @@
#
# This file is part of the Antares project.

import math

from andromede.libs.standard import DEMAND_MODEL, GENERATOR_MODEL, NODE_BALANCE_MODEL
from andromede.libs.standard_sc import (
CONVERTOR_MODEL,
CONVERTOR_MODEL_MOD,
DECOMPOSE_1_FLOW_INTO_2_FLOW,
TWO_INPUTS_CONVERTOR_MODEL,
)
from andromede.simulation import TimeBlock, build_problem
from andromede.simulation import OutputValues, TimeBlock, build_problem
from andromede.study import (
ConstantData,
DataBase,
Expand All @@ -27,10 +29,32 @@
create_component,
)

"""
for every following test we have two electrical production with an electrolyzer converting to a gaz flow
we always have:
first electric production:
- p_max = 70
- cost = 10
second electric production:
- p_max = 80
- cost = 20
gaz production:
- p_max = 30
- cost = 40
first production conversion:
- alpha = 0.7
second production conversion:
- alpha = 0.5
for a gaz demand of 100
"""


def test_electrolyzer_n_inputs_1():
"""
Test with an electrolyzer for each inputs
flow_ep1 * alpha_ez1 + flow_ep2 * alpha_ez2 + flow_gp
"""
elec_node_1 = Node(model=NODE_BALANCE_MODEL, id="e1")
electric_prod_1 = create_component(model=GENERATOR_MODEL, id="ep1")
Expand All @@ -46,17 +70,17 @@ def test_electrolyzer_n_inputs_1():

database = DataBase()

database.add_data("ep1", "p_max", ConstantData(100))
database.add_data("ep1", "cost", ConstantData(30))
database.add_data("ep1", "p_max", ConstantData(70))
database.add_data("ep1", "cost", ConstantData(10))
database.add_data("ez1", "alpha", ConstantData(0.7))

database.add_data("ep2", "p_max", ConstantData(100))
database.add_data("ep2", "cost", ConstantData(30))
database.add_data("ez2", "alpha", ConstantData(0.7))
database.add_data("ep2", "p_max", ConstantData(80))
database.add_data("ep2", "cost", ConstantData(20))
database.add_data("ez2", "alpha", ConstantData(0.5))

database.add_data("gd", "demand", ConstantData(70))
database.add_data("gp", "p_max", ConstantData(10))
database.add_data("gp", "cost", ConstantData(40))
database.add_data("gd", "demand", ConstantData(100))
database.add_data("gp", "p_max", ConstantData(30))
database.add_data("gp", "cost", ConstantData(15))

network = Network("test")
network.add_node(elec_node_1)
Expand Down Expand Up @@ -94,13 +118,27 @@ def test_electrolyzer_n_inputs_1():
problem = build_problem(network, database, TimeBlock(1, [0]), scenarios)
status = problem.solver.Solve()

output = OutputValues(problem)
ep1_gen = output.component("ep1").var("generation").value
ep2_gen = output.component("ep2").var("generation").value
gp_gen = output.component("gp").var("generation").value
print(ep1_gen)
print(ep2_gen)
print(gp_gen)

assert math.isclose(ep1_gen, 70)
assert math.isclose(ep2_gen, 42)
assert math.isclose(gp_gen, 30)

assert status == problem.solver.OPTIMAL
assert problem.solver.Objective().Value() == 2971.4285714285716
assert math.isclose(problem.solver.Objective().Value(), 1990)


def test_electrolyzer_n_inputs_2():
"""
Test with one electrolyzer that has two inputs
flow_ep1 * alpha1_ez + flow_ep2 * alpha2_ez + flow_gp
"""

elec_node_1 = Node(model=NODE_BALANCE_MODEL, id="e1")
Expand All @@ -116,15 +154,19 @@ def test_electrolyzer_n_inputs_2():
electrolyzer = create_component(model=TWO_INPUTS_CONVERTOR_MODEL, id="ez")

database = DataBase()
database.add_data("gd", "demand", ConstantData(70))
database.add_data("ep1", "p_max", ConstantData(100))
database.add_data("ep2", "p_max", ConstantData(100))
database.add_data("ep1", "cost", ConstantData(30))
database.add_data("ep2", "cost", ConstantData(30))

database.add_data("ez", "alpha1", ConstantData(0.7))
database.add_data("ez", "alpha2", ConstantData(0.5))
database.add_data("gp", "p_max", ConstantData(10))
database.add_data("gp", "cost", ConstantData(40))

database.add_data("ep1", "p_max", ConstantData(70))
database.add_data("ep1", "cost", ConstantData(10))

database.add_data("ep2", "p_max", ConstantData(80))
database.add_data("ep2", "cost", ConstantData(20))

database.add_data("gd", "demand", ConstantData(100))
database.add_data("gp", "p_max", ConstantData(30))
database.add_data("gp", "cost", ConstantData(15))

network = Network("test")
network.add_node(elec_node_1)
Expand Down Expand Up @@ -160,13 +202,29 @@ def test_electrolyzer_n_inputs_2():
problem = build_problem(network, database, TimeBlock(1, [0]), scenarios)
status = problem.solver.Solve()

output = OutputValues(problem)
ep1_gen = output.component("ep1").var("generation").value
ep2_gen = output.component("ep2").var("generation").value
gp_gen = output.component("gp").var("generation").value
print(ep1_gen)
print(ep2_gen)
print(gp_gen)

assert math.isclose(ep1_gen, 70)
assert math.isclose(ep2_gen, 42)
assert math.isclose(gp_gen, 30)

assert status == problem.solver.OPTIMAL
assert problem.solver.Objective().Value() == 2971.4285714285716
assert math.isclose(problem.solver.Objective().Value(), 1990)


def test_electrolyzer_n_inputs_3():
"""
Test with a consumption_electrolyzer with two inputs
(flow_ep1 + flow_ep2) * alpha_ez + flow_gp
The result is different since we only have one alpha at 0.7
"""
elec_node_1 = Node(model=NODE_BALANCE_MODEL, id="e1")
elec_node_2 = Node(model=NODE_BALANCE_MODEL, id="e2")
Expand All @@ -184,14 +242,18 @@ def test_electrolyzer_n_inputs_3():
)

database = DataBase()
database.add_data("gd", "demand", ConstantData(70))
database.add_data("ep1", "p_max", ConstantData(100))
database.add_data("ep2", "p_max", ConstantData(100))
database.add_data("ep1", "cost", ConstantData(30))
database.add_data("ep2", "cost", ConstantData(30))

database.add_data("ez", "alpha", ConstantData(0.7))
database.add_data("gp", "p_max", ConstantData(10))
database.add_data("gp", "cost", ConstantData(40))

database.add_data("ep1", "p_max", ConstantData(70))
database.add_data("ep1", "cost", ConstantData(10))

database.add_data("ep2", "p_max", ConstantData(80))
database.add_data("ep2", "cost", ConstantData(20))

database.add_data("gd", "demand", ConstantData(100))
database.add_data("gp", "p_max", ConstantData(30))
database.add_data("gp", "cost", ConstantData(15))

network = Network("test")
network.add_node(elec_node_1)
Expand Down Expand Up @@ -233,13 +295,26 @@ def test_electrolyzer_n_inputs_3():
problem = build_problem(network, database, TimeBlock(1, [0]), scenarios)
status = problem.solver.Solve()

output = OutputValues(problem)
ep1_gen = output.component("ep1").var("generation").value
ep2_gen = output.component("ep2").var("generation").value
gp_gen = output.component("gp").var("generation").value

assert math.isclose(ep1_gen, 70)
assert math.isclose(ep2_gen, 30)
assert math.isclose(gp_gen, 30)

assert status == problem.solver.OPTIMAL
assert problem.solver.Objective().Value() == 2971.4285714285716
assert math.isclose(problem.solver.Objective().Value(), 1750)


def test_electrolyzer_n_inputs_4():
"""
Test with one electrolyzer with one input that takes every inputs
flow_ep1 + flow_ep2 * alpha_ez + flow_gp
with the same values this one is infeasible since the electrolyzer model as in the drawio doesn't work with two inputs
"""
elec_node_1 = Node(model=NODE_BALANCE_MODEL, id="e1")
elec_node_2 = Node(model=NODE_BALANCE_MODEL, id="e2")
Expand All @@ -254,14 +329,18 @@ def test_electrolyzer_n_inputs_4():
electrolyzer = create_component(model=CONVERTOR_MODEL, id="ez")

database = DataBase()
database.add_data("gd", "demand", ConstantData(70))
database.add_data("ep1", "p_max", ConstantData(100))
database.add_data("ep2", "p_max", ConstantData(100))
database.add_data("ep1", "cost", ConstantData(30))
database.add_data("ep2", "cost", ConstantData(30))

database.add_data("ez", "alpha", ConstantData(0.7))
database.add_data("gp", "p_max", ConstantData(10))
database.add_data("gp", "cost", ConstantData(40))

database.add_data("ep1", "p_max", ConstantData(70))
database.add_data("ep1", "cost", ConstantData(10))

database.add_data("ep2", "p_max", ConstantData(80))
database.add_data("ep2", "cost", ConstantData(20))

database.add_data("gd", "demand", ConstantData(100))
database.add_data("gp", "p_max", ConstantData(30))
database.add_data("gp", "cost", ConstantData(15))

network = Network("test")
network.add_node(elec_node_1)
Expand Down Expand Up @@ -297,5 +376,4 @@ def test_electrolyzer_n_inputs_4():
problem = build_problem(network, database, TimeBlock(1, [0]), scenarios)
status = problem.solver.Solve()

assert status == problem.solver.OPTIMAL
assert problem.solver.Objective().Value() == 2971.4285714285716
assert status == problem.solver.INFEASIBLE
22 changes: 6 additions & 16 deletions tests/unittests/test_quota_co2.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -10,20 +10,10 @@
#
# This file is part of the Antares project.

from andromede.expression import literal, param, var
from andromede.expression.expression import port_field
import math

from andromede.libs.standard import DEMAND_MODEL, LINK_MODEL, NODE_BALANCE_MODEL
from andromede.libs.standard_sc import C02_POWER_MODEL, QUOTA_CO2_MODEL
from andromede.model import (
Constraint,
ModelPort,
PortField,
PortType,
float_parameter,
float_variable,
model,
)
from andromede.model.model import PortFieldDefinition, PortFieldId
from andromede.simulation import OutputValues, TimeBlock, build_problem
from andromede.study import (
ConstantData,
Expand Down Expand Up @@ -98,7 +88,7 @@ def test_quota_co2():
l12_flow = output.component("L12").var("flow").value

assert status == problem.solver.OPTIMAL
assert problem.solver.Objective().Value() == 5500
assert oil1_p == 50
assert coal1_p == 50
assert l12_flow == -50
assert math.isclose(problem.solver.Objective().Value(), 5500)
assert math.isclose(oil1_p, 50)
assert math.isclose(coal1_p, 50)
assert math.isclose(l12_flow, -50)

0 comments on commit d0e8e26

Please sign in to comment.