add benchmark cases in a separate file

Author: hexaeder

benchmark/benchmark_cases.jl

@@ -0,0 +1,82 @@
+"""
+    BenchmarkCase
+
+Represents a single benchmark case for NetworkDynamics.
+
+Fields:
+- name: String identifier for the benchmark
+- constructor: Function that creates the Network for a given N
+- Ns: Vector of network sizes to benchmark
+- description: Documentation string
+"""
+struct BenchmarkCase
+    name::String
+    constructor::Function
+    Ns::Vector{Int}
+    description::String
+end
+
+# Helper functions to reduce code duplication
+function make_ws_graph(N, k, β=0.0; rng=StableRNG(1))
+    watts_strogatz(N, k, β; rng)
+end
+
+function make_mixed_vertices(v1, v2, N; rng=StableRNG(1))
+    vertex = [v1, v2]
+    vertex[shuffle(rng, vcat([1 for _ in 1:N÷2], [2 for _ in 1:N÷2]))]
+end
+
+# Standard benchmark cases
+const BENCHMARK_CASES = [
+    BenchmarkCase(
+        "diffusion_static_edge",
+        N -> begin
+            g = make_ws_graph(N, N ÷ 2)
+            Network(g, diffusion_vertex(), diffusion_edge())
+        end,
+        [100, 300, 1000, 3000],
+        "Diffusion network with static edges"
+    ),
+
+    BenchmarkCase(
+        "diffusion_ode_edge",
+        N -> begin
+            g = make_ws_graph(N, N ÷ 2)
+            Network(g, diffusion_vertex(), diffusion_dedge())
+        end,
+        [100, 300, 1000, 3000],
+        "Diffusion network with ODE edges"
+    ),
+
+    BenchmarkCase(
+        "kuramoto_homogeneous",
+        N -> begin
+            g = make_ws_graph(N, 3, 0.8)
+            Network(g, kuramoto_vertex_2d(), static_kuramoto_edge())
+        end,
+        [100, 1_000, 10_000, 100_000],
+        "Homogeneous Kuramoto oscillators"
+    ),
+
+    BenchmarkCase(
+        "kuramoto_heterogeneous",
+        N -> begin
+            g = make_ws_graph(N, 3, 0.8)
+            vertices = make_mixed_vertices(kuramoto_vertex_1d(), kuramoto_vertex_2d(), N)
+            Network(g, vertices, static_kuramoto_edge())
+        end,
+        [100, 1_000, 10_000, 100_000],
+        "Heterogeneous Kuramoto oscillators"
+    ),
+
+    BenchmarkCase(
+        "powergrid",
+        N -> begin
+            g = make_ws_graph(N, 3, 0.8)
+            vertices = make_mixed_vertices(pqnode(), generator(), N)
+            Network(g, vertices, piline())
+        end,
+        [100, 1_000, 10_000, 100_000],
+        "Power grid network with PQ nodes and generators"
+    )
+]