neural/inf/rs/cozo/cozo-core/tests/air_routes.rs

/*
 * Copyright 2022, The Cozo Project Authors.
 *
 * This Source Code Form is subject to the terms of the Mozilla Public License, v. 2.0.
 * If a copy of the MPL was not distributed with this file,
 * You can obtain one at https://mozilla.org/MPL/2.0/.
 */
#![cfg(feature = "graph-algo")]

use std::env;
use std::str::FromStr;
use std::time::Instant;

use approx::AbsDiffEq;
use env_logger::Env;
use lazy_static::{initialize, lazy_static};
use serde_json::json;

use cozo::DbInstance;

lazy_static! {
    static ref TEST_DB: DbInstance = {
        env_logger::Builder::from_env(Env::default().default_filter_or("info")).init();
        let creation = Instant::now();
        let path = "_test_air_routes";
        _ = std::fs::remove_file(path);
        _ = std::fs::remove_dir_all(path);
        let db_kind = env::var("COZO_TEST_DB_ENGINE").unwrap_or("mem".to_string());
        println!("Using {} engine", db_kind);
        let db = DbInstance::new(&db_kind, path, Default::default()).unwrap();
        dbg!(creation.elapsed());

        let init = Instant::now();
        db.run_default(r##"
            res[idx, label, typ, code, icao, desc, region, runways, longest, elev, country, city, lat, lon] <~
                CsvReader(types: ['Int', 'Any', 'Any', 'Any', 'Any', 'Any', 'Any', 'Int?', 'Float?', 'Float?', 'Any', 'Any', 'Float?', 'Float?'],
                          url: 'file://./tests/air-routes-latest-nodes.csv',
                          has_headers: true)

            ?[code, icao, desc, region, runways, longest, elev, country, city, lat, lon] :=
                res[idx, label, typ, code, icao, desc, region, runways, longest, elev, country, city, lat, lon],
                label == 'airport'

            :replace airport {
                code: String
                =>
                icao: String,
                desc: String,
                region: String,
                runways: Int,
                longest: Float,
                elev: Float,
                country: String,
                city: String,
                lat: Float,
                lon: Float
            }
        "##).unwrap();

        db.run_default(
            r##"
            res[idx, label, typ, code, icao, desc] <~
                CsvReader(types: ['Int', 'Any', 'Any', 'Any', 'Any', 'Any'],
                          url: 'file://./tests/air-routes-latest-nodes.csv',
                          has_headers: true)
            ?[code, desc] :=
                res[idx, label, typ, code, icao, desc],
                label == 'country'

            :replace country {
                code: String
                =>
                desc: String
            }
        "##,
        )
        .unwrap();

        db.run_default(
            r##"
            res[idx, label, typ, code, icao, desc] <~
                CsvReader(types: ['Int', 'Any', 'Any', 'Any', 'Any', 'Any'],
                          url: 'file://./tests/air-routes-latest-nodes.csv',
                          has_headers: true)
            ?[idx, code, desc] :=
                res[idx, label, typ, code, icao, desc],
                label == 'continent'

            :replace continent {
                code: String
                =>
                desc: String
            }
        "##,
        )
        .unwrap();

        db.run_default(
            r##"
            res[idx, label, typ, code] <~
                CsvReader(types: ['Int', 'Any', 'Any', 'Any'],
                          url: 'file://./tests/air-routes-latest-nodes.csv',
                          has_headers: true)
            ?[idx, code] :=
                res[idx, label, typ, code],

            :replace idx2code { idx: Int => code: String }
        "##,
        )
        .unwrap();

        db.run_default(
            r##"
            res[] <~
                CsvReader(types: ['Int', 'Int', 'Int', 'String', 'Float?'],
                          url: 'file://./tests/air-routes-latest-edges.csv',
                          has_headers: true)
            ?[fr, to, dist] :=
                res[idx, fr_i, to_i, typ, dist],
                typ == 'route',
                *idx2code[fr_i, fr],
                *idx2code[to_i, to]

            :replace route { fr: String, to: String => dist: Float }
        "##,
        )
        .unwrap();

        db.run_default(
            r##"
            res[] <~
                CsvReader(types: ['Int', 'Int', 'Int', 'String'],
                          url: 'file://./tests/air-routes-latest-edges.csv',
                          has_headers: true)
            ?[entity, contained] :=
                res[idx, fr_i, to_i, typ],
                typ == 'contains',
                *idx2code[fr_i, entity],
                *idx2code[to_i, contained]


            :replace contain { entity: String, contained: String }
        "##,
        )
        .unwrap();

        db.run_default("::remove idx2code").unwrap();

        dbg!(init.elapsed());
        db
    };
}

#[test]
fn dfs() {
    initialize(&TEST_DB);
    let dfs = Instant::now();
    let rows = TEST_DB
        .run_default(
            r#"
        starting[] <- 'PEK'
        ?[] <~ DFS(*route[], *airport[code], starting[], condition: (code == 'LHR'))
    "#,
        )
        .unwrap()
        .rows;
    assert_eq!(rows.len(), 1);
    let row = rows.get(0).unwrap();
    assert_eq!(row.get(0).unwrap().get_str().unwrap(), "PEK");
    assert_eq!(row.get(1).unwrap().get_str().unwrap(), "LHR");
    let path = row.get(2).unwrap().get_slice().unwrap();
    assert_eq!(path.first().unwrap().get_str().unwrap(), "PEK");
    assert_eq!(path.last().unwrap().get_str().unwrap(), "LHR");
    dbg!(dfs.elapsed());
}

#[test]
fn empty() {
    initialize(&TEST_DB);
    let res = TEST_DB.run_default(
        r#"
        ?[id, name] <- [[]]
    "#,
    );
    assert!(res.is_err());
}

#[test]
fn parallel_counts() {
    initialize(&TEST_DB);
    let res = TEST_DB
        .run_default(
            r#"
        a[count(fr)] := *route{fr}
        b[count(fr)] := *route{fr}
        c[count(fr)] := *route{fr}
        d[count(fr)] := *route{fr}
        e[count(fr)] := *route{fr}
        ?[x] := a[a], b[b], c[c], d[d], e[e], x = a + b + c + d + e
        "#,
        )
        .unwrap()
        .rows
        .remove(0)
        .remove(0)
        .get_int()
        .unwrap();
    assert_eq!(res, 50637 * 5);
}

#[test]
fn bfs() {
    initialize(&TEST_DB);
    let bfs = Instant::now();
    let rows = TEST_DB
        .run_default(
            r#"
        starting[] <- 'PEK'
        ?[] <~ BFS(*route[], *airport[code], starting[], condition: (code == 'LHR'))
    "#,
        )
        .unwrap()
        .rows;

    assert_eq!(rows.len(), 1);
    let row = rows.get(0).unwrap();
    assert_eq!(row.get(0).unwrap().get_str().unwrap(), "PEK");
    assert_eq!(row.get(1).unwrap().get_str().unwrap(), "LHR");
    let path = row.get(2).unwrap().get_slice().unwrap();
    assert_eq!(path.first().unwrap().get_str().unwrap(), "PEK");
    assert_eq!(path.last().unwrap().get_str().unwrap(), "LHR");
    dbg!(bfs.elapsed());
}

#[test]
fn scc() {
    initialize(&TEST_DB);
    let scc = Instant::now();
    let _ = TEST_DB
        .run_default(
            r#"
        res[] <~ StronglyConnectedComponents(*route[], *airport[code]);
        ?[grp, code] := res[code, grp], grp != 0;
    "#,
        )
        .unwrap()
        .rows;
    dbg!(scc.elapsed());
}

#[test]
fn cc() {
    initialize(&TEST_DB);
    let cc = Instant::now();
    let _ = TEST_DB
        .run_default(
            r#"
        res[] <~ ConnectedComponents(*route[], *airport[code]);
        ?[grp, code] := res[code, grp], grp != 0;
    "#,
        )
        .unwrap()
        .rows;
    dbg!(cc.elapsed());
}

#[test]
fn astar() {
    initialize(&TEST_DB);
    let astar = Instant::now();
    let _ = TEST_DB.run_default(r#"
        code_lat_lon[code, lat, lon] := *airport{code, lat, lon}
        starting[code, lat, lon] := code = 'HFE', *airport{code, lat, lon};
        goal[code, lat, lon] := code = 'LHR', *airport{code, lat, lon};
        ?[] <~ ShortestPathAStar(*route[], code_lat_lon[node, lat1, lon1], starting[], goal[goal, lat2, lon2], heuristic: haversine_deg_input(lat1, lon1, lat2, lon2) * 3963);
    "#).unwrap().rows;
    dbg!(astar.elapsed());
}

#[test]
fn deg_centrality() {
    initialize(&TEST_DB);
    let deg_centrality = Instant::now();
    TEST_DB
        .run_default(
            r#"
        deg_centrality[] <~ DegreeCentrality(*route[a, b]);
        ?[total, out, in] := deg_centrality[node, total, out, in];
        :order -total;
        :limit 10;
    "#,
        )
        .unwrap()
        .rows;
    dbg!(deg_centrality.elapsed());
}

#[test]
fn dijkstra() {
    initialize(&TEST_DB);
    let dijkstra = Instant::now();

    TEST_DB
        .run_default(
            r#"
        starting[] <- 'JFK';
        ending[] <- 'KUL';
        res[] <~ ShortestPathDijkstra(*route[], starting[], ending[]);
        ?[path] := res[src, dst, cost, path];
    "#,
        )
        .unwrap()
        .rows;

    dbg!(dijkstra.elapsed());
}

#[test]
fn yen() {
    initialize(&TEST_DB);
    let yen = Instant::now();

    TEST_DB
        .run_default(
            r#"
        starting[] <- 'PEK';
        ending[] <- 'SIN';
        ?[] <~ KShortestPathYen(*route[], starting[], ending[], k: 5);
    "#,
        )
        .unwrap()
        .rows;

    dbg!(yen.elapsed());
}

#[test]
fn starts_with() {
    initialize(&TEST_DB);
    let starts_with = Instant::now();
    let rows = TEST_DB
        .run_default(
            r#"
         ?[code] := *airport{code}, starts_with(code, 'US');
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        json!([
            ["USA"],
            ["USH"],
            ["USJ"],
            ["USK"],
            ["USM"],
            ["USN"],
            ["USQ"],
            ["UST"],
            ["USU"]
        ])
    );

    dbg!(starts_with.elapsed());
}

#[test]
fn range_check() {
    initialize(&TEST_DB);
    let range_check = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        r[code, dist] := *airport{code}, *route{fr: code, dist};
        ?[dist] := r['PEK', dist], dist > 7000, dist <= 7722;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        json!([[7176.0], [7270.0], [7311.0], [7722.0]])
    );
    dbg!(range_check.elapsed());
}

#[test]
fn no_airports() {
    initialize(&TEST_DB);
    let no_airports = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[desc] := *country{code, desc}, not *airport{country: code};
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        json!([
            ["Andorra"],
            ["Liechtenstein"],
            ["Monaco"],
            ["Pitcairn"],
            ["San Marino"]
        ])
    );
    dbg!(no_airports.elapsed());
}

