<gh_stars>0
import pathlib
import re

import numpy as np
import pytest
from scipy.io import loadmat

from pysoc.bnb import pb3wc

BNB_MAT_FILES_FOLDERPATH = pathlib.Path(__file__).parent.parent / 'mat_files'
PB3WC_FOLDERPATH = BNB_MAT_FILES_FOLDERPATH / 'pb3wc'
PB3WC_INPUT_FOLDERPATH = PB3WC_FOLDERPATH / 'input'
PB3WC_RESULTS_FOLDERPATH = PB3WC_FOLDERPATH / 'results'
CAO_INPUT_FILENAME = 'ny40_nu15_nd5_input.mat'

_matfiles_list = [mat.name for mat in PB3WC_RESULTS_FOLDERPATH.glob('*.mat')]
_param_tuple_list = [(CAO_INPUT_FILENAME, res_file)
                     for res_file in _matfiles_list]

# pytest routine to check results against mat files generated by matlab routine


@pytest.mark.parametrize("input_filename,result_filename", _param_tuple_list)
def test_same_result(input_filename, result_filename):
    """Test routine that verifies that results from python routine are in
    conformity with matlab results.
    """
    input_contents = loadmat(PB3WC_INPUT_FOLDERPATH / input_filename)
    results_contents = loadmat(PB3WC_RESULTS_FOLDERPATH / result_filename)

    # input variables
    Gy = input_contents["G"]
    Gyd = input_contents["Gd"]
    Juu = input_contents["Juu"]
    Jud = input_contents["Jud"]
    Wd = input_contents["Wd"]
    Wd = np.diag(Wd)
    Wny = input_contents["Wn"]
    Wny = np.diag(Wny)

    # matlab results
    res_to_load = ['B', 'sset', 'ops', 'flag']
    B_mat, sset_mat, ops_mat, flag_mat = tuple(
        results_contents[var] for var in res_to_load)

    # generate python results
    n = int(re.search(r'_n(\d+)_', result_filename).group(1))
    nc = int(re.search(r'_nc(\d+)', result_filename).group(1))
    B_py, sset_py, ops_py, _, flag_py = pb3wc(Gy, Gyd, Wd, Wny, Juu, Jud, n,
                                              nc=nc)

    assert np.allclose(B_py, B_mat.flatten()), "B failed"
    assert np.allclose(sset_py, sset_mat), "sset failed"
    assert np.allclose(ops_py, ops_mat.flatten()), "ops failed"
    assert flag_py == flag_mat, 'flag failed'


# if __name__ == "__main__":
#     SOC_FOLDERPATH = pathlib.Path(
#         __file__).parent.parent.parent / 'soc' / 'mat_files'
#     SOC_INPUT = SOC_FOLDERPATH / 'input'
#     SOC_FILENAME = 'matrices_self_opt_100pts.mat'

#     mat_contents = loadmat(SOC_INPUT / SOC_FILENAME)

#     # input variables
#     Gy = mat_contents["Gy"]
#     Gyd = mat_contents["Gyd"]
#     Juu = mat_contents["Juu"]
#     Jud = mat_contents["Jud"]
#     Wd = mat_contents["Wd"]
#     Wny = mat_contents["Wny"]
#     Wd = np.diag(Wd)
#     Wny = np.diag(Wny)

#     n = 2
#     nc = 45

#     B_py, sset_py, ops_py, _, flag_py = pb3wc(Gy, Gyd, Wd, Wny, Juu, Jud, n,
#                                               nc=nc)
#     a = 1

import pathlib
import re

import numpy as np
import pytest
from scipy.io import loadmat

from pysoc.bnb import pb3wc

BNB_MAT_FILES_FOLDERPATH = pathlib.Path(__file__).parent.parent / 'mat_files'
PB3WC_FOLDERPATH = BNB_MAT_FILES_FOLDERPATH / 'pb3wc'
PB3WC_INPUT_FOLDERPATH = PB3WC_FOLDERPATH / 'input'
PB3WC_RESULTS_FOLDERPATH = PB3WC_FOLDERPATH / 'results'
CAO_INPUT_FILENAME = 'ny40_nu15_nd5_input.mat'

