satellite-constellation-test.cc

Go to the documentation of this file.

/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */
/*
 * Copyright (c) 2014 Magister Solutions Ltd
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License version 2 as
 * published by the Free Software Foundation;
 *
 * 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, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 * Author: Bastien Tauran <bastien.tauran@viveris.fr>
 */

 
#include "ns3/cbr-application.h"
#include "ns3/cbr-helper.h"
#include "ns3/config.h"
#include "ns3/log.h"
#include "ns3/packet-sink-helper.h"
#include "ns3/packet-sink.h"
#include "ns3/satellite-env-variables.h"
#include "ns3/satellite-gw-mac.h"
#include "ns3/satellite-helper.h"
#include "ns3/satellite-id-mapper.h"
#include "ns3/satellite-isl-arbiter-unicast.h"
#include "ns3/satellite-isl-arbiter.h"
#include "ns3/satellite-orbiter-feeder-phy.h"
#include "ns3/satellite-orbiter-net-device.h"
#include "ns3/satellite-orbiter-user-phy.h"
#include "ns3/satellite-phy-rx-carrier.h"
#include "ns3/satellite-topology.h"
#include "ns3/satellite-ut-mac-state.h"
#include "ns3/simulation-helper.h"
#include "ns3/simulator.h"
#include "ns3/singleton.h"
#include "ns3/string.h"
#include "ns3/test.h"
 
#include <fstream>
#include <iostream>
 
using namespace ns3;

class SatConstellationTest1 : public TestCase
{
  public:
    SatConstellationTest1();
    virtual ~SatConstellationTest1();
 
  private:
    virtual void DoRun(void);
};
 
// Add some help text to this case to describe what it is intended to test
SatConstellationTest1::SatConstellationTest1()
    : TestCase("This case tests that a topology is correctly loaded.")
{
}
 
// This destructor does nothing but we include it as a reminder that
// the test case should clean up after itself
SatConstellationTest1::~SatConstellationTest1()
{
}
 