#[test]
fn no_routes_airport() {
    initialize(&TEST_DB);
    let no_routes_airports = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[code] := *airport{code}, not *route{fr: code}, not *route{to: code}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["AFW"],["APA"],["APK"],["BID"],["BVS"],["BWU"],["CRC"],["CVT"],["EKA"],["GYZ"],
            ["HFN"],["HZK"],["ILG"],["INT"],["ISL"],["KGG"],["NBW"],["NFO"],["PSY"],["RIG"],
            ["SFD"],["SFH"],["SXF"],["TUA"],["TWB"],["TXL"],["VCV"],["YEI"]
        ]"#
        )
        .unwrap()
    );
    dbg!(no_routes_airports.elapsed());
}

#[test]
fn runway_distribution() {
    initialize(&TEST_DB);
    let no_routes_airports = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[runways, count(code)] := *airport{code, runways}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        json!([
            [1, 2429],
            [2, 775],
            [3, 227],
            [4, 53],
            [5, 14],
            [6, 4],
            [7, 2]
        ])
    );
    dbg!(no_routes_airports.elapsed());
}

#[test]
fn most_out_routes() {
    initialize(&TEST_DB);
    let most_out_routes = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[fr, count(fr)] := *route{fr};
        ?[code, n] := route_count[code, n], n > 180;
        :sort -n;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["FRA",310],["IST",309],["CDG",293],["AMS",283],["MUC",270],["ORD",265],["DFW",253],
            ["DXB",248],["PEK",248],["ATL",242],["DME",232],["LGW",232],["LHR",221],["DEN",217],
            ["MAN",216],["LAX",214],["PVG",213],["STN",211],["MAD",206],["VIE",206],["JFK",204],
            ["BCN",203],["EWR",203],["BER",202],["FCO",201],["DUS",199],["IAH",199],["MIA",196],
            ["YYZ",195],["BRU",194],["CPH",194],["DOH",187],["DUB",185],["CLT",184],["SVO",181]
            ]"#
        )
        .unwrap()
    );
    dbg!(most_out_routes.elapsed());
}

#[test]
fn most_out_routes_again() {
    initialize(&TEST_DB);
    let most_out_routes_again = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[count(fr), fr] := *route{fr};
        ?[code, n] := route_count[n, code], n > 180;
        :sort -n;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["FRA",310],["IST",309],["CDG",293],["AMS",283],["MUC",270],["ORD",265],["DFW",253],
            ["DXB",248],["PEK",248],["ATL",242],["DME",232],["LGW",232],["LHR",221],["DEN",217],
            ["MAN",216],["LAX",214],["PVG",213],["STN",211],["MAD",206],["VIE",206],["JFK",204],
            ["BCN",203],["EWR",203],["BER",202],["FCO",201],["DUS",199],["IAH",199],["MIA",196],
            ["YYZ",195],["BRU",194],["CPH",194],["DOH",187],["DUB",185],["CLT",184],["SVO",181]
            ]"#
        )
        .unwrap()
    );
    dbg!(most_out_routes_again.elapsed());
}