_matfiles_list = [mat.name for mat in PB3WC_RESULTS_FOLDERPATH.glob('*.mat')]
_param_tuple_list = [(CAO_INPUT_FILENAME, res_file)
                     for res_file in _matfiles_list]

# pytest routine to check results against mat files generated by matlab routine


@pytest.mark.parametrize("input_filename,result_filename", _param_tuple_list)
def test_same_result(input_filename, result_filename):
    """Test routine that verifies that results from python routine are in
    conformity with matlab results.
    """
    input_contents = loadmat(PB3WC_INPUT_FOLDERPATH / input_filename)
    results_contents = loadmat(PB3WC_RESULTS_FOLDERPATH / result_filename)

    # input variables
    Gy = input_contents["G"]
    Gyd = input_contents["Gd"]
    Juu = input_contents["Juu"]
    Jud = input_contents["Jud"]
    Wd = input_contents["Wd"]
    Wd = np.diag(Wd)
    Wny = input_contents["Wn"]
    Wny = np.diag(Wny)

    # matlab results
    res_to_load = ['B', 'sset', 'ops', 'flag']
    B_mat, sset_mat, ops_mat, flag_mat = tuple(
        results_contents[var] for var in res_to_load)

    # generate python results
    n = int(re.search(r'_n(\d+)_', result_filename).group(1))
    nc = int(re.search(r'_nc(\d+)', result_filename).group(1))
    B_py, sset_py, ops_py, _, flag_py = pb3wc(Gy, Gyd, Wd, Wny, Juu, Jud, n,
                                              nc=nc)

    assert np.allclose(B_py, B_mat.flatten()), "B failed"
    assert np.allclose(sset_py, sset_mat), "sset failed"
    assert np.allclose(ops_py, ops_mat.flatten()), "ops failed"
    assert flag_py == flag_mat, 'flag failed'


# if __name__ == "__main__":
#     SOC_FOLDERPATH = pathlib.Path(
#         __file__).parent.parent.parent / 'soc' / 'mat_files'
#     SOC_INPUT = SOC_FOLDERPATH / 'input'
#     SOC_FILENAME = 'matrices_self_opt_100pts.mat'

#     mat_contents = loadmat(SOC_INPUT / SOC_FILENAME)

#     # input variables
#     Gy = mat_contents["Gy"]
#     Gyd = mat_contents["Gyd"]
#     Juu = mat_contents["Juu"]
#     Jud = mat_contents["Jud"]
#     Wd = mat_contents["Wd"]
#     Wny = mat_contents["Wny"]
#     Wd = np.diag(Wd)
#     Wny = np.diag(Wny)

#     n = 2
#     nc = 45

#     B_py, sset_py, ops_py, _, flag_py = pb3wc(Gy, Gyd, Wd, Wny, Juu, Jud, n,
#                                               nc=nc)
#     a = 1
# pytest routine to check results against mat files generated by matlab routine Test routine that verifies that results from python routine are in     conformity with matlab results. # input variables # matlab results # generate python results # if __name__ == "__main__": #     SOC_FOLDERPATH = pathlib.Path( #         __file__).parent.parent.parent / 'soc' / 'mat_files' #     SOC_INPUT = SOC_FOLDERPATH / 'input' #     SOC_FILENAME = 'matrices_self_opt_100pts.mat' #     mat_contents = loadmat(SOC_INPUT / SOC_FILENAME) #     # input variables #     Gy = mat_contents["Gy"] #     Gyd = mat_contents["Gyd"] #     Juu = mat_contents["Juu"] #     Jud = mat_contents["Jud"] #     Wd = mat_contents["Wd"] #     Wny = mat_contents["Wny"] #     Wd = np.diag(Wd) #     Wny = np.diag(Wny) #     n = 2 #     nc = 45 #     B_py, sset_py, ops_py, _, flag_py = pb3wc(Gy, Gyd, Wd, Wny, Juu, Jud, n, #                                               nc=nc) #     a = 1