Esempio n. 1
0
def simulator_factory(
        config_dict,
        benchmark_os_path,
        property_checker,
        plt,
        plant_pvt_init_data,
        parallel=False,
        test_params=None,
        sim_args=None
        ):

    # ##!!##logger.debug('requested simulator creation')

    sim_type = config_dict['plant_description']
    if sim_type == 'matlab':
        logger.info('creating matlab simulator')

        # get the m_file's path

        m_file_path = config_dict['plant_path']
        abs_m_file_path = fp.construct_path(m_file_path, benchmark_os_path)
        if fp.validate_file_names([abs_m_file_path]):

            # return MatlabSim(m_file_path, benchmark_os_path, parallel)

            return MEngPy(m_file_path, benchmark_os_path, parallel, sim_args)
        else:
            raise err.FileNotFound('file does not exist: ' + m_file_path)
    elif sim_type == 'simulink':
        return SimulinkSim()
    elif sim_type == 'python':
        logger.info('creating Native Python simulator')

        # get the file path

        python_file_path = config_dict['plant_path']
        #(module_name, file_ext) = python_file_path.split('.')
        (module_name, file_ext) = fp.split_filename_ext(python_file_path)
        if file_ext != 'py':
            raise err.Fatal('Python file extension py expected, found: {}'.format(file_ext))
        module_path = fp.construct_path(python_file_path, benchmark_os_path)
        if fp.validate_file_names([module_path]):
            return NativeSim(module_name, module_path,
                    property_checker, plt, plant_pvt_init_data, parallel)
        else:
            raise err.FileNotFound('file does not exist: ' + python_file_path)
    elif sim_type == 'test':
        return TestSim(test_params)
    else:
        raise err.Fatal('unknown sim type : {}'.format(sim_type))
Esempio n. 2
0
def simulator_factory(
        config_dict,
        benchmark_os_path,
        property_checker,
        plt,
        plant_pvt_init_data,
        parallel=False,
        test_params=None,
        sim_args=None
        ):

    # ##!!##logger.debug('requested simulator creation')

    sim_type = config_dict['plant_description']
    if sim_type == 'matlab':
        logger.info('creating matlab simulator')

        # get the m_file's path

        m_file_path = config_dict['plant_path']
        abs_m_file_path = fp.construct_path(m_file_path, benchmark_os_path)
        if fp.validate_file_names([abs_m_file_path]):

            # return MatlabSim(m_file_path, benchmark_os_path, parallel)

            return MEngPy(m_file_path, benchmark_os_path, parallel, sim_args)
        else:
            raise err.FileNotFound('file does not exist: ' + m_file_path)
    elif sim_type == 'simulink':
        return SimulinkSim()
    elif sim_type == 'python':
        logger.info('creating Native Python simulator')

        # get the file path

        python_file_path = config_dict['plant_path']
        #(module_name, file_ext) = python_file_path.split('.')
        (module_name, file_ext) = fp.split_filename_ext(python_file_path)
        if file_ext != 'py':
            raise err.Fatal('Python file extension py expected, found: {}'.format(file_ext))
        module_path = fp.construct_path(python_file_path, benchmark_os_path)
        if fp.validate_file_names([module_path]):
            return NativeSim(module_name, module_path,
                    property_checker, plt, plant_pvt_init_data, parallel)
        else:
            raise err.FileNotFound('file does not exist: ' + python_file_path)
    elif sim_type == 'test':
        return TestSim(test_params)
    else:
        raise err.Fatal('unknown sim type : {}'.format(sim_type))
