powerio-matrix 0.2.1

Sparse matrices and graph views for power system case files, built on powerio.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32

//! Square 2-D lattice / grid topology. `n` is rounded up to the nearest
//! perfect square.

use rand::SeedableRng;
use rand_chacha::ChaCha8Rng;

use crate::network::Network;

use super::SynthSpec;
use super::tree::{make_branch, make_buses, net};

pub fn generate_lattice(spec: &SynthSpec) -> Network {
    let side = ((spec.n as f64).sqrt().ceil() as usize).max(2);
    let n = side * side;
    let mut rng = ChaCha8Rng::seed_from_u64(spec.seed);

    let buses = make_buses(n);
    let mut branches = Vec::with_capacity(2 * side * (side - 1));
    for r in 0..side {
        for c in 0..side {
            let idx = r * side + c;
            if c + 1 < side {
                branches.push(make_branch(idx + 1, idx + 2, spec, &mut rng));
            }
            if r + 1 < side {
                branches.push(make_branch(idx + 1, idx + side + 1, spec, &mut rng));
            }
        }
    }

    net(format!("synth_lattice_{side}x{side}"), buses, branches)
}