//
// SatConstellationTest1 TestCase implementation
//
void
SatConstellationTest1::DoRun(void)
{
    Config::Reset();
 
    // Set simulation output details
    SatEnvVariables::GetInstance()->DoInitialize();
    SatEnvVariables::GetInstance()->SetOutputVariables("test-sat-constellation", "test1", true);

    Config::SetDefault("ns3::SatConf::ForwardLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));
    Config::SetDefault("ns3::SatConf::ReturnLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));

    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionEachRequest", BooleanValue(false));
    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionPeriod",
                       TimeValue(MilliSeconds(10)));
    Config::SetDefault("ns3::SatHelper::GwUsers", UintegerValue(3));
 
    Ptr<SimulationHelper> simulationHelper =
        CreateObject<SimulationHelper>("test-sat-constellation");
 
    // Creating the reference system.
    simulationHelper->SetBeamSet({43, 30});
    simulationHelper->SetUserCountPerUt(5);
 
    simulationHelper->LoadScenario("constellation-eutelsat-geo-2-sats-no-isls");
 
    simulationHelper->CreateSatScenario();
 
    NodeContainer sats = Singleton<SatTopology>::Get()->GetOrbiterNodes();
    NodeContainer gws = Singleton<SatTopology>::Get()->GetGwNodes();
    NodeContainer uts = Singleton<SatTopology>::Get()->GetUtNodes();
    NodeContainer gwUsers = Singleton<SatTopology>::Get()->GetGwUserNodes();
    NodeContainer utUsers = Singleton<SatTopology>::Get()->GetUtUserNodes();
 
    uint32_t countNetDevices = 0;
    for (uint32_t i = 0; i < sats.GetN(); i++)
    {
        for (uint32_t j = 0; j < sats.Get(i)->GetNDevices(); j++)
        {
            if (DynamicCast<SatOrbiterNetDevice>(sats.Get(i)->GetDevice(j)) != nullptr)
            {
                countNetDevices += 1;
            }
        }
    }
 
    NS_TEST_ASSERT_MSG_EQ(sats.GetN(), 2, "Topology must contain 2 satellites");
    NS_TEST_ASSERT_MSG_EQ(countNetDevices,
                          2,
                          "Topology must contain 2 satellite Orbiter Net Devices");
    NS_TEST_ASSERT_MSG_EQ(gws.GetN(), 2, "Topology must contain 2 GWs");
    NS_TEST_ASSERT_MSG_EQ(uts.GetN(), 3, "Topology must contain 3 UTs");
    NS_TEST_ASSERT_MSG_EQ(gwUsers.GetN(), 3, "Topology must contain 3 GW users");
    NS_TEST_ASSERT_MSG_EQ(utUsers.GetN(), 15, "Topology must contain 15 UT users");
 
    GeoCoordinate sat1 = GeoCoordinate(sats.Get(0)->GetObject<SatMobilityModel>()->GetPosition());
    GeoCoordinate sat2 = GeoCoordinate(sats.Get(1)->GetObject<SatMobilityModel>()->GetPosition());
 
    GeoCoordinate gw1 = GeoCoordinate(gws.Get(0)->GetObject<SatMobilityModel>()->GetPosition());
    GeoCoordinate gw2 = GeoCoordinate(gws.Get(1)->GetObject<SatMobilityModel>()->GetPosition());
 
    GeoCoordinate ut1 = GeoCoordinate(uts.Get(0)->GetObject<SatMobilityModel>()->GetPosition());
    GeoCoordinate ut2 = GeoCoordinate(uts.Get(1)->GetObject<SatMobilityModel>()->GetPosition());
    GeoCoordinate ut3 = GeoCoordinate(uts.Get(2)->GetObject<SatMobilityModel>()->GetPosition());
 
    NS_TEST_ASSERT_MSG_EQ_TOL(sat1.GetLatitude(),
                              0.013043,
                              0.001,
                              "Incorrect latitude for satellite 1");
    NS_TEST_ASSERT_MSG_EQ_TOL(sat1.GetLongitude(),
                              7.20835,
                              0.001,
                              "Incorrect longitude for satellite 1");
    NS_TEST_ASSERT_MSG_EQ_TOL(sat1.GetAltitude(),
                              35786400,
                              100,
                              "Incorrect altitude for satellite 1");
 
    NS_TEST_ASSERT_MSG_EQ_TOL(sat2.GetLatitude(),
                              -0.0219759,
                              0.001,
                              "Incorrect latitude for satellite 2");
    NS_TEST_ASSERT_MSG_EQ_TOL(sat2.GetLongitude(),
                              47.7047,
                              0.001,
                              "Incorrect longitude for satellite 2");
    NS_TEST_ASSERT_MSG_EQ_TOL(sat2.GetAltitude(),
                              35779000,
                              100,
                              "Incorrect altitude for satellite 2");
 
    NS_TEST_ASSERT_MSG_EQ_TOL(gw1.GetLatitude(), 48.85, 0.001, "Incorrect latitude for GW 1");
    NS_TEST_ASSERT_MSG_EQ_TOL(gw1.GetLongitude(), 2.34, 0.001, "Incorrect longitude for GW 1");
    NS_TEST_ASSERT_MSG_EQ(gw1.GetAltitude(), 0, "Incorrect altitude for GW 1");
 
    NS_TEST_ASSERT_MSG_EQ_TOL(gw2.GetLatitude(), 55.75, 0.001, "Incorrect latitude for GW 2");
    NS_TEST_ASSERT_MSG_EQ_TOL(gw2.GetLongitude(), 37.62, 0.001, "Incorrect longitude for GW 2");
    NS_TEST_ASSERT_MSG_EQ(gw2.GetAltitude(), 0, "Incorrect altitude for GW 2");
 
    NS_TEST_ASSERT_MSG_EQ_TOL(ut1.GetLatitude(), 48.8534, 0.001, "Incorrect latitude for UT 1");
    NS_TEST_ASSERT_MSG_EQ_TOL(ut1.GetLongitude(), 2.3488, 0.001, "Incorrect longitude for UT 1");
    NS_TEST_ASSERT_MSG_EQ(ut1.GetAltitude(), 0, "Incorrect altitude for UT 1");
 
    NS_TEST_ASSERT_MSG_EQ_TOL(ut2.GetLatitude(), 55.755, 0.001, "Incorrect latitude for UT 2");
    NS_TEST_ASSERT_MSG_EQ_TOL(ut2.GetLongitude(), 37.6218, 0.001, "Incorrect longitude for UT 2");
    NS_TEST_ASSERT_MSG_EQ(ut2.GetAltitude(), 0, "Incorrect altitude for UT 2");
 
    NS_TEST_ASSERT_MSG_EQ_TOL(ut3.GetLatitude(), 55.754, 0.001, "Incorrect latitude for UT 3");
    NS_TEST_ASSERT_MSG_EQ_TOL(ut3.GetLongitude(), 37.621, 0.001, "Incorrect longitude for UT 3");
    NS_TEST_ASSERT_MSG_EQ(ut3.GetAltitude(), 0, "Incorrect altitude for UT 3");
 
    Simulator::Destroy();
}

class SatConstellationTest2 : public TestCase
{
  public:
    SatConstellationTest2();
    virtual ~SatConstellationTest2();
 
  private:
    virtual void DoRun(void);
 
    std::vector<std::string> Split(std::string s, char del);
    void TestFileValue(std::string path, uint32_t time, uint32_t expectedValue);
};
 
// Add some help text to this case to describe what it is intended to test
SatConstellationTest2::SatConstellationTest2()
    : TestCase("This case tests that a application throughputs PerEntity are correct.")
{
}
 
// This destructor does nothing but we include it as a reminder that
// the test case should clean up after itself
SatConstellationTest2::~SatConstellationTest2()
{
}
 
std::vector<std::string>
SatConstellationTest2::Split(std::string s, char del)
{
    std::stringstream ss(s);
    std::string word;
    std::vector<std::string> tokens;
    while (!ss.eof())
    {
        std::getline(ss, word, del);
        tokens.push_back(word);
    }
    return tokens;
}
 
void
SatConstellationTest2::TestFileValue(std::string path, uint32_t time, uint32_t expectedValue)
{
    std::string line;
    std::ifstream myfile(path);
    std::string delimiter = " ";
    std::string token;
    bool valueFound;
 
    if (myfile.is_open())
    {
        valueFound = false;
        while (std::getline(myfile, line))
        {
            std::vector<std::string> tokens = Split(line, ' ');
            if (tokens.size() == 2)
            {
                if ((uint32_t)std::stoi(tokens[0]) == time)
                {
                    NS_TEST_ASSERT_MSG_EQ_TOL(std::stof(tokens[1]),
                                              expectedValue,
                                              expectedValue / 10,
                                              "Incorrect throughput for statistic " << path);
                    valueFound = true;
                    break;
                }
            }
        }
        myfile.close();
        if (!valueFound)
        {
            NS_TEST_ASSERT_MSG_EQ(0, 1, "Cannot find time " << time << " for trace file " << path);
        }
    }
    else
    {
        NS_TEST_ASSERT_MSG_EQ(0, 1, "Cannot read trace file " << path);
    }
}
 
//
// SatConstellationTest2 TestCase implementation
//
void
SatConstellationTest2::DoRun(void)
{
    Config::Reset();
    Singleton<SatIdMapper>::Get()->Reset();
 
    // Set simulation output details
    SatEnvVariables::GetInstance()->DoInitialize();
    SatEnvVariables::GetInstance()->SetOutputVariables("test-sat-constellation", "test2", true);

    Config::SetDefault("ns3::SatConf::ForwardLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));
    Config::SetDefault("ns3::SatConf::ReturnLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));

    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionEachRequest", BooleanValue(false));
    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionPeriod",
                       TimeValue(MilliSeconds(10)));
    Config::SetDefault("ns3::SatHelper::GwUsers", UintegerValue(3));
 
    Ptr<SimulationHelper> simulationHelper =
        CreateObject<SimulationHelper>("test-sat-constellation/test2");
 
    // Creating the reference system.
    simulationHelper->SetBeamSet({43, 30});
    simulationHelper->SetUserCountPerUt(5);
 
    simulationHelper->LoadScenario("constellation-eutelsat-geo-2-sats-no-isls");
 
    simulationHelper->CreateSatScenario();
 
    simulationHelper->GetTrafficHelper()->AddCbrTraffic(
        SatTrafficHelper::FWD_LINK,
        SatTrafficHelper::UDP,
        MilliSeconds(10),
        512,
        NodeContainer(Singleton<SatTopology>::Get()->GetGwUserNode(0)),
        Singleton<SatTopology>::Get()->GetUtUserNodes(),
        Seconds(1.0),
        Seconds(29.0),
        Seconds(0));
 
    simulationHelper->GetTrafficHelper()->AddCbrTraffic(
        SatTrafficHelper::RTN_LINK,
        SatTrafficHelper::UDP,
        MilliSeconds(10),
        512,
        NodeContainer(Singleton<SatTopology>::Get()->GetGwUserNode(0)),
        Singleton<SatTopology>::Get()->GetUtUserNodes(),
        Seconds(1.0),
        Seconds(29.0),
        Seconds(0));
 
    Ptr<SatStatsHelperContainer> s = simulationHelper->GetStatisticsContainer();
 
    s->AddGlobalFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddGlobalRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerGwFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerGwRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerSatFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerSatRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerBeamFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerBeamRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerUtFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerUtRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
 
    simulationHelper->SetSimulationTime(Seconds(10));
    simulationHelper->RunSimulation();
 
    Simulator::Destroy();
 
    // Global throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-global-fwd-app-throughput-scatter-0.txt",
                  5,
                  6144);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-global-rtn-app-throughput-scatter-0.txt",
                  5,
                  6144);
 
    // Per satellite throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-sat-fwd-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-sat-fwd-app-throughput-scatter-2.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-sat-rtn-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-sat-rtn-app-throughput-scatter-2.txt",
                  5,
                  4096);
 
    // Per beam throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-fwd-app-throughput-scatter-1-30.txt",
                  5,
                  0);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-fwd-app-throughput-scatter-1-43.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-fwd-app-throughput-scatter-2-30.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-fwd-app-throughput-scatter-2-43.txt",
                  5,
                  0);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-rtn-app-throughput-scatter-1-30.txt",
                  5,
                  0);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-rtn-app-throughput-scatter-1-43.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-rtn-app-throughput-scatter-2-30.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-beam-rtn-app-throughput-scatter-2-43.txt",
                  5,
                  0);
 
    // Per GW throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-gw-fwd-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-gw-fwd-app-throughput-scatter-2.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-gw-rtn-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-gw-rtn-app-throughput-scatter-2.txt",
                  5,
                  4096);
 
    // Per UT throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-ut-fwd-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-ut-fwd-app-throughput-scatter-2.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-ut-fwd-app-throughput-scatter-3.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-ut-rtn-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-ut-rtn-app-throughput-scatter-2.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test2/"
                  "stat-per-ut-rtn-app-throughput-scatter-3.txt",
                  5,
                  2048);
}