Esempio n. 3
0
def parse(file_path, plant_pvt_init_data):
    test_des_file = file_path
    path = fp.get_abs_base_path(file_path)
    sys_name = fp.split_filename_ext(fp.get_file_name_from_path(file_path))[0]

    test_dir = fp.get_abs_base_path(test_des_file)
    SYS.path.append(test_dir)

    try:
        sut = imp.load_source('test_des', test_des_file)
    except IOError as e:
        print('File not found: {}'.format(test_des_file))
        raise e

    sut.plant_pvt_init_data = plant_pvt_init_data

    try:
        assert (len(sut.initial_set) == 2)
        assert (len(sut.error_set) == 2)
        assert (len(sut.ci) == 2)
        assert (len(sut.pi) == 2)

        assert (len(sut.ci[0]) == len(sut.ci[1]))
        assert (len(sut.pi[0]) == len(sut.pi[1]))

        assert (len(sut.initial_set[0]) == len(sut.initial_set[1]))
        assert (len(sut.error_set[0]) == len(sut.error_set[1]))

        assert (len(sut.initial_set[0]) == len(sut.error_set[0]))

        assert (len(sut.initial_set[0]) == len(sut.grid_eps))

        num_dims = NumDims(
            si=len(sut.initial_controller_integer_state),
            sf=len(sut.initial_controller_float_state),
            s=len(sut.initial_controller_integer_state) +
            len(sut.initial_controller_float_state),
            x=len(sut.initial_set[0]),
            u=sut.num_control_inputs,
            ci=len(sut.ci[0]),
            pi=len(sut.pi[0]),
            d=len(sut.initial_discrete_state),
            pvt=0,
        )

        if not hasattr(sut, 'ci_grid_eps') or len(sut.ci_grid_eps) == 0:
            sut.ci_grid_eps = []

#         if hasattr(sut, 'ci_grid_eps'):
#             if len(sut.ci_grid_eps) == 0:
#                 sut.ci_grid_eps = None
#         else:
#             sut.ci_grid_eps = None

        if not hasattr(sut, 'pi_grid_eps') or len(sut.pi_grid_eps) == 0:
            sut.pi_grid_eps = []

#         if hasattr(sut, 'pi_grid_eps'):
#             if len(sut.pi_grid_eps) == 0:
#                 sut.pi_grid_eps = None
#         else:
#             sut.pi_grid_eps = None

        if num_dims.ci == 0:
            ci = None
        else:
            ci = cns.IntervalCons(np.array(sut.ci[0]), np.array(sut.ci[1]))

        if num_dims.pi == 0:
            pi = None
        else:
            pi = cns.IntervalCons(np.array(sut.pi[0]), np.array(sut.pi[1]))

        assert ((pi is None and sut.pi_grid_eps == [])
                or (pi is not None and sut.pi_grid_eps != []))
        assert ((ci is None and sut.ci_grid_eps == [])
                or (ci is not None and sut.ci_grid_eps != []))

        if pi is not None:
            assert (len(sut.pi[0]) == len(sut.pi_grid_eps))

        if ci is not None:
            assert (len(sut.ci[0]) == len(sut.ci_grid_eps))

        plant_config_dict = {
            'plant_description': sut.plant_description,
            'plant_path': sut.plant_path,
            #'refinement_factor': sut.refinement_factor,
            'eps': np.array(sut.grid_eps, dtype=float),
            #'pi_eps': np.array(sut.pi_grid_eps, dtype=float),
            'num_samples': int(sut.num_samples),
            'delta_t': float(sut.delta_t),
            #'type': 'string',
            'type': 'value',
        }
        #print(plant_config_dict)
        #for k, d in plant_config_dict.iteritems():
        #    plant_config_dict[k] = str(d)

        # TODO: fix the list hack, once the controller ifc matures

        init_cons = cns.IntervalCons(np.array(sut.initial_set[0]),
                                     np.array(sut.initial_set[1]))

        init_cons_list = [init_cons]

        final_cons = cns.IntervalCons(np.array(sut.error_set[0]),
                                      np.array(sut.error_set[1]))
        T = sut.T
        delta_t = sut.delta_t

        if sut.controller_path is not None:
            controller_object_str = fp.split_filename_ext(
                sut.controller_path)[0]
            controller_path_dir_path = fp.construct_path(
                sut.controller_path_dir_path, path)
        else:
            controller_object_str = None
            controller_path_dir_path = None

        initial_discrete_state = sut.initial_discrete_state
        initial_controller_state = (sut.initial_controller_integer_state +
                                    sut.initial_controller_float_state)
        MAX_ITER = sut.MAX_ITER

        num_segments = int(np.ceil(T / delta_t))

        controller_path = sut.controller_path
        plant_pvt_init_data = sut.plant_pvt_init_data

        min_smt_sample_dist = sut.min_smt_sample_dist
        ci_grid_eps = sut.ci_grid_eps
        pi_grid_eps = sut.pi_grid_eps

        #num_sim_samples = sut.num_sim_samples
        #METHOD = sut.METHOD
        #symbolic_analyzer = sut.symbolic_analyzer
        #plot = sut.plot
        #MODE = sut.MODE
        #opts = Options(plot, MODE, num_sim_samples, METHOD, symbolic_analyzer)
    except AttributeError as e:
        raise MissingSystemDefError(e)

    # ROI is not mandatory and defaults to xi \in (-inf, inf)
    try:
        ROI = cns.IntervalCons(*sut.ROI)
        assert (ROI.dim == num_dims.x)
    except AttributeError:
        ROI = cns.top2ic(num_dims.x)

    sys = System(sys_name, controller_path, num_dims, plant_config_dict,
                 delta_t, controller_path_dir_path, controller_object_str,
                 path, plant_pvt_init_data, min_smt_sample_dist, ci_grid_eps,
                 pi_grid_eps)
    prop = Property(T, init_cons_list, init_cons, final_cons, ci, pi,
                    initial_discrete_state, initial_controller_state, MAX_ITER,
                    num_segments, ROI)

    print('system loaded...')
    return sys, prop