#[test]
fn most_routes() {
    initialize(&TEST_DB);
    let most_routes = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[a, count(a)] := *route{fr: a}
        route_count[a, count(a)] := *route{to: a}
        ?[code, n] := route_count[code, n], n > 400
        :sort -n;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["FRA",620],["IST",618],["CDG",587],["AMS",568],["MUC",541],["ORD",529],["DFW",506],
            ["PEK",497],["DXB",496],["ATL",484],["DME",465],["LGW",464],["LHR",442],["DEN",434],
            ["MAN",431],["LAX",428],["PVG",426],["STN",423],["MAD",412],["VIE",412],["JFK",407],
            ["BCN",406],["EWR",406],["BER",404],["FCO",402]]"#
        )
        .unwrap()
    );
    dbg!(most_routes.elapsed());
}

#[test]
fn airport_with_one_route() {
    initialize(&TEST_DB);
    let airport_with_one_route = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[fr, count(fr)] := *route{fr}
        ?[count(a)] := route_count[a, n], n == 1;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[777]]));
    dbg!(airport_with_one_route.elapsed());
}

#[test]
fn single_runway_with_most_routes() {
    initialize(&TEST_DB);
    let single_runway_with_most_routes = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        single_or_lgw[code] := code = 'LGW'
        single_or_lgw[code] := *airport{code, runways}, runways == 1
        out_counts[a, count(a)] := single_or_lgw[a], *route{fr: a}
        ?[code, city, out_n] := out_counts[code, out_n], *airport{code, city}

        :order -out_n;
        :limit 10;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        ["LGW","London",232],["STN","London",211],["CTU","Chengdu",139],["LIS","Lisbon",139],
        ["BHX","Birmingham",130],["LTN","London",130],["SZX","Shenzhen",129],
        ["CKG","Chongqing",122],["STR","Stuttgart",121],["CRL","Brussels",117]]"#
        )
        .unwrap()
    );
    dbg!(single_runway_with_most_routes.elapsed());
}

#[test]
fn most_routes_in_canada() {
    initialize(&TEST_DB);
    let most_routes_in_canada = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ca_airports[code, count(code)] := *airport{code, country: 'CA'}, *route{fr: code}
        ?[code, city, n_routes] := ca_airports[code, n_routes], *airport{code, city}

        :order -n_routes;
        :limit 10;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        json!([
            ["YYZ", "Toronto", 195],
            ["YUL", "Montreal", 123],
            ["YVR", "Vancouver", 106],
            ["YYC", "Calgary", 75],
            ["YEG", "Edmonton", 48],
            ["YHZ", "Halifax", 45],
            ["YWG", "Winnipeg", 38],
            ["YOW", "Ottawa", 36],
            ["YZF", "Yellowknife", 21],
            ["YQB", "Quebec City", 20]
        ])
    );
    dbg!(most_routes_in_canada.elapsed());
}

#[test]
fn uk_count() {
    initialize(&TEST_DB);
    let uk_count = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
       ?[region, count(region)] := *airport{country: 'UK', region}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        json!([["GB-ENG", 27], ["GB-NIR", 3], ["GB-SCT", 25], ["GB-WLS", 3]])
    );
    dbg!(uk_count.elapsed());
}

#[test]
fn airports_by_country() {
    initialize(&TEST_DB);
    let airports_by_country = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        airports_by_country[country, count(code)] := *airport{code, country}
        ?[country, count] := airports_by_country[country, count];
        ?[country, count] := *country{code: country}, not airports_by_country[country, _], count = 0

        :order count
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
    ["AD",0],["LI",0],["MC",0],["PN",0],["SM",0],["AG",1],["AI",1],["AL",1],["AS",1],["AW",1],
    ["BB",1],["BH",1],["BI",1],["BJ",1],["BL",1],["BM",1],["BN",1],["BT",1],["CC",1],["CF",1],
    ["CW",1],["CX",1],["DJ",1],["DM",1],["ER",1],["FO",1],["GD",1],["GF",1],["GI",1],["GM",1],
    ["GN",1],["GP",1],["GU",1],["GW",1],["HK",1],["IM",1],["JE",1],["KM",1],["KP",1],["KS",1],
    ["KW",1],["LB",1],["LS",1],["LU",1],["LV",1],["MD",1],["MF",1],["ML",1],["MO",1],["MQ",1],
    ["MS",1],["MT",1],["NC",1],["NE",1],["NF",1],["NI",1],["NR",1],["PM",1],["PW",1],["QA",1],
    ["SL",1],["SR",1],["SS",1],["ST",1],["SV",1],["SX",1],["SZ",1],["TG",1],["TL",1],["TM",1],
    ["TV",1],["VC",1],["WS",1],["YT",1],["AM",2],["BF",2],["CI",2],["EH",2],["FK",2],["GA",2],
    ["GG",2],["GQ",2],["GT",2],["GY",2],["HT",2],["HU",2],["JM",2],["JO",2],["KG",2],["KI",2],
    ["KN",2],["LC",2],["LR",2],["ME",2],["MH",2],["MK",2],["MP",2],["MU",2],["PY",2],["RE",2],
    ["RW",2],["SC",2],["SG",2],["SH",2],["SI",2],["SK",2],["SY",2],["TT",2],["UY",2],["VG",2],
    ["VI",2],["WF",2],["BQ",3],["BY",3],["CG",3],["CY",3],["EE",3],["GE",3],["KH",3],["KY",3],
    ["LT",3],["MR",3],["RS",3],["ZW",3],["BA",4],["BG",4],["BW",4],["FM",4],["OM",4],["SN",4],
    ["TC",4],["TJ",4],["UG",4],["AF",5],["AZ",5],["BE",5],["CM",5],["CZ",5],["NA",5],["NL",5],
    ["PA",5],["SD",5],["TD",5],["TO",5],["AT",6],["CH",6],["CK",6],["GH",6],["HN",6],["IL",6],
    ["IQ",6],["LK",6],["SO",6],["BD",7],["CV",7],["DO",7],["IE",7],["IS",7],["MW",7],["PR",7],
    ["DK",8],["HR",8],["LA",8],["MV",8],["TN",8],["TW",9],["YE",9],["ZM",9],["AE",10],["FJ",10],
    ["MN",10],["CD",11],["EG",11],["LY",11],["MZ",11],["NP",11],["TZ",11],["UZ",11],["CU",12],
    ["BZ",13],["CR",13],["MG",13],["PL",13],["AO",14],["GL",14],["KE",14],["RO",14],["BO",15],
    ["EC",15],["KR",15],["UA",15],["ET",16],["MA",16],["CL",17],["MM",17],["SB",17],["BS",18],
    ["NG",19],["PT",19],["FI",20],["ZA",20],["KZ",21],["PK",21],["PE",22],["VN",22],["NZ",25],
    ["PG",26],["SA",26],["VU",26],["VE",27],["DZ",30],["TH",33],["DE",34],["MY",35],["AR",38],
    ["IT",38],["GR",39],["PF",39],["SE",39],["PH",40],["ES",43],["IR",45],["NO",49],["CO",51],
    ["TR",52],["UK",58],["FR",59],["MX",60],["JP",65],["ID",70],["IN",77],["BR",117],["RU",129],
    ["AU",132],["CA",205],["CN",217],["US",586]]"#
        )
        .unwrap()
    );
    dbg!(airports_by_country.elapsed());
}