class SatConstellationTest3 : public TestCase
{
  public:
    SatConstellationTest3();
    virtual ~SatConstellationTest3();
 
  private:
    virtual void DoRun(void);
 
    std::vector<std::string> Split(std::string s, char del);
    void TestFileValue(std::string path, uint32_t time, uint32_t expectedValue);
};
 
// Add some help text to this case to describe what it is intended to test
SatConstellationTest3::SatConstellationTest3()
    : TestCase("This case tests that a application throughputs PerEntity are correct.")
{
}
 
// This destructor does nothing but we include it as a reminder that
// the test case should clean up after itself
SatConstellationTest3::~SatConstellationTest3()
{
}
 
std::vector<std::string>
SatConstellationTest3::Split(std::string s, char del)
{
    std::stringstream ss(s);
    std::string word;
    std::vector<std::string> tokens;
    while (!ss.eof())
    {
        std::getline(ss, word, del);
        tokens.push_back(word);
    }
    return tokens;
}
 
void
SatConstellationTest3::TestFileValue(std::string path, uint32_t time, uint32_t expectedValue)
{
    std::string line;
    std::ifstream myfile(path);
    std::string delimiter = " ";
    std::string token;
    bool valueFound;
 
    if (myfile.is_open())
    {
        valueFound = false;
        while (std::getline(myfile, line))
        {
            std::vector<std::string> tokens = Split(line, ' ');
            if (tokens.size() == 2)
            {
                if ((uint32_t)std::stoi(tokens[0]) == time)
                {
                    NS_TEST_ASSERT_MSG_EQ_TOL(std::stof(tokens[1]),
                                              expectedValue,
                                              expectedValue / 10,
                                              "Incorrect throughput for statistic " << path);
                    valueFound = true;
                    break;
                }
            }
        }
        myfile.close();
        if (!valueFound)
        {
            NS_TEST_ASSERT_MSG_EQ(0, 1, "Cannot find time " << time << " for trace file " << path);
        }
    }
    else
    {
        NS_TEST_ASSERT_MSG_EQ(0, 1, "Cannot read trace file " << path);
    }
}
 