Esempio n. 4
0
def parse(file_path, plant_pvt_init_data):
    test_des_file = file_path
    path = fp.get_abs_base_path(file_path)
    sys_name = fp.split_filename_ext(
                fp.get_file_name_from_path(
                    file_path)
                )[0]

    test_dir = fp.get_abs_base_path(test_des_file)
    SYS.path.append(test_dir)

    try:
        sut = imp.load_source('test_des', test_des_file)
    except IOError as e:
        print('File not found: {}'.format(test_des_file))
        raise e

    sut.plant_pvt_init_data = plant_pvt_init_data

    try:
        assert(len(sut.initial_set) == 2)
        assert(len(sut.error_set) == 2)
        assert(len(sut.ci) == 2)
        assert(len(sut.pi) == 2)

        assert(len(sut.ci[0]) == len(sut.ci[1]))
        assert(len(sut.pi[0]) == len(sut.pi[1]))

        assert(len(sut.initial_set[0]) == len(sut.initial_set[1]))
        assert(len(sut.error_set[0]) == len(sut.error_set[1]))

        assert(len(sut.initial_set[0]) == len(sut.error_set[0]))

        assert(len(sut.initial_set[0]) == len(sut.grid_eps))

        num_dims = NumDims(
            si=len(sut.initial_controller_integer_state),
            sf=len(sut.initial_controller_float_state),
            s=len(sut.initial_controller_integer_state)+len(sut.initial_controller_float_state),
            x=len(sut.initial_set[0]),
            u=sut.num_control_inputs,
            ci=len(sut.ci[0]),
            pi=len(sut.pi[0]),
            d=len(sut.initial_discrete_state),
            pvt=0,
            )

        if not hasattr(sut, 'ci_grid_eps') or len(sut.ci_grid_eps) == 0:
            sut.ci_grid_eps = []