#[test]
fn n_airports_by_continent() {
    initialize(&TEST_DB);
    let n_airports_by_continent = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        airports_by_continent[cont, count(code)] := *airport{code}, *contain[cont, code]
        ?[cont, max(count)] := *continent{code: cont}, airports_by_continent[cont, count]
        ?[cont, max(count)] := *continent{code: cont}, count = 0
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[["AF",326],["AN",0],["AS",972],["EU",605],["NA",994],["OC",305],["SA",339]]"#
        )
        .unwrap()
    );
    dbg!(n_airports_by_continent.elapsed());
}

#[test]
fn routes_per_airport() {
    initialize(&TEST_DB);
    let routes_per_airport = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        given[] <- [['A' ++ 'U' ++ 'S'],['AMS'],['JFK'],['DUB'],['MEX']]
        ?[code, count(code)] := given[code], *route{fr: code}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[["AMS",283],["AUS",98],["DUB",185],["JFK",204],["MEX",116]]"#
        )
        .unwrap()
    );
    dbg!(routes_per_airport.elapsed());
}

#[test]
fn airports_by_route_number() {
    initialize(&TEST_DB);
    let airports_by_route_number = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[fr, count(fr)] := *route{fr}
        ?[n, collect(code)] := route_count[code, n], n = 106;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[106, ["TFS", "YVR"]]]));
    dbg!(airports_by_route_number.elapsed());
}

#[test]
fn out_from_aus() {
    initialize(&TEST_DB);
    let out_from_aus = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        out_by_runways[runways, count(code)] := *route{fr: 'AUS', to: code}, *airport{code, runways}
        two_hops[count(a)] := *route{fr: 'AUS', to: a}, *route{fr: a}
        ?[max(total), collect(coll)] := two_hops[total], out_by_runways[n, ct], coll = [n, ct];
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"[[8354,[[1,9],[2,24],[3,30],[4,24],[5,5],[6,4],[7,2]]]]"#)
            .unwrap()
    );
    dbg!(out_from_aus.elapsed());
}

#[test]
fn const_return() {
    initialize(&TEST_DB);
    let const_return = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[name, count(code)] := *airport{code, region: 'US-OK'}, name = 'OK';
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([["OK", 4]]));
    dbg!(const_return.elapsed());
}

#[test]
fn multi_res() {
    initialize(&TEST_DB);
    let multi_res = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        total[count(code)] := *airport{code}
        high[count(code)] := *airport{code, runways}, runways >= 6
        low[count(code)] := *airport{code, runways}, runways <= 2
        four[count(code)] := *airport{code, runways}, runways == 4
        france[count(code)] := *airport{code, country: 'FR'}

        ?[total, high, low, four, france] := total[total], high[high], low[low],
                                                  four[four], france[france];
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"3504,6,3204,53,59"#).unwrap()
    );
    dbg!(multi_res.elapsed());
}

#[test]
fn multi_unification() {
    initialize(&TEST_DB);
    let multi_unification = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        target_airports[collect(code, 5)] := *airport{code}
        ?[a, count(a)] := target_airports[targets], a in targets, *route{fr: a}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"[["AAA",4],["AAE",8],["AAL",17],["AAN",5],["AAQ",11]]"#)
            .unwrap()
    );
    dbg!(multi_unification.elapsed());
}

#[test]
fn num_routes_from_eu_to_us() {
    initialize(&TEST_DB);
    let num_routes_from_eu_to_us = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        routes[unique(r)] := *contain['EU', fr],
                             *route{fr, to},
                             *airport{code: to, country: 'US'},
                             r = [fr, to]
        ?[n] := routes[rs], n = length(rs);
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[435]]));
    dbg!(num_routes_from_eu_to_us.elapsed());
}

#[test]
fn num_airports_in_us_with_routes_from_eu() {
    initialize(&TEST_DB);
    let num_airports_in_us_with_routes_from_eu = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[count_unique(to)] := *contain['EU', fr],
                               *route{fr, to},
                               *airport{code: to, country: 'US'}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[45]]));
    dbg!(num_airports_in_us_with_routes_from_eu.elapsed());
}

#[test]
fn num_routes_in_us_airports_from_eu() {
    initialize(&TEST_DB);
    let num_routes_in_us_airports_from_eu = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[to, count(to)] := *contain['EU', fr], *route{fr, to}, *airport{code: to, country: 'US'}
        :order count(to);
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["ANC",1],["BNA",1],["CHS",1],["CLE",1],["IND",1],["MCI",1],["BDL",2],["BWI",2],
            ["CVG",2],["MSY",2],["PHX",2],["SJC",2],["STL",2],["PDX",3],["RDU",3],["SAN",3],
            ["AUS",4],["PIT",4],["RSW",4],["SLC",4],["SFB",5],["SWF",5],["TPA",5],["DTW",6],
            ["MSP",6],["OAK",6],["DEN",7],["FLL",7],["PVD",7],["CLT",8],["IAH",8],["LAS",11],
            ["DFW",12],["SEA",12],["MCO",14],["ATL",15],["SFO",20],["IAD",22],["PHL",22],["BOS",26],
            ["LAX",26],["ORD",27],["MIA",28],["JFK",42],["EWR",43]]"#
        )
        .unwrap()
    );
    dbg!(num_routes_in_us_airports_from_eu.elapsed());
}

#[test]
fn routes_from_eu_to_us_starting_with_l() {
    initialize(&TEST_DB);
    let routes_from_eu_to_us_starting_with_l = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[eu_code, us_code] := *contain['EU', eu_code],
                               starts_with(eu_code, 'L'),
                               *route{fr: eu_code, to: us_code},
                               *airport{code: us_code, country: 'US'}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
    ["LGW","AUS"],["LGW","BOS"],["LGW","DEN"],["LGW","FLL"],["LGW","JFK"],["LGW","LAS"],
    ["LGW","LAX"],["LGW","MCO"],["LGW","MIA"],["LGW","OAK"],["LGW","ORD"],["LGW","SEA"],
    ["LGW","SFO"],["LGW","TPA"],["LHR","ATL"],["LHR","AUS"],["LHR","BNA"],["LHR","BOS"],
    ["LHR","BWI"],["LHR","CHS"],["LHR","CLT"],["LHR","DEN"],["LHR","DFW"],["LHR","DTW"],
    ["LHR","EWR"],["LHR","IAD"],["LHR","IAH"],["LHR","JFK"],["LHR","LAS"],["LHR","LAX"],
    ["LHR","MIA"],["LHR","MSP"],["LHR","MSY"],["LHR","ORD"],["LHR","PDX"],["LHR","PHL"],
    ["LHR","PHX"],["LHR","PIT"],["LHR","RDU"],["LHR","SAN"],["LHR","SEA"],["LHR","SFO"],
    ["LHR","SJC"],["LHR","SLC"],["LIS","ATL"],["LIS","BOS"],["LIS","EWR"],["LIS","IAD"],
    ["LIS","JFK"],["LIS","MIA"],["LIS","ORD"],["LIS","PHL"],["LIS","SFO"]
            ]"#
        )
        .unwrap()
    );
    dbg!(routes_from_eu_to_us_starting_with_l.elapsed());
}

