new version, see release notes

Author: mcosovic

Project.toml

@@ -1,7 +1,7 @@
 name = "JuliaGrid"
 uuid = "040ba96c-7d3d-4d4a-931c-37903f53a565"
 authors = ["Mirsad Cosovic <mirsad.cosovic@gmail.com>"]
-version = "0.5.0"
+version = "0.5.1"
 
 [deps]
 HDF5 = "f67ccb44-e63f-5c2f-98bd-6dc0ccc4ba2f"

docs/src/background/releasenotes.md

@@ -1,5 +1,18 @@
 # [Release Notes](@id ReleaseNotes)
 
+
+---
+
+## Version 0.5.1
+
+Release Date: July 12, 2025
+
+#### Added
+  * Introduced the `optimal` keyword in the function [`powerSystem`](@ref powerSystem) to skip importing data related to optimal power flow analyses, see [Partial Load for Faster Import](@ref PartialLoadFasterImportManual).
+
+#### Other
+  * Removed unused parameters from the `PowerSystem` type.
+
 ---
 
 ## Version 0.5.0

docs/src/manual/powerSystemModel.md

@@ -67,6 +67,15 @@ system = powerSystem("C:/hdf5/case14.h5")
 
 ---
 
+##### [Partial Load for Faster Import](@id PartialLoadFasterImportManual)
+A significant portion of data in input files may relate to optimal power flow analysis. If AC or DC optimal power flow is not needed, users can skip importing this additional data by setting the `optimal` keyword to `false`. For example:
+```julia
+system = powerSystem("C:/matpower/case14.m"; optimal = false)
+```
+This reduces memory usage and speeds up data loading, which is especially useful for large-scale systems.
+
+---
+
 ##### Model from Scratch
 Alternatively, the model can be built from scratch using built-in functions, for example:
 ```@example buildModelScratch

src/backend/internal.jl

@@ -56,7 +56,7 @@ Changing the unit of `active` power is reflected in the following quantities:
   [`@branch`](@ref @branch): if `type = 2`: `minFromBus`, `maxFromBus`, `minToBus`, `maxToBus`;
 * [`addGenerator!`](@ref addGenerator!), [`updateGenerator!`](@ref updateGenerator!),
   [`@generator`](@ref @generator): `active`, `minActive`, `maxActive`, `lowActive`,
-  `upActive`, `loadFollowing`, `reserve10min`, `reserve30min`;
+  `upActive`;
 * [`cost!`](@ref cost!): if `active`: `piecewise`, `polynomial`;
 * [`addWattmeter!`](@ref addWattmeter!), [`updateWattmeter!`](@ref updateWattmeter!):
   `active`, `variance`;
@@ -67,7 +67,7 @@ Changing the unit of `reactive` power unit is reflected in the following quantit
   `reactive`, `susceptance`;
 * [`addGenerator!`](@ref addGenerator!), [`updateGenerator!`](@ref updateGenerator!),
   [`@generator`](@ref @generator): `reactive`, `minReactive`, `maxReactive`,
-  `minLowReactive`, `maxLowReactive`, `minUpReactive`, `maxUpReactive`, `reactiveRamp`;
+  `minLowReactive`, `maxLowReactive`, `minUpReactive`, `maxUpReactive`;
 * [`cost!`](@ref cost!): if `reactive`: `piecewise`, `polynomial`;
 * [`addVarmeter!`](@ref addVarmeter!), [`updateVarmeter!`](@ref updateVarmeter!):
   `reactive`, `variance`;
@@ -482,17 +482,7 @@ macro default(mode::Symbol)
             template.generator.maxUpReactive.value = 0.0
             template.generator.maxUpReactive.pu = true
 
-            template.generator.loadFollowing.value = 0.0
-            template.generator.loadFollowing.pu = true
-            template.generator.reactiveRamp.value = 0.0
-            template.generator.reactiveRamp.pu = true
-            template.generator.reserve10min.value = 0.0
-            template.generator.reserve10min.pu = true
-            template.generator.reserve30min.value = 0.0
-            template.generator.reserve30min.pu = true
-
             template.generator.status = Int8(1)
-            template.generator.area = 0
             template.generator.label = "?"
             template.generator.key = String
         end

src/backend/precompile.jl

@@ -58,5 +58,8 @@ PrecompileTools.@setup_workload begin
         ########## Observability Analysis ###########
         islands = islandTopologicalFlow(mt)
         restorationGram!(mt, mt, islands)
+
+        ########## Default ###########
+        @default(template)
     end
 end

src/backend/utility.jl

@@ -762,6 +762,15 @@ function errorCovariance(pmu::PMU, varianceIm::Float64, L2::Float64, idx::Int64)
     end
 end
 
+function errorOptimal(system::PowerSystem)
+    if !system.bus.layout.optimal
+        throw(ErrorException(
+            "The power system object lacks optimal power flow data, likely due to " *
+            "optimal = false in the powerSystem function call.")
+        )
+    end
+end
+
 ##### Info Messages #####
 function infoObjective(label::IntStr)
     @info(

src/data/case14.h5

@@ -87,20 +87,20 @@
 
 %% bus names
 mpc.bus_name = {
-	'Bus 1     HV';
-	'Bus 2     HV';
-	'Bus 3     HV';
-	'Bus 4     HV';
-	'Bus 5     HV';
-	'Bus 6     LV';
-	'Bus 7     ZV';
-	'Bus 8     TV';
-	'Bus 9     LV';
-	'Bus 10    LV';
-	'Bus 11    LV';
-	'Bus 12    LV';
-	'Bus 13    LV';
-	'Bus 14    LV';
+	'Bus 1 HV';
+	'Bus 2 HV';
+	'Bus 3 HV';
+	'Bus 4 HV';
+	'Bus 5 HV';
+	'Bus 6 LV';
+	'Bus 7 ZV';
+	'Bus 8 TV';
+	'Bus 9 LV';
+	'Bus 10 LV';
+	'Bus 11 LV';
+	'Bus 12 LV';
+	'Bus 13 LV';
+	'Bus 14 LV';
 };
 
 % Warnings from cdf2matp conversion:

src/data/case14.m

@@ -165,12 +165,7 @@ Base.@kwdef mutable struct GeneratorTemplate
     upActive::ContainerTemplate = ContainerTemplate()
     minUpReactive::ContainerTemplate = ContainerTemplate()
     maxUpReactive::ContainerTemplate = ContainerTemplate()
-    loadFollowing::ContainerTemplate = ContainerTemplate()
-    reactiveRamp::ContainerTemplate = ContainerTemplate()
-    reserve10min::ContainerTemplate = ContainerTemplate()
-    reserve30min::ContainerTemplate = ContainerTemplate()
     status::Int8 = Int8(1)
-    area::Int64 = 0
     label::String = "?"
     key::DataType = String
 end
@@ -369,7 +364,6 @@ end
 
 Base.@kwdef struct GeneratorKey
     label::IntStrMiss = missing
-    area::IntMiss = missing
     status::IntMiss = missing
     active::FltIntMiss = missing
     reactive::FltIntMiss = missing
@@ -384,10 +378,6 @@ Base.@kwdef struct GeneratorKey
     upActive::FltIntMiss = missing
     minUpReactive::FltIntMiss = missing
     maxUpReactive::FltIntMiss = missing
-    loadFollowing::FltIntMiss = missing
-    reserve10min::FltIntMiss = missing
-    reserve30min::FltIntMiss = missing
-    reactiveRamp::FltIntMiss = missing
 end
 
 Base.@kwdef struct CostKey

src/definition/internal.jl

@@ -26,6 +26,7 @@ mutable struct BusLayout
     slack::Int64
     label::Int64
     pattern::Int64
+    const optimal::Bool
 end
 
 mutable struct BusSupply
@@ -130,13 +131,6 @@ mutable struct GeneratorCapability
     maxUpReactive::Vector{Float64}
 end
 
-mutable struct GeneratorRamping
-    loadFollowing::Vector{Float64}
-    reserve10min::Vector{Float64}
-    reserve30min::Vector{Float64}
-    reactiveRamp::Vector{Float64}
-end
-
 mutable struct Cost
     model::Vector{Int8}
     polynomial::OrderedDict{Int64, Vector{Float64}}
@@ -154,7 +148,6 @@ end
 
 mutable struct GeneratorLayout
     bus::Vector{Int64}
-    area::Vector{Float64}
     status::Vector{Int8}
     inservice::Int64
     label::Int64
@@ -169,7 +162,6 @@ A composite type used in [`PowerSystem`](@ref PowerSystem) to store generator-re
 - `label::LabelDict`: Generator labels.
 - `output::GeneratorOutput`: Active and reactive power outputs.
 - `capability::GeneratorCapability`: Power output constraints.
-- `ramping::GeneratorRamping`: Ramp rate limits.
 - `voltage::GeneratorVoltage`: Voltage magnitude setpoints.
 - `cost::GeneratorCost`: Costs associated with active and reactive power outputs.
 - `layout::GeneratorLayout`: Generator layout, including operational statuses.
@@ -179,7 +171,6 @@ mutable struct Generator
     label::LabelDict
     output::GeneratorOutput
     capability::GeneratorCapability
-    ramping::GeneratorRamping
     voltage::GeneratorVoltage
     cost::GeneratorCost
     layout::GeneratorLayout