#         if hasattr(sut, 'ci_grid_eps'):
#             if len(sut.ci_grid_eps) == 0:
#                 sut.ci_grid_eps = None
#         else:
#             sut.ci_grid_eps = None

        if not hasattr(sut, 'pi_grid_eps') or len(sut.pi_grid_eps) == 0:
            sut.pi_grid_eps = []

#         if hasattr(sut, 'pi_grid_eps'):
#             if len(sut.pi_grid_eps) == 0:
#                 sut.pi_grid_eps = None
#         else:
#             sut.pi_grid_eps = None

        if num_dims.ci == 0:
            ci = None
        else:
            ci = cns.IntervalCons(np.array(sut.ci[0]), np.array(sut.ci[1]))

        if num_dims.pi == 0:
            pi = None
        else:
            pi = cns.IntervalCons(np.array(sut.pi[0]), np.array(sut.pi[1]))

        assert((pi is None and sut.pi_grid_eps == []) or
               (pi is not None and sut.pi_grid_eps != []))
        assert((ci is None and sut.ci_grid_eps == []) or
               (ci is not None and sut.ci_grid_eps != []))

        if pi is not None:
            assert(len(sut.pi[0]) == len(sut.pi_grid_eps))

        if ci is not None:
            assert(len(sut.ci[0]) == len(sut.ci_grid_eps))

        plant_config_dict = {'plant_description': sut.plant_description,
                             'plant_path': sut.plant_path,
                             #'refinement_factor': sut.refinement_factor,
                             'eps': np.array(sut.grid_eps, dtype=float),
                             #'pi_eps': np.array(sut.pi_grid_eps, dtype=float),
                             'num_samples': int(sut.num_samples),
                             'delta_t': float(sut.delta_t),
                             #'type': 'string',
                             'type': 'value',
                             }
        #print(plant_config_dict)
        #for k, d in plant_config_dict.iteritems():
        #    plant_config_dict[k] = str(d)

        # TODO: fix the list hack, once the controller ifc matures

        init_cons = cns.IntervalCons(np.array(sut.initial_set[0]),
                                     np.array(sut.initial_set[1]))

        init_cons_list = [init_cons]

        final_cons = cns.IntervalCons(np.array(sut.error_set[0]),
                                      np.array(sut.error_set[1]))
        T = sut.T
        delta_t = sut.delta_t

        if sut.controller_path is not None:
            controller_object_str = fp.split_filename_ext(sut.controller_path)[0]
            controller_path_dir_path = fp.construct_path(sut.controller_path_dir_path, path)
        else:
            controller_object_str = None
            controller_path_dir_path = None

        initial_discrete_state = sut.initial_discrete_state
        initial_controller_state = (sut.initial_controller_integer_state +
                                    sut.initial_controller_float_state)
        MAX_ITER = sut.MAX_ITER

        num_segments = int(np.ceil(T / delta_t))

        controller_path = sut.controller_path
        plant_pvt_init_data = sut.plant_pvt_init_data

        min_smt_sample_dist = sut.min_smt_sample_dist
        ci_grid_eps = sut.ci_grid_eps
        pi_grid_eps = sut.pi_grid_eps

        #num_sim_samples = sut.num_sim_samples
        #METHOD = sut.METHOD
        #symbolic_analyzer = sut.symbolic_analyzer
        #plot = sut.plot
        #MODE = sut.MODE
        #opts = Options(plot, MODE, num_sim_samples, METHOD, symbolic_analyzer)
    except AttributeError as e:
        raise MissingSystemDefError(e)

    sys = System(sys_name, controller_path, num_dims, plant_config_dict,
                 delta_t, controller_path_dir_path,
                 controller_object_str, path, plant_pvt_init_data,
                 min_smt_sample_dist, ci_grid_eps,
                 pi_grid_eps)
    prop = Property(T, init_cons_list, init_cons, final_cons, ci,
                    pi, initial_discrete_state, initial_controller_state,
                    MAX_ITER, num_segments)

    print('system loaded...')
    return sys, prop