#[test]
fn len_of_names_count() {
    initialize(&TEST_DB);
    let len_of_names_count = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[sum(n)] := *route{fr: 'AUS', to},
                     *airport{code: to, city},
                     n = length(city)
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"891.0"#).unwrap()
    );
    dbg!(len_of_names_count.elapsed());
}

#[test]
fn group_count_by_out() {
    initialize(&TEST_DB);
    let group_count_by_out = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[count(fr), fr] := *route{fr}
        rc[max(n), a] := route_count[n, a]
        rc[max(n), a] := *airport{code: a}, n = 0
        ?[n, count(a)] := rc[n, a]
        :order n;
        :limit 10;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[[0,29],[1,777],[2,649],[3,357],[4,234],[5,149],[6,140],[7,100],[8,73],[9,64]]"#
        )
        .unwrap()
    );
    dbg!(group_count_by_out.elapsed());
}

#[test]
fn mean_group_count() {
    initialize(&TEST_DB);
    let mean_group_count = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        route_count[count(fr), fr] := *route{fr};
        rc[max(n), a] := route_count[n, a] or (*airport{code: a}, n = 0);
        ?[mean(n)] := rc[n, _];
    "#,
        )
        .unwrap()
        .rows;

    let v = rows.get(0).unwrap().get(0).unwrap().get_float().unwrap();

    assert!(v.abs_diff_eq(&14.451198630136986, 1e-8));
    dbg!(mean_group_count.elapsed());
}

#[test]
fn n_routes_from_london_uk() {
    initialize(&TEST_DB);
    let n_routes_from_london_uk = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[code, count(code)] := *airport{code, city: 'London', region: 'GB-ENG'}, *route{fr: code}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[["LCY",51],["LGW",232],["LHR",221],["LTN",130],["STN",211]]"#
        )
        .unwrap()
    );
    dbg!(n_routes_from_london_uk.elapsed());
}

#[test]
fn reachable_from_london_uk_in_two_hops() {
    initialize(&TEST_DB);
    let reachable_from_london_uk_in_two_hops = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        lon_uk_airports[code] := *airport{code, city: 'London', region: 'GB-ENG'}
        one_hop[to] := lon_uk_airports[fr], *route{fr, to}, not lon_uk_airports[to];
        ?[count_unique(a3)] := one_hop[a2], *route{fr: a2, to: a3}, not lon_uk_airports[a3];
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[2353]]));
    dbg!(reachable_from_london_uk_in_two_hops.elapsed());
}

#[test]
fn routes_within_england() {
    initialize(&TEST_DB);
    let routes_within_england = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        eng_aps[code] := *airport{code, region: 'GB-ENG'}
        ?[fr, to] := eng_aps[fr], *route{fr, to}, eng_aps[to],
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
    ["BHX","NCL"],["BRS","NCL"],["EMA","SOU"],["EXT","ISC"],["EXT","MAN"],["EXT","NQY"],
    ["HUY","NWI"],["ISC","EXT"],["ISC","LEQ"],["ISC","NQY"],["LBA","LHR"],["LBA","NQY"],
    ["LBA","SOU"],["LCY","MAN"],["LCY","NCL"],["LEQ","ISC"],["LGW","NCL"],["LGW","NQY"],
    ["LHR","LBA"],["LHR","MAN"],["LHR","NCL"],["LHR","NQY"],["LPL","NQY"],["MAN","EXT"],
    ["MAN","LCY"],["MAN","LHR"],["MAN","NQY"],["MAN","NWI"],["MAN","SEN"],["MAN","SOU"],
    ["MME","NWI"],["NCL","BHX"],["NCL","BRS"],["NCL","LCY"],["NCL","LGW"],["NCL","LHR"],
    ["NCL","SOU"],["NQY","EXT"],["NQY","ISC"],["NQY","LBA"],["NQY","LGW"],["NQY","LHR"],
    ["NQY","LPL"],["NQY","MAN"],["NQY","SEN"],["NWI","HUY"],["NWI","MAN"],["NWI","MME"],
    ["SEN","MAN"],["SEN","NQY"],["SOU","EMA"],["SOU","LBA"],["SOU","MAN"],["SOU","NCL"]]"#
        )
        .unwrap()
    );
    dbg!(routes_within_england.elapsed());
}

#[test]
fn routes_within_england_time_no_dup() {
    initialize(&TEST_DB);
    let routes_within_england_time_no_dup = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        eng_aps[code] := *airport{code, region: 'GB-ENG'}
        ?[pair] := eng_aps[fr], *route{fr, to}, eng_aps[to], pair = sorted([fr, to]);
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
    "BHX","NCL","BRS","NCL","EMA","SOU","EXT","ISC","EXT","MAN","EXT","NQY",
    "HUY","NWI","ISC","LEQ","ISC","NQY","LBA","LHR","LBA","NQY","LBA","SOU",
    "LCY","MAN","LCY","NCL","LGW","NCL","LGW","NQY","LHR","MAN","LHR","NCL",
    "LHR","NQY","LPL","NQY","MAN","NQY","MAN","NWI","MAN","SEN","MAN","SOU",
    "MME","NWI","NCL","SOU","NQY","SEN"
    ]"#
        )
        .unwrap()
    );
    dbg!(routes_within_england_time_no_dup.elapsed());
}

#[test]
fn hard_route_finding() {
    initialize(&TEST_DB);
    let hard_route_finding = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        reachable[to, choice(p)] := *route{fr: 'AUS', to}, to != 'YYZ', p = ['AUS', to];
        reachable[to, choice(p)] := reachable[b, prev], *route{fr: b, to},
                                    to != 'YYZ', p = append(prev, to)
        ?[p] := reachable['YPO', p]

        :limit 1;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"["AUS","YYC","YQT","YTS","YMO","YFA","ZKE","YAT","YPO"]"#
        )
        .unwrap()
    );
    dbg!(hard_route_finding.elapsed());
}

#[test]
fn na_from_india() {
    initialize(&TEST_DB);
    let na_from_india = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[ind_a, na_a] := *airport{code: ind_a, country: 'IN'},
                          *route{fr: ind_a, to: na_a},
                          *airport{code: na_a, country},
                          country in ['US', 'CA']
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
     ["BOM","EWR"],["BOM","JFK"],["BOM","YYZ"],["DEL","EWR"],["DEL","IAD"],["DEL","JFK"],
     ["DEL","ORD"],["DEL","SFO"],["DEL","YVR"],["DEL","YYZ"]
     ]"#
        )
        .unwrap()
    );
    dbg!(na_from_india.elapsed());
}

