Struct net_ensembles::sw::SwEnsemble

pub struct SwEnsemble<T: Node, R: Rng> where
    T: Node,
    R: Rng,  { /* fields omitted */ }

Implements small-world graph ensemble

• for more details look at documentation of module sw

Sampling

• for markov steps look at MarkovChain trait
• for simple sampling look at SimpleSample trait

Other

• for topology functions look at GenericGraph
• to access underlying topology or manipulate additional data look at WithGraph trait
• to use or swap the random number generator, look at HasRng trait

Minimal example

use net_ensembles::{SwEnsemble, EmptyNode};
use net_ensembles::traits::*; // I recommend always using this
use rand_pcg::Pcg64; //or whatever you want to use as rng
use net_ensembles::rand::SeedableRng; // I use this to seed my rng, but you can use whatever

let rng = Pcg64::seed_from_u64(12);

// now create small-world ensemble with 200 nodes
// and a rewiring probability of 0.3 for each edge
let sw_ensemble = SwEnsemble::<EmptyNode, Pcg64>::new(200, 0.3, rng);

Simple sampling example

use net_ensembles::{SwEnsemble, EmptyNode};
use net_ensembles::traits::*; // I recommend always using this
use rand_pcg::Pcg64; //or whatever you want to use as rng
use net_ensembles::rand::SeedableRng; // I use this to seed my rng, but you can use whatever
use std::fs::File;
use std::io::{BufWriter, Write};

let rng = Pcg64::seed_from_u64(122);

// now create small-world ensemble with 100 nodes
// and a rewiring probability of 0.3 for each edge
let mut sw_ensemble = SwEnsemble::<EmptyNode, Pcg64>::new(100, 0.3, rng);

// setup file for writing
let f = File::create("simple_sample_sw_example.dat")
    .expect("Unable to create file");
let mut f = BufWriter::new(f);
f.write_all(b"#diameter bi_connect_max average_degree\n")
    .unwrap();

// simple sample for 10 steps
sw_ensemble.simple_sample(10,
    |ensemble|
    {
        let diameter = ensemble.graph()
            .diameter()
            .unwrap();

        let bi_connect_max = ensemble.graph()
            .clone()
            .vertex_biconnected_components(false)[0];

        let average_degree = ensemble.graph()
            .average_degree();

        write!(f, "{} {} {}\n", diameter, bi_connect_max, average_degree)
            .unwrap();
    }
);

// or just collect this into a vector to print or do whatever
let vec = sw_ensemble.simple_sample_vec(10,
    |ensemble|
    {
        let diameter = ensemble.graph()
            .diameter()
            .unwrap();

        let transitivity = ensemble.graph()
            .transitivity();
        (diameter, transitivity)
    }
);
println!("{:?}", vec);

Save and load example

• only works if feature "serde_support" is enabled
• Note: "serde_support" is enabled by default
• I need the #[cfg(feature = "serde_support")] to ensure the example does compile if you opt out of the default feature
• you can do not have to use serde_json, look here for more info

use net_ensembles::traits::*; // I recommend always using this
use serde_json;
use rand_pcg::Pcg64;
use net_ensembles::{SwEnsemble, EmptyNode, rand::SeedableRng};
use std::fs::File;

let rng = Pcg64::seed_from_u64(95);
// create small-world ensemble
let sw_ensemble = SwEnsemble::<EmptyNode, Pcg64>::new(200, 0.3, rng);

#[cfg(feature = "serde_support")]
{
    // storing the ensemble in a file:

    let sw_file = File::create("store_SW.dat")
          .expect("Unable to create file");

    // or serde_json::to_writer(sw_file, &sw_ensemble);
    serde_json::to_writer_pretty(sw_file, &sw_ensemble);

    // loading ensemble from file:

    let mut read = File::open("store_SW.dat")
        .expect("Unable to open file");

    let sw: SwEnsemble::<EmptyNode, Pcg64> = serde_json::from_reader(read).unwrap();
}

Implementations

impl<T, R> SwEnsemble<T, R> where
    T: Node + SerdeStateConform,
    R: Rng, 

pub fn new(n: usize, r_prob: f64, rng: R) -> Self

Initialize

• create new SwEnsemble graph with n vertices
• r_prob is probability of rewiring for each edge
• rng is consumed and used as random number generator in the following
• internally uses SwGraph<T>::new(n)

pub fn make_connected(&mut self)

Experimental! Connect the connected components

• resets edges, to connect the connected components
• intended as starting point for a markov chain, if you require connected graphs
• do not use this to independently (simple-) sample connected networks, as this will skew the statistics
• This is still experimental, this member might change the internal functionallity resulting in different connected networks, without prior notice
• This member might be removed in braking releases

pub fn draw_edge(&mut self) -> (usize, usize)

• draws random edge (i0, i1)
• edge rooted at i0
• uniform probability
• result dependent on order of adjecency lists
• mut because it uses the rng

pub fn sort_adj(&mut self)

Sort adjecency lists

If you depend on the order of the adjecency lists, you can sort them

Performance

• internally uses pattern-defeating quicksort as long as that is the standard
• sorts an adjecency list with length d in worst-case: O(d log(d))
• is called for each adjecency list, i.e., self.vertex_count() times

pub fn r_prob(&self) -> f64

• returns rewiring probability

pub fn set_r_prob(&mut self, r_prob: f64)

• set new value for rewiring probability

Note

• will only set the value, which will be used from now on
• if you also want to create a new sample, call randomize afterwards

pub fn contained_iter_neighbors_mut(
    &mut self,
    index: usize
) -> NContainedIterMut<T, SwContainer<T>>

• retuns GenericGraph::contained_iter_neighbors_mut
• otherwise you would not have access to this function, since no mut access to the graph is allowed