//
// SatConstellationTest3 TestCase implementation
//
void
SatConstellationTest3::DoRun(void)
{
    Config::Reset();
    Singleton<SatIdMapper>::Get()->Reset();
 
    // Set simulation output details
    SatEnvVariables::GetInstance()->DoInitialize();
    SatEnvVariables::GetInstance()->SetOutputVariables("test-sat-constellation", "test3", true);

    Config::SetDefault("ns3::SatConf::ForwardLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));
    Config::SetDefault("ns3::SatConf::ReturnLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));

    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionEachRequest", BooleanValue(false));
    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionPeriod",
                       TimeValue(MilliSeconds(10)));
    Config::SetDefault("ns3::SatHelper::GwUsers", UintegerValue(3));
 
    Ptr<SimulationHelper> simulationHelper =
        CreateObject<SimulationHelper>("test-sat-constellation/test3");
 
    // Creating the reference system.
    simulationHelper->SetBeamSet({43, 30});
    simulationHelper->SetUserCountPerUt(5);
 
    simulationHelper->LoadScenario("constellation-eutelsat-geo-2-sats-isls");
 
    simulationHelper->CreateSatScenario();
 
    simulationHelper->GetTrafficHelper()->AddCbrTraffic(
        SatTrafficHelper::FWD_LINK,
        SatTrafficHelper::UDP,
        MilliSeconds(10),
        512,
        NodeContainer(Singleton<SatTopology>::Get()->GetGwUserNode(0)),
        Singleton<SatTopology>::Get()->GetUtUserNodes(),
        Seconds(1.0),
        Seconds(29.0),
        Seconds(0));
 
    simulationHelper->GetTrafficHelper()->AddCbrTraffic(
        SatTrafficHelper::RTN_LINK,
        SatTrafficHelper::UDP,
        MilliSeconds(10),
        512,
        NodeContainer(Singleton<SatTopology>::Get()->GetGwUserNode(0)),
        Singleton<SatTopology>::Get()->GetUtUserNodes(),
        Seconds(1.0),
        Seconds(29.0),
        Seconds(0));
 
    Ptr<SatStatsHelperContainer> s = simulationHelper->GetStatisticsContainer();
 
    s->AddGlobalFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddGlobalRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerGwFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerGwRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerSatFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerSatRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerBeamFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerBeamRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerUtFwdAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
    s->AddPerUtRtnAppThroughput(SatStatsHelper::OUTPUT_SCATTER_FILE);
 
    simulationHelper->SetSimulationTime(Seconds(10));
    simulationHelper->RunSimulation();
 
    Simulator::Destroy();
 
    // Global throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-global-fwd-app-throughput-scatter-0.txt",
                  5,
                  6144);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-global-rtn-app-throughput-scatter-0.txt",
                  5,
                  6144);
 
    // Per satellite throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-sat-fwd-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-sat-fwd-app-throughput-scatter-2.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-sat-rtn-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-sat-rtn-app-throughput-scatter-2.txt",
                  5,
                  4096);
 
    // Per beam throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-fwd-app-throughput-scatter-1-30.txt",
                  5,
                  0);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-fwd-app-throughput-scatter-1-43.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-fwd-app-throughput-scatter-2-30.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-fwd-app-throughput-scatter-2-43.txt",
                  5,
                  0);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-rtn-app-throughput-scatter-1-30.txt",
                  5,
                  0);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-rtn-app-throughput-scatter-1-43.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-rtn-app-throughput-scatter-2-30.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-beam-rtn-app-throughput-scatter-2-43.txt",
                  5,
                  0);
 
    // Per GW throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-gw-fwd-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-gw-fwd-app-throughput-scatter-2.txt",
                  5,
                  4096);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-gw-rtn-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-gw-rtn-app-throughput-scatter-2.txt",
                  5,
                  4096);
 
    // Per UT throughput
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-ut-fwd-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-ut-fwd-app-throughput-scatter-2.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-ut-fwd-app-throughput-scatter-3.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-ut-rtn-app-throughput-scatter-1.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-ut-rtn-app-throughput-scatter-2.txt",
                  5,
                  2048);
    TestFileValue("contrib/satellite/data/sims/test-sat-constellation/test3/"
                  "stat-per-ut-rtn-app-throughput-scatter-3.txt",
                  5,
                  2048);
}

