import numpy as np
from SimPleAC_save import save_obj, load_obj
from SimPleAC_pof_simulate import pof_parameters
import pickle as pickle
from gpkit.small_scripts import mag

if __name__ == "__main__":
    # Retrieving pof parameters
    [
        model, methods, gammas, number_of_iterations,
        min_num_of_linear_sections, max_num_of_linear_sections, verbosity,
        linearization_tolerance, number_of_time_average_solves,
        uncertainty_sets, nominal_solution, directly_uncertain_vars_subs,
        parallel, nominal_number_of_constraints, nominal_solve_time
    ] = pof_parameters()

    # # Loading margins results
    margin = {}
    nmargins = len(gammas)
    margins = np.linspace(0., 1., nmargins)

    margin['solutions'] = {}
    for i in range(nmargins):
        margin['solutions'][margins[i]] = pickle.load(
            open("marginResults/" + str(margins[i]), 'rb'))
    margin['number_of_constraints'] = load_obj('marginnumber_of_constraints',
                                               'marginResults')
    margin['simulation_results'] = load_obj('marginsimulation_results',
                                            'marginResults')
Пример #2
0
from builtins import range
import numpy as np
from SimPleAC_setup import SimPleAC_setup
from SimPleAC_save import save_obj
from gpkit.small_scripts import mag
from gpkit import Model

from SimPleAC_pof_simulate import pof_parameters
from robust.simulations.simulate import RobustGPTools
import pickle as pickle

if __name__ == '__main__':
    [model, methods, gammas, number_of_iterations,
    min_num_of_linear_sections, max_num_of_linear_sections, verbosity, linearization_tolerance,
    number_of_time_average_solves, uncertainty_sets, nominal_solution, directly_uncertain_vars_subs, parallel,
    nominal_number_of_constraints, nominal_solve_time] = pof_parameters()

    # Loading directly_uncertain_vars_subs
    try:
        pickle_in = open('directly_uncertain_dict.pickle', 'rb')
        directly_uncertain_vars_subs = pickle.load(pickle_in)
        pickle_in.close()
    except:
        print('Warning: Please run pof_simulate first for consistent MC results.')


    nGammas = len(gammas)
    solutions = {}
    simulation_results = {}
    number_of_constraints = {}
    for i in range(nGammas):