#[test]
fn eu_cities_reachable_from_fll() {
    initialize(&TEST_DB);
    let eu_cities_reachable_from_fll = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[city] := *route{fr: 'FLL', to}, *contain['EU', to], *airport{code: to, city}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
     ["Barcelona"],["Copenhagen"],["London"],["Madrid"],["Oslo"],["Paris"],["Stockholm"]
     ]"#
        )
        .unwrap()
    );
    dbg!(eu_cities_reachable_from_fll.elapsed());
}

#[test]
fn clt_to_eu_or_sa() {
    initialize(&TEST_DB);
    let clt_to_eu_or_sa = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[to] := *route{fr: 'CLT', to}, c_name in ['EU', 'SA'], *contain[c_name, to]
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
     ["BCN"],["CDG"],["DUB"],["FCO"],["FRA"],["GIG"],["GRU"],["LHR"],["MAD"],["MUC"]
     ]"#
        )
        .unwrap()
    );
    dbg!(clt_to_eu_or_sa.elapsed());
}

#[test]
fn london_to_us() {
    initialize(&TEST_DB);
    let london_to_us = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[fr, to] := fr in ['LHR', 'LCY', 'LGW', 'LTN', 'STN'],
                     *route{fr, to}, *airport{code: to, country: 'US'}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
     ["LGW","AUS"],["LGW","BOS"],["LGW","DEN"],["LGW","FLL"],["LGW","JFK"],["LGW","LAS"],
     ["LGW","LAX"],["LGW","MCO"],["LGW","MIA"],["LGW","OAK"],["LGW","ORD"],["LGW","SEA"],
     ["LGW","SFO"],["LGW","TPA"],["LHR","ATL"],["LHR","AUS"],["LHR","BNA"],["LHR","BOS"],
     ["LHR","BWI"],["LHR","CHS"],["LHR","CLT"],["LHR","DEN"],["LHR","DFW"],["LHR","DTW"],
     ["LHR","EWR"],["LHR","IAD"],["LHR","IAH"],["LHR","JFK"],["LHR","LAS"],["LHR","LAX"],
     ["LHR","MIA"],["LHR","MSP"],["LHR","MSY"],["LHR","ORD"],["LHR","PDX"],["LHR","PHL"],
     ["LHR","PHX"],["LHR","PIT"],["LHR","RDU"],["LHR","SAN"],["LHR","SEA"],["LHR","SFO"],
     ["LHR","SJC"],["LHR","SLC"],["STN","BOS"],["STN","EWR"],["STN","IAD"],["STN","SFB"]
     ]"#
        )
        .unwrap()
    );
    dbg!(london_to_us.elapsed());
}

#[test]
fn tx_to_ny() {
    initialize(&TEST_DB);
    let tx_to_ny = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[fr, to] := *airport{code: fr, region: 'US-TX'},
                     *route{fr, to}, *airport{code: to, region: 'US-NY'}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        ["AUS","BUF"],["AUS","EWR"],["AUS","JFK"],["DAL","LGA"],["DFW","BUF"],["DFW","EWR"],
        ["DFW","JFK"],["DFW","LGA"],["HOU","EWR"],["HOU","JFK"],["HOU","LGA"],["IAH","EWR"],
        ["IAH","JFK"],["IAH","LGA"],["SAT","EWR"],["SAT","JFK"]
     ]"#
        )
        .unwrap()
    );
    dbg!(tx_to_ny.elapsed());
}

#[test]
fn denver_to_mexico() {
    initialize(&TEST_DB);
    let denver_to_mexico = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[city] := *route{fr: 'DEN', to}, *airport{code: to, country: 'MX', city}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        ["Cancun"],["Cozumel"],["Guadalajara"],["Mexico City"],["Monterrey"],
        ["Puerto Vallarta"],["San Josรฉ del Cabo"]
        ]"#
        )
        .unwrap()
    );
    dbg!(denver_to_mexico.elapsed());
}

#[test]
fn three_cities() {
    initialize(&TEST_DB);
    let three_cities = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        three[code] := city in ['London', 'Munich', 'Paris'], *airport{code, city}
        ?[s, d] := three[s], *route{fr: s, to: d}, three[d]
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        ["CDG","LCY"],["CDG","LGW"],["CDG","LHR"],["CDG","LTN"],["CDG","MUC"],["LCY","CDG"],
        ["LCY","MUC"],["LCY","ORY"],["LGW","CDG"],["LGW","MUC"],["LHR","CDG"],["LHR","MUC"],
        ["LHR","ORY"],["LTN","CDG"],["LTN","MUC"],["LTN","ORY"],["MUC","CDG"],["MUC","LCY"],
        ["MUC","LGW"],["MUC","LHR"],["MUC","LTN"],["MUC","ORY"],["MUC","STN"],["ORY","LCY"],
        ["ORY","LHR"],["ORY","MUC"],["STN","MUC"]
        ]"#
        )
        .unwrap()
    );
    dbg!(three_cities.elapsed());
}

#[test]
fn long_distance_from_lgw() {
    initialize(&TEST_DB);
    let long_distance_from_lgw = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[city, dist] := *route{fr: 'LGW', to, dist},
                         dist > 4000, *airport{code: to, city}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["Austin",4921.0],["Beijing",5070.0],["Bridgetown",4197.0],["Buenos Aires",6908.0],["Calgary",4380.0],
            ["Cancun",4953.0],["Cape Town",5987.0],["Chengdu",5156.0],["Chongqing",5303.0],["Colombo",5399.0],
            ["Denver",4678.0],["Duong Dong",6264.0],["Fort Lauderdale",4410.0],["Havana",4662.0],["Hong Kong",5982.0],
            ["Kigali",4077.0],["Kingston",4680.0],["Langkawi",6299.0],["Las Vegas",5236.0],["Los Angeles",5463.0],
            ["Malรฉ",5287.0],["Miami",4429.0],["Montego Bay",4699.0],["Oakland",5364.0],["Orlando",4341.0],
            ["Port Louis",6053.0],["Port of Spain",4408.0],["Punta Cana",4283.0],["Rayong",6008.0],
            ["Rio de Janeiro",5736.0],["San Francisco",5374.0],["San Jose",5419.0],["Seattle",4807.0],
            ["Shanghai",5745.0],["Singapore",6751.0],["St. George",4076.0],["Taipei",6080.0],["Tampa",4416.0],
            ["Tianjin",5147.0],["Vancouver",4731.0],["Varadero",4618.0],["Vieux Fort",4222.0]
        ]"#
        )
            .unwrap()
    );
    dbg!(long_distance_from_lgw.elapsed());
}

