def METHOD_NAME(self):
    log = get_logger(level='warning')
    for folder, prefix, freqs in CASES:
        bdf = BDF(debug=False, log=log)
        basepath = os.path.join(TEST_PATH, folder)
        bdf_filename = os.path.join(basepath, f'{prefix}.bdf')
        op2_filename = os.path.join(basepath, f'{prefix}.op2')
        bdf.read_bdf(bdf_filename)
        op2 = read_op2(
            op2_filename,
            debug=False, log=log,
            exclude_results=['element_forces', 'strain'],
        )
        op2_new = data_in_material_coord(bdf, op2)
        for freq in freqs:
            for vecname in stress_vectors:
                vector = getattr(op2_new, vecname).get(1)
                if vector is None:
                    continue
                if 'center' in prefix:
                    name = os.path.join(basepath, f'{vecname}_center_freq_{freq:1.1f}.txt')
                else:
                    name = os.path.join(basepath, f'{vecname}_corner_freq_{freq:1.1f}.txt')
                if not os.path.isfile(name):
                    raise AssertionError(f'Not found reference result {name}')
                ref_stress = np.loadtxt(name)
                mag = ref_stress[:, 1::2]
                phase = ref_stress[:, 2::2]
                if freq == 1.0:
                    data = vector.data[0]
                elif freq == 9.5:
                    data = vector.data[17]
                eids = get_eids_from_op2_vector(vector)
                check = eids != 0
                if 'cquad8' in vecname:
                    assert np.allclose(np.abs(data[check][0::10, :]), mag[0::10], rtol=RTOL)
                    assert np.allclose(calc_phasedeg(data[check][1::10, :]), phase[1::10], rtol=RTOL)
                else:
                    assert np.allclose(np.abs(data[check]), mag, rtol=RTOL)
                    assert np.allclose(calc_phasedeg(data[check]), phase, rtol=RTOL)