Circularities dev branch

Cargo.toml

@@ -78,7 +78,6 @@ missing_panics_doc = "allow"
 implicit_hasher = "allow"
 cast_lossless = "allow"
 return_self_not_must_use = "allow"
-from_iter_instead_of_collect = "allow"
 
 [[bench]]
 name = "examples"

benches/assets.rs

@@ -95,13 +95,15 @@ fn calculate_seed_prices(
     candidates: &[AssetRef],
     writer: &mut DataWriter,
 ) -> Prices {
-    let solution_existing = DispatchRun::new(model, base_year_assets, BASE_YEAR)
+    let market_demands = collect_preset_demands_for_year(&model.commodities, BASE_YEAR);
+    let solution_existing = DispatchRun::new(model, base_year_assets, BASE_YEAR, &market_demands)
         .run("bench setup: without candidates", writer)
         .expect("Dispatch without candidates failed");
-    let solution_with_candidates = DispatchRun::new(model, base_year_assets, BASE_YEAR)
-        .with_candidates(candidates)
-        .run("bench setup: with candidates", writer)
-        .expect("Dispatch with candidates failed");
+    let solution_with_candidates =
+        DispatchRun::new(model, base_year_assets, BASE_YEAR, &market_demands)
+            .with_candidates(candidates)
+            .run("bench setup: with candidates", writer)
+            .expect("Dispatch with candidates failed");
 
     calculate_prices(
         model,

examples/circularity/processes.csv

@@ -1,19 +1,19 @@
-id,description,regions,primary_output,start_year,end_year,capacity_to_activity
-GASDRV,Dry gas extraction,GBR,GASPRD,2020,,1.0
-OAGRSV,Oil and associated gas extraction,GBR,OILCRD,2020,,1.0
-GASPRC,Gas processing,GBR,GASNAT,2020,,1.0
-OILREF,Oil refinery,GBR,DIESEL,2020,,1.0
-OILRF2,O2G refinery,GBR,GASOLI,2020,,1.0
-WNDFRM,Wind farm,GBR,ELCTRI,2020,,31.54
-GASCGT,Gas combined cycle turbine,GBR,ELCTRI,2020,,31.54
-H2YGEN,H2 turbine,GBR,ELCTRI,2020,,31.54
-H2YPRO,H2 electrolyer,GBR,H2YPRD,2020,,1.0
-TPETCR,Petrol car,GBR,TPASKM,2020,,1.0
-TDIECR,Diesel car,GBR,TPASKM,2020,,1.0
-TELCCR,Electric car,GBR,TPASKM,2020,,1.0
-THYBCR,Plug-in hybrid car,GBR,TPASKM,2020,,1.0
-RGASBR,Gas boiler,GBR,RSHEAT,2020,,1.0
-RELCHP,Heat pump,GBR,RSHEAT,2020,,1.0
-BIOPRO,Biomass production,all,BIOPRD,2020,,1.0
-BIOPLL,Biomass pelletiser,all,BIOPEL,2020,,1.0
-RBIOBL,Biomass boiler,all,RSHEAT,2020,,1.0
+id,description,regions,primary_output,start_year,end_year,capacity_to_activity,feedback_process
+GASDRV,Dry gas extraction,GBR,GASPRD,2020,,1.0,
+OAGRSV,Oil and associated gas extraction,GBR,OILCRD,2020,,1.0,
+GASPRC,Gas processing,GBR,GASNAT,2020,,1.0,
+OILREF,Oil refinery,GBR,DIESEL,2020,,1.0,
+OILRF2,O2G refinery,GBR,GASOLI,2020,,1.0,
+WNDFRM,Wind farm,GBR,ELCTRI,2020,,31.54,
+GASCGT,Gas combined cycle turbine,GBR,ELCTRI,2020,,31.54,
+H2YGEN,H2 turbine,GBR,ELCTRI,2020,,31.54,true
+H2YPRO,H2 electrolyer,GBR,H2YPRD,2020,,1.0,true
+TPETCR,Petrol car,GBR,TPASKM,2020,,1.0,
+TDIECR,Diesel car,GBR,TPASKM,2020,,1.0,
+TELCCR,Electric car,GBR,TPASKM,2020,,1.0,
+THYBCR,Plug-in hybrid car,GBR,TPASKM,2020,,1.0,
+RGASBR,Gas boiler,GBR,RSHEAT,2020,,1.0,
+RELCHP,Heat pump,GBR,RSHEAT,2020,,1.0,
+BIOPRO,Biomass production,all,BIOPRD,2020,,1.0,
+BIOPLL,Biomass pelletiser,all,BIOPEL,2020,,1.0,
+RBIOBL,Biomass boiler,all,RSHEAT,2020,,1.0,

schemas/input/model.yaml

@@ -95,10 +95,6 @@ properties:
     type: integer
     description: Number of iterations to perform when calculating prices for cyclically-dependent markets
     default: 1
-  capacity_margin:
-    type: number
-    description: Slack proportion for assets selected during cycle balancing to absorb small demand shifts
-    default: 0.2
   mothball_years:
     type: integer
     default: 0

src/fixture.rs

@@ -325,6 +325,7 @@ pub fn process(
         capacity_to_activity: ActivityPerCapacity(1.0),
         investment_constraints: process_investment_constraints,
         tranche_size: None,
+        feedback_process: false,
     }
 }
 

src/graph.rs

@@ -14,6 +14,7 @@ use std::io::Write as IoWrite;
 use std::path::Path;
 use std::sync::Arc;
 
+pub mod feedback_suggest;
 pub mod investment;
 pub mod validate;
 
@@ -41,11 +42,21 @@ pub enum GraphNode {
 /// An edge in the commodity graph
 pub enum GraphEdge {
     /// An edge representing a primary flow of a process
-    #[display("{_0}")]
-    Primary(ProcessID),
+    #[display("{process_id}")]
+    Primary {
+        /// The process responsible for the flow.
+        process_id: ProcessID,
+        /// Whether the process participates in a feedback conversion.
+        feedback: bool,
+    },
     /// An edge representing a secondary (non-primary) flow of a process
-    #[display("{_0}")]
-    Secondary(ProcessID),
+    #[display("{process_id}")]
+    Secondary {
+        /// The process responsible for the flow.
+        process_id: ProcessID,
+        /// Whether the process participates in a feedback conversion.
+        feedback: bool,
+    },
     /// An edge representing a service demand
     #[display("DEMAND")]
     Demand,
@@ -150,9 +161,15 @@ fn create_commodities_graph_for_region_year(
                 source_node_index,
                 target_node_index,
                 if is_primary {
-                    GraphEdge::Primary(process.id.clone())
+                    GraphEdge::Primary {
+                        process_id: process.id.clone(),
+                        feedback: process.feedback_process,
+                    }
                 } else {
-                    GraphEdge::Secondary(process.id.clone())
+                    GraphEdge::Secondary {
+                        process_id: process.id.clone(),
+                        feedback: process.feedback_process,
+                    }
                 },
             );
         }
@@ -181,7 +198,7 @@ pub fn build_commodity_graphs_for_model(
 fn get_edge_attributes(_: &CommoditiesGraph, edge_ref: EdgeReference<GraphEdge>) -> String {
     match edge_ref.weight() {
         // Use dashed lines for secondary flows
-        GraphEdge::Secondary(_) => "style=dashed".to_string(),
+        GraphEdge::Secondary { .. } => "style=dashed".to_string(),
         // Other edges use default attributes
         _ => String::new(),
     }