def test_api_optim(cleanup): results_folder_path = pth.join(RESULTS_FOLDER_PATH, "api_optim") configuration_file_path = pth.join(results_folder_path, "oad_process.yml") # Generation of configuration file ---------------------------------------- api.generate_configuration_file(configuration_file_path, True) # Generation of inputs ---------------------------------------------------- # We get the same inputs as in tutorial notebook source_xml = pth.join( root_folder_path, "src", "fastoad", "notebooks", "01_tutorial", "data", "CeRAS01_baseline.xml", ) api.generate_inputs(configuration_file_path, source_xml, overwrite=True) # Run optim --------------------------------------------------------------- problem = api.optimize_problem(configuration_file_path, True) assert not problem.optim_failed # Check that weight-performances loop correctly converged _check_weight_performance_loop(problem) # Design Variable assert_allclose(problem["data:geometry:wing:aspect_ratio"], 14.52, atol=1e-2) # Constraint assert_allclose(problem["data:geometry:wing:span"], 44.88, atol=1e-2) # Objective assert_allclose(problem["data:mission:sizing:needed_block_fuel"], 18900.0, atol=1)
def test_api_eval_mission(cleanup): results_folder_path = pth.join(RESULTS_FOLDER_PATH, "api_eval_mission") configuration_file_path = pth.join(results_folder_path, "oad_process.yml") api._PROBLEM_CONFIGURATOR = MissionConfigurator() # Generation of configuration file ---------------------------------------- api.generate_configuration_file(configuration_file_path, True) # Generation of inputs ---------------------------------------------------- # We get the same inputs as in tutorial notebook source_xml = pth.join( root_folder_path, "src", "fastoad", "notebooks", "01_tutorial", "data", "CeRAS01_baseline.xml", ) api.generate_inputs(configuration_file_path, source_xml, overwrite=True) # Run model --------------------------------------------------------------- problem = api.evaluate_problem(configuration_file_path, True) api._PROBLEM_CONFIGURATOR = None # Check that weight-performances loop correctly converged _check_weight_performance_loop(problem) assert_allclose(problem["data:handling_qualities:static_margin"], 0.05, atol=1e-2) assert_allclose(problem["data:geometry:wing:MAC:at25percent:x"], 17.149, atol=1e-2) assert_allclose(problem["data:weight:aircraft:MTOW"], 74846, atol=1) assert_allclose(problem["data:geometry:wing:area"], 126.581, atol=1e-2) assert_allclose(problem["data:geometry:vertical_tail:area"], 27.535, atol=1e-2) assert_allclose(problem["data:geometry:horizontal_tail:area"], 35.848, atol=1e-2) assert_allclose(problem["data:mission:sizing:needed_block_fuel"], 19495, atol=1)
def test_api_eval_breguet(cleanup): results_folder_path = pth.join(RESULTS_FOLDER_PATH, "api_eval_breguet") configuration_file_path = pth.join(results_folder_path, "oad_process.yml") # Generation of configuration file ---------------------------------------- api.generate_configuration_file(configuration_file_path, True) # Generation of inputs ---------------------------------------------------- # We get the same inputs as in tutorial notebook source_xml = pth.join(DATA_FOLDER_PATH, "CeRAS01_notebooks.xml") api.generate_inputs(configuration_file_path, source_xml, overwrite=True) # Run model --------------------------------------------------------------- problem = api.evaluate_problem(configuration_file_path, True) # Check that weight-performances loop correctly converged _check_weight_performance_loop(problem) assert_allclose(problem["data:handling_qualities:static_margin"], 0.05, atol=1e-2) assert_allclose(problem["data:geometry:wing:MAC:at25percent:x"], 17.149, atol=1e-2) assert_allclose(problem["data:weight:aircraft:MTOW"], 74892, atol=1) assert_allclose(problem["data:geometry:wing:area"], 126.732, atol=1e-2) assert_allclose(problem["data:geometry:vertical_tail:area"], 27.565, atol=1e-2) assert_allclose(problem["data:geometry:horizontal_tail:area"], 35.884, atol=1e-2) assert_allclose(problem["data:mission:sizing:needed_block_fuel"], 19527, atol=1)
def _generate_conf_file(args): """Generates a sample TOML file.""" try: api.generate_configuration_file(args.conf_file, args.force) except FastFileExistsError: if _query_yes_no( 'Configuration file "%s" already exists. Do you want to overwrite it?' % args.conf_file): api.generate_configuration_file(args.conf_file, True) else: print("No file written.")
def run_non_regression_test( conf_file, legacy_result_file, result_dir, use_xfoil=False, global_tolerance=1e-2, vars_to_check=None, specific_tolerance=5.0e-3, check_weight_perfo_loop=True, ): """ Convenience function for non regression tests :param conf_file: FAST-OAD configuration file :param legacy_result_file: reference data for inputs and outputs :param result_dir: relative name, folder will be in RESULTS_FOLDER_PATH :param use_xfoil: if True, XFOIL computation will be activated :param vars_to_check: variables that will be concerned by specific_tolerance :param specific_tolerance: test will fail if absolute relative error between computed and reference values is beyond this value for variables in vars_to_check :param global_tolerance: test will fail if absolute relative error between computed and reference values is beyond this value for ANY variable :param check_weight_perfo_loop: if True, consistency of weights will be checked """ results_folder_path = pth.join(RESULTS_FOLDER_PATH, result_dir) configuration_file_path = pth.join(results_folder_path, conf_file) # Copy of configuration file and generation of problem instance ------------------ api.generate_configuration_file(configuration_file_path) # just ensure folders are created... shutil.copy(pth.join(DATA_FOLDER_PATH, conf_file), configuration_file_path) configurator = FASTOADProblemConfigurator(configuration_file_path) configurator._set_configuration_modifier(XFOILConfigurator(use_xfoil)) # Generation of inputs ---------------------------------------- ref_inputs = pth.join(DATA_FOLDER_PATH, legacy_result_file) configurator.write_needed_inputs(ref_inputs) # Get problem with inputs ------------------------------------- problem = configurator.get_problem(read_inputs=True) problem.setup() # Run model --------------------------------------------------------------- problem.run_model() problem.write_outputs() om.view_connections( problem, outfile=pth.join(results_folder_path, "connections.html"), show_browser=False ) if check_weight_perfo_loop: _check_weight_performance_loop(problem) ref_data = DataFile(pth.join(DATA_FOLDER_PATH, legacy_result_file)) row_list = [] for ref_var in ref_data: try: value = problem.get_val(ref_var.name, units=ref_var.units)[0] except KeyError: continue row_list.append( { "name": ref_var.name, "units": ref_var.units, "ref_value": ref_var.value[0], "value": value, } ) df = pd.DataFrame(row_list) df["rel_delta"] = (df.value - df.ref_value) / df.ref_value df["rel_delta"][(df.ref_value == 0) & (abs(df.value) <= 1e-10)] = 0.0 df["abs_rel_delta"] = np.abs(df.rel_delta) pd.set_option("display.max_rows", None) pd.set_option("display.max_columns", None) pd.set_option("display.width", 1000) pd.set_option("display.max_colwidth", 120) print(df.sort_values(by=["abs_rel_delta"])) if vars_to_check is not None: for name in vars_to_check: assert_allclose(df.ref_value, df.value, rtol=global_tolerance) row = df.loc[df.name == name] assert_allclose(row.ref_value, row.value, rtol=specific_tolerance) else: assert np.all(df.abs_rel_delta < specific_tolerance)