#[test]
fn long_routes_one_dir() {
    initialize(&TEST_DB);
    let long_routes_one_dir = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[fr, dist, to] := *route{fr, to, dist}, dist > 8000, fr < to;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        ["AKL",8186.0,"ORD"],["AKL",8818.0,"DXB"],["AKL",9025.0,"DOH"],["ATL",8434.0,"JNB"],
        ["AUH",8053.0,"DFW"],["AUH",8139.0,"SFO"],["AUH",8372.0,"LAX"],["CAN",8754.0,"MEX"],
        ["DFW",8022.0,"DXB"],["DFW",8105.0,"HKG"],["DFW",8574.0,"SYD"],["DOH",8030.0,"IAH"],
        ["DOH",8287.0,"LAX"],["DXB",8085.0,"SFO"],["DXB",8150.0,"IAH"],["DXB",8321.0,"LAX"],
        ["EWR",8047.0,"HKG"],["EWR",9523.0,"SIN"],["HKG",8054.0,"JFK"],["HKG",8135.0,"IAD"],
        ["IAH",8591.0,"SYD"],["JED",8314.0,"LAX"],["JFK",8504.0,"MNL"],["JFK",9526.0,"SIN"],
        ["LAX",8246.0,"RUH"],["LAX",8756.0,"SIN"],["LHR",9009.0,"PER"],["MEL",8197.0,"YVR"],
        ["PEK",8884.0,"PTY"],["SCL",8208.0,"TLV"],["SEA",8059.0,"SIN"],["SFO",8433.0,"SIN"]]"#
        )
        .unwrap()
    );
    dbg!(long_routes_one_dir.elapsed());
}

#[test]
fn longest_routes() {
    initialize(&TEST_DB);
    let longest_routes = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[fr, dist, to] := *route{fr, to, dist}, dist > 4000, fr < to;
        :sort -dist;
        :limit 20;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["JFK",9526.0,"SIN"],["EWR",9523.0,"SIN"],["AKL",9025.0,"DOH"],["LHR",9009.0,"PER"],
            ["PEK",8884.0,"PTY"],["AKL",8818.0,"DXB"],["LAX",8756.0,"SIN"],["CAN",8754.0,"MEX"],
            ["IAH",8591.0,"SYD"],["DFW",8574.0,"SYD"],["JFK",8504.0,"MNL"],["ATL",8434.0,"JNB"],
            ["SFO",8433.0,"SIN"],["AUH",8372.0,"LAX"],["DXB",8321.0,"LAX"],["JED",8314.0,"LAX"],
            ["DOH",8287.0,"LAX"],["LAX",8246.0,"RUH"],["SCL",8208.0,"TLV"],["MEL",8197.0,"YVR"]]"#
        )
        .unwrap()
    );
    dbg!(longest_routes.elapsed());
}

#[test]
fn longest_routes_from_each_airports() {
    initialize(&TEST_DB);
    let longest_routes_from_each_airports = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[fr, max(dist), choice(to)] := *route{fr, dist, to}
        :limit 10;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        ["AAA",968.0,"FAC"],["AAE",1161.0,"ALG"],["AAL",1693.0,"AAR"],["AAN",1613.0,"CAI"],
        ["AAQ",2122.0,"BAX"],["AAR",1585.0,"AAL"],["AAT",267.0,"URC"],["AAX",69.0,"POJ"],
        ["AAY",531.0,"SAH"],["ABA",2096.0,"DME"]]"#
        )
        .unwrap()
    );
    dbg!(longest_routes_from_each_airports.elapsed());
}

#[test]
fn total_distance_from_three_cities() {
    initialize(&TEST_DB);
    let total_distance_from_three_cities = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        three[code] := city in ['London', 'Munich', 'Paris'], *airport{code, city}
        ?[sum(dist)] := three[a], *route{fr: a, dist}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"2739039.0"#).unwrap()
    );
    dbg!(total_distance_from_three_cities.elapsed());
}

#[test]
fn total_distance_within_three_cities() {
    initialize(&TEST_DB);
    let total_distance_within_three_cities = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        three[code] := city in ['London', 'Munich', 'Paris'], *airport{code, city}
        ?[sum(dist)] := three[a], *route{fr: a, dist, to}, three[to]
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"10282.0"#).unwrap()
    );
    dbg!(total_distance_within_three_cities.elapsed());
}

#[test]
fn specific_distance() {
    initialize(&TEST_DB);
    let specific_distance = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[dist] := *route{fr: 'AUS', to: 'MEX', dist}
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"748.0"#).unwrap()
    );
    dbg!(specific_distance.elapsed());
}

#[test]
fn n_routes_between() {
    initialize(&TEST_DB);
    let n_routes_between = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        us_a[a] := *contain['US', a]
        ?[count(fr)] := *route{fr, to, dist}, dist >= 100, dist <= 200,
                        us_a[fr], us_a[to]
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"597"#).unwrap()
    );
    dbg!(n_routes_between.elapsed());
}

#[test]
fn one_stop_distance() {
    initialize(&TEST_DB);
    let one_stop_distance = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[code, dist] := *route{fr: 'AUS', to: code, dist: dis1},
                         *route{fr: code, to: 'LHR', dist: dis2},
                         dist = dis1 + dis2
        :order dist;
        :limit 10;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["DTW",4893.0],["YYZ",4901.0],["ORD",4912.0],["PIT",4916.0],["BNA",4923.0],
            ["DFW",4926.0],["BOS",4944.0],["EWR",4953.0],["IAD",4959.0],["JFK",4960.0]]"#
        )
        .unwrap()
    );
    dbg!(one_stop_distance.elapsed());
}

#[test]
fn airport_most_routes() {
    initialize(&TEST_DB);
    let airport_most_routes = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[fr, count(fr)] := *route{fr}
        :order -count(fr);
        :limit 10;
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            ["FRA",310],["IST",309],["CDG",293],["AMS",283],["MUC",270],["ORD",265],["DFW",253],
            ["DXB",248],["PEK",248],["ATL",242]]"#
        )
        .unwrap()
    );
    dbg!(airport_most_routes.elapsed());
}

#[test]
fn north_of_77() {
    initialize(&TEST_DB);
    let north_of_77 = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[city, latitude] := *airport{lat, city}, lat > 77, latitude = round(lat)
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"[["Longyearbyen",78.0],["Qaanaaq",77.0]]"#).unwrap()
    );
    dbg!(north_of_77.elapsed());
}

#[test]
fn greenwich_meridian() {
    initialize(&TEST_DB);
    let greenwich_meridian = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[code] := *airport{lon, code}, lon > -0.1, lon < 0.1
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"[["CDT"], ["LCY"], ["LDE"], ["LEH"]]"#).unwrap()
    );
    dbg!(greenwich_meridian.elapsed());
}

#[test]
fn box_around_heathrow() {
    initialize(&TEST_DB);
    let box_around_heathrow = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        h_box[lon, lat] := *airport{code: 'LHR', lon, lat}
        ?[code] := h_box[lhr_lon, lhr_lat], *airport{code, lon, lat},
                    abs(lhr_lon - lon) < 1, abs(lhr_lat - lat) < 1
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"[["LCY"], ["LGW"], ["LHR"], ["LTN"], ["SOU"], ["STN"]]"#)
            .unwrap()
    );
    dbg!(box_around_heathrow.elapsed());
}

