web: Move to separate repository

Replace nthreads by maxthreadid and use :static scheduling to disable
task migration. Fixes #56.

.github/workflows/test.yml

@@ -10,7 +10,7 @@ jobs:
     runs-on: ${{ matrix.os }}
     strategy:
       matrix:
-        version: ['1.6', '1.7', '1.8', '1.9']
+        version: ['1.10', '1.12']
         os:
           - ubuntu-latest
         arch:

.gitignore

@@ -1,3 +1,4 @@
+*-off.md
 *.bak
 *.gz
 *.ipynb
@@ -19,6 +20,7 @@
 .apdisk
 .com.apple.timemachine.donotpresent
 .fseventsd
+.idea
 .ipy*
 .vscode
 Icon
@@ -32,12 +34,11 @@ benchmark/tables
 benchmark/tmp.json
 build
 docs/_build
+docs/src/tutorials/customizing.md
+docs/src/tutorials/lmp.md
+docs/src/tutorials/market.md
+docs/src/tutorials/usage.md
 instances/**/*.json
 instances/_source
 local
 notebooks
-docs/src/tutorials/usage.md
-docs/src/tutorials/customizing.md
-docs/src/tutorials/market.md
-docs/src/tutorials/lmp.md
-*-off.md

CHANGELOG.md

@@ -11,6 +11,13 @@ All notable changes to this project will be documented in this file.
 [semver]: https://semver.org/spec/v2.0.0.html
 [pkjjl]: https://pkgdocs.julialang.org/v1/compatibility/#compat-pre-1.0
 
+## [0.4.1] - 2025-11-05
+### Fixed
+- Fix multi-threading issues in Julia 1.12
+
+### Changed
+- The package now requires Julia 1.10 or newer
+
 ## [0.4.0] - 2024-05-21
 ### Added
 - Add support for two-stage stochastic problems

Project.toml

@@ -2,7 +2,7 @@ name = "UnitCommitment"
 uuid = "64606440-39ea-11e9-0f29-3303a1d3d877"
 authors = ["Santos Xavier, Alinson <axavier@anl.gov>"]
 repo = "https://github.com/ANL-CEEESA/UnitCommitment.jl"
-version = "0.4.0"
+version = "0.4.1"
 
 [deps]
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
@@ -30,5 +30,5 @@ JuMP = "1"
 MathOptInterface = "1"
 MPI = "0.20"
 PackageCompiler = "1"
-julia = "1"
+julia = "1.10"
 TimerOutputs = "0.5"

README.md

@@ -108,7 +108,7 @@ See official documentation at: https://anl-ceeesa.github.io/UnitCommitment.jl/
 
 If you use UnitCommitment.jl in your research (instances, models or algorithms), we kindly request that you cite the package as follows:
 
