simulation.py

# -*- coding: utf-8 -*-
"""
    Pricing_aware_MEC_offloading.simulation
    ~~~~~~~~~~~~~~~~~~~~~~~~~

    Simulation for the Pricing_aware_MEC_offloading

    :copyright: (c) 2019 by Giorgos Mitsis.
    :license: MIT License, see LICENSE for more details.
"""

import numpy as np
import matplotlib.pyplot as plt

from parameters import *
from helper_functions import *
from pricing_aware_MEC_offloading import *

import time
import dill

# Keep only three decimal places when printing numbers
np.set_printoptions(formatter={'float': lambda x: "{0:0.3f}".format(x)})

case = {"users": "homo"}
# Ns = [1,2,5,10,25,50,75,100]
Ns = [25]
# cpars = np.linspace(0.1,0.9,33)
cpars = [0.5]

ans = [0.2]
# ans = np.linspace(0.1,1,10)

kns = [1.2]
# kns = [0.31]
# kns = np.linspace(0.2,2,10)

for N in Ns:
    # Set random parameter in order to generate the same parameters
    print("Generating new parameters")
    np.random.seed()
    params = set_parameters(case, N)

    print("Number of users: " + str(params["N"]))

    for an in ans:

        params["an"] = an * np.ones(N)

        print("Sensitivity an: " + str(params["an"][0]))

        for kn in kns:

            params["kn"] = kn * np.ones(N)

            print("Sensitivity kn: " + str(params["kn"][0]))

            for cpar in cpars:

                params["cpar"] = cpar
                params["c"] = cpar * params["bn"]/params["dn"] * (1 - (1/(params["tn"]*params["en"])))

                print("Cost parameter: "+ str(params["cpar"]))

                # if bpars is -1 we don't want constant offloading
                bpars = [-1,0,0.5,1] if params["CONSTANT_OFFLOADING"] else [-1]

                while bpars:
                    bpar = bpars.pop()
                    params["bpar"] = bpar

                    # if bpar is -1 then we don't want constant offloading
                    if bpar == -1:
                        params["CONSTANT_OFFLOADING"] = False
                        params["bpar"] = 0

                    for repetition in range(1):
                        print("Repetition no: " + str(repetition+1))

                        results = {}

                        start = time.time()

                        # Run main simulation
                        results = main(params)
                        # check_all_parameters(**params)
                        # check_best_parameters(**params)

                        end = time.time()
                        running_time = end - start
                        print("Time of simulation:")
                        print(running_time)


                        results["N"] = N
                        results["time"] = running_time
                        results["repetition"] = repetition

                        if params["SAVE_RESULTS"] == True:
                            if params["CONSTANT_OFFLOADING"]:
                                constant_str = "_b_constant_" + str(round(bpar,3))
                            else:
                                constant_str = ""
                            outfile = 'saved_runs/results/individual/' + case["users"] + constant_str + "_N_" + str(N) + "_an_" + str(round(an,3)) + "_kn_" + str(round(kn,3)) + "_c_" + str(round(cpar,3)) + "_" + str(repetition)

                            with open(outfile, 'wb') as fp:
                                dill.dump(results, fp)

if params["GENERATE_FIGURES"] and not params["SAVE_FIGS"]:
    plt.show()