class SatConstellationTest4 : public TestCase
{
  public:
    SatConstellationTest4();
    virtual ~SatConstellationTest4();
 
  private:
    virtual void DoRun(void);
};
 
// Add some help text to this case to describe what it is intended to test
SatConstellationTest4::SatConstellationTest4()
    : TestCase(
          "This case tests that a topology with hundreds of LEO satellites is correctly loaded.")
{
}
 
// This destructor does nothing but we include it as a reminder that
// the test case should clean up after itself
SatConstellationTest4::~SatConstellationTest4()
{
}
 
//
// SatConstellationTest4 TestCase implementation
//
void
SatConstellationTest4::DoRun(void)
{
    Config::Reset();
 
    // Set simulation output details
    SatEnvVariables::GetInstance()->DoInitialize();
    SatEnvVariables::GetInstance()->SetOutputVariables("test-sat-constellation", "test4", true);

    Config::SetDefault("ns3::SatConf::ForwardLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));
    Config::SetDefault("ns3::SatConf::ReturnLinkRegenerationMode",
                       EnumValue(SatEnums::REGENERATION_NETWORK));

    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionEachRequest", BooleanValue(false));
    Config::SetDefault("ns3::SatSGP4MobilityModel::UpdatePositionPeriod",
                       TimeValue(MilliSeconds(10)));
    Config::SetDefault("ns3::SatOrbiterHelper::IslArbiterType", EnumValue(SatEnums::UNICAST));
    Config::SetDefault("ns3::SatHelper::GwUsers", UintegerValue(3));

    Config::SetDefault("ns3::SatHelper::BeamNetworkAddress", Ipv4AddressValue("20.1.0.0"));
    Config::SetDefault("ns3::SatHelper::GwNetworkAddress", Ipv4AddressValue("10.1.0.0"));
    Config::SetDefault("ns3::SatHelper::UtNetworkAddress", Ipv4AddressValue("250.1.0.0"));
 
    Ptr<SimulationHelper> simulationHelper =
        CreateObject<SimulationHelper>("test-sat-constellation");
 
    // Creating the reference system.
    simulationHelper->SetBeamSet({1, 43, 60, 64});
    simulationHelper->SetUserCountPerUt(5);
 
    simulationHelper->LoadScenario("constellation-telesat-351-sats");
 
    simulationHelper->CreateSatScenario();
 
    NodeContainer sats = Singleton<SatTopology>::Get()->GetOrbiterNodes();
    NodeContainer gws = Singleton<SatTopology>::Get()->GetGwNodes();
    NodeContainer uts = Singleton<SatTopology>::Get()->GetUtNodes();
    NodeContainer gwUsers = Singleton<SatTopology>::Get()->GetGwUserNodes();
    NodeContainer utUsers = Singleton<SatTopology>::Get()->GetUtUserNodes();
 
    uint32_t countNetDevices = 0;
    uint32_t countIslNetDevice = 0;
    for (uint32_t i = 0; i < sats.GetN(); i++)
    {
        for (uint32_t j = 0; j < sats.Get(i)->GetNDevices(); j++)
        {
            if (DynamicCast<SatOrbiterNetDevice>(sats.Get(i)->GetDevice(j)) != nullptr)
            {
                countNetDevices += 1;
                countIslNetDevice += DynamicCast<SatOrbiterNetDevice>(sats.Get(i)->GetDevice(j))
                                         ->GetIslsNetDevices()
                                         .size();
            }
        }
    }
 
    NS_TEST_ASSERT_MSG_EQ(sats.GetN(), 351, "Topology must contain 351 satellites");
    NS_TEST_ASSERT_MSG_EQ(countNetDevices,
                          351,
                          "Topology must contain 351 satellite Orbiter Net Devices");
    NS_TEST_ASSERT_MSG_EQ(countIslNetDevice,
                          1404,
                          "Topology must contain 1404 (351*4) satellite ISL Net Devices");
    NS_TEST_ASSERT_MSG_EQ(gws.GetN(), 2, "Topology must contain 2 GWs");
    NS_TEST_ASSERT_MSG_EQ(uts.GetN(), 3, "Topology must contain 3 UTs");
    NS_TEST_ASSERT_MSG_EQ(gwUsers.GetN(), 3, "Topology must contain 3 GW users");
    NS_TEST_ASSERT_MSG_EQ(utUsers.GetN(), 15, "Topology must contain 15 UT users");
 
    for (uint32_t i = 0; i < sats.GetN(); i++)
    {
        Ptr<SatIslArbiter> arbiter =
            DynamicCast<SatOrbiterNetDevice>(sats.Get(i)->GetDevice(0))->GetArbiter();
        Ptr<SatIslArbiterUnicast> arbiterUnicast = DynamicCast<SatIslArbiterUnicast>(arbiter);
        NS_TEST_ASSERT_MSG_NE(arbiterUnicast,
                              nullptr,
                              "Arbiter of satellite " << i
                                                      << " is not of type SatIslArbiterUnicast");
        for (uint32_t j = 0; j < sats.GetN(); j++)
        {
            if (i != j)
            {
                NS_TEST_ASSERT_MSG_NE(arbiterUnicast->Decide(i, j, nullptr),
                                      -1,
                                      "No entry in arbiter of satellite " << i << " to satellite "
                                                                          << j);
            }
        }
    }
 
    Simulator::Destroy();
}
 
// The TestSuite class names the TestSuite as sat-constellation-test, identifies what type of
// TestSuite (Type::SYSTEM), and enables the TestCases to be run. Typically, only the constructor
// for this class must be defined
//
class SatConstellationTestSuite : public TestSuite
{
  public:
    SatConstellationTestSuite();
};
 
SatConstellationTestSuite::SatConstellationTestSuite()
    : TestSuite("sat-constellation-test", Type::SYSTEM)
{
    AddTestCase(new SatConstellationTest1, TestCase::Duration::QUICK); // Test topology loading
    AddTestCase(new SatConstellationTest2,
                TestCase::Duration::QUICK); // Test good throughputs without ISLs
    AddTestCase(new SatConstellationTest3,
                TestCase::Duration::QUICK); // Test good throughputs with ISLs
    AddTestCase(new SatConstellationTest4,
                TestCase::Duration::QUICK); // Test topology loading with huge constellation
}
 
// Allocate an instance of this TestSuite
static SatConstellationTestSuite satConstellationTestSuite;