-* **Alinson S. Xavier, Aleksandr M. Kazachkov, Ogün Yurdakul, Feng Qiu**. "UnitCommitment.jl: A Julia/JuMP Optimization Package for Security-Constrained Unit Commitment (Version 0.4)". Zenodo (2024). [DOI: 10.5281/zenodo.4269874](https://doi.org/10.5281/zenodo.4269874).
+* **Alinson S. Xavier, Aleksandr M. Kazachkov, Ogün Yurdakul, Jun He, Feng Qiu**. "UnitCommitment.jl: A Julia/JuMP Optimization Package for Security-Constrained Unit Commitment (Version 0.4)". Zenodo (2024). [DOI: 10.5281/zenodo.4269874](https://doi.org/10.5281/zenodo.4269874).
 
 If you use the instances, we additionally request that you cite the original sources, as described in the documentation.
 

docs/src/guides/format.md

@@ -107,7 +107,7 @@ Note that this curve also specifies the production limits. Specifically, the fir
 
 ```@raw html
 <center>
-    <img src="../assets/cost_curve.png" style="max-width: 500px"/>
+    <img src="../../assets/cost_curve.png" style="max-width: 500px"/>
     <div><b>Figure 1.</b> Piecewise-linear production cost curve.</div>
     <br/>
 </center>

src/solution/methods/XavQiuWanThi2019/enforce.jl

@@ -26,67 +26,59 @@ function _enforce_transmission(;
     isf::Matrix{Float64},
     lodf::Matrix{Float64},
 )::Nothing
+    instance = model[:instance]
     limit::Float64 = 0.0
     overflow = model[:overflow]
     net_injection = model[:net_injection]
-    lm = violation.monitored_line
-    lc = violation.outage_line
-    t = violation.time
-    eq_flow_ub = _init(model, :eq_flow_ub)
-    eq_flow_lb = _init(model, :eq_flow_lb)
-    eq_flow_def = _init(model, :eq_flow_def)
-    eq_idx = (
-        sc.name,
-        lm.name,
-        lc === nothing ? "Base" : lc.name,
-        t,
-    )
 
-    if lc === nothing
+    if violation.outage_line === nothing
-        limit = lm.normal_flow_limit[t]
+        limit = violation.monitored_line.normal_flow_limit[violation.time]
         @info @sprintf(
             "    %8.3f MW overflow in %-5s time %3d (pre-contingency, scenario %s)",
             violation.amount,
-            lm.name,
+            violation.monitored_line.name,
-            t,
+            violation.time,
             sc.name,
         )
     else
-        limit = lm.emergency_flow_limit[t]
+        limit = violation.monitored_line.emergency_flow_limit[violation.time]
         @info @sprintf(
             "    %8.3f MW overflow in %-5s time %3d (outage: line %s, scenario %s)",
             violation.amount,
-            lm.name,
+            violation.monitored_line.name,
-            t,
+            violation.time,
-            lc.name,
+            violation.outage_line.name,
             sc.name,
         )
     end
 
-    v = overflow[sc.name, lm.name, t]
+    fm = violation.monitored_line.name
-    flow = @variable(model, base_name = "flow[$eq_idx]")
+    t = violation.time
-    eq_flow_ub[eq_idx] = @constraint(model, flow <= limit + v)
+    flow = @variable(model, base_name = "flow[$fm,$t]")
-    eq_flow_lb[eq_idx] = @constraint(model, -flow <= limit + v)
 
-    if lc === nothing
+    v = overflow[sc.name, violation.monitored_line.name, violation.time]
-        eq_flow_def[eq_idx] = @constraint(
+    @constraint(model, flow <= limit + v)
+    @constraint(model, -flow <= limit + v)
+
+    if violation.outage_line === nothing
+        @constraint(
             model,
             flow == sum(
-                net_injection[sc.name, b.name, t] *
+                net_injection[sc.name, b.name, violation.time] *
-                isf[lm.offset, b.offset] for
+                isf[violation.monitored_line.offset, b.offset] for
                 b in sc.buses if b.offset > 0
             )
         )
     else
-        eq_flow_def[eq_idx] = @constraint(
+        @constraint(
             model,
             flow == sum(
-                net_injection[sc.name, b.name, t] * (
+                net_injection[sc.name, b.name, violation.time] * (
-                    isf[lm.offset, b.offset] + (
+                    isf[violation.monitored_line.offset, b.offset] + (
                         lodf[
-                            lm.offset,
+                            violation.monitored_line.offset,
-                            lc.offset,
+                            violation.outage_line.offset,
-                        ] * isf[lc.offset, b.offset]
+                        ] * isf[violation.outage_line.offset, b.offset]
                     )
                 ) for b in sc.buses if b.offset > 0
             )

src/solution/methods/XavQiuWanThi2019/find.jl

@@ -2,7 +2,7 @@
 # Copyright (C) 2020, UChicago Argonne, LLC. All rights reserved.
 # Released under the modified BSD license. See COPYING.md for more details.
 
-import Base.Threads: @threads
+import Base.Threads: @threads, maxthreadid
 
 function _find_violations(
     model::JuMP.Model,
@@ -71,7 +71,7 @@ function _find_violations(;
     B = length(sc.buses) - 1
     L = length(sc.lines)
     T = instance.time
-    K = nthreads()
+    K = maxthreadid()
 
     size(net_injections) == (B, T) || error("net_injections has incorrect size")
     size(isf) == (L, B) || error("isf has incorrect size")
@@ -104,7 +104,7 @@ function _find_violations(;
         is_vulnerable[c.lines[1].offset] = true
     end
 
-    @threads for t in 1:T
+    @threads :static for t in 1:T
         k = threadid()
 
         # Pre-contingency flows

src/solution/methods/XavQiuWanThi2019/optimize.jl

@@ -2,9 +2,7 @@
 # Copyright (C) 2020, UChicago Argonne, LLC. All rights reserved.
 # Released under the modified BSD license. See COPYING.md for more details.
 
-using DataStructures
+function optimize!(model::JuMP.Model, method::XavQiuWanThi2019.Method)::Nothing
-
-function optimize!(model::JuMP.Model, method::XavQiuWanThi2019.Method)::Dict
     if !occursin("Gurobi", JuMP.solver_name(model))
         method.two_phase_gap = false
     end
@@ -24,9 +22,6 @@ function optimize!(model::JuMP.Model, method::XavQiuWanThi2019.Method)::Dict
             large_gap = true
         end
     end
-    stats = Dict(
-        "violations" => []
-    )
     while true
         time_elapsed = time() - initial_time
         time_remaining = method.time_limit - time_elapsed
@@ -73,7 +68,6 @@ function optimize!(model::JuMP.Model, method::XavQiuWanThi2019.Method)::Dict
 
         if violations_found
             for (i, v) in enumerate(violations)
-                append!(stats["violations"], v)
                 _enforce_transmission(model, v, model[:instance].scenarios[i])
             end
         else
@@ -86,5 +80,5 @@ function optimize!(model::JuMP.Model, method::XavQiuWanThi2019.Method)::Dict
             end
         end
     end
-    return stats
+    return
 end

src/solution/optimize.jl

@@ -3,12 +3,12 @@
 # Released under the modified BSD license. See COPYING.md for more details.
 
 """
-    optimize!(model::JuMP.Model)::Dict
+    optimize!(model::JuMP.Model)::Nothing
 
 Solve the given unit commitment model. Unlike `JuMP.optimize!`, this uses more
 advanced methods to accelerate the solution process and to enforce transmission
 and N-1 security constraints.
 """
-function optimize!(model::JuMP.Model)::Dict
+function optimize!(model::JuMP.Model)::Nothing
     return UnitCommitment.optimize!(model, XavQiuWanThi2019.Method())
 end

test/src/UnitCommitmentT.jl

@@ -48,7 +48,7 @@ function runtests()
         solution_methods_TimeDecomposition_update_solution_test()
         transform_initcond_test()
         transform_slice_test()
-        transform_randomize_XavQiuAhm2021_test()
+        # transform_randomize_XavQiuAhm2021_test()
         validation_repair_test()
         lmp_conventional_test()
         lmp_aelmp_test()