Trait Implementations

impl<T, R> AsRef<GenericGraph<T, SwContainer<T>>> for SwEnsemble<T, R> where
    T: Node,
    R: Rng, 

fn as_ref(&self) -> &SwGraph<T>

Performs the conversion.

impl<T, R> Borrow<GenericGraph<T, SwContainer<T>>> for SwEnsemble<T, R> where
    T: Node,
    R: Rng, 

fn borrow(&self) -> &SwGraph<T>

Immutably borrows from an owned value.

impl<T: Clone + Node, R: Clone + Rng> Clone for SwEnsemble<T, R> where
    T: Node,
    R: Rng, 

fn clone(&self) -> SwEnsemble<T, R>

Returns a copy of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<T: Debug + Node, R: Debug + Rng> Debug for SwEnsemble<T, R> where
    T: Node,
    R: Rng, 

fn fmt(&self, f: &mut Formatter) -> Result

Formats the value using the given formatter.

impl<'de, T: Node, R: Rng> Deserialize<'de> for SwEnsemble<T, R> where
    T: Node,
    R: Rng,
    T: Deserialize<'de>,
    R: Deserialize<'de>, 

fn deserialize<__D>(__deserializer: __D) -> Result<Self, __D::Error> where
    __D: Deserializer<'de>, 

Deserialize this value from the given Serde deserializer.

impl<T, R> GraphIteratorsMut<T, GenericGraph<T, SwContainer<T>>, SwContainer<T>> for SwEnsemble<T, R> where
    T: Node + SerdeStateConform,
    R: Rng, 

fn contained_iter_neighbors_mut(
    &mut self,
    index: usize
) -> NContainedIterMut<T, SwContainer<T>>

index``&mut T``sort_adj

fn contained_iter_neighbors_mut_with_index(
    &mut self,
    index: usize
) -> INContainedIterMut<T, SwContainer<T>>

index``(index_neighbor: usize, neighbor: &mut T)``sort_adj

fn contained_iter_mut(&mut self) -> ContainedIterMut<T, SwContainer<T>>

Node``EmptyNode

impl<T, R> HasRng<R> for SwEnsemble<T, R> where
    T: Node + SerdeStateConform,
    R: Rng, 

fn rng(&mut self) -> &mut R

Access RNG

If, for some reason, you want access to the internal random number generator: Here you go

fn swap_rng(&mut self, rng: R) -> R

Swap random number generator

• returns old internal rng

impl<T, R> MarkovChain<SwChangeState, SwChangeState> for SwEnsemble<T, R> where
    T: Node + SerdeStateConform,
    R: Rng, 

fn m_step(&mut self) -> SwChangeState

Markov step

• use this to perform a markov step
• keep in mind, that it is not unlikely for a step to do Nothing as it works by drawing an edge and then reseting it with r_prob, else the edge is rewired
• result SwChangeState can be used to undo the step with self.undo_step(result)
• result should never be InvalidAdjecency or GError if used on a valid graph

fn undo_step(&mut self, step: SwChangeState) -> SwChangeState

Undo a markov step

• rewires edge to old state
• Note: cannot undo InvalidAdjecency or GError, will just return InvalidAdjecency or GError respectively
• returns result of rewire

Important:

Restored graph is the same as before the random step except the order of nodes in the adjacency list might be shuffled!

fn undo_step_quiet(&mut self, step: SwChangeState)

Undo a Monte Carlo step

• rewires edge to old state
• panics if you try to undo InvalidAdjecency or GError
• panics if rewire result (SwChangeState) is invalid (i.e. !result.is_valid())

Important:

Restored graph is the same as before the random step except the order of nodes in the adjacency list might be shuffled!

fn m_steps(&mut self, count: usize) -> Vec<S>

Markov steps * use this to perform multiple markov steps at once * result Vec<S> can be used to undo the steps with self.undo_steps(result)

fn undo_steps(&mut self, steps: Vec<S>) -> Vec<Res>

Undo markov steps * Note: uses undo_step in correct order and returns result ## Important: * look at specific implementation of undo_step, every thing mentioned there applies to each step

fn undo_steps_quiet(&mut self, steps: Vec<S>)

Undo markov steps * Note: uses undo_step_quiet in correct order ## Important: * look at specific implementation of undo_step_quiet, every thing mentioned there applies to each step

impl<T: Node, R: Rng> Serialize for SwEnsemble<T, R> where
    T: Node,
    R: Rng,
    T: Serialize,
    R: Serialize, 

fn serialize<__S>(&self, __serializer: __S) -> Result<__S::Ok, __S::Error> where
    __S: Serializer, 

Serialize this value into the given Serde serializer.

impl<T, R> SimpleSample for SwEnsemble<T, R> where
    T: Node + SerdeStateConform,
    R: Rng, 

fn randomize(&mut self)

Randomizes the edges according to small-world model

• this is used by SwEnsemble::new to create the initial topology
• you can use this for sampling the ensemble
• runs in O(vertices)

fn simple_sample<F>(&mut self, times: usize, f: F) where
    F: FnMut(&Self), 

do the following times times: 1) f(self) 2) self.randomize()

fn simple_sample_vec<F, G>(&mut self, times: usize, f: F) -> Vec<G> where
    F: FnMut(&Self) -> G, 

do the following times times: 1) res = f(self) 2) self.randomize() ## res is collected into Vector

impl<T, R> WithGraph<T, GenericGraph<T, SwContainer<T>>> for SwEnsemble<T, R> where
    T: Node + SerdeStateConform,
    R: Rng, 

fn at(&self, index: usize) -> &T

fn at_mut(&mut self, index: usize) -> &mut T

fn graph(&self) -> &SwGraph<T>

GenericGraph