#[test]
fn dfw_by_region() {
    initialize(&TEST_DB);
    let dfw_by_region = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[region, collect(to)] := *route{fr: 'DFW', to},
                                  *airport{code: to, country: 'US', region},
                                  region in ['US-CA', 'US-TX', 'US-FL', 'US-CO', 'US-IL']
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"
    "US-CA",["BFL","BUR","FAT","LAX","MRY","OAK","ONT","PSP","SAN","SBA","SFO","SJC","SMF","SNA",
    ["US-CO",["ASE","COS","DEN","DRO","EGE","GJT","GUC","HDN","MTJ"]],
    ["US-FL",["ECP","EYW","FLL","GNV","JAX","MCO","MIA","PBI","PNS","RSW","SRQ","TLH","TPA","VPS"]],
    ["US-IL",["BMI","CMI","MLI","ORD","PIA","SPI"]],
    ["US-TX",["ABI","ACT","AMA","AUS","BPT","BRO","CLL","CRP","DRT","ELP","GGG","GRK","HOU","HRL",
              "IAH","LBB","LRD","MAF","MFE","SAT","SJT","SPS","TYR"]]]
        "#)
            .unwrap()
    );
    dbg!(dfw_by_region.elapsed());
}

#[test]
fn great_circle_distance() {
    initialize(&TEST_DB);
    let great_circle_distance = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        ?[deg_diff] := *airport{code: 'SFO', lat: a_lat, lon: a_lon},
                       *airport{code: 'NRT', lat: b_lat, lon: b_lon},
                        deg_diff = round(haversine_deg_input(a_lat, a_lon, b_lat, b_lon));
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"1.0"#).unwrap()
    );
    dbg!(great_circle_distance.elapsed());
}

#[test]
fn aus_to_edi() {
    initialize(&TEST_DB);
    let aus_to_edi = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        us_uk_airports[code] := *airport{code, country: 'UK'}
        us_uk_airports[code] := *airport{code, country: 'US'}
        routes[to, shortest(path)] := *route{fr: 'AUS', to}, us_uk_airports[to],
                                        path = ['AUS', to];
        routes[to, shortest(path)] := routes[a, prev], *route{fr: a, to},
                                        us_uk_airports[to],
                                        path = append(prev, to);
        ?[path] := routes['EDI', path];
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(r#"["AUS", "BOS", "EDI"]"#).unwrap()
    );
    dbg!(aus_to_edi.elapsed());
}

#[test]
fn reachable_from_lhr() {
    initialize(&TEST_DB);
    let reachable_from_lhr = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        routes[to, shortest(path)] := *route{fr: 'LHR', to},
                                      path = ['LHR', to];
        routes[to, shortest(path)] := routes[a, prev], *route{fr: a, to},
                                      path = append(prev, to);
        ?[len, path] := routes[_, path], len = length(path);

        :order -len;
        :limit 10;
        "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
            [8,["LHR","YYZ","YTS","YMO","YFA","ZKE","YAT","YPO"]],
            [7,["LHR","AUH","BNE","ISA","BQL","BEU","BVI"]],
            [7,["LHR","AUH","BNE","WTB","SGO","CMA","XTG"]],
            [7,["LHR","CAN","ADL","AYQ","MEB","WMB","PTJ"]],
            [7,["LHR","DEN","ANC","AKN","PIP","UGB","PTH"]],
            [7,["LHR","DEN","ANC","ANI","CHU","CKD","RDV"]],
            [7,["LHR","DEN","ANC","ANI","CHU","CKD","SLQ"]],
            [7,["LHR","DEN","ANC","BET","NME","TNK","WWT"]],
            [7,["LHR","KEF","GOH","JAV","JUV","NAQ","THU"]],
            [7,["LHR","YUL","YGL","YPX","AKV","YIK","YZG"]]]
        "#
        )
        .unwrap()
    );
    dbg!(reachable_from_lhr.elapsed());
}

#[test]
fn furthest_from_lhr() {
    initialize(&TEST_DB);
    let furthest_from_lhr = Instant::now();

    let rows = TEST_DB
        .run_default(
            r#"
        routes[to, min_cost(cost_pair)] := *route{fr: 'LHR', to, dist},
                                             path = ['LHR', to],
                                             cost_pair = [path, dist];
        routes[to, min_cost(cost_pair)] := routes[a, prev], *route{fr: a, to, dist},
                                           path = append(first(prev), to),
                                           cost_pair = [path, last(prev) + dist];
        ?[cost, path] := routes[dst, cost_pair], cost = last(cost_pair), path = first(cost_pair);

        :order -cost;
        :limit 10;
        "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(
        rows["rows"],
        serde_json::Value::from_str(
            r#"[
        [12922.0,["LHR","JNB","HLE","ASI","BZZ"]],[12093.0,["LHR","PVG","CHC","IVC"]],
        [12015.0,["LHR","NRT","AKL","WLG","TIU"]],[12009.0,["LHR","PVG","CHC","DUD"]],
        [11910.0,["LHR","NRT","AKL","WLG","WSZ"]],[11900.0,["LHR","PVG","CHC","HKK"]],
        [11805.0,["LHR","PVG","CHC"]],[11766.0,["LHR","PVG","BNE","ZQN"]],
        [11758.0,["LHR","NRT","AKL","BHE"]],[11751.0,["LHR","NRT","AKL","NSN"]]]
        "#
        )
        .unwrap()
    );
    dbg!(furthest_from_lhr.elapsed());
}

#[test]
fn skip_limit() {
    initialize(&TEST_DB);
    let rows = TEST_DB
        .run_default(
            r#"
        ?[a] := a in [9, 9, 8, 9, 8, 7, 7, 6, 5, 9, 4, 4, 3]
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[3], [4], [5], [6], [7], [8], [9]]));

    let rows = TEST_DB
        .run_default(
            r#"
        ?[a] := a in [9, 9, 8, 9, 8, 7, 7, 6, 5, 9, 4, 4, 3]
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[3], [4], [5], [6], [7], [8], [9]]));

    let rows = TEST_DB
        .run_default(
            r#"
        ?[a] := a in [9, 9, 8, 9, 8, 7, 7, 6, 5, 9, 4, 4, 3]
        :limit 2
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[8], [9]]));

    let rows = TEST_DB
        .run_default(
            r#"
        ?[a] := a in [9, 9, 8, 9, 8, 7, 7, 6, 5, 9, 4, 4, 3]
        :limit 2
        :offset 1
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[7], [8]]));

    let rows = TEST_DB
        .run_default(
            r#"
        ?[a] := a in [9, 9, 8, 9, 8, 7, 7, 6, 5, 9, 4, 4, 3]
        :limit 100
        :offset 1
    "#,
        )
        .unwrap()
        .into_json();

    assert_eq!(rows["rows"], json!([[3], [4], [5], [6], [7], [8]]));
}

Graph