def METHOD_NAME():
    test_optimise_storage_size()
    energysystem = solph.EnergySystem()
    energysystem.restore()
    # Results
    results = energysystem.results["main"]
    meta = energysystem.results["meta"]
    electricity_bus = views.node(results, "electricity")
    my_results = electricity_bus["sequences"].sum(axis=0).to_dict()
    storage = energysystem.groups["storage"]
    my_results["storage_invest"] = results[(storage, None)]["scalars"][
        "invest"
    ]
    stor_invest_dict = {
        "storage_invest": 2040000,
        (("electricity", "None"), "duals"): 10800000000321,
        (("electricity", "demand"), "flow"): 105867395,
        (("electricity", "excess_bel"), "flow"): 211771291,
        (("electricity", "storage"), "flow"): 2350931,
        (("pp_gas", "electricity"), "flow"): 5148414,
        (("pv", "electricity"), "flow"): 7488607,
        (("storage", "electricity"), "flow"): 1880745,
        (("wind", "electricity"), "flow"): 305471851,
    }
    for key in stor_invest_dict.keys():
        assert int(round(my_results[key])) == int(round(stor_invest_dict[key]))
    # Solver results
    assert str(meta["solver"]["Termination condition"]) == "optimal"
    assert meta["solver"]["Error rc"] == 0
    assert str(meta["solver"]["Status"]) == "ok"
    # Problem results
    assert meta["problem"]["Lower bound"] == 4.231675777e17
    assert meta["problem"]["Upper bound"], 4.231675777e17
    assert meta["problem"]["Number of variables"] == 2807
    assert meta["problem"]["Number of constraints"] == 2808
    assert meta["problem"]["Number of nonzeros"] == 1197
    assert meta["problem"]["Number of objectives"] == 1
    assert str(meta["problem"]["Sense"]) == "minimize"
    # Objective function
    assert round(meta["objective"]) == 423167578261115584