rs-graph 0.20.1

A library for graph algorithms and combinatorial optimization
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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74

/*
 * Copyright (c) 2021 Frank Fischer <frank-fischer@shadow-soft.de>
 *
 * This program is free software: you can redistribute it and/or
 * modify it under the terms of the GNU General Public License as
 * published by the Free Software Foundation, either version 3 of the
 * License, or (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see  <http://www.gnu.org/licenses/>
 */

use std::collections::HashMap;
use std::error::Error;
use std::fs::read_dir;
use std::path::Path;

use rs_graph::dimacs;
use rs_graph::mcf::{MinCostFlow, NetworkSimplex, SolutionState};
use rs_graph::Net;
use rs_graph::{EdgeVec, NodeVec};

#[test]
fn test_network_simplex() -> Result<(), Box<dyn Error>> {
    let mut values = HashMap::new();

    for entry in read_dir(Path::new("tests/mcf"))? {
        let entry = entry?;
        if entry.path().extension().map(|ext| ext == "sol").unwrap_or(false) {
            let (value, _) = dimacs::min::read_solution_from_file(&entry.path().to_string_lossy())?;
            if let Some(file_stem) = entry.path().file_stem().map(|s| s.to_string_lossy().to_string()) {
                values.insert(file_stem, value);
            }
        }
    }

    for entry in read_dir(Path::new("tests/mcf"))? {
        let entry = entry?;
        if entry.path().extension().map(|ext| ext == "min").unwrap_or(false) {
            let instance = dimacs::min::read_from_file(&entry.path().to_string_lossy())?;
            let g: Net = instance.graph;

            let balances = NodeVec::new_from_vec(&g, instance.balances);
            let lower = EdgeVec::new_from_vec(&g, instance.lower);
            let upper = EdgeVec::new_from_vec(&g, instance.upper);
            let costs = EdgeVec::new_from_vec(&g, instance.costs);

            let mut spx = NetworkSimplex::new(&g);
            spx.set_balances(|u| balances[u]);
            spx.set_lowers(|e| lower[e]);
            spx.set_uppers(|e| upper[e]);
            spx.set_costs(|e| costs[e]);

            let state = spx.solve();
            assert_eq!(state, SolutionState::Optimal);
            if let Some(value) = entry
                .path()
                .file_name()
                .and_then(|s| values.get(s.to_string_lossy().as_ref()))
            {
                assert_eq!(*value, spx.value());
            } else {
                panic!("Can't find solution file for {:?}", entry.path());
            }
        }
    }